• With a generic lattice model for electrons occupying a semi-infinite 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 time-reversal symmetry, with this connection, we rigorously prove the correspondence between the change of Kramers degeneracy of the surface states and the bulk time-reversal $Z_2$ invariant. The theory is applicable to systems of (topological) insulators, superconductors, and semi-metals. Examples for solving the edge states of electrons with/without the spin-orbit interactions in graphene with a hard zigzag edge and that in a two-dimensional $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 first-principles calculations. We calculate the electronic structure, surface states, and Z$_{2}$ topological invariant. The band structure of ZrRuP without the spin-orbit coupling (SOC) exhibits a nodal line located at $\Gamma$-K-M 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 particle-hole 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 mean-field theory. The results for the orbital magnetization of interacting Dirac fermions are compared with that of non-interacting 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 spin-rotation symmetry, similar to that of the A-phase of $^{3}$He as suggested by recent nuclear-magnetic 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 fully-gapped 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 zero-energy SABSs which are Majorana fermions, together with the TSSs not gapped by the pairing, produce a zero-energy peak in the surface density of states (SDOS). The SABSs are expected to be stable against short-range 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 Carbon-Lithium(C4Li) compound using the first-principles 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 type-II Dirac cone with higher anisotropy along the {\Gamma} to X line. The obtained Fermi surface sheets of the hole-pocket and the electron-pocket near the type-II Dirac cone are separated from each other, and they would touch each other when the Fermi level is doped to cross the type-II Dirac cone. The evolution of Fermi surface with doping is also discussed. The bands crossing from T to W make a line-node at the intersection of kx={\pi} and ky={\pi} mirror planes. The C4Li is a novel material with both nonsymmorphic protected Dirac cone and type-II 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 Larkin-Ovchinnikov 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 partly-gapped (fully-gapped) 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 iron-based superconductors. Such a vacancy can produce a local hopping correction ranging from $-0.22$ eV to 0.12 eV and always induce two in-gap 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 energy-dependent 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.
  • Two-dimensional Weyl superconductor is the most elusive member of a group of materials with Weyl fermions as low-energy excitations. Here, we propose to realize this state in a heterostructure consisting of thin films of half-metal and spin-singlet superconductor. In particular, for the $d$-wave case, a very robust two-dimensional Weyl superconductor (dWSC) is realized independent of the orientation of the spontaneous magnetization of the half-metal. The quasiparticle spectra of the dWSC show interesting evolution with the direction of the magnetization, featured by a series of Lifshitz transitions in the zero-energy contour of the quasiparticle spectrum. For a general magnetization orientation of the half-metal, 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 type-I and type-II Weyl nodes. This is also the first example of a type-II 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 half-metal replaced by a ferromagnetic metal with a partially spin-polarized 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 mean-field theory. At the fillings $\nu = 4n+2$ with $n = 0, 1, \cdots$, the system is in the high-symmetry state with the Landau levels four-fold 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 one-dimensional ($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 number-conserving system consisting of two fermionic chains coupled only by pair-hopping processes. It is commonly believed that significant interchain single-particle 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 leg-interchange 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 single-particle tunnelings, comparable to the intrachain hopping strength. Our work highlights the unusual impacts of different topologically trivial band structures on these interaction-driven 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 half-metal/superconductor heterostructure, by tuning the direction of the magnetization of the half-metal 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 counter-propagating 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 two-dimensional 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 group-based self-consistent mean-field approach, and calculating the self-energy correction due to the topological insulator, we find that the spin texture and the spin-momentum 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 anti-spinons, rather than electrons and holes. They can participate in the transport of spinons, leading to nontrivial properties of the spin transport.
  • Using a four-band Hamiltonian, we study the phase boundary of spin-polarized-current state (SPCS) of interacting electrons in bilayer graphene. The model of spin-polarized-current 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 first-principles 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 narrow-gap 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 spin-polarized 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 triplet-pairing superconductivity as well as two-dimensional spintronics.
  • Recent experiments on Zn-substituted 122-type iron-based 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 first-principle calculations of the electronic structure reveal that the Zn 3d orbitals are far below the Fermi level and chemically inactive, while the Zn 4s-orbital is partially occupied and its wave function overlapping with those 3d-orbitals of neighboring Fe-ions. This suggests that the impurity effect is originating in the Zn 4s-orbital, not its 3d-orbitals. Employing a phenomenological two-orbital lattice model for 122-FeSCs and the self-consistent Bogoliubov-de Gennes equations, we study how the Zn-impurities 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 parity-like 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 two-leg fermionic ladder model with the number-conserving pair hopping that hosting the Majorana bound states. By utilizing density matrix renormalization group and exact diagonalization, we identify a novel type-II topological Meissner $($\emph{topo}-Meissner$)$ phase $($as distinguished from the low-field type-I \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 ground-state 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 low-dimensional quantum systems.
  • We have investigated the complex two-dimensional electron system with unequal Rashba and Dresselhaus spin-orbit 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 two-orbital model originally designed for iron-based 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 low-energy in-gap 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 eight-band 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 spin-polarized-current phase for the ground-state 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$_{1-x}$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 two-orbital model for both BaFe$_{2}$As$_{2}$ and BaFe$_{2}$P$_{2}$, together with the self-consistent lattice Bogoliubov-de Gennes (BdG) equation, we calculate the phase diagram against the P content $x$ in which the doped isovalent P-atoms are treated as impurities. We show that our numerical results can qualitatively compare with the experimental measurements.
  • The pairing symmetry is examined in highly electron-doped Ba(Fe$_{1-x}$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 two-orbital model, we consider nearest-neighbor and next-nearest-neighbor 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 time-reversal 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$_{1-x}$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$_{1-x}$Co$_x$As)$_2$. A real-space 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 optimally-doped Eu(Fe$_{1-x}$Ir$_{x}$)$_{2}$As$_{2}$ ($\mathit{x}$ = 0.12) with superconducting transition temperature $\mathit{T_{SC}}$ = 22 K was investigated by single-crystal 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 tetragonal-to-orthorhombic structural phase transition was found in this compound within the experimental uncertainty, in which the spin-density-wave (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 first-principles 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 Fe-SDW 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 quasi-two-dimensional topological crystalline insulators. In this report, we explore interaction-induced 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 inter-orbital 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 Z2-like invariant that preserves parity. More specifically, the interaction-induced novel phase of the quasi-two-dimensional topological crystalline insulator is protected by mirror reflection symmetries and therefore may be realized in pnictides.
  • Topological phases originating from spin-orbit coupling have attracted great attention recently. In this work, we use cellular dynamical mean field theory with the continuous-time quantum Monte Carlo solver to study the Kane-Mele-Hubbard model supplemented with an additional third-neighbor hopping term. For weak interactions, the third-neighbor hopping term drives a topological phase transition between a topological insulator and a trivial insulator, consistent with previous fermion sign-free 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 third-neighbor hopping. Furthermore, we extend the range of interactions into the strong coupling regime and find an easy-plane anti-ferromagnetic insulating state across a wide range of third-neighbor hopping. In contrast to the weak coupling regime, no topological phase transition occurs at strong coupling, and the ground state is a trivial anti-ferromagnetic insulating state. A comprehensive finite temperature phase diagram in the interaction-third-neighbor hopping plane is provided.