Quantum metrology is being gradually studied for weak measurement systems. For weak measurement systems with thermal state pointer, we find that in the displacement space corresponding to imaginary weak values, the maximal QFI after successful postselection can attain the level of thermal fluctuations, without surpassing total QFI, and that QFI which increases with increasing temperature can constantly improve the measurement precision. These results are much better than that of weak measurement with pure state (i.e., Gaussian state) pointer. On the other hand, in Kerr nonlinear interaction systems with weak measurement, and by using thermal state pointer, we obtain in the phase space successful postselection and postselected measurements both achieve the Heisenberg limit of quantum metrology, and show weak measurement with thermal states only obtain classical Fisher information (CFI) which increases with increasing temperature and achieves classical enhanced scaling of N^2. Moreover, weak measurement with thermal states has an advantage over that with coherent states or mixed states of the light because generating these states with more large uncertainty are limited under the current technology, but thermal states with more large uncertainly are very easy to achieve with increasing temperature in nature, regardless of thermal states of the light or the matter.

The magnetic field response of the Mott-insulating honeycomb iridate Na$_{2}$IrO$_{3}$ is investigated using torque magnetometry measurements in magnetic fields up to 60 tesla. A peak-dip structure is observed in the torque response at magnetic fields corresponding to an energy scale close to the zigzag ordering ($\approx 15~K$) temperature. Using exact diagonalization calculations, we show that such a distinctive signature in the torque response constrains the effective spin models for these classes of Kitaev materials to ones with dominant ferromagnetic Kitaev interactions, while alternative models with dominant antiferromagnetic Kitaev interactions are excluded. We further show that at high magnetic fields, long range spin correlation functions decay rapidly, signaling a transition to a long-sought-after field-induced quantum spin liquid beyond the peak-dip structure. Kitaev systems are thus revealed to be excellent candidates for field-induced quantum spin liquids, similar physics having been suggested in another Kitaev material $\alpha-$RuCl$_{3}$.

We propose a scheme to entangle multiple material qubits through interaction with single photons via non-exciting processes associated with strongly coupling systems. The basic idea is based on the material state dependent reflection and transmission for the input photons. Thus, the material qubits in several systems can be entangled when one photon interacts with each system in cascade and the photon paths are mixed by the photon detection. The character of non-exciting of the material qubits does not change the state of material qubit and thus ensures the possibility of purifying entangled states by using more photons under realistic imperfect parameters. It also guarantees directly scaling up the scheme to entangle more qubits. Detailed analysis of fidelity and success probability of the scheme in the frame of an optical Fabry-P\'{e}rot (FP) cavity based strongly coupling system is presented. It is shown that a two-qubit entangled state with fidelity above 0.99 is promised with only two photons by using currently feasible experimental parameters. Our scheme can also be directly implemented on other strongly coupled system.

For a generic Higgs boson, measuring the relative sign and magnitude of its couplings with the $W$ and $Z$ bosons is essential in determining its origin. Such a test is also indispensable for the 125-GeV Higgs boson. We propose that the ratio of the $HWW$ and $HZZ$ couplings $\lambda_{WZ}$ can be directly determined through the $W^+W^-H$ production, where $H$ denotes a generic Higgs boson, owing to the tree-level interference effect. While this is impractical at the LHC due to the limited sensitivity, it can be done at future $e^+e^-$ colliders, such as a 500-GeV ILC with the beam polarization $P(e^-,e^+)=(-0.8,+0.3)$ in the $jj\ell^{\pm}bb$ and $\ell^{\pm}\ell^{\pm}\ell^{\mp}jj$ channels. The discovery potential of a general ratio and the power to discriminate it from the SM value are studied in detail. Combining the cross section of $e^+e^-\to W^+ W^- H$ with the measurements of $HZZ$ coupling at the HL-LHC, one can further improve the sensitivity of $\lambda_{WZ}$.

Based on the maximum likelihood estimation principle, we derive a collaborative estimation framework that fuses several different estimators and yields a better estimate. Applying it to compressive sensing (CS), we propose a collaborative CS (CCS) scheme consisting of a bank of $K$ CS systems that share the same sensing matrix but have different sparsifying dictionaries. This CCS system is expected to yield better performance than each individual CS system, while requiring the same time as that needed for each individual CS system when a parallel computing strategy is used. We then provide an approach to designing optimal CCS systems by utilizing a measure that involves both the sensing matrix and dictionaries and hence allows us to simultaneously optimize the sensing matrix and all the $K$ dictionaries. An alternating minimization-based algorithm is derived for solving the corresponding optimal design problem. With a rigorous convergence analysis, we show that the proposed algorithm is convergent. Experiments are carried out to confirm the theoretical results and show that the proposed CCS system yields significant improvements over the existing CS systems in terms of the signal recovery accuracy.

We consider the existence of solutions of the following $p(x)$-Laplacian Dirichlet problem without the Ambrosetti-Rabinowitz condition: $-\mbox{div}(|\nabla u|^{p(x)-2}\nabla u)=f(x,u) \text{ in }\Omega,$ and $u=0,\text{ on }\partial \Omega.$ We give a new growth condition and we point out its importance for checking the Cerami compactness condition. We prove the existence of solutions of the above problem via the critical point theory, and also provide some multiplicity properties. Our results extend previous work by Q. Zhang and C. Zhao, Existence of strong solutions of a $p(x)$-Laplacian Dirichlet problem without the Ambrosetti-Rabinowitz condition, Comp. Math. Appl. 69 (2015), 1-12, and we establish the existence of solutions under weaker hypotheses on the nonlinear term.

Due to their large bulk band gap, bismuthene, antimonene, and arsenene on a SiC substrate offer intriguing new opportunities for room-temperature quantum spin Hall (QSH) applications. Although edge states have been observed in the local density of states (LDOS) of bismuthene/SiC, there has been no experimental evidence until now that they are spin-polarized and topologically protected. Here, we show that for experimentally relevant armchair nanoribbons, we find a distinct behavior of the gap induced in the QSH edge states for out-of-plane magnetic fields (a few meV) versus in-plane magnetic fields (negligible gap) which is the hallmark of their topological origin. Further, we predict experimentally testable fingerprints of this behavior in the LDOS and in ballistic magnetotransport. While we focus on bismuthene/SiC, our main findings are also applicable to other honeycomb-based QSH systems.

Deep neural networks have evolved remarkably over the past few years and they are currently the fundamental tools of many intelligent systems. At the same time, the computational complexity and resource consumption of these networks also continue to increase. This will pose a significant challenge to the deployment of such networks, especially in real-time applications or on resource-limited devices. Thus, network acceleration has become a hot topic within the deep learning community. As for hardware implementation of deep neural networks, a batch of accelerators based on FPGA/ASIC have been proposed in recent years. In this paper, we provide a comprehensive survey of recent advances in network acceleration, compression and accelerator design from both algorithm and hardware points of view. Specifically, we provide a thorough analysis of each of the following topics: network pruning, low-rank approximation, network quantization, teacher-student networks, compact network design and hardware accelerators. Finally, we will introduce and discuss a few possible future directions.

In this paper we show that for an $\text{Sp}(k+1)$ invariant metric $\hat{g}$ on $\mathbb{S}^{4k+3}$ $(k\geq 1)$ close to the round metric, the conformally compact Einstein (CCE) manifold $(M, g)$ with $(\mathbb{S}^{4k+3}, [\hat{g}])$ as its conformal infinity is unique up to isometries. Moreover, by the result in [LiQingShi], $g$ is the Graham-Lee metric on the unit ball $B_1\subset \mathbb{R}^{4k+4}$. We also give an a priori estimate on the Einstein metric $g$. Based on the estimate and Graham-Lee and Lee's seminal perturbation result, we use the continuity method directly to obtain an existence result of the non-positively curved CCE metric with prescribed conformal infinity $(\mathbb{S}^{4k+3}, [\hat{g}])$ when the metric $\hat{g}$ is $\text{Sp}(k+1)$-invariant.

A complete understanding of the onset and subsequent evolution of confined flares has not been achieved. Earlier studies mainly analyzed disk events so as to reveal their magnetic topology and cause of confinement. In this study, taking advantage of a tandem of instruments working at different wavelengths of X-rays, EUVs, and microwaves, we present dynamic details of a confined flare observed on the northwestern limb of the solar disk on July 24th, 2016. The entire dynamic evolutionary process starting from its onset is consistent with a loop-loop interaction scenario. The X-ray profiles manifest an intriguing double-peak feature. From spectral fitting, it is found that the first peak is non-thermally dominated while the second peak is mostly multi-thermal with a hot (~10 MK) and a super-hot (~30 MK) component. This double-peak feature is unique in that the two peaks are clearly separated by 4 minutes, and the second peak reaches up to 25-50 keV; in addition, at energy bands above 3 keV the X-ray fluxes decline significantly between the two peaks. This, together with other available imaging and spectral data, manifest a two-stage energy release process. A comprehensive analysis is carried out to investigate the nature of this two-stage process. We conclude that the second stage with the hot and super-hot sources mainly involves direct heating through loop-loop reconnection at a relatively high altitude in the corona. The uniqueness of the event characteristics and complete data set make the study a nice addition to present literature on solar flares.

We measure complete and continuous Wigner functions of a two-level cesium atom in both a nearly pure state and highly mixed states. We apply the method [T. Tilma et al., Phys. Rev. Lett. 117, 180401 (2016)] of strictly constructing continuous Wigner functions for qubit or spin systems. We find that the Wigner function of all pure states of a qubit has negative regions and the negativity completely vanishes when the purity of an arbitrary mixed state is less than $\frac{2}{3}$. We experimentally demonstrate these findings using a single cesium atom confined in an optical dipole trap, which undergoes a nearly pure dephasing process. Our method can be applied straightforwardly to multi-atom systems for measuring the Wigner function of their collective spin state.

Recently, Belle Collaboration has reported the measurement of the spin-flipping transition $\Upsilon(4S) \to h_b(1P)\eta$ with an unexpectedly large branching ratio: $\mathcal{B}(\Upsilon(4S) \to h_b(1P)\eta) =(2.18\pm 0.11\pm 0.18)\times 10^{-3}$. Such a large branching fraction contradicts with the anticipated suppression for the spin flip. In this work, we examine the effects induced by intermediate bottomed meson loops and point out that these effects are significantly important. Using the effective Lagrangian approach (ELA), we find the experimental data on $\Upsilon(4S) \to h_b(1P)\eta$ can be accommodated with the reasonable inputs. We then explore the decays $\Upsilon(5S,6S)\to h_b(1P)\eta$ and find that these two channels also have sizable branching fractions. We also calculate these these processes in the framework of nonrelativistic effective field theory (NREFT). For the decays $\Upsilon(4S) \to h_b(1P) \eta$, the NREFT results are at the same order of magnitude but smaller than the ELA results by a factor of $2$ to $5$. For the decays $\Upsilon(5S, 6S) \to h_b(1P) \eta$ the NREFT results are smaller than the ELA results by approximately one order of magnitude. We suggest future experiment Belle-II to search for the $\Upsilon(5S, 6S)\to h_b(1P) \eta$ decays which will be helpful to understand the transition mechanism.

We investigate the decays of $\bar{B}^0_s$, $\bar{B}^0$ and $B^-$ into $\eta_c$ plus a scalar meson in a theoretical framework by taking into account the dominant process for the weak decay of $\bar{B}$ meson into $\eta_c$ and a $q\bar{q}$ pair. After hadronization of this $q\bar{q}$ component into pairs of pseudoscalar mesons we obtain certain weights for the pseudoscalar meson-pseudoscalar meson components. The calculation is based on the postulation that the scalar mesons $f_0(500)$, $f_0(980)$ and $a_0(980)$ are dynamically generated states from the pseudoscalar meson-pseudoscalar meson interactions in $S$-wave. Up to a global normalization factor, the $\pi \pi$, $K \bar{K}$ and $\pi \eta$ invariant mass distributions for the decays of $\bar B^0_s \to \eta_c \pi^+ \pi^-$, $\bar{B}^0_s \to \eta_c K^+ K^-$, $\bar B^0 \to \eta_c \pi^+ \pi^-$, $\bar{B}^0 \to \eta_c K^+ K^-$, $\bar{B}^0 \to \eta_c \pi^0 \eta$, $B^- \to \eta_c K^0 K^-$ and $B^- \to \eta_c \pi^- \eta$ are predicted. Comparison is made with the limited experimental information available and other theoretical calcualtions. Further comparison of these results with coming LHCb measurements will be very valuable to make progress in our understanding of the nature of the low lying scalar mesons, $f_0(500), f_0(980)$ and $a_0(980)$.

Fast simulation tools are highly appreciated in particle physics phenomenology studies, especially in the exploration of the physics potential of future experimental facilities. The Circular Electron Positron Collider is a proposed Higgs and Z factory that can precisely measure the Higgs boson properties and the electroweak precision observables. A fast-simulation toolkit dedicated to the CEPC detector has been developed using Delphes. The comparison shows that this fast simulation tool is highly consistent with the full simulation, on a set of benchmark distributions. Therefore, we recommend this fast simulation toolkit for CEPC phenomenological investigations.

A coupling of a dark photon $A'$ from a $U(1)_{A'}$ with the standard model (SM) particles can be generated through kinetic mixing represented by a parameter $\epsilon$. A non-zero $\epsilon$ also induces a mixing between $A'$ and $Z$ if dark photon mass $m_{A'}$ is not zero. This mixing can be large when $m_{A'}$ is close to $m_Z$ even if the parameter $\epsilon$ is small. Many efforts have been made to constrain the parameter $\epsilon$ for a low dark photon mass $m_{A'}$ compared with the $Z$ boson mass $m_Z$. We study the search for dark photon in $e^+e^- \to \gamma A' \to \gamma \mu^+ \mu^-$ for a dark photon mass $m_{A'}$ as large as kinematically allowed at future $e^+e^-$ colliders. For large $m_{A'}$, care should be taken to properly treat possible large mixing between $A'$ and $Z$. We obtain sensitivities to the parameter $\epsilon$ for a wide range of dark photon mass at planed $e^+\;e^-$ colliders, such as Circular Electron Positron Collider (CEPC), International Linear Collider (ILC) and Future Circular Collider (FCC-ee). For the dark photon mass $20~\text{GeV}\lesssim m_{A^{\prime}}\lesssim 330~\text{GeV}$, the $2\sigma$ exclusion limits on the mixing parameter are $\epsilon\lesssim 10^{-3}-10^{-2}$. The CEPC with $\sqrt{s}=240~\text{GeV}$ and FCC-ee with $\sqrt{s}=160~\text{GeV}$ are more sensitive than the constraint from current LHCb measurement once the dark photon mass $m_{A^{\prime}}\gtrsim 50~\text{GeV}$. For $m_{A^{\prime}}\gtrsim 220~\text{GeV}$, the sensitivity at the FCC-ee with $\sqrt{s}=350~\text{GeV}$ and $1.5~\text{ab}^{-1}$ is better than that at the 13~TeV LHC with $300~\text{fb}^{-1}$, while the sensitivity at the CEPC with $\sqrt{s}=240~\text{GeV}$ and $5~\text{ab}^{-1}$ can be even better than that at 13~TeV LHC with $3~\text{ab}^{-1}$ for $m_{A^{\prime}}\gtrsim 180~\text{GeV}$.

The variability of young stellar objects is mostly driven by star-disk interactions. In long-term photometric monitoring of the accreting T Tauri star GI Tau, we detect extinction events with typical depths of $\Delta V \sim 2.5$ mag that last for days-to-months and often appear to occur stochastically. In 2014 - 2015, extinctions that repeated with a quasi-period of 21 days over several months is the first empirical evidence of slow warps predicted from MHD simulations to form at a few stellar radii away from the central star. The reddening is consistent with $R_V=3.85\pm0.5$ and, along with an absence of diffuse interstellar bands, indicates that some dust processing has occurred in the disk. The 2015 -- 2016 multi-band lightcurve includes variations in spot coverage, extinction, and accretion, each of which results in different traces in color-magnitude diagrams. This lightcurve is initially dominated by a month-long extinction event and return to the unocculted brightness. The subsequent light-curve then features spot modulation with a 7.03 day period, punctuated by brief, randomly-spaced extinction events. The accretion rate measured from $U$-band photometry ranges from $1.3\times10^{-8}$ to $1.1\times10^{-10}$ M$_\odot$ yr$^{-1}$ (excluding the highest and lowest 5% of high- and low- accretion rate outliers), with an average of $4.7 \times 10^{-9}$ M$_\odot$ yr$^{-1}$. A total of 50% of the mass is accreted during bursts of $>12.8\times10^{-9}$ M$_\odot$ yr${^{-1}}$, which indicates limitations on analyses of disk evolution using single-epoch accretion rates.

In this paper, we present some generalized monogamy inequalities based on negativity and convex-roof extended negativity (CREN). These monogamy relations are satisfied by the negativity of $N$-qubit quantum systems $ABC_1\cdots C_{N-2}$, under the partitions $AB|C_1\cdots C_{N-2}$ and $ABC_1|C_2\cdots C_{N-2}$. Furthermore, the $W$-class states are used to test these generalized monogamy inequalities.

Finding the electromagnetic (EM) counterpart of binary compact star merger, especially the binary neutron star (BNS) merger, is critically important for gravitational wave (GW) astronomy, cosmology and fundamental physics. On Aug. 17, 2017, Advanced LIGO and \textit{Fermi}/GBM independently triggered the first BNS merger, GW170817, and its high energy EM counterpart, GRB 170817A, respectively, resulting in a global observation campaign covering gamma-ray, X-ray, UV, optical, IR, radio as well as neutrinos. The High Energy X-ray telescope (HE) onboard \textit{Insight}-HXMT (Hard X-ray Modulation Telescope) is the unique high-energy gamma-ray telescope that monitored the entire GW localization area and especially the optical counterpart (SSS17a/AT2017gfo) with very large collection area ($\sim$1000 cm$^2$) and microsecond time resolution in 0.2-5 MeV. In addition, \textit{Insight}-HXMT quickly implemented a Target of Opportunity (ToO) observation to scan the GW localization area for potential X-ray emission from the GW source. Although it did not detect any significant high energy (0.2-5 MeV) radiation from GW170817, its observation helped to confirm the unexpected weak and soft nature of GRB 170817A. Meanwhile, \textit{Insight}-HXMT/HE provides one of the most stringent constraints (~10$^{-7}$ to 10$^{-6}$ erg/cm$^2$/s) for both GRB170817A and any other possible precursor or extended emissions in 0.2-5 MeV, which help us to better understand the properties of EM radiation from this BNS merger. Therefore the observation of \textit{Insight}-HXMT constitutes an important chapter in the full context of multi-wavelength and multi-messenger observation of this historical GW event.

An asymptotically exact theory of spectrum and transverse distribution of magnetization in long-wave magnons is presented. It is based on exact analytical solution of linearized Landau-Lifshitz equation in a film. The quantization of the transverse wave parameters and wave vectors is studied.

With the development of Big Data and cloud data sharing, privacy preserving data publishing becomes one of the most important topics in the past decade. As one of the most influential privacy definitions, differential privacy provides a rigorous and provable privacy guarantee for data publishing. Differentially private interactive publishing achieves good performance in many applications; however, the curator has to release a large number of queries in a batch or a synthetic dataset in the Big Data era. To provide accurate non-interactive publishing results in the constraint of differential privacy, two challenges need to be tackled: one is how to decrease the correlation between large sets of queries, while the other is how to predict on fresh queries. Neither is easy to solve by the traditional differential privacy mechanism. This paper transfers the data publishing problem to a machine learning problem, in which queries are considered as training samples and a prediction model will be released rather than query results or synthetic datasets. When the model is published, it can be used to answer current submitted queries and predict results for fresh queries from the public. Compared with the traditional method, the proposed prediction model enhances the accuracy of query results for non-interactive publishing. Experimental results show that the proposed solution outperforms traditional differential privacy in terms of Mean Absolute Value on a large group of queries. This also suggests the learning model can successfully retain the utility of published queries while preserving privacy.

Deliberate control of magnon transportation will lead to an energy-efficient technology for information transmission and processing. Y3Fe5O12(YIG), exhibiting extremely large magnon diffusion length due to the low magnetic damping constant, has been intensively investigated for decades. While most of the previous works focused on the determination of magnon diffusion length by various techniques, herein we demonstrated how to tune magnon diffusion by light excitation. We found that the diffusion length of thermal magnons is strongly dependent on light wavelength when the magnon is generated by exposing YIG directly to laser beam. The diffusion length, determined by a nonlocal geometry at room temperature, is ~30 um for the magnons produced by visible light (400-650 nm), and ~136-156 um for the laser between 808 nm and 980 nm. The diffusion distance is much longer than the reported value. In addition to thermal gradient, we found that light illumination affected the electron configuration of the Fe3+ ion in YIG. Long wavelength laser triggers a high spin to low spin state transition of the Fe3+ ions in FeO6 octahedron. This in turn causes a substantial softening of the magnon thus a dramatic increase in diffusion distance. The present work paves the way towards an efficient tuning of magnon transport behavior which is crucially important for magnon spintronics.

Observations for the $\delta$ Scuti star AN Lyn have been made between 2008 and 2016 with the 85-cm telescope at Xinglong station of National Astronomical Observatories of China, the 84-cm telescope at SPM Observatory of Mexico and the Nanshan One meter Wide field Telescope of Xinjiang Observatory of China. Data in $V$ in 50 nights and in $R$ in 34 nights are obtained in total. The bi-site observations from both Xinglong Station and SPM Observatory in 2014 are analyzed with Fourier Decomposition to detect pulsation frequencies. Two independent frequencies are resolved, including one non-radial mode. A number of stellar model tracks are constructed with the MESA code and the fit of frequencies leads to the best-fit model with the stellar mass of $M = 1.70 \pm 0.05~\mathrm{M_{\odot}}$ , the metallicity abundance of $Z = 0.020 \pm 0.001$, the age of $1.33 \pm 0.01$ billion years and the period change rate $1/P\cdot \mathrm{d}P/\mathrm{d}t =1.06 \times 10^{-9} ~\mathrm{yr^{-1}} $, locating the star at the evolutionary stage close to the terminal age main sequence (TAMS). The O-C diagram provides the period change rate of $1/P \cdot \mathrm{d}P /\mathrm{d}t =4.5(8)\times10^{-7}~\mathrm{yr^{-1}}$. However, the period change rate calculated from the models is smaller in two orders than the one derived from the O-C diagram. Together with the sinusoidal function signature, the period variations are regarded to be dominated by the light-travel time effect of the orbital motion of a three-body system with two low-luminosity components, rather than the stellar evolutionary effect.

Quasiparticle excitations of free electrons in condensed-matter physics, characterized by the dimensionality of the band crossing, can find their elementary-particle analogs in high-energy physics, such as Majorana, Weyl, and Dirac fermions. While crystalline symmetry allows more quasiparticle excitations and exotic fermions to emerge. Using symmetry analysis and {\it ab-initio} calculations, we propose that the 3D honeycomb crystal NaCu$_3$Te$_2$ hosts triply degenerate nodal points (TDNPs) which are perfectly separated from the bulk states. We find a tunable phase transition between TDNPs and a weak TI triggered by a symmetry-allowed perturbation, and we further reveal the crucial role played by the spin-orbital coupling (SOC) for the emergence of the TDNPs in this system. Such topological non-trivial ternary compound not only serves as a perfect candidate for studying three-component fermions, but also provides a beautiful playground for understanding the topological phase transitions between TDNPs, TIs and trivial insulators, which distinguishes this system from other TDNP candidates.

{\it Victory}, i.e. \underline{vi}enna \underline{c}omputational \underline{to}ol deposito\underline{ry}, is a collection of numerical tools for solving the parquet equations for the Hubbard model and similar many body problems. The parquet formalism is a self-consistent theory at both the single- and two-particle levels, and can thus describe individual fermions as well as their collective behavior on equal footing. This is essential for the understanding of various emergent phases and their transitions in many-body systems, in particular for cases in which a single-particle description fails. Our implementation of {\it victory} is in modern Fortran and it fully respects the structure of various vertex functions in both momentum and Matsubara frequency space. We found the latter to be crucial for the convergence of the parquet equations, as well as for the correct determination of various physical observables. In this release, we thoroughly explain the program structure and the controlled approximations to efficiently solve the parquet equations, i.e. the two-level kernel approximation and the high-frequency regulation.

SrRuO$_3$ heterostructures grown in the (111) direction are a rare example of thin film ferromagnets. By means of density functional theory plus dynamical mean field theory we show that the half-metallic ferromagnetic state with an ordered magnetic moment of 2$\mu_{B}$/Ru survives the ultimate dimensional confinement down to a bilayer, even at elevated temperatures of 500$\,$K. In the minority channel, the spin-orbit coupling opens a gap at the linear band crossing corresponding to $\frac34$ filling of the $t_{2g}$ shell. We demonstrate that the respective state is Haldane's quantum anomalous Hall state with Chern number $C$=1, without an external magnetic field or magnetic impurities.