
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
hightemperature superconductivity. While a twodimensional 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 nonmagnetically ordered
Cu$_3$Zn(OD)$_6$FBr are similarly gapped although the content of
interkagomelayer 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 Diracnodalline state in a wellknown
conventional hightemperature superconductor, MgB2. Firstprinciples
calculations show that the Diracnodalline structure exhibits a unique
onedimensional dispersive Diracnodal line, protected by both
spatialinversion and timereversal symmetry, which connects the electron and
hole Dirac states. Most importantly, we show that the topological
superconducting phase can be realized with a conventional swave
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) longrange
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$_{4x}$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
hexagonaltoorthorhombic structural change with decreasing temperature in the
$x$ = 0 sample. While confirming the 3D antiferromagnetic nature of
lowtemperature 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
hightemperature 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$_{4x}$Zn$_x$(OH)$_6$FBr may be the longthought experimental playground for
the systematic investigations of the quantum phase transition from a longrange
antiferromagnet to a topologically ordered quantum spin liquid.

We describe electron spin resonance in a quantum spin liquid with significant
spinorbit 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 meanfield Hamiltonians, obtained with the help of
projective symmetry group analysis, which model a putative quantum spin liquid
state of the triangular rareearth 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
firstprinciples study, we provide a comprehensive theoretical characterization
of this insulating compound. The Fe$^{3+} $spin degree of freedom is quantified
as a quasi1D $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
orbitaldependent Hund's coupling between the spin and the hole is revealed.
Several important material properties are analyzed, including (a) factors
affecting the small $pd$ chargetransfer gap; (b) role of the extra interchain
Fe; and (c) the quasi1D 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 CurieWeiss 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 spin1/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
spin1/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 hightemperature 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 holelike 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
largesurface 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 smallpocket 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,3benzenedicarboxylate) [Cu(1,3bdc)] 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,3bdc) just a trivial
weak ferromagnet?". Using the the density functional theory, we have
systematically studied the electronic and magnetic properties of Cu(1,3bdc),
putting forth a theoretical basis to clarify this novel material. We present
numerical evidences of a dominating antiferromagnetic (AFM) exchange between
nearestneighbor (NN) Cu$^{2+}$ as experimentally extracted from the
hightemperature susceptibility data. We further show that beyond the NN AFM
exchange, the additional interactions in Cu(1,3bdc) have similar strength as
those in the wellstudied 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 spinwave continuum (SWC) observed in recent neutron scattering
experiments in Herbertsmithite, we use Gutzwillerprojected wave functions to
study dynamic spin structure factor $S(\mathbf{q},\omega)$ of spin liquid
states on the kagome lattice. Spin1 excited states in spin liquids are
represented by Gutzwillerprojected twospinon excited wave functions. We
investigate three different spin liquid candidates, spinon Fermisurface spin
liquid (FSL), Dirac spin liquid (DSL) and randomflux 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.
spin1/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
densitymatrixrenormalizationgroup 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
Schwingerboson parton wavefunctions, 12 more states can be realized in
Abrikosovfermion parton construction. Within a larger gauge group than SU(2),
the rest 44 states are also found in a spin3/2 system. Among 20 types of
Abrikosovfermion 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 wellknown 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 firstprinciples 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 antisite
disorder in Mg/Zn doped Barlowite is estimated to be one order of magnitude
lower than that in Herbertsmithite. The singleelectron 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
YangRiceZhang 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 holelike 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 electronlike
Fermi pocket ($\alpha$ orbit) and a new holelike 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 manybody 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
Gutzwillerprojected wave functions (GWF) carries a nontrivial 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 shortrange correlations, (2) a trivial topological order via
nondegenerate ground state, and zero topological entanglement entropy, (3) a
nontrivial $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 twodimensional bosonic lattice systems.

We implement the universal wave function overlap (UWFO) method to extract
modular $S$ and $T$ matrices for topological orders in Gutzwillerprojected
parton wave functions (GPWFs). The modular $S$ and $T$ matrices generate a
projective representation of $SL(2,\mathbb{Z})$ on the degenerategroundstate
Hilbert space on a torus and may fully characterize the 2+1D topological
orders, i.e. the quasiparticle 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 nonuniversal 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
Gutzwillerprojected 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 spin1/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 Gutzwillerprojected spin Hall
(GSH) states. We compute the twobody reduceddensitymatrices (2RDM) of GPWFs
for those target states, and compare them to the 2RDM of the DSL. This allows
us to design local spin models with only twobody 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 spin1 systems
on kagome lattice, and design local spin models that may realize CSL, VBS and
symmetryprotected 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 Odopants 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
Gutzwillerprojected spin Hall state. The inplane 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 longrange 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 hightemperature 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 firstprinciples design of a twodimensional (2D)
IndiumPhenylene Organometallic Framework (IPOF) that realizes a nearly FCB
right around the Fermi level by combining lattice geometry, spinorbit coupling
and ferromagnetism. An effective fourband model is constructed to reproduce
the firstprinciples 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 pseudospin1/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 caxis 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 bondoperator mean field calculation to describe
the magnetic properties of the bilayer twisted pseudospin1/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$_{2x}$La$_x$CuO$_{6+\delta}$
(LaBi2201) in strong magnetic fields, we identify the nonsuperconducting
pseudogap phase in the cuprates as a Luttingervolume 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=1p$. Both the phase string theory and the dopon
theory are introduced to describe the LvvFL. For the dopon theory, we can
obtain a semiquantitative 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 welldefined quasiparticles near the Fermi surface
and can be detected in the angleresolved photoemission spectroscopy. The DKSL
has the topological Kitaev spin liquid (KSL) as its neutral background. It
violates the WiedemannFranz 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, spinorbit 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 spinorbit coupling and ferromagnetism 
overall much higher than the energy scale of FQH states in semiconductors.