
With a generic lattice model for electrons occupying a semiinfinite crystal
with a hard surface, we study the eigenstates of the system with a bulk band
gap (or with nodal points in the gap). The exact solution to the wave functions
of scattering states is obtained. From the scattering states, we derive the
criterion for the existence of surface states. The wave functions and the
energy of the surface states are then determined. We obtain a connection
between the wave functions of the bulk states and the surface states. For
electrons in a system with timereversal symmetry, with this connection, we
rigorously prove the correspondence between the change of Kramers degeneracy of
the surface states and the bulk timereversal $Z_2$ invariant. The theory is
applicable to systems of (topological) insulators, superconductors, and
semimetals. Examples for solving the edge states of electrons with/without the
spinorbit interactions in graphene with a hard zigzag edge and that in a
twodimensional $d$wave superconductor with a (1,1) edge are given in
appendices.

In this paper, we investigate the topological properties of ZrRuP and
ScPd$_{3}$ using the firstprinciples calculations. We calculate the electronic
structure, surface states, and Z$_{2}$ topological invariant. The band
structure of ZrRuP without the spinorbit coupling (SOC) exhibits a nodal line
located at $\Gamma$KM high symmetry line in the Brillouin zone(BZ). The nodal
line is protected by the mirror symmetry corresponding to the plane $k_{z}=0$.
The geometry of Fermi surface is a connected torus centered at K point with
open orbits at the boundary of BZ, when the Fermi level is shifted up 0.25 eV.
The inclusion of SOC leads to the gap opening in the band structure which
corresponds adiabatically to transition from nodal line semimetal to
topological insulator. The corresponding surface states were obtained from thin
film calculations, and the result shows surface states emerging along the gap
region. We also identify the ScPd$_{3}$ compound as a Dirac nodal line
semimetal. This structure has nontrivial surface states and the corresponding
Z$_{2}$ topological invariant is (1;000). Both the ScPd$_{3}$ and ZrRuP have
the Fermi surfaces of the bands which make nodal lines with open orbits.

We present a formalism to calculate the orbital magnetization of interacting
Dirac fermions under a magnetic field. In this approach, the divergence
difficulty is overcome with a special limit of the derivative of the
thermodynamic potential with respect to the magnetic field. The formalism
satisfies the particlehole symmetry of the Dirac fermions system. We apply the
formalism to the interacting Dirac fermions in graphene. The charge and spin
orderings and the exchange interactions between all the Landau levels are taken
into account by the meanfield theory. The results for the orbital
magnetization of interacting Dirac fermions are compared with that of
noninteracting cases.

We study theoretically the topological surface states (TSSs) and the possible
surface Andreev bound states (SABSs) of Cu$_{x}$Bi$_{2}$Se$_{3}$ which is known
to be a topological insulator at $x=0$. The superconductivity (SC) pairing of
this compound is assumed to have the broken spinrotation symmetry, similar to
that of the Aphase of $^{3}$He as suggested by recent nuclearmagnetic
resonance experiments. For both spheroidal and corrugated cylindrical Fermi
surfaces with the hexagonal warping terms, we show that the bulk SC gap is
rather anisotropic; the minimum of the gap is negligibly small as comparing to
the maximum of the gap. This would make the fullygapped pairing effectively
nodal. For a clean system, our results indicate the bulk of this compound to be
a topological superconductor with the SABSs appearing inside the bulk SC gap.
The zeroenergy SABSs which are Majorana fermions, together with the TSSs not
gapped by the pairing, produce a zeroenergy peak in the surface density of
states (SDOS). The SABSs are expected to be stable against shortrange
nonmagnetic impurities, and the local SDOS is calculated around a nonmagnetic
impurity. The relevance of our results to experiments is discussed.

In this work, we predict a novel band structure for CarbonLithium(C4Li)
compound using the firstprinciples method. We show that it exhibits two Dirac
points near the Fermi level; one located at W point originating from the
nonsymmophic symmetry of the compound, and the other one behaves like a typeII
Dirac cone with higher anisotropy along the {\Gamma} to X line. The obtained
Fermi surface sheets of the holepocket and the electronpocket near the
typeII Dirac cone are separated from each other, and they would touch each
other when the Fermi level is doped to cross the typeII Dirac cone. The
evolution of Fermi surface with doping is also discussed. The bands crossing
from T to W make a linenode at the intersection of kx={\pi} and ky={\pi}
mirror planes. The C4Li is a novel material with both nonsymmorphic protected
Dirac cone and typeII Dirac cone near the Fermi level which may exhibit
exceptional topological property for electronic applications.

Two common approaches of studying theoretically the property of a
superconductor are shown to have significant differences, when they are applied
to the LarkinOvchinnikov state of Weyl metals. In the first approach the
pairing term is restricted by a cutoff energy to the neighborhood of the Fermi
surface, whereas in the second approach the pairing term is extended to the
whole Brillouin zone. We explore their difference by considering two minimal
models for the Weyl metal. For a model giving a single pair of Weyl pockets,
both two approaches give a partlygapped (fullygapped) bulk spectrum for small
(large) pairing amplitude. However, for very small cutoff energy, a portion of
the Fermi surface can be completely unaffected by the pairing term in the first
approach. For the other model giving two pairs of Weyl pockets, while the bulk
spectrum for the first approach can be fully gapped, the one from the second
approach has a robust line node, and the surface states are also changed
qualitatively by the pairing. We elucidate the above differences by topological
arguments and analytical analyses. A factor common to both of the two models is
the tilting of the Weyl cones which leads to asymmetric normal state band
structure with respect to the Weyl nodes. For the Weyl metal with two pairs of
Weyl pockets, the band folding leads to a double degeneracy in the effective
model, which distinguishes the pairing of the second approach from all others.

The local electronic states around a single As (Te, Se) vacancy are
investigated in order to shed light on the role of ligands in a series of
ironbased superconductors. Such a vacancy can produce a local hopping
correction ranging from $0.22$ eV to 0.12 eV and always induce two ingap
resonance peaks in the local density of states (LDOS) at the fixed symmetrical
bias voltages, which are rather robust and irrelevant to the phase of
superconducting order parameter. The LDOS images near the defect predominantly
possess $0^o$ and $45^o$ stripes. These energydependent charge modulations
created by quasiparticle interference are originated in the nesting effect
between the inner (outer) hole Fermi surface around $\Gamma$ point and the
inner (outer) electron Fermi surface around $M$ point.

Twodimensional Weyl superconductor is the most elusive member of a group of
materials with Weyl fermions as lowenergy excitations. Here, we propose to
realize this state in a heterostructure consisting of thin films of halfmetal
and spinsinglet superconductor. In particular, for the $d$wave case, a very
robust twodimensional Weyl superconductor (dWSC) is realized independent of
the orientation of the spontaneous magnetization of the halfmetal. The
quasiparticle spectra of the dWSC show interesting evolution with the direction
of the magnetization, featured by a series of Lifshitz transitions in the
zeroenergy contour of the quasiparticle spectrum. For a general magnetization
orientation of the halfmetal, the state is a combination of a superconducting
component and a normal fluid component and is different from all known forms of
pairings. In addition, we find a transition between typeI and typeII Weyl
nodes. This is also the first example of a typeII Weyl node in the presence of
superconducting correlation. The symmetries and topological properties of the
system are analyzed. We also study the phases in the heterostructure with the
halfmetal replaced by a ferromagnetic metal with a partially spinpolarized
Fermi surface.

By taking into account the charge and spin orderings and the exchange
interactions between all the Landau levels, we investigate the integer quantum
Hall effect of electrons in graphene using the meanfield theory. At the
fillings $\nu = 4n+2$ with $n = 0, 1, \cdots$, the system is in the
highsymmetry state with the Landau levels fourfold degenerated. We show that
with doping the degenerated lowest empty levels can be sequentially filled one
level by one level, the filled level is lower than the empty ones because of
the symmetry breaking. This result explains the step $\Delta\nu$ = 1 in the
integer quantized Hall conductivity of the experimental observations. We also
present in the supplemental material a high efficient method for dealing with
huge number of the Coulomb couplings between all the levels.

One promising avenue to study onedimensional ($1$D) topological phases is to
realize them in synthetic materials such as cold atomic gases. Intriguingly, it
is possible to realize Majorana boundary modes in a $1$D numberconserving
system consisting of two fermionic chains coupled only by pairhopping
processes. It is commonly believed that significant interchain singleparticle
tunneling necessarily destroys these Majorana modes, as it spoils the
$\mathbb{Z}_2$ fermion parity symmetry that protects them. In this Letter, we
present a new mechanism to overcome this obstacle, by piercing a (synthetic)
magnetic $\pi$flux through each plaquette of the Fermi ladder. Using
bosonization, we show that in this case there exists an exact leginterchange
symmetry that is robust to interchain hopping, and acts as fermion parity at
long wavelengths. We utilize density matrix renormalization group and exact
diagonalization to verify that the resulting model exhibits Majorana boundary
modes up to large singleparticle tunnelings, comparable to the intrachain
hopping strength. Our work highlights the unusual impacts of different
topologically trivial band structures on these interactiondriven topological
phases, and identifies a distinct route to stabilizing Majorana boundary modes
in $1$D fermionic ladders.

We find a series of topological phase transitions in a
halfmetal/superconductor heterostructure, by tuning the direction of the
magnetization of the halfmetal film. These include transitions between a
topological superconducting phase with a bulk gap and another phase without a
bulk gap but has a ubiquitous local gap. At the same time, the edge states
change from counterpropagating Majorana edge modes to unidirectional Majorana
edge modes. In addition, we find transitions between the second phase and a
nodal phase which turns out to be a twodimensional Weyl superconductor with
Fermi line edge states. We identify the topological invariants relevant to each
phase and the symmetry that protects the Weyl superconductivity phase.

A realization of Majorana fermions is proposed in the frustrated magnets via
the topological proximity effect. Specifically, we consider a theoretical
model, where a topological insulator is coupled to a frustrated magnetic
material through the spin exchange interaction. Using the renormalization
groupbased selfconsistent meanfield approach, and calculating the
selfenergy correction due to the topological insulator, we find that the spin
texture and the spinmomentum locking of the Dirac cone will be inherited by
the spinons in the nearby frustrated magnets. This leads to a particular
topological state of matter that supports the Majorana excitations. Unlike the
conventional realization in SC systems, these Majorana fermions are the
combination of spinons and antispinons, rather than electrons and holes. They
can participate in the transport of spinons, leading to nontrivial properties
of the spin transport.

Using a fourband Hamiltonian, we study the phase boundary of
spinpolarizedcurrent state (SPCS) of interacting electrons in bilayer
graphene. The model of spinpolarizedcurrent state has previously been shown
to resolve a number of experimental puzzles in bilayer graphene. The phase
boundaries of the SPCS with and without the external voltage between the two
layers are obtained in this work. An unusual phase boundary where there are two
transition temperatures for a given carrier concentration is found at finite
external voltage. The physics of this phenomenon is explained.

By means of firstprinciples calculations, we predict two new types of
partially hydrogenated graphene systems: C$_{6}$H$_{1}$ and C$_{6}$H$_{5}$,
which are shown to be ferromagnetic (FM) semimetal and FM narrowgap
semiconductor, respectively. When properly doped, the Fermi surface of the two
systems consists of an electron pocket or six hole patches in the first
Brillouin zone with completely spinpolarized charge carries. If
superconductivity exists in these systems, the stable pairing symmetries are
shown to be $p+ip$ for electron doped case, and anisotropic $p+ip$ for hole
doped case. The predicted systems may provide fascinating platforms for
studying the novel properties of ferromagnetism and tripletpairing
superconductivity as well as twodimensional spintronics.

Recent experiments on Znsubstituted 122type ironbased superconductors
(FeSCs) at electron and hole doped region provide us with a testing ground
for understanding the effect of Zn impurities in these systems. Our
firstprinciple calculations of the electronic structure reveal that the Zn 3d
orbitals are far below the Fermi level and chemically inactive, while the Zn
4sorbital is partially occupied and its wave function overlapping with those
3dorbitals of neighboring Feions. This suggests that the impurity effect is
originating in the Zn 4sorbital, not its 3dorbitals. Employing a
phenomenological twoorbital lattice model for 122FeSCs and the
selfconsistent Bogoliubovde Gennes equations, we study how the Znimpurities
suppress the superconductivity in electron and hole doped compounds. Our
obtained results qualitatively agree with the experimental measurements.

We theoretically study the quantum anomalies in the superconducting Weyl
metals based on the topological field theory. It is demonstrated that the Fermi
arc and the surface Andreev bound state, characteristic of the superconducting
Weyl metals, are the manifestations of two underlying phenomenon, namely the
chiral anomaly and the paritylike anomaly, respectively. The first anomaly is
inherited from the Berry curvature around the original Weyl points, while the
second is the result of the superconductivity. We show that, all the
fascinating topological behavior of the superconducting Weyl metals, either
intranode FFLO or the internode BCS pairing state, can be satisfactorily
described and predicted by our topological field theory.

We study the \emph{orbital effects} of the synthetic magnetic fields in an
interacting square lattice twoleg fermionic ladder model with the
numberconserving pair hopping that hosting the Majorana bound states. By
utilizing density matrix renormalization group and exact diagonalization, we
identify a novel typeII topological Meissner $($\emph{topo}Meissner$)$ phase
$($as distinguished from the lowfield typeI \emph{topo}Meissner state$)$
when threading a \emph{high} magnetic flux through the plaquette of the ladder,
which not only exhibits a large uniformly circulating chiral current along the
legs, the characteristics of the celebrated Meissner state, but also
accommodates a topologically protected groundstate manifold due to the
\emph{reentrant} emergence of the Majorana end modes. Our work reveals some
interesting interference effects resulting from the interplay between the gauge
fields and the strong interactions in establishing the intrinsic topological
states of matter in lowdimensional quantum systems.

We have investigated the complex twodimensional electron system with unequal
Rashba and Dresselhaus spinorbit interactions in the presence of a
perpendicular magnetic field. The spin polarizations are obtained in a wide
range of magnetic fields. It is shown that such a system is hard to be
magnetized. We also find that the resonant charge and spin Hall conductances
occurs simultaneously at a certain magnetic field, at which two (nearly)
degenerate Landau levels are filled partly. The resonant Hall effects are
universal in this type of semiconductor materials, and could have potential
application for semiconductor spintronics.

Recently topological superconducting states has attracted a lot of interest.
In this work, we consider a topo logical superconductor with $Z_2$ topological
mirror order [1] and s$\pm$wave superconducting pairing symmetry, within a
twoorbital model originally designed for ironbased superconductivity [2]. We
predict the existence of gapless edge states. We also study the local
electronic structure around an adsorbed interstitial magnetic impurity in the
system, and find the existence of lowenergy ingap bound states even with a
weak spin polar ization on the impurity. We also discuss the relevance of our
results to the recent STM experiment on Fe(Te,Se) compound with adsorbed Fe
impurity [3], for which our density functional calculations show the Fe
impurity is spin polarized.

By taking into account the possibility of all the intralayer as well as the
interlayer current orderings, we derive an eightband model for interacting
electrons in bilayer graphene. With the numerical solution to the model, we
show that only the current orderings between the same sublattice sites can
exist within the range of the physical interacting strength. This result
confirms our previous model of spinpolarizedcurrent phase for the
groundstate of interacting electrons in bilayer graphene that resolves a
number of experimental puzzles.

Recent experiments demonstrated that isovalent doping system gives the
similar phase diagram as the heterovalent doped cases. For example, with the
phosphorous (P)doping, the magnetic order in
BaFe$_{2}$(As$_{1x}$P$_{x}$)$_{2}$ compound is first suppressed, then the
superconductivity dome emerges to an extended doping region but eventually it
disappears at large $x$. With the help of a minimal twoorbital model for both
BaFe$_{2}$As$_{2}$ and BaFe$_{2}$P$_{2}$, together with the selfconsistent
lattice Bogoliubovde Gennes (BdG) equation, we calculate the phase diagram
against the P content $x$ in which the doped isovalent Patoms are treated as
impurities. We show that our numerical results can qualitatively compare with
the experimental measurements.

The pairing symmetry is examined in highly electrondoped
Ba(Fe$_{1x}$Co$_x$As)$_2$ and A$_y$Fe$_2$Se$_2$ (with A=K, Cs) compounds, with
similar crystallographic and electronic band structures. Starting from a
phenomenological twoorbital model, we consider nearestneighbor and
nextnearestneighbor intraorbital pairing interactions on the Fe square
lattice. In this model, we find a unified description of the evolution from
$s_\pm$wave pairing ($2.0 < n \lesssim 2.4$) to $d$wave pairing ($2.4
\lesssim n \lesssim 2.5$) as a function of electron filling. In the crossover
region a novel timereversal symmetry breaking state with $s_\pm+id$ pairing
symmetry emerges. This minimal model offers an overall picture of the evolution
of superconductivity with electron doping for both $s_\pm$wave
[Ba(Fe$_{1x}$Co$_x$As)$_2$] and $d$wave [A$_y$Fe$_2$Se$_2$] pairing, as long
as the dopants only play the role of a charge reservoir. However, the situation
is more complicated for Ba(Fe$_{1x}$Co$_x$As)$_2$. A realspace study further
shows that when the impurity scattering effects of Co dopants are taken into
account, the superconductivity is completely suppressed for $n > 2.4$. This
preempts any observation of $d$wave pairing in this compound, in contrast to
A$_y$Fe$_2$Se$_2$.

The magnetic order of the localized Eu$^{2+}$ spins in optimallydoped
Eu(Fe$_{1x}$Ir$_{x}$)$_{2}$As$_{2}$ ($\mathit{x}$ = 0.12) with superconducting
transition temperature $\mathit{T_{SC}}$ = 22 K was investigated by
singlecrystal neutron diffraction. The Eu$^{2+}$ moments were found to be
ferromagnetically aligned along the $\mathit{c}$direction with an ordered
moment of 7.0(1) $\mu_{B}$ well below the magnetic phase transition temperature
$\mathit{T_{C}}$ = 17 K. No evidence of the tetragonaltoorthorhombic
structural phase transition was found in this compound within the experimental
uncertainty, in which the spindensitywave (SDW) order of the Fe sublattice is
supposed to be completely suppressed and the superconductivity gets fully
developed. The ferromagnetic groud state of the Eu$^{2+}$ spins in
Eu(Fe$_{0.88}$Ir$_{0.12}$)$_{2}$As$_{2}$ was supported by the firstprinciples
density functional calculation. In addition, comparison of the electronic
structure calculations between Eu(Fe$_{0.875}$Ir$_{0.125}$)$_{2}$As$_{2}$ and
the parent compound EuFe$_{2}$As$_{2}$ indicates stronger hybridization and
more expanded bandwith due to the Ir substitution, which together with the
introduction of electrons might work against the FeSDW in favor of the
superconductivity.

The electronic band structure of iron pnictides exhibits four Dirac cones,
which are due to crystal symmetry and orbital bonding orientation. This
hallmark signature presents the pnictide family as an ideal candidate in the
search for quasitwodimensional topological crystalline insulators. In this
report, we explore interactioninduced topological phases which cannot be
described by conventional local order parameters. Based on a model Hamiltonian
our symmetry analysis shows that sponta neous novel topological phases may be
realized in compounds with tetragonal crystal field symmetry, where the
electrons occupy the two degenerate t2g energy levels at low temperature. We
identify two stable topological phases in the ground state, which emerge from
spontaneous orbital current order. These currents are driven by electronic
correlations caused by interorbital Coulomb interactions. The first
topological phase is an anomalous orbital Hall phase, characterized by a
nonzero Chern number, while the second topological phase has a vanishing Chern
number, though with an extra Z2like invariant that preserves parity. More
specifically, the interactioninduced novel phase of the quasitwodimensional
topological crystalline insulator is protected by mirror reflection symmetries
and therefore may be realized in pnictides.

Topological phases originating from spinorbit coupling have attracted great
attention recently. In this work, we use cellular dynamical mean field theory
with the continuoustime quantum Monte Carlo solver to study the
KaneMeleHubbard model supplemented with an additional thirdneighbor hopping
term. For weak interactions, the thirdneighbor hopping term drives a
topological phase transition between a topological insulator and a trivial
insulator, consistent with previous fermion signfree quantum Monte Carlo
results [H.Hung et al., Phys. Rev. B 89, 235104 (2014)]. At finite
temperatures, the Dirac cones of the zero temperature topological phase
boundary give rise to a metallic regime of finite width in the thirdneighbor
hopping. Furthermore, we extend the range of interactions into the strong
coupling regime and find an easyplane antiferromagnetic insulating state
across a wide range of thirdneighbor hopping. In contrast to the weak coupling
regime, no topological phase transition occurs at strong coupling, and the
ground state is a trivial antiferromagnetic insulating state. A comprehensive
finite temperature phase diagram in the interactionthirdneighbor hopping
plane is provided.