• A quantum spin liquid with a $Z_2$ topological order has long been thought to be important for the application of quantum computing and may be related to high-temperature superconductivity. While a two-dimensional kagome antiferromagnet may host such a state, strong experimental evidences are still lacking. Here we show that the spin excitations from the kagome planes in magnetically ordered Cu$_4$(OD)$_6$FBr and non-magnetically ordered Cu$_3$Zn(OD)$_6$FBr are similarly gapped although the content of inter-kagome-layer Cu$^{2+}$ ions changes dramatically. This suggests that the spin triplet gap and continuum of the intrinsic kagome antiferromagnet are robust against the interlayer magnetic impurities. Our results show that the ground state of Cu$_3$Zn(OD)$_6$FBr is a gapped quantum spin liquid with Z$_2$ topological order.
  • Topological superconductors are an intriguing and elusive quantum phase, characterized by topologically protected gapless surface/edge states residing in a bulk superconducting gap, which hosts Majorana fermions. Unfortunately, all currently known topological superconductors have a very low transition temperature, limiting experimental measurements of Majorana fermions. Here we discover the existence of a topological Dirac-nodal-line state in a well-known conventional high-temperature superconductor, MgB2. First-principles calculations show that the Dirac-nodal-line structure exhibits a unique one-dimensional dispersive Dirac-nodal line, protected by both spatial-inversion and time-reversal symmetry, which connects the electron and hole Dirac states. Most importantly, we show that the topological superconducting phase can be realized with a conventional s-wave superconducting gap, evidenced by the topological edge mode of the MgB2 thin films showing chiral edge states. Our discovery may enable the experimental measurement of Majorana fermions at high temperature.
  • Barlowite Cu$_4$(OH)$_6$FBr shows three dimensional (3D) long-range antiferromagnetism, which is fully suppressed in Cu$_3$Zn(OH)$_6$FBr with a kagome quantum spin liquid ground state. Here we report systematic studies on the evolution of magnetism in the Cu$_{4-x}$Zn$_x$(OH)$_{6}$FBr system, as a function of $x$, to bridge the two limits of Cu$_4$(OH)$_6$FBr ($x=0$) and Cu$_3$Zn(OH)$_6$FBr ($x=1$). Neutron diffraction measurements reveal a hexagonal-to-orthorhombic structural change with decreasing temperature in the $x$ = 0 sample. While confirming the 3D antiferromagnetic nature of low-temperature magnetism, the magnetic moments on some Cu$^{2+}$ sites in the kagome planes are found to be vanishingly small, suggesting strong frustration already exist in barlowite. Substitution of interlayer Cu$^{2+}$ with Zn$^{2+}$, with gradually increasing $x$, completely suppresses the bulk magnetic order at around $x$ = 0.4, but leaves a local, secondary magnetic order up to $x\sim 0.8$ with slight decrease of its transition temperature. The high-temperature magnetic susceptibility and specific heat measurements further suggest that the intrinsic magnetic properties of kagome spin liquid planes may already appear from $x>0.3$ samples. Our results therefore reveal that the Cu$_{4-x}$Zn$_x$(OH)$_6$FBr may be the long-thought experimental playground for the systematic investigations of the quantum phase transition from a long-range antiferromagnet to a topologically ordered quantum spin liquid.
  • We describe electron spin resonance in a quantum spin liquid with significant spin-orbit coupling. We find that the resonance directly probes spinon continuum which makes it an efficient and informative probe of exotic excitations of the spin liquid. Specifically, we consider spinon resonance of three different spinon mean-field Hamiltonians, obtained with the help of projective symmetry group analysis, which model a putative quantum spin liquid state of the triangular rare-earth antiferromagnet YbMgGaO4. The band of absorption is found to be very broad and exhibit strong van Hove singularities of single spinon spectrum as well as pronounced polarization dependence.
  • NaFe$_{0.5}$Cu$_{0.5}$As represents a rare exception in the metallic iron pnictide family, in which a small insulating gap is opened. Based on first-principles study, we provide a comprehensive theoretical characterization of this insulating compound. The Fe$^{3+} $spin degree of freedom is quantified as a quasi-1D $S=\frac{5}{2}$ Heisenberg model. The itinerant As hole state is downfolded to a $p_{xy}$-orbital hopping model on a square lattice. A unique orbital-dependent Hund's coupling between the spin and the hole is revealed. Several important material properties are analyzed, including (a) factors affecting the small $p-d$ charge-transfer gap; (b) role of the extra interchain Fe; and (c) the quasi-1D spin excitation in the Fe chains. The experimental manifestations of these properties are discussed.
  • We report a new kagome quantum spin liquid candidate Cu$_3$Zn(OH)$_6$FBr, which does not experience any phase transition down to 50 mK, more than three orders lower than the antiferromagnetic Curie-Weiss temperature ($\sim$ 200 K). A clear gap opening at low temperature is observed in the uniform spin susceptibility obtained from $^{19}$F nuclear magnetic resonance measurements. We observe the characteristic magnetic field dependence of the gap as expected for fractionalized spin-1/2 spinon excitations. Our experimental results provide firm evidence for spin fractionalization in a topologically ordered spin system, resembling charge fractionalization in the fractional quantum Hall state.
  • We apply symmetric tensor network state (TNS) to study the nearest neighbor spin-1/2 antiferromagnetic Heisenberg model on Kagome lattice. Our method keeps track of the global and gauge symmetries in TNS update procedure and in tensor renormalization group (TRG) calculation. We also introduce a very sensitive probe for the gap of the ground state -- the modular matrices, which can also determine the topological order if the ground state is gapped. We find that the ground state of Heisenberg model on Kagome lattice is a gapped spin liquid with the $\mathbb{Z}_2$-topological order (or toric code type), which has a long correlation length $\xi\sim 10$ unit cell length. We justify that the TRG method can handle very large systems with over thousands of spins. Such a long $\xi$ explains the gapless behaviors observed in simulations on smaller systems with less than 300 spins or shorter than 10 unit cell length. We also discuss experimental implications of the topological excitations encoded in our symmetric tensors.
  • The discovery of high-temperature superconductivity in FeSe/STO has trigged great research interest to reveal a range of exotic physical phenomena in this novel material. Here we present a temperature dependent magnetotransport measurement for ultrathin FeSe/STO films with different thickness and protection layers. Remarkably, a surprising linear magnetoresistance (LMR) is observed around the superconducting transition temperatures but absent otherwise. The experimental LMR can be reproduced by magnetotransport calculations based on a model of magnetic field dependent disorder induced by spin fluctuation. Thus, the observed LMR in coexistence with superconductivity provides the first magnetotransport signature for spin fluctuation around the superconducting transition region in ultrathin FeSe/STO films.
  • The Fermi surface in underdoped cuprates is reconstructed by the charge density wave (CDW) order in the pseudogap phase. Theoretical proposals can be divided into two classes: one assumes the underlying Fermi surface without CDW as a conventional large surface; the other assumes small hole-like Fermi pockets. In both scenarios, we theoretically study the quantum oscillation and find three evenly spaced peaks in the oscillation spectra. The central dominant peak is induced by the CDW order. Its effective mass is strongly enhanced as the CDW vanishes in agreement with experiments. But the two scenarios have different understandings of the subdominant satellite peaks. In the large-surface scenario they are induced by the interlayer tunneling between the bilayer CuO$_{2}$ planes. Their effective masses are also enhanced with descreasing CDW. In the small-pocket scenario one of the subdominant peaks comes from the original small Fermi pockets of the pseudogap state. Its effective mass is nearly independent of the CDW strength and increases monotonically with the doping. We propose future quantum oscillation experiments to test these different predictions and thus to clarify the underlying Fermi surface structure of the pseudogap state.
  • Cu(1,3-benzenedicarboxylate) [Cu(1,3-bdc)] contains structurally perfect kagom\'{e} planes formed by Cu$^{2+}$ ions without the presence of diamagnetic defects. This organometallic compound should have served as a precious platform to explore quantum frustrated magnetism, yet the experimental results so far are mysterious, leading to questions such as "Is Cu(1,3-bdc) just a trivial weak ferromagnet?". Using the the density functional theory, we have systematically studied the electronic and magnetic properties of Cu(1,3-bdc), putting forth a theoretical basis to clarify this novel material. We present numerical evidences of a dominating antiferromagnetic (AFM) exchange between nearest-neighbor (NN) Cu$^{2+}$ as experimentally extracted from the high-temperature susceptibility data. We further show that beyond the NN AFM exchange, the additional interactions in Cu(1,3-bdc) have similar strength as those in the well-studied kagom\'{e} antiferromagnet, Herbertsmithite, by designing a comparative study. In the end, the origin of the phase transition and FM signals observed under low temperature is discussed based on the DFT results.
  • Motivated by spin-wave continuum (SWC) observed in recent neutron scattering experiments in Herbertsmithite, we use Gutzwiller-projected wave functions to study dynamic spin structure factor $S(\mathbf{q},\omega)$ of spin liquid states on the kagome lattice. Spin-1 excited states in spin liquids are represented by Gutzwiller-projected two-spinon excited wave functions. We investigate three different spin liquid candidates, spinon Fermi-surface spin liquid (FSL), Dirac spin liquid (DSL) and random-flux spin liquid (RSL). FSL and RSL have low energy peaks in $S(\mathbf{q},\omega)$ at $K$ points in the extended magnetic Brillouin zone, in contrast to experiments where low energy peaks are found at $M$ points. There is no obviuos contradiction between DSL and neutron scattering measurements. Besides a fractionalized spin (i.e. spin-1/2), spinons in DSL carry a fractionalized crystal momentum which is potentially detectable in SWC in the neutron scattering measurements.
  • We study different ways of symmetry fractionalization in $Z_{2}$ spin liquids on the triangular lattice. Our classification can be used to identify the symmetry fractionalization in the $Z_{2}$ spin liquid reported in recent density-matrix-renormalization-group simulations for $J_{1}-J_{2}$ spin model on the triangular lattice. We find 64 types of symmetry enriched $Z_{2}$ spin liquid states on triangular lattice. Besides 8 states constructed in Schwinger-boson parton wavefunctions, 12 more states can be realized in Abrikosov-fermion parton construction. Within a larger gauge group than SU(2), the rest 44 states are also found in a spin-3/2 system. Among 20 types of Abrikosov-fermion Gutzwillerprojected wavefunctions, No.B5 state is a promising candidate for the $Z_{2}$ spin liquid for $J_{1}-J_{2}$ spin model on the triangular lattice. No.B5 lies close to Dirac spin liquid (DSL). However, variational Monte Carlo simulation find that DSL has a good variational energy and $J_{1}-J_{2}$ spin model cannot open a gap for spinons on top of DSL to stabilize No.B5 state.
  • Barlowite $Cu_4(OH)_6FBr$ is a newly found mineral containing $Cu^{2+}$ kagome planes. Despite similarities in many aspects to Herbertsmithite $Cu_3Zn(OH)_6Cl_2$, the well-known quantum spin liquid (QSL) candidate, intrinsic Barlowite turns out not to be a QSL, possibly due to the presence of $Cu^{2+}$ ions in between kagome planes that induce interkagome magnetic interaction [PRL, 113, 227203 (2014)]. Using first-principles calculation, we systematically study the feasibility of selective substitution of the interkagome Cu ions with isovalent nonmagnetic ions. Unlike previous speculation of using larger dopants, such as $Cd^{2+}$ and $Ca^{2+}$, we identify the most ideal stoichiometric doping elements to be Mg and Zn in forming $Cu_3Mg(OH)_6FBr$ and $Cu_3Zn(OH)_6FBr$ with the highest site selectivity and smallest lattice distortion. The equilibirium anti-site disorder in Mg/Zn- doped Barlowite is estimated to be one order of magnitude lower than that in Herbertsmithite. The single-electron band structure and orbital component analysis show that the proposed selective doping effectively mitigates the difference between Barlowite and Herbertsmithite.
  • Recent experiments have revealed incommensurate charge density wave (CDW) in the pseudogap regime in underdoped cuprates, e.g. YBa$_2$Cu$_3$O$_{6+\delta}$ and HgBa$_2$CuO$_{4+\delta}$. However, its relationship with the pseudogap is still controversial. In this work, we take a phenomenological synthesis of the doped resonating valence bond (RVB) state and the CDW order. Starting from the Yang-Rice-Zhang Green's function ansatz for the doped RVB state [Phys. Rev. B {\bf 73}, 174501 (2006)], in which the Fermi surface is partially truncated into four nodal hole-like Fermi pockets by the antinodal RVB gap, we show that the CDW order at the wavevectors connecting the tips of the Fermi arcs (the hotspots) induces Fermi surface reconstruction, giving rise to an electron-like Fermi pocket ($\alpha$ orbit) and a new hole-like Fermi pocket ($\beta$ orbit). The $\alpha$ orbit is formed by joining the Fermi arcs at the hotspots and it dominates the quantum oscillation Fourier spectrum, while the $\beta$ orbit is formed by joining the outer patches of the original hole pockets, which has vanishingly small spectral weight. The areas enclosed by these orbits are extracted from the density of states oscillation in magnetic field and quantitatively agree with the experiments.
  • Gutzwiller projection is a way to construct many-body wave functions that could carry topological order or symmetry protected topological (SPT) order. However, an important issue is to determine whether or not a given Gutzwiller-projected wave functions (GWF) carries a non-trivial SPT order, and which SPT order is carried by the wavefunction. In this paper, we numerically study the SPT order in a spin $S = 1$ GWF on the Kagome lattice. Using the standard Monte Carlo method, we directly confirm that the GWF has (1) gapped bulk with short-range correlations, (2) a trivial topological order via nondegenerate ground state, and zero topological entanglement entropy, (3) a non-trivial $U(1)\times U(1)$ SPT order via the Hall conductances of the protecting $U(1)\times U(1)$ symmetry, and (4) symmetry protected gapless boundary. This represents numerical evidence of continuous symmetry protected topological order in two-dimensional bosonic lattice systems.
  • We implement the universal wave function overlap (UWFO) method to extract modular $S$ and $T$ matrices for topological orders in Gutzwiller-projected parton wave functions (GPWFs). The modular $S$ and $T$ matrices generate a projective representation of $SL(2,\mathbb{Z})$ on the degenerate-ground-state Hilbert space on a torus and may fully characterize the 2+1D topological orders, i.e. the quasi-particle statistics and chiral central charge (up to $E_8$ bosonic quantum Hall states). We used the variational Monte Carlo method to computed the $S$ and $T$ matrices of the chiral spin liquid (CSL) constructed by the GPWF on the square lattice, and confirm that the CSL carries the same topological order as the $\nu=\frac{1}{2}$ bosonic Laughlin state. We find that the non-universal exponents in UWFO can be small and direct numerical computation is able to be applied on relatively large systems. We also discuss the UWFO method for GPWFs on other Bravais lattices in two and three dimensions by using the Monte Carlo method. UWFO may be a powerful method to calculate the topological order in GPWFs.
  • We introduce a method to design a local spin Hamiltonian to realize a Gutzwiller-projected parton wave functions (GPWF) as its ground state. For example, the Dirac spin liquid (DSL) state is quite close to the true ground state of the spin-1/2 Heisenberg model on kagome lattice. We examine what kind of perturbations we should add in order to drive the DSL to more stable chiral spin liquid (CSL), valence bond solid (VBS) or Gutzwiller-projected spin Hall (GSH) states. We compute the two-body reduced-density-matrices (2-RDM) of GPWFs for those target states, and compare them to the 2-RDM of the DSL. This allows us to design local spin models with only two-body interactions that may realize those interesting target states. Our results agree very well with recent numerical calculations for CSL on kagome lattice. We also study spin-1 systems on kagome lattice, and design local spin models that may realize CSL, VBS and symmetry-protected topological (SPT) states. Our work establishes a directional guide for further numerical simulations.
  • Generally strong charge and magnetic inhomogeneities are observed in NQR/NMR experiments on underdoped cuprates. It is not the case for the underdoped HgBa$_2$CuO$_{4+\delta}$, the most symmetric and highest T$_c$ single layer cuprate, whose magnetic inhomogeneity is strongly suppressed. Also neutron scattering experiments reveal a unique pair of weakly dispersive magnetic modes in this material. We propose that these special properties stem from the symmetric positioning of the O-dopants between adjacent CuO$_2$ layers that lead to a strong superexchange interaction between a pair of hole spins. In this Letter we present a theoretical model, which gives a consistent explanation to the anomalous magnetic properties of this material.
  • A planar Mott insulator with easy plane N\'eel order can be mapped onto a Gutzwiller-projected spin Hall state. The in-plane substitution of Zn$^{2+}$ and Li$^{+}$ for Cu$^{2+}$ brings zero modes around the spin vacancy. The Gutzwiller projection precludes single particle excitations, however, preserves the zero modes which have a local spin moment for Zn$^{2+}$ and a local charged hole for Li$^{+}$, respectively. While the local spin moment for Zn$^{2+}$ is screened by the long-range spin correlations, the active charge degree of freedom for Li$^{+}$ impurity twists the spin background with long ranged disturbances. This proposal explains the very different roles of the Zn$^{2+}$ and Li$^{+}$ impurities on the magnetic perturbations in the CuO$_2$ plane in the undoped cuprate compounds.
  • By combining exotic band dispersion with nontrivial band topology, an interesting type of band structure, namely the flat chern band (FCB), has recently been proposed to spawn high-temperature fractional quantum hall states. Despite the proposal of several theoretical lattice models, however, it remains a doubt whether such a "romance of flatland" could exist in a real material. Here, we present a first-principles design of a two-dimensional (2D) Indium-Phenylene Organometallic Framework (IPOF) that realizes a nearly FCB right around the Fermi level by combining lattice geometry, spin-orbit coupling and ferromagnetism. An effective four-band model is constructed to reproduce the first-principles results. Our design in addition provides a general strategy to synthesize topologically nontrivial materials in virtue of organic chemistry and nanotechnology.
  • We use the bilayer twisted Hubbard model to study the low energy electronic structure of Sr3Ir2O7. Sr3Ir2O7 is suggested as a Mott insulator and described by the bilayer twisted pseudospin-1/2 model. The bilayer twists in Sr3Ir2O7 bring the intrinsic anisotropy and the system is no longer SU(2) invariant. The anisotropy selects an easy c-axis collinear antiferromagnetic magnetic structure at low temperatures which breaks the discrete Ising symmetry. The transverse spin excitations are gapped and the longitudinal one is weakly dispersive. We implement the bond-operator mean field calculation to describe the magnetic properties of the bilayer twisted pseudospin-1/2 model in good agreement with experimental measurements in Sr3Ir2O7. The bilayer twisted Hubbard model is believed to be the good zeroth approximation for the low energy electronic structure of Sr3Ir2O7.
  • We have realized that our results have a strong overlap with those previously published in J Chovan, N Papanicolaou, and S Komineas, Phys.Rev.B65,064433(2002).
  • Based on the NMR measurements on Bi$_2$Sr$_{2-x}$La$_x$CuO$_{6+\delta}$ (La-Bi2201) in strong magnetic fields, we identify the non-superconducting pseudogap phase in the cuprates as a Luttinger-volume violating Fermi liquid (LvvFL). This state is a zero temperature quantum liquid that does not break translational symmetry, and yet, the Fermi surface encloses a volume smaller than the large one given by the Luttinger theorem. The particle number enclosed by the small Fermi surface in the LvvFL equals the doping level $p$, not the total electron number $n_e=1-p$. Both the phase string theory and the dopon theory are introduced to describe the LvvFL. For the dopon theory, we can obtain a semi-quantitative agreement with the NMR experiments.
  • In this paper, we study the lightly doped Kitaev spin liquid (DKSL) and find it to be a Fermi liquid state. The DKSL has a closed Fermi surface around the center of the Brillouin zone with the quantized volume determined by the Luttinger's theorem. It has well-defined quasiparticles near the Fermi surface and can be detected in the angle-resolved photoemission spectroscopy. The DKSL has the topological Kitaev spin liquid (KSL) as its neutral background. It violates the Wiedemann-Franz law and has a large Wilson ration. These results has the potential experimental test in the iridates upon doping.
  • We study quantum spin systems with a proper combination of geometric frustration, spin-orbit coupling and ferromagnetism. We argue that such a system is likely to be in a chiral spin state, a fractional quantum Hall (FQH) state for bosonic spin degrees of freedom. The energy scale of the bosonic FQH state is of the same order as the spin-orbit coupling and ferromagnetism --- overall much higher than the energy scale of FQH states in semiconductors.