Source author record

Zhong-Yi Lu

Zhong-Yi Lu appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

62works
14topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

62 published item(s)

preprint2026arXiv

Strain-triggered high-temperature superconducting transition in two-dimensional carbon allotrope

Driving non-superconducting materials into a superconducting state through specific modulation is a key focus in the field of superconductivity. Pressure is a powerful method that can switch a three-dimensional (3D) material between non-superconducting and superconducting states. In the two-dimensional (2D) case, strain engineering plays a similar role to pressure. However, purely strain-induced superconductivity in 2D systems remains exceedingly scarce. Using first-principles calculations, we demonstrate that a superconducting transition can be induced solely by applying biaxial tensile strain in a 2D carbon allotrope, THO-graphene, which is composed of triangles, hexagons, and octagons. Free-standing THO-graphene is non-superconducting. Surprisingly, the electron-phonon coupling in strained THO-graphene is enhanced strong enough to pair electrons and realize superconductivity, with the highest superconducting transition temperature reaching 45 K. This work not only provides a notable example of controlling metal-superconductor transition in 2D system just via strain, but also sets a new record of superconducting transition temperature for 2D elemental superconductors.

preprint2026arXiv

Strategic Over-Parameterization for Generalizable Low-Rank Adaptation

Adapting large language models (LLMs) to downstream tasks via full fine-tuning is increasingly impractical due to its computational and memory demands. Parameter-efficient fine-tuning (PEFT) approaches such as Low-Rank Adaptation (LoRA) mitigate this by confining updates to a compact set of trainable parameters, but this aggressive reduction often sacrifices generalization, especially under transfer across heterogeneous tasks and domains. We revisit the tension between parameter efficiency and adaptation capacity, and ask whether the two are truly at odds. We answer in the negative by introducing LoRA-Over, a framework grounded in a simple principle: enrich the optimization landscape during training, then collapse the enrichment at inference. LoRA-Over injects auxiliary parameters into the low-rank adapters during training to broaden the effective hypothesis space, and through a decomposition-based reformulation folds them back into a standard low-rank structure with negligible reconstruction error, keeping inference cost identical to vanilla LoRA. Since not all weight matrices benefit equally from added capacity, we further propose two scheduling strategies, one statically predefined and one dynamically determined at runtime, that direct extra capacity where most needed. We evaluate LoRA-Over on language understanding (GLUE, T5-Base), dialogue (MT-Bench), arithmetic reasoning (GSM8K), and code generation (HumanEval), using LLaMA 2-7B and LLaMA 3.1-8B. Across all benchmarks and scales, LoRA-Over consistently outperforms vanilla LoRA, showing that principled over-parameterization designed to vanish at inference is an effective lever for improving PEFT generalization. Code will be released upon acceptance.

preprint2025arXiv

Correlated electronic structure of high-temperature superconductor Ba$_2$CuO$_{3+δ}$

Cuprate superconductors have attracted extensive attention due to high critical temperatures. Conventional cuprates typically contain perfect CuO$_2$ planes which are considered as a key factor to superconductivity since the superconductivity takes place in them. However, in Ba$_2$CuO$_{3+δ}$ with $δ=0.2$ and O-depleted CuO$_2$ planes, superconductivity still arises even with a transition temperature as high as 73 K. Using combined density functional theory and dynamical mean-field theory (DFT+DMFT) calculations, we investigated the electronic correlation and electronic structure of Ba$_2$CuO$_{3.25}$ with alternating quasi-one-dimensional (1D) CuO planes and O-depleted CuO$_2$ planes. We find that although different from the usual cuprates, the Cu atoms are still dominated by a 3$d^9$ configuration and the system is of a new kind of correlated single-orbital physics. The quasi-1D CuO planes, composed of parallel Cu-O chains, are slightly hole-doped quasi-1D Mott insulator, while the O-depleted CuO$_2$ planes are more hole doped, with a 2D correlated electronic structure, and may host superconductivity.

preprint2025arXiv

Kinetically accessible 1D magnetic chains of transition-metal chalcogenides and halides on van der Waals surfaces

One-dimensional (1D) chains offer unique opportunities for nanoelectronics and spintronics, yet their experimental realization remains challenging because 1D motifs are often thermodynamically disfavored relative to higher-dimensional phases. Here we present a high-throughput first-principles exploration of 1D single-atomic transition-metal chalcogenide and halide chains, screening 6,832 candidates constructed from binary combinations of 28 metals and 8 non-metals. To assess kinetic accessibility, we compare the formation energetics of 1D chains with competing two-dimensional polymorphs at the nucleation stage across relevant chemical-potential windows, using nucleation-stage thermodynamic selectivity as a proxy. This workflow identifies 183 kinetically accessible 1D chains. Interpretable machine-learning analysis reveals two simple stability descriptors as key drivers of 1D stabilization. The accessible chains exhibit diverse magnetic configurations with different magnetic characters. We further uncover their pronounced magnetoelastic couplings, exemplified by CrTe with giant magnetostriction reaching 5.93%. Finally, we show that selected metallic ferromagnetic chains retain robust edge magnetism on superconducting substrates, laying the groundwork for proximity-induced topological superconductivity and Majorana zero modes.

preprint2023arXiv

Emergent Electronic Kagome Lattice in Correlated Charge-Density-Wave State of 1T-TaS$_2$

Quantum materials with tunable correlated and/or topological electronic states, such as the electronic Kagome lattice, provide an ideal platform to study the exotic quantum properties. However, the real-space investigations on the correlated electronic Kagome lattice have been rarely reported. Herein, we report on the electronic Kagome lattice emerging in the correlated charge-density-wave (CDW) state of 1T-TaS$_2$ at ~200 K via variable-temperature scanning tunneling microscopy (VT-STM). This emergent Kagome lattice can be considered a fractional electron-filling superstructure with reduced translational and rotational symmetries, confirmed by STM measurements and density functional theory simulations. The characteristic band structure and density of states of this electronic Kagome lattice are further explored based on theoretical calculations. Our results demonstrate a self-organized electronic Kagome lattice from the correlated CDW state via the effective tuning parameter of temperature and provide a platform to directly explore the interplay of correlated electrons and topological physics.

preprint2022arXiv

An Anderson impurity interacting with the helical edge states in a quantum spin Hall insulator

Using the natural orbitals renormalization group (NORG) method, we have investigated the screening of the local spin of an Anderson impurity interacting with the helical edge states in a quantum spin Hall insulator. We find that there is a local spin formed at the impurity site and the local spin is completely screened by electrons in the quantum spin Hall insulator. Meanwhile, the local spin is screened dominantly by a single active natural orbital. We then show that the Kondo screening mechanism becomes transparent and simple in the framework of natural orbitals formalism. We project the active natural orbital respectively into real space and momentum space to characterize its structure. And we confirm the spin-momentum locking property of the edge states based on the occupancy of a Bloch state in the edge to which the impurity couples. Furthermore, we study the dynamical property of the active natural orbital represented by the local density of states, from which we observe the Kondo resonance peak.

preprint2022arXiv

Compressing LSTM Networks by Matrix Product Operators

Long Short Term Memory(LSTM) models are the building blocks of many state-of-the-art natural language processing(NLP) and speech enhancement(SE) algorithms. However, there are a large number of parameters in an LSTM model. This usually consumes a large number of resources to train the LSTM model. Also, LSTM models suffer from computational inefficiency in the inference phase. Existing model compression methods (e.g., model pruning) can only discriminate based on the magnitude of model parameters, ignoring the issue of importance distribution based on the model information. Here we introduce the MPO decomposition, which describes the local correlation of quantum states in quantum many-body physics and is used to represent the large model parameter matrix in a neural network, which can compress the neural network by truncating the unimportant information in the weight matrix. In this paper, we propose a matrix product operator(MPO) based neural network architecture to replace the LSTM model. The effective representation of neural networks by MPO can effectively reduce the computational consumption of training LSTM models on the one hand, and speed up the computation in the inference phase of the model on the other hand. We compare the MPO-LSTM model-based compression model with the traditional LSTM model with pruning methods on sequence classification, sequence prediction, and speech enhancement tasks in our experiments. The experimental results show that our proposed neural network architecture based on the MPO approach significantly outperforms the pruning approach.

preprint2022arXiv

Hydrogenation induced magnetic and electronic transitions in monolayer electride Gd$_2$C: A first-principles study

The recently synthesized two-dimensional electride Gd$_2$C was proposed to be a ferromagnetic metal that possesses multiple pairs of Weyl points and may display a large anomalous Hall conductivity [Liu \textit{et al.}, Phys. Rev. Lett. \textbf{125}, 187203 (2020)]. In view of its layered structure, here we carry out first-principles studies on the magnetic and electronic properties of Gd$_2$C in the ultrathin monolayer limit. We find that monolayer Gd$_2$C remains ferromagnetic like the bulk form and the hydrogenation can effectively tune its magnetism and electronic structure. With one-sided coverage of hydrogen atoms, monolayer Gd$_2$C becomes a half-metal with one spin channel around the Fermi level. For two-sided hydrogenation, monolayer Gd$_2$C transforms to an antiferromagnetic insulator with a band gap of 0.8 eV. Our studies show that monolayer electride Gd$_2$C can perform multiple magnetic and electronic transitions with different levels of hydrogenation and may be also adopted to construct a planar heterojunction with selective area adsorption of hydrogen atoms, which has promising applications in future electronic and spintronic devices.

preprint2022arXiv

Magnetic correlation between two local spins in a quantum spin Hall insulator

Two spins located at the edge of a quantum spin Hall insulator may interact with each other via indirect spin-exchange interaction mediated by the helical edge states, namely the RKKY interaction, which can be measured by the magnetic correlation between the two spins. By means of the newly developed natural orbitals renormalization group (NORG) method, we investigated the magnetic correlation between two Kondo impurities interacting with the helical edge states, based on the Kane-Mele model defined in a finite zigzag graphene nanoribbon with spin-orbital coupling (SOC). We find that the SOC effect breaks the symmetry in spatial distribution of the magnetic correlation, leading to anisotropy in the RKKY interaction. Specifically, the total correlation is always ferromagnetic (FM) when the two impurities are located at the same sublattice, while it is always antiferromagnetic (AFM) when at the different sublattices. Meanwhile, the behavior of the in-plane correlation is consistent with that of the total correlation. However, the out-of-plane correlation can be tuned from FM to AFM by manipulating either the Kondo coupling or the interimpurity distance. Furthermore, the magnetic correlation is tunable by the SOC, especially that the out-of-plane correlation can be adjusted from FM to AFM by increasing the strength of SOC. Dynamic properties of the system, represented by the spin-staggered excitation spectrum and the spin-staggered susceptibility at the two impurity sites, are finally explored. It is shown that the spin-staggered susceptibility is larger when the two impurities are located at the different sublattices than at the same sublattice, which is consistent with the behavior of the out-of-plane correlation. On the other hand, our study further demonstrates that the NORG is an effective numerical method for studying the quantum impurity systems.

preprint2022arXiv

Order parameter for the multichannel Kondo model at quantum criticality

A multichannel Kondo model, where two or more equivalent but independent channels of electrons compete to screen a spin-1/2 impurity, shows overcompensation of the impurity spin, leading to the non-Fermi-liquid behavior in various thermodynamic and transport properties. However, when the channel symmetry is broken, an impurity quantum phase transition can occur at zero temperature. Identification of an order parameter describing the impurity quantum phase transition is very difficult since it is beyond the conventional Landau-Ginzburg-Wilson theory. By employing the natural orbitals renormalization group method, we study both two-channel and threechannel Kondo models, from the perspective of spin correlation between the impurity and electrons in electronic channels. Here we demonstrate that by introducing the spin-correlation ratio as an order parameter we can characterize impurity quantum phase transitions driven by channel asymmetry. In particular, the universal critical exponents $β$ of the spin-correlation ratio and $ν$ of the correlation length are explicitly determined by finite-sizescaling analysis, namely, $β= 0.10(1), ν= 2.0(1)$, and $β= 0.10(1), ν= 2.5(1)$ for the two-channel and three-channel Kondo models, respectively.

preprint2022arXiv

Superconductivity in monolayer Ba$_2$N electride: a first-principles study

The exploration of superconductivity in low-dimensional materials has attracted intensive attention for decades. Based on first-principles electronic structure calculations, we have systematically investigated the electronic and superconducting properties of the two-dimensional electride Ba$_2$N in the monolayer limit. Our results show that monolayer Ba$_2$N has a low work function of 3.0 eV and a predicted superconducting transition temperature ($T_c$) of 3.4 K. The superconductivity can be further improved with the tensile strain, which results from the increase of density of states at the Fermi level as well as the enhanced coupling between inner-layer electrons and phonons. Remarkably, at the 4$\%$ tensile strain, the acoustic branches have noticeable softening at the K point of Brillouin zone and the superconducting $T_c$ can reach 10.8 K. The effect of lattice strain on the electron transfer from the superficial region to the inner-layer region of monolayer Ba$_2$N may also apply to other electride materials and influence their physical properties.

preprint2022arXiv

Two-dimensional anisotropic Dirac materials PtN4C2 and Pt2N8C6 with quantum spin and valley Hall effects

We propose two novel two-dimensional topological Dirac materials, planar PtN4C2 and Pt2N8C6, which exhibit graphene-like electronic structures with linearly dispersive Dirac-cone states exactly at the Fermi level. Moreover, the Dirac cone is anisotropic, resulting in anisotropic Fermi velocities and making it possible to realize orientation-dependent quantum devices. Using the first-principles electronic structure calculations, we have systemically studied the structural, electronic, and topological properties. We find that spin-orbit coupling opens a sizable topological band gap so that the materials can be classified as quantum spin Hall insulators as well as quantum valley Hall insulators. Helical edge states that reside in the insulating band gap connecting the bulk conduction and valence bands are observed. Our work not only expands the Dirac cone material family, but also provides a new avenue to searching for more two-dimensional topological quantum spin and valley Hall insulators.

preprint2022arXiv

Two-dimensional quadratic double Weyl semimetal

Unconventional Weyl semimetals have attracted intensive research interest in condensed matter physics and materials science, but they are very rare in two dimensions. In this work, based on symmetry analysis and the first-principles electronic structure calculations, we predict that the Si/Bi van der Waals heterostructure is a two-dimensional unconventional quadratic double Weyl semimetal with strong spin-orbit coupling (SOC). Although unprotected by the C3v double group symmetry of the heterostructure, the two-dimensional quadratic double Weyl semimetal is stable for compressive strains up to 6.64%. The system transforms into a trivial semimetal with further increasing strain, where the phase boundary is a two-dimensional triple degenerate semimetal state. Furthermore, the Kane-Mele tight-binding model calculations show that the quadratic double Weyl phase is derived from the competition between the Rashba SOC and the proximity-effect-enhanced intrinsic SOC. On the other hand, by breaking mirror symmetry, the quadratic double Weyl semimetal transforms into a quantum spin Hall insulator as well as a quantum valley Hall insulator phase. Thus, the Si/Bi heterostructure is an excellent platform for studying the exotic physics of two-dimensional double Weyl semimetal and other novel topological phases.

preprint2021arXiv

Coexistence of topological Weyl and nodal-ring states in ferromagnetic and ferrimagnetic double perovskites

Magnetic topological quantum materials have attracted great attention due to their exotic topological quantum physics induced by the interplay among crystalology, magnetism, and topology, which is of profound importance to fundamental research and technology applications. However, limited materials are experimentally available, most of whom are realized by magnetic impurity doping or heterostructural constructions. In this work, based on the first-principles calculations, we predict that double perovskite Ba2CdReO6 is an intrinsic ferromagnetic topological semi-half-metal, while the ferrimagnetic double perovskite with space group symmetry Fm-3m, such as Ba2FeMoO6, belongs to a topological half-metal. One pair of Weyl points and fully spin-polarized nodal-ring states are found in the vicinity of the Fermi level in Ba2CdReO6. Its two-dimensional nearly flat drumhead surface states are fully spin-polarized. In Ba2FeMoO6, however, there exist four pairs of Weyl points and two fully spin-polarized nodal-rings near the Fermi level. These topological properties are stable in the presence of spin-orbit coupling. This makes these materials be an appropriate platform for studying the emerging intriguing properties, especially for the applications in spintronics, information technology, and topological superconductivity.

preprint2021arXiv

Intrinsic ferromagnetic and antiferromagnetic axion insulators in van der Waals materials Mn\emph{X}$_{2}$\emph{B}$_{2}$\emph{T}$_{6}$ family

The MnBi$_{2}$Te$_{4}$ family has attracted significant attention due to its rich topological states such as the quantum anomalous Hall (QAH) insulator state, the axion insulator state, and the magnetic Weyl semimetal state. Nevertheless, the intrinsic antiferromagnetic (AFM) interlayer coupling in MnBi$_{2}$Te$_{4}$ partly hinders the realization of "high-temperature" QAH effect. Here, by using first-principles electronic structure calculations, we design a new class of materials Mn\emph{X}$_{2}$\emph{B}$_{2}$\emph{T}$_{6}$ (\emph{X}=Ge, Sn, or Pb; \emph{B}=Sb or Bi; \emph{T}=Se or Te) based on the \emph{X}$_{2}$\emph{B}$_{2}$\emph{T}$_{5}$ structures rather than the Bi$_{2}$Te$_{3}$ family. We find that each septuple-layer Mn\emph{B}$_{2}$\emph{T}$_{4}$ is sandwiched by two [\emph{X}\emph{T}] layers, which may turn the AFM interlayer coupling into a ferromagnetic (FM) coupling. The calculations specifically demonstrate that \emph{MnGe}$_{2}$\emph{Sb}$_{2}$\emph{Te}$_{6}$, \emph{MnGe}$_{2}$\emph{Bi}$_{2}$\emph{Te}$_{6}$, and \emph{MnPb}$_{2}$\emph{Bi}$_{2}$\emph{Te}$_{6}$ are FM axion insulators, while MnGe$_{2}$Sb$_{2}$Se$_{6}$, MnGe$_{2}$Bi$_{2}$Se$_{6}$, MnSn$_{2}$Sb$_{2}$Te$_{6}$, and MnSn$_{2}$Bi$_{2}$Te$_{6}$ are A-type AFM axion insulators. These seven materials all have an out-of-plane easy axis of magnetization. The Mn\emph{X}$_{2}$\emph{B}$_{2}$\emph{T}$_{6}$ family thus offers a promising platform beyond the MnBi$_{2}$Te$_{4}$ family for the realization of quantized magnetoelectric effect and "high-temperature" QAH effect in future experiments.

preprint2021arXiv

Pressure induced superconductivity in WB2 and ReB2 through modifying the B layers

The recent discovery of superconductivity up to 32 K in the pressurized MoB2 reignites the interests in exploring high-Tc superconductors in transition-metal diborides. Inspired by that work, we turn our attention to the 5d transition-metal diborides. Here we systematically investigate the responses of both structural and physical properties of WB2 and ReB2 to external pressure, which possess different types of boron layers. Similar to MoB2, the pressure-induced superconductivity was also observed in WB2 above 60 GPa with a maximum Tc of 15 K at 100 GPa, while no superconductivity was detected in ReB2 in this pressure range. Interestingly, the structures at ambient pressure for both WB2 and ReB2 persist to high pressure without structural phase transitions. Theoretical calculations suggest that the ratio of flat boron layers in this class of transition-metal diborides may be crucial for the appearance of high Tc. The combined theoretical and experimental results highlight the effect of geometry of boron layers on superconductivity and shed light on the exploration of novel high-Tc superconductors in borides.

preprint2021arXiv

Two-dimensional Dirac nodal-line semimetal protected by symmetry

Dirac nodal line semimetals (DNLSs) host relativistic quasiparticles in their one-dimensional (1D) Dirac nodal line (DNL) bands that are protected by certain crystalline symmetries. Their novel low-energy fermion quasiparticle excitations and transport properties invite studies of relativistic physics in the solid state where their linearly dispersing Dirac bands cross at continuous lines with four-fold degeneracy. In materials studied up to now, the four-fold degeneracy, however, has been vulnerable to suppression by the ubiquitous spin-orbit coupling (SOC). Despite the current effort to discover 3D DNLSs that are robust to SOC by theory, positive experimental evidence is yet to emerge. In 2D DNLSs, because of the decreased total density of states as compared with their 3D counterparts, it is anticipated that their physical properties would be dominated by the electronic states defined by the DNL. It has been even more challenging, however, to discover robust 2D DNLSs against SOC because of their lowered symmetry; no such materials have yet been predicted by theory. By combining molecular beam epitaxy growth, STM, nc-AFM characterisation, with DFT calculations and space group theory analysis, here we reveal a novel class of 2D crystalline DNLSs that host the exact symmetry that protects them against SOC. The discovered quantum material is a brick phase 3-AL Bi(110), whose symmetry protection and thermal stability are imparted by the compressive vdW epitaxial growth on black phosphorus substrates. The BP substrate templates the growth of 3-AL Bi(110) nano-islands in a non-symmorphic space group structure. This crystalline symmetry protects the DNL electronic phase against SOC independent of any orbital or elemental factors. We theoretically establish that this intrinsic symmetry imparts a general, robust protection of DNL in a series of isostructural 2D quantum materials.

preprint2020arXiv

AFeSe2 (A=Tl, K, Rb, or Cs): Iron-based superconducting analog of the cuprates

It has long been a challenging task to find compounds with similar crystal and electronic structures as cuprate superconductors with low dimensionality and strong antiferromagnetic fluctuations. The parent compounds of cuprate superconductors are Mott insulators with strong in-plane antiferromagnetic exchange interactions between Cu moments. Here we show, based on first-principles density functional calculations, that AFeSe2 (A=Tl, K, Rb, or Cs) exhibit many of the physical properties common to the cuprate parent compounds: (1) the FeSe2 layer in AFeSe2 is similar in crystalline and electronic structures to the CuO2 plane in cuprates, although Se atoms are not coplanar to the square Fe-lattice; (2) they are antiferromagnetic insulators, but with relatively small charge excitation gaps; (3) their ground states are Neel antiferromagnetic ordered, similar as in cuprates; and (4) the antiferromagnetic exchange interactions between Fe moments are larger than in other iron-based superconducting materials, but comparable to those in cuprates. Like cuprates, these compounds may become high-Tc superconductors upon doping of charge carriers either by chemical substitution or intercalation or by liquid or solid gating.

preprint2020arXiv

Bethe-Slater-curve-like behavior and interlayer spin-exchange coupling mechanisms in two-dimensional magnetic bilayers

Layered magnets have recently received tremendous attention, however, spin-exchange coupling mechanism across their interlayer regions is yet to be revealed. Here, we report a Bethe-Slater-curve (BSC) like behavior in nine transition metal dichalcogenide bilayers (MX2, M=V, Cr, Mn; X=S, Se, Te) and established interlayer spin-exchange coupling mechanisms at their van der Waals gaps using first-principle calculations. The BSC-like behavior offers a distance-dependent interlayer anti-ferromagnetic (AFM) to ferromagnetic (FM) transition. This phenomenon is explained with the spin-exchange coupling mechanisms established using bilayer CrSe2 as a prototype in this work. The Se pz wavefunctions from two adjacent interfacial Se sublayers overlap at the interlayer region. The spin alignment of the region determines interlayer magnetic coupling. At a shorter interlayer distance, Pauli repulsion at the overlapped region dominates and thus favors anti-parallel oriented spins leading to interlayer AFM. For a longer distance, kinetic energy gain of polarized electrons across the bilayer balances the Pauli repulsion and the bilayer thus prefers an interlayer FM state. In light of this, the AFM-FM transition is a result of competition between Pauli and Coulomb repulsion and kinetic energy gain. All these results open a new route to tune interlayer magnetism and the revealed spin-exchange coupling mechanisms are paramount additions to those previously established ones.

preprint2020arXiv

Compressing deep neural networks by matrix product operators

A deep neural network is a parametrization of a multilayer mapping of signals in terms of many alternatively arranged linear and nonlinear transformations. The linear transformations, which are generally used in the fully connected as well as convolutional layers, contain most of the variational parameters that are trained and stored. Compressing a deep neural network to reduce its number of variational parameters but not its prediction power is an important but challenging problem toward the establishment of an optimized scheme in training efficiently these parameters and in lowering the risk of overfitting. Here we show that this problem can be effectively solved by representing linear transformations with matrix product operators (MPOs), which is a tensor network originally proposed in physics to characterize the short-range entanglement in one-dimensional quantum states. We have tested this approach in five typical neural networks, including FC2, LeNet-5, VGG, ResNet, and DenseNet on two widely used data sets, namely, MNIST and CIFAR-10, and found that this MPO representation indeed sets up a faithful and efficient mapping between input and output signals, which can keep or even improve the prediction accuracy with a dramatically reduced number of parameters. Our method greatly simplifies the representations in deep learning, and opens a possible route toward establishing a framework of modern neural networks which might be simpler and cheaper, but more efficient.

preprint2019arXiv

Combined spontaneous symmetry-breaking and symmetry-protected topological order from cluster charge interaction

The study of symmetry-protected topological states in presence of electron correlations has recently aroused great interest as rich and exotic phenomena can emerge. Here, we report a concrete example by employing large-scale unbiased quantum Monte Carlo study of the Kane-Mele model with cluster charge interactions. The ground-state phase diagram for the model at half filling is established. Our simulation identifies the coexistence of a symmetry-protected topological order with a symmetry-breaking Kekul$\acute{e}$ valence bond order and shows that the spontaneous symmetry-breaking is accompanied by an interaction-driven topological phase transition (TPT). This TPT features appearance of zeros of single-particle Green's function and gap closing in spin channel rather than single-particle excitation spectrum, and thus has no mean-field correspondence.

preprint2019arXiv

Correlation Effects in Quadrupole Insulators: a Quantum Monte Carlo Study

The quadrupole insulator, a high-order topological insulator, with on-site Hubbard interaction is numerically studied by large-scale projector quantum Monte Carlo (PQMC) simulations. The Green's function formalism is successfully used to characterize topological properties in interacting quadrupole insulators for the first time. We find that the topological quadrupole insulator is stable against weak interactions and turns into a trivial antiferromagnetic (AFM) insulator by a continuous topological phase transition (TPT) for strong interactions. The critical exponents related to the TPT are estimated to be $ν=0.67(4)$, $β=0.40(2)$, which are distinct from those of the known AFM transitions and suggest a new universality class.

preprint2019arXiv

Interlayer quantum transport in Dirac semimetal BaGa$_2$

Quantum limit is quite easy to achieve once the band crossing exists exactly at the Fermi level ($E_F$) in topological semimetals. In multilayered Dirac fermion system, the density of Dirac fermions on the zeroth Landau levels (LLs) increases in proportion to the magnetic field, resulting in intriguing angle- and field-dependent interlayer tunneling conductivity near the quantum limit. BaGa$_2$ is an example of multilayered Dirac semimetal with anisotropic Dirac cone close to $E_F$, providing a good platform to study its interlayer transport properties. In this paper, we report the negative interlayer magnetoresistance (NIMR, I//c and B//c) induced by the tunneling of Dirac fermions on the zeroth LLs of neighbouring Ga layers in BaGa$_2$. When the field deviates from the c-axis, the interlayer resistivity $ρ_{zz}(θ)$ increases and finally results in a peak with the field perpendicular to the c-axis. These unusual interlayer transport properties (NIMR and resistivity peak with B$\perp$c) are observed together for the first time in Dirac semimetal under ambient pressure and are well explained by the model of tunneling between Dirac fermions in the quantum limit.

preprint2019arXiv

Quantum spin Hall effect in monolayer and bilayer TaIrTe$_{4}$

Generally, stacking two quantum spin Hall insulators gives rise to a trivial insulator. Here, based on first-principles electronic structure calculations, we confirm that monolayer TaIrTe$_{4}$ is a quantum spin Hall insulator and remarkably find that bilayer TaIrTe$_{4}$ is still a quantum spin Hall insulator. Theoretical analysis indicates that the covalent-like interlayer interaction in combination with the small bandgap at time-reversal invariant $Γ$ point results in new band inversion in bilayer TaIrTe$_{4}$, namely, the emergence of quantum spin Hall phase. Meanwhile, a topological phase transition can be observed by increasing the interlayer distance in bilayer TaIrTe$_{4}$. Considering that bulk TaIrTe$_{4}$ is a type-II Weyl semimetal, layered TaIrTe$_{4}$ thus provides an ideal platform to realize different topological phases at different dimensions.

preprint2019arXiv

Strong coupling superconductivity in trilayer film LiB$_2$C$_2$

Coupling between $σ$-bonding electrons and phonons is generally very strong. To metallize $σ$-electrons provides a promising route to hunt for new high-T$_c$ superconductors. Based on this picture and first-principles density functional calculation with Wannier interpolation for electronic structure and lattice dynamics, we predict that trilayer film LiB$_2$C$_2$ is a good candidate to realize this kind of high-T$_c$ superconductivity. By solving the anisotropic Eliashberg equations, we find that free-standing trilayer LiB$_2$C$_2$ is a phonon-mediated superconductor with T$_c$ exceeding the liquid-nitrogen temperature at ambient pressure. The transition temperature can be further raised to 125 K by applying a biaxial tensile strain.

preprint2018arXiv

Natural orbitals renormalization group approach to a Kondo singlet

A magnetic impurity embedded in a metal host is collectively screened by a cloud of conduction electrons to form a Kondo singlet below a characteristic energy scale $T_K$, the Kondo temperature, through the mechanism of the Kondo effect. We have reinvestigated the Kondo singlet by means of the newly developed natural orbitals renormalization group (NORG) method. We find that, in the framework of natural orbitals formalism, the Kondo screening mechanism becomes transparent and simple, while the intrinsic structure of Kondo singlet is clearly resolved. For a single impurity Kondo system, there exits a single active natural orbital which screens the magnetic impurity dominantly. In the perspective of entanglement, the magnetic impurity is entangled dominantly with the active natural orbital, i.e., the subsystem formed by the active natural orbital and the magnetic impurity basically disentangles from the remaining system. We have also studied the structures of the active natural orbital respectively projected into real space and momentum space. Moreover, the dynamical properties, represented by one-particle Green's functions defined at impurity site with active natural orbital, were obtained by using correction vector method. In order to clarify the spatial extension of the Kondo screening cloud, the concept of Kondo correlation energy was introduced. With this concept we obtain a characteristic length scale beyond which the Kondo screening cloud is hardly detected in experiment. Our numerical results indicate that this characteristic length scale usually is just a few nanometers, which interprets why it is difficult to detect the Kondo screening cloud experimentally in a metal host.

preprint2016arXiv

Bona fide interaction-driven topological phase transition in correlated SPT states

It is expected that the interplay between non-trivial band topology and strong electron correlation will lead to very rich physics. Thus a controlled study of the competition between topology and correlation is of great interest. Here, employing large-scale quantum Monte Carlo (QMC) simulations, we provide a concrete example of the Kane-Mele-Hubbard (KMH) model on an AA stacking bilayer honeycomb lattice with inter-layer antiferromagnetic interaction. Our simulation identified several different phases: a quantum spin-Hall insulator (QSH), a $xy$-plane antiferromagnetic Mott insulator ($xy$-AFM) and an inter-layer dimer-singlet insulator (dimer-singlet). Most importantly, a bona fide topological phase transition between the QSH and the dimer-singlet insulators, purely driven by the inter-layer antiferromagnetic interaction is found. At the transition, the spin and charge gap of the system close while the single-particle excitations remain gapped, which means that this transition has no mean field analogue and it can be viewed as a transition between bosonic SPT states. At one special point, this transition is described by a $(2+1)d$ $O(4)$ nonlinear sigma model (NLSM) with {\it exact} $SO(4)$ symmetry, and a topological term at {\it exactly} $Θ= π$. Relevance of this work towards more general interacting SPT states is discussed.

preprint2016arXiv

Diagnosis of interaction-driven topological phase via exact diagonalization

We propose a general scheme for diagnosing interaction-driven topological phases in the weak interaction regime using exact diagonalization (ED). The scheme comprises the analysis of eigenvalues of the point-group operators for the many-body eigenstates and the correlation functions for physical observables to extract the symmetries of the order parameters and the topological numbers of the underlying ground states at the thermodynamic limit from a relatively small size system afforded by ED. As a concrete example, we investigate the interaction effects on the half-filled spinless fermions on the checkerboard lattice with a quadratic band crossing point. Numerical results support the existence of a spontaneous quantum anomalous Hall phase purely driven by a nearest-neighbor weak repulsive interaction, separated from a nematic Mott insulator phase at strong repulsive interaction by a first-order phase transition.

preprint2016arXiv

Perfect charge compensation in extremely large magnetoresistance materials LaSb and LaBi revealed by the first-principles calculations

By the first-principles electronic structure calculations, we have systematically studied the electronic structures of recently discovered extremely large magnetoresistance (XMR) materials LaSb and LaBi. We find that both LaSb and LaBi are semimetals with the electron and hole carriers in perfect balance. The calculated carrier densities in the order of $10^{20}$ cm$^{-3}$ are in good agreement with the experimental values, implying long mean free time of carriers and thus high carrier mobilities. With a semiclassical two-band model, the perfect charge compensation and high carrier mobilities naturally explain (i) the XMR observed in LaSb and LaBi; (ii) the non-saturating quadratic dependence of XMR on external magnetic field; and (iii) the resistivity plateau in the turn-on temperature behavior at very low temperatures. The explanation of these features without resorting to the topological effect indicates that they should be the common characteristics of all perfectly electron-hole compensated semimetals.

preprint2016arXiv

Quantum critical point of Dirac fermion mass generation without spontaneous symmetry breaking

We study a lattice model of interacting Dirac fermions in $(2+1)$ dimension space-time with an SU(4) symmetry. While increasing interaction strength, this model undergoes a {\it continuous} quantum phase transition from the weakly interacting Dirac semimetal to a fully gapped and nondegenerate phase without condensing any Dirac fermion bilinear mass operator. This unusual mechanism for mass generation is consistent with recent studies of interacting topological insulators/superconductors, and also consistent with recent progresses in lattice QCD community.

preprint2016arXiv

Visualizing a Bosonic Symmetry Protected Topological Phase in an Interacting Fermion Model

Symmetry protected topological (SPT) phases in free fermion and interacting bosonic systems have been classified, but the physical phenomena of interacting fermionic SPT phases have not been fully explored. Here, employing large-scale quantum Monte Carlo simulation, we investigate the edge physics of a bilayer Kane-Mele-Hubbard model with zigzag ribbon geometry. Our unbiased numerical results show that the fermion edge modes are gapped out by interaction, while the bosonic edge modes remain gapless at the $(1+1)d$ boundary, before the bulk quantum phase transition to a topologically trivial phase. Therefore, finite fermion gaps both in the bulk and on the edge, together with the robust gapless bosonic edge modes, prove that our system becomes an emergent bosonic SPT phase at low energy, which is, for the first time, directly observed in an interacting fermion lattice model.

preprint2015arXiv

Correlated Dirac semimetal by periodized cluster dynamical mean-field theory

The periodized cluster dynamical mean-field theory (PCDMFT) combined with exact diagonalization as impurity solver has been applied to the half-filled standard Hubbard model on the honeycomb lattice. A correlated Dirac semimetal is found for weak interactions and it transforms into an antiferromagnetic insulating phase for strong interactions via a first-order quantum phase transition, not intervened by a spin liquid phase in between. In this application, the PCDMFT introduces the partial translation symmetry, but cures well the problem due to the translation symmetry breaking in the cluster dynamical mean-field theory studies for the same model, which give rise to a spurious insulating phase in the weakly interacting region.

preprint2015arXiv

First-principles study of FeSe epitaxial films on SrTiO3

The discovery of high temperature superconductivity in FeSe films on SrTiO3 substrate has inspired great experimental and theoretical interests. First-principles density functional theory calculations, which have played an important role in the study of bulk iron-based superconductors, also participate in the investigation of interfacial superconductivity. In this article, we review the calculation results on the electronic and magnetic structures of FeSe epitaxial films, emphasizing on the interplay between different degrees of freedom, such as charge, spin, and lattice vibrations. Furthermore, the comparison between FeSe monolayer and bilayer films on SrTiO3 is discussed.

preprint2015arXiv

First-principles study of magnetic frustration in FeSe epitaxial films on SrTiO$_3$

The effects of electron doping and phonon vibrations on the magnetic properties of monolayer and bilayer FeSe epitaxial films on SrTiO$_3$ have been studied, respectively, using first-principles calculations with van der Waals correction. For monolayer FeSe epitaxial film, the combined effect of electron doping and phonon vibrations readily leads to magnetic frustration between the collinear antiferromagnetic state and the checkerboard antiferromagnetic Néel state. For bilayer FeSe epitaxial film, such magnetic frustration is much more easily induced by electron doping in its bottom layer than its top layer. The underlying physics is that the doped electrons are accumulated at the interface between the FeSe layers and the substrate. These results are consistent with existing experimental studies.

preprint2015arXiv

Natural Orbitals Renormalization Group Approach to the Two-Impurity Kondo Critical Point

The problem of two magnetic impurities in a normal metal exposes the two opposite tendencies in the formation of a singlet ground state, driven respectively by the single-ion Kondo effect with conduction electrons to screen impurity spins or the Ruderman-Kittel-Kasuya-Yosida interaction between the two impurities to directly form impurity spin singlet. However, whether the competition between these two tendencies can lead to a quantum critical point has been debated over more than two decades. Here, we study this problem by applying the newly proposed natural orbitals renormalization group method to a lattice version of the two-impurity Kondo model with a direct exchange $K$ between the two impurity spins. The method allows for unbiased accessing the ground state wave functions and low-lying excitations for sufficiently large system sizes. We demonstrate the existence of a quantum critical point, characterized by the power-law divergence of impurity staggered susceptibility with critical exponent $γ= 0.60(1)$, on the antiferromagnetic side of $K$ when the interimpurity distance $R$ is even lattice spacing, while a crossover behavior is recovered when $R$ is odd lattice spacing. These results have ultimately resolved the long-standing discrepancy between the numerical renormalization group and quantum Monte Carlo studies, confirming a link of this two-impurity Kondo critical point to a hidden particle-hole symmetry predicted by the local Fermi liquid theory.

preprint2015arXiv

Nematic antiferromagnetic states in bulk FeSe

We revisit bulk FeSe through the systematic first-principles electronic structure calculations. We find that there are a series of staggered $n$-mer antiferromagnetic (AFM) states with corresponding energies below that of the collinear AFM state which is the ground state for the parent compounds of most iron-based superconductors. Here the staggered $n$-mer ($n$ any integer $>1$) means that a set of $n$ adjacent spins parallel on a line along $b$-axis with spins in antiparallel between $n$-mers and along $a$-axis. Among them, the lowest energy states are quasi-degenerate staggered dimer and staggered trimer AFM states as well as their any staggered combinations. Thus, to have the largest entropy to minimize the free energy at low temperature, the most favorable state is such a quasi-one-dimensional antiferromagnet in which along $b$-axis a variety of $n$-mers, mostly dimers and trimers, are randomly antiparallel aligned while along $a$-axis spins are antiparallel aligned, i.e. actually a nematic paramagnet. This finding accounts well for the absence of long-range magnetic order in bulk FeSe and meanwhile indicates the dominant stripe spin fluctuation and the nematicity as spin-driven.

preprint2015arXiv

Quantum Monte Carlo study of strange correlator in interacting topological insulators

Distinguishing the nontrivial symmetry-protected topological (SPT) phase from the trivial insulator phase in the presence of electron-electron interaction is an urgent question to the study of topological insulators, due to the fact that most of the topological indices defined for free electron systems are very likely unsuitable for interacting cases. In this work, we demonstrate that the strange correlator is a sensitive diagnosis to detect SPT states in interacting systems. Employing large-scale quantum Monte Carlo (QMC) simulations, we investigate the interaction-driven quantum phase transition in the Kane-Mele-Hubbard model. The transition from the quantum spin Hall insulator at weak interaction to an antiferromagnetic Mott insulator at strong interaction can be readily detected by the momentum space behavior of the strange correlator in single-particle, spin, and pairing sectors. The interaction effects on the symmetry-protected edge states in various sectors, i.e., the helical Luttinger liquid behavior, are well captured in the QMC measurements of strange correlators. Moreover, we demonstrate that the strange correlator is technically easier to implement in QMC and more robust in performance than other proposed numerical diagnoses for interacting topological states, as only static correlations are needed. The attempt in this work paves the way for using the strange correlator to study interaction-driven topological phase transitions in fermionic as well as bosonic systems.

preprint2015arXiv

Strong charge density wave fluctuation and sliding state in PdTeI with quasi-1D PdTe chains

In quasi-one-dimensional (quasi-1D) system, the charge density wave (CDW) transition temperature TCDW is usually lower than the mean-field-theory predicted TMF and a CDW fluctuation region exists between them. Here, we investigate the physical properties of PdTeI single crystal containing quasi-1D PdTe chains. Surprisingly, we find that the carrier concentration decreases gradually before the long-range CDW ordering state occurring at T1 ~ 110 K, reflecting the existence of strong CDW fluctuation with possible pseudogap state at T >> T1 because of dynamic charge separation of Pd ions (Pd3+ -> Pd2+ + Pd4+). Moreover, the sliding CDW state appears below T2 ~ 6 K. Combined such low T2 with the feature of multiple quasi-1D bands, PdTeI exhibits exotic crossover behavior from negative to huge positive magnetoresistance under magnetic field and field-induced localization. Thus, PdTeI provides a novel platform for studying the CDW fluctuation and the interplay between magnetic field and CDW state.

preprint2015arXiv

Topological invariants for interacting topological insulators: I. Efficient numerical evaluation scheme and implementations

The aim of this series of two papers is to discuss topological invariants for interacting topological insulators (TIs). In the first paper (I), we provide a paradigm of efficient numerical evaluation scheme for topological invariants, in which we demystify the procedures and techniques employed in calculating Z2 invariant and spin Chern number via zero-frequency single-particle Green's function in quantum Monte Carlo (QMC) simulations. Here we introduce a periodization process to overcome the ubiquitous finite-size effect, so that the calculated spin Chern number shows ideally quantized values. We also show that making use of symmetry properties of the underlying systems can greatly reduce the computational effort. To demonstrate the effectiveness of our numerical evaluation scheme, especially the periodization process, of topological invariants, we apply it on two independent two-dimensional models of interacting topological insulators. In the subsequent paper (II), we apply the scheme developed here to wider classes of models of interacting topological insulators, for which certain limitation of constructing topological invariant via single-particle Green's functions will be presented.

preprint2015arXiv

Topological invariants for interacting topological insulators: II. Breakdown of the Green's function formalism

Topological phase transitions in free fermion systems can be characterized by closing of single-particle gap and change in topological invariants. However, in the presence of electronic interactions, topological phase transitions are more complicated. In paper I of this series (arXiv:1510.07816), we have developed an efficient scheme to evaluate the topological invariants based on Green's function formalism. Here, in paper II, we demonstrate four interaction-drive topological phase transitions (TPTs) in two-dimensional (2D) interacting topological insulators (TIs) via large-scale quantum Monte Carlo (QMC) simulations, based on the scheme of evaluating topological invariants presented in paper I. Across these transitions, the defining symmetries of the TIs have been neither explicitly nor spontaneously broken. In the first two models, the topological invariants calculated from Green's function formalism succeed in characterizing interaction-driven TPTs. However, in the second two models, we find single-particle gap does not close and the topological invariants constructed from single-particle Green's function acquire no change across the TPTs. Unexpected breakdown of the Green's function formalism in constructing topological invariants is thus discovered. We thence classify the TPTs in interacting TIs into two categories: those have noninteracting correspondence can be characterized successfully by the topological invariants constructed from Green's functions, while for the others that do not have noninteracting correspondence, the Green's function formalism experiences a breakdown but more interesting and exciting phenomena, such as emergent collective critical modes at the transition, arise. Discussion on the success and breakdown of topological invariants constructed from the Green's function formalism in the context of symmetry protected topological (SPT) states is presented.

preprint2014arXiv

Direct observation of fragile Mott insulators on plaquette Hubbard lattices

Employing extensive cellular dynamical mean-field theory (CDMFT) calculations with exact diagonalization impurity solver, we investigate the ground state phase diagrams and non-magnetic metal-insulator transitions of the half-filled Hubbard model on two plaquette -- the 1/5 depleted and checkerboard -- square lattices. We identify three different insulators in the phase diagrams: dimer insulator, antiferromagnetic insulator, and plaquette insulator. And we demonstrate that the plaquette insulator is a novel fragile Mott insulator (FMI) which features a nontrivial one-dimensional irreducible representation of the $C_{4v}$ crystalline point-group and cannot be adiabatically connected to any band insulator with time-reversal symmetry. Furthermore, we study the non-magnetic quantum phase transitions from the metal to the FMI and find that this Mott metal-insulator transition is characterized by the splitting of the non-interacting bands due to interaction effects.

preprint2014arXiv

Electronic transport of a large scale system studied by renormalized transfer matrix method: application to armchair graphene nanoribbons between quantum wires

Study on the electronic transport of a large scale two dimensional system by the transfer matrix method (TMM) based on the Schördinger equation suffers from the numerical instability. To address this problem, we propose a renormalized transfer matrix method (RTMM) by setting up a set of linear equations from U times of multiplication of traditional transfer matrix (U=N/S}with N and S being the atom number of length and the transfer step), and smaller S is required for wider systems. Then we solve the above linear equations by Gauss elimination method and further optimize to reduce the computational complexity from O(U^3M^3) to O(UM^3), in which M is the atom number of the width. Applying RTMM, we study transport properties of large scale pure and long-range correlated disordered armchair graphene nanoribbon (AGR) (carbon atoms up to 10^6 for pure case) between quantum wire contacts. As for pure AGR, the conductance is superlinear with the Fermi energy and the conductance is linear with the width while independent of the length, showing characteristics of ballistic transport. As for disordered AGR with long-range correlation, there is metal-insulator transition induced by the correlation strength of disorder. It is straightforward to extend RTMM to investigate transport in large scale system with irregular structure.

preprint2014arXiv

Neutron Scattering Measurements of Spatially Anisotropic Magnetic Exchange Interactions in Semiconducting K0.85Fe1.54Se2 (TN=280 K)

We use neutron scattering to study the spin excitations associated with the stripe antiferromagnetic (AFM) order in semiconducting K$_{0.85}$Fe$_{1.54}$Se$_2$ ($T_N$=$280$ K). We show that the spin wave spectra can be accurately described by an effective Heisenberg Hamiltonian with highly anisotropic in-plane couplings at $T$= $5$ K. At high temperature ($T$= $300$ K) above $T_N$, short range magnetic correlation with anisotropic correlation lengths are observed. Our results suggest that, despite the dramatic difference in the Fermi surface topology, the in-plane anisotropic magnetic couplings are a fundamental property of the iron based compounds; this implies that their antiferromagnetism may originate from local strong correlation effects rather than weak coupling Fermi surface nesting.

preprint2014arXiv

Quantum Renormalization Groups Based on Natural Orbitals

We propose a new concept upon the renormalization group (RG) procedure for an interacting many-electron correlated system in the framework of natural orbitals, and formulate an algorithm for this RG approach. To demonstrate its effectiveness, we apply this new approach on a quantum cluster-impurity model with four impurities in comparison with the exact diagonalization method. We also find a shortcut to dramatically improving this RG algorithm. Further discussion is presented with the cluster dynamical mean-field theory and multi-impurity/orbital Kondo problems.

preprint2013arXiv

Investigation of atomic and electronic structures of MgOFeSe studied by the first-principles calculations

In order to assist the search of new superconductors in iron selenide materials by intercalation, we calculate the crystal and electronic structures of MgOFeSe using the first-principles density functional theory. MgOFeSe is isotructural to the parent compound of iron pnictide superconductor LaOFeAs. In LaOFeAs, the anion O$^{2-}$ is located at the center of each LaO tetrahedra. But for MgOFeSe, we find that the crystal structure with the cation Mg$^{2+}$ as the tetrahedral center in the MgO layer is energetically more stable. The low energy band structures around the Fermi surfaces of MgOFeSe are contributed mainly by Fe 3$d$ orbitals. The ground state of MgOFeSe is collinearly antiferromagnetically ordered. The height of Se atoms above the Fe-Fe layer is about 1.38 Å, which is close to the height of As from the Fe-Fe layer in the iron pnictide superconductors with optimal superconducting transition temperatures.

preprint2013arXiv

Prediction of phonon-mediated high temperature superconductivity in stoichiometric Li$_2$B$_3$C

The discovery of superconductivity in Magnesium Diborate (MgB$_2$) has stimulated great interest in the search of new superconductors with similar lattice structures. Unlike cuprate or iron-based superconductors, MgB$_2$ is indisputably a phonon-mediated high temperature superconductor. The emergence of high temperature superconductivity in this material results from the strong coupling between the boron $σ$-bonding electrons around the Fermi level and the bond-stretching optical phonon modes. Here we show, based on the first-principles calculations, that Li$_2$B$_3$C is such a good candidate of superconductor whose superconducting transition temperature (T$_c$) might be even higher than MgB$_2$. Li$_2$B$_3$C consists of alternating graphene-like boron-carbon layers and boron-boron layers with intercalated lithium atoms between them. Similar to MgB$_2$, Li$_2$B$_3$C is inherently metallic and possesses two $σ$- and two $π$-electron bands around the Fermi energy. The superconducting pairs are glued predominately by the strong interaction between boron $σ$-bonding electrons and various optical phonon modes.

preprint2013arXiv

The influence of spin-phonon coupling on antiferromagnetic spin fluctuations in FeSe under pressure: the First-principles calculations with van der Waals corrections

The electronic structures, lattice dynamics, and magnetic properties of crystal β-FeSe under hydrostatic pressure have been studied by using the first-principles electronic structure calculations with van der Waals (vdW) corrections. With applied pressures, the energy bands around the Fermi energy level consisting mainly of Fe-3d orbitals show obvious energy shifts and occupation variations, and meanwhile the frequencies of all optical phonon modes increase. Among these phonon modes, the A1g mode, which relates to the Se height from the Fe-Fe plane, shows a clear frequency jump in the range between 5 and 6 GPa. This is also the pressure range within which the highest superconducting transition temperature Tc of FeSe is reached in experiments. In comparison with the other phonon modes, the zero-point vibration of the A1g mode induces the strongest fluctuation of local magnetic moment on Fe under a pressure from 0 to 9 GPa, and the induced fluctuation reaches the maximum around 5 GPa. These results suggest that the effect of phonon via spin-phonon coupling need to be considered when exploring the superconducting mechanism in iron-based superconductors.

preprint2012arXiv

Atomic and electronic structures of FeSe monolayer and bilayer thin films on SrTiO$_3$ (001): a first-principles study

By the first-principles electronic structure calculations, we have studied electronic structures of FeSe monolayer and bilayer thin films on SrTiO$_3$ (001) with SrO-termination or TiO$_2$ termination. We find that both FeSe monolayer and bilayer on either termination behave like a slightly doped semiconductor and a collinear antiferromagnetic order on Fe ions. There is no substantial charge transfer between the FeSe layers and the substrate. FeSe is adhered to the SrTiO$_3$ surface by a dipole-dipole interaction. The Fermi surface is mainly the contribution of Fe-3d orbitals. A valence band contributed mainly by the O-$2p$ orbitals in the TiO$_2$ layer is located slightly below the Fermi level, which can become conducting upon a small doping of holes.

preprint2012arXiv

Layered pnictide-oxide Na$_2$Ti$_2$Pn$_2$O (Pn=As, Sb): a paradigm for spin density waves

From the first-principles calculations, we have studied the electronic and magnetic structures of compound Na$_2$Ti$_2$Pn$_2$O (Pn = As or Sb). We find that in the ground state Na$_2$Ti$_2$As$_2$O is a blocked checkerboard antiferromagnetic semiconductor with a small band gap of about 0.15 eV, in contrast, Na$_2$Ti$_2$Sb$_2$O is a bi-collinear antiferromagnetic semimetal, both with a small moment of about 0.5$μ_B$ around each Ti atom. We show that there is a strong Fermi surface nesting in Na$_2$Ti$_2$Pn$_2$O. And we verify that the blocked checkerboard and bi-collinear antiferromagnetic states both are the spin density waves induced by the Fermi surface nesting. A tetramer structural distortion is found in company with the formation of a blocked checkerboard antiferromagnetic order, in good agreement with the experimentally observed commensurate structural distortion but with space group symmetry retained after the anomaly happening. Further analysis and discussion in connection with experimental observations are given as well.

preprint2011arXiv

Cluster dynamical mean field theory of quantum phases on a honeycomb lattice

We have studied the ground state of the half-filled Hubbard model on a honeycomb lattice by performing the cluster dynamical mean field theory calculations with exact diagonalization on the cluster-impurity solver. Through using elaborate numerical analytic continuation, we identify the existence of a `spin liquid' from the on-site interaction U=0 to $U_c$ (between $4.6t$ and $4.85t$) with a smooth crossover correspondingly from the charge fluctuation dominating phase into the charge correlation dominating phase. The semi-metallic state exits only at U=0. We further find that the magnetic phase transition at $U_c$ from the `spin liquid' to the Néel antiferromagnetic Mott insulating phase is a first-order quantum phase transition. We also show that the charge fluctuation plays a substantial role on keeping the `spin liquid' phase against the emergence of a magnetic order.

preprint2011arXiv

Electronic structures and magnetic orders of Fe-vacancies ordered ternary iron selenides TlFe$_{1.5}$Se$_2$ and AFe$_{1.5}$Se$_2$ (A=K, Rb, or Cs)

By the first-principles electronic structure calculations, we find that the ground state of the Fe-vacancies ordered TlFe$_{1.5}$Se$_2$ is a quasi-two-dimensional collinear antiferromagnetic semiconductor with an energy gap of 94 meV, in agreement with experimental measurements. This antiferromagnetic order is driven by the Se-bridged antiferromagnetic superexchange interactions between Fe moments. Similarly, we find that crystals AFe$_{1.5}$Se$_2$ (A=K, Rb, or Cs) are also antiferromagnetic semiconductors but with a zero-gap semiconducting state or semimetallic state nearly degenerated with the ground states. Thus rich physical properties and phase diagrams are expected.

preprint2011arXiv

Exact results for intrinsic electronic transport in graphene

We present exact results for the electronic transport properties of graphene sheets connected to two metallic electrodes. Our results, obtained by transfer-matrix methods, are valid for all sheet widths and lengths. In the limit of large width-to-length ratio relevant to recent experiments, we find a Dirac-point conductivity of $2e^2/\sqrt{3}h$ and a sub-Poissonian Fano factor of $2 - 3\sqrt{3}/π\simeq 0.346$ for armchair graphene; for the zigzag geometry these are respectively 0 and 1. Our results reflect essential effects from both the topology of graphene and the electronic structure of the leads, giving a complete microscopic understanding of the unique intrinsic transport in graphene.

preprint2011arXiv

First-principles study of pressure-induced magnetic phase transitions in ternary iron selenide K$_{0.8}$Fe$_{1.6}$Se$_2$

We have studied the pressure effect on electronic structures and magnetic orders of ternary iron selenide K$_{0.8}$Fe$_{1.6}$Se$_2$ by the first-principles electronic structure calculations. At low pressure, the compound is in the blocked checkerboard antiferromagnetic (AFM) semiconducting phase, as observed by the neutron scatting measurements. Applying pressure induces two phase transitions, first from the blocked checkerboard AFM semiconducting phase to a collinear AFM metallic phase around 12 GPa, and then to a non-magnetic metallic phase around 25 GPa, respectively. Our results help to clarify the recent experimental measurements under pressure.

preprint2011arXiv

Ternary iron selenide K$_{0.8}$Fe$_{1.6}$Se$_2$ is an antiferromagnetic semiconductor

We have studied electronic and magnetic structures of K$_{0.8+x}$Fe$_{1.6}$Se$_2$ by performing the first-principles electronic structure calculations. The ground state of the Fe-vacancies ordered K$_{0.8}$Fe$_{1.6}$Se$_2$ is found to be a quasi-two-dimensional blocked checkerboard antiferromagnetic (AFM) semiconductor with an energy gap of 594 meV and a large ordering magnetic moment of 3.37 $μ_B$ for each Fe atom, in excellent agreement with the neutron scattering measurement. The underlying mechanism is the chemical-bonding-driven tetramer lattice distortion. K$_{0.8+x}$Fe$_{1.6}$Se$_2$ with finite $x$ is a doped AFM semiconductor with low conducting carrier concentration which is approximately proportional to the excess potassium content, consistent qualitatively with the infrared observation. Our study reveals the importance of the interplay between antiferromagnetism and superconductivity in these materials. This suggests that K$_{0.8}$Fe$_{1.6}$Se$_2$, instead of KFe$_2$Se$_2$, should be regarded as a parent compound from which the superconductivity emerges upon electron or hole doping.

preprint2010arXiv

Electronic and magnetic structures of ternary iron selenides AFe$_2$Se$_2$ (A=K, Cs, or Tl)

By the first-principles electronic structure calculations, we find that the ground state of ternary iron selenides AFe$_2$Se$_2$ (A=K, Cs, or Tl) is in a bi-collinear antiferromagnetic order, in which the Fe local moments ($\sim2.8μ_B$) align ferromagnetically along a diagonal direction and antiferromagnetically along the other diagonal direction on the Fe-Fe square lattice. This bi-collinear antiferromagnetic order results from the interplay among the nearest, the next nearest, and the next next nearest neighbor superexchange interactions, mediated by Se $4p$-orbitals.

preprint2010arXiv

Electronic structure of Fe1.04(Te0.66Se0.34)

We report the electronic structure of the iron-chalcogenide superconductor, Fe1.04(Te0.66Se0.34), obtained with high resolution angle-resolved photoemission spectroscopy and density functional calculations. In photoemission measurements, various photon energies and polarizations are exploited to study the Fermi surface topology and symmetry properties of the bands. The measured band structure and their symmetry characters qualitatively agree with our density function theory calculations of Fe(Te0.66Se0.34), although the band structure is renormalized by about a factor of three. We find that the electronic structures of this iron-chalcogenides and the iron-pnictides have many aspects in common, however, significant differences exist near the Gamma-point. For Fe1.04(Te0.66Se0.34), there are clearly separated three bands with distinct even or odd symmetry that cross the Fermi energy (EF) near the zone center, which contribute to three hole-like Fermi surfaces. Especially, both experiments and calculations show a hole-like elliptical Fermi surface at the zone center. Moreover, no sign of spin density wave was observed in the electronic structure and susceptibility measurements of this compound.

preprint2010arXiv

Electronic structures of ternary iron arsenides AFe$_2$As$_2$ (A=Ba, Ca, or Sr)

We have studied the electronic and magnetic structures of the ternary iron arsenides AFe$_2$As$_2$ (A = Ba, Ca, or Sr) using the first-principles density functional theory. The ground states of these compounds are in a collinear antiferromagnetic order, resulting from the interplay between the nearest and the next-nearest neighbor superexchange antiferromagnetic interactions bridged by As $4p$ orbitals. The correction from the spin-orbit interaction to the band structure is small. The pressure can reduce dramatically the magnetic moment and diminish the collinear antiferromagnetic order. Based on the calculations, we propose that the low energy dynamics of these materials is described effectively by a $t-J_H-J_1-J_2$-type model.

preprint2010arXiv

Pressure- and temperature-induced structural phase transitions of CaFe$_2$As$_2$ and BaFe$_2$As$_2$ studied in the Hund's rule correlation picture

With the proposed Hund's rule correlation picture, i.e. the fluctuating Fe local moments with the As-bridged antiferromagnetic superexchange interactions, the exceptional collapsed tetragonal phase and related phase transitions observed in CaFe$_2$As$_2$ are well understood. With the same framework, a pressure-temperature phase diagram is predicted for BaFe$_2$As$_2$ as well, in which a paramagnetic tetragonal and a collinear antiferromagnetic orthorhombic structures to nonmagnetic tetragonal structure transitions take place around 4-8 GPa and 10-15 GPa respectively, and a nonmagnetic tetragonal to a nonmagnetic collapsed tetragonal structure transition takes place over 26 GPa. Our study helps better understand the complex correlation among crystal structure, magnetism, and electronic structure in pnictides, a precondition to understand the superconductivity in pnictides.

preprint2010arXiv

The layered compound CaClFeP is an Arsenic-free high $T_c$ iron-pnictide

We first analyze why the iron pnictides with high $T_c$ superconductivity so far are As-based, by the Hund's rule correlation picture, then examine the P-based and Sb-based cases, respectively. Consequently, we propose that CaClFeP with ZrCuSiAs-type structure is an As-free high $T_c$ iron-pnictide. The subsequent density functional theory calculations show that the ground state of CaClFeP is of a collinearly antiferromagnetic order on Fe moments with structural distortion, resulting from the interplay between the strong nearest and next-nearest neighbor antiferromagnetic superexchange interactions bridged by P atoms, similar as the As-based pnictides. The other P-based pnictides are either nonmagnetic or magnetic but with weak exchange interactions. The Sb-based pnictides unlikely show high $T_c$ superconductivity because of the existence of robust ferromagnetic order.

preprint2009arXiv

Atomic and electronic structures of ternary iron arsenides $A$Fe$_2$As$_2$(001) surfaces ($A$=Ba, Sr, or Ca)

By the first-principles electronic structure calculations, we find that energetically the most favorable cleaved $A$Fe$_2$As$_2$(001) surface ($A$=Ba, Sr, or Ca) is $A$-terminated with a $(\sqrt{2}\times \sqrt{2})R45^{\circ}$ or $(1\times 2)$ order. The $(1\times 2)$ ordered structure yields a $(1\times 2)$ dimerized STM image, in agreement with the experimental observation. The $A$ atoms are found to diffuse on the surface with a small energy barrier so that the cleaving process may destroy the $A$ atoms ordering. At the very low temperatures this may result in an As-terminated surface with the $A$ atoms in randomly assembling. The As-terminated BaFe$_2$As$_2$ surface in orthorhombic phase is $(\sqrt{2}\times\sqrt{2})R45^{\circ}$ buckled, giving rise to a switchable $(\sqrt{2}\times \sqrt{2})R45^{\circ}$ STM pattern upon varying the applied bias. No any reconstruction is found for the other As-terminated surfaces. There are surface states crossing or nearby the Fermi energy in the As-terminated and $(1\times 2)$ $A$-terminated surfaces. A unified physical picture is thus established to help understand the cleaved $A$Fe$_2$As$_2$(001) surfaces.

preprint1995arXiv

Impurity Energy Level Within The Haldane Gap

An impurity bond $J{'}$ in a periodic 1D antiferromagnetic, spin 1 chain with exchange $J$ is considered. Using the numerical density matrix renormalization group method, we find an impurity energy level in the Haldane gap, corresponding to a bound state near the impurity bond. When $J{'}<J$ the level changes gradually from the edge of the Haldane gap to the ground state energy as the deviation $dev=(J-J{'})/J$ changes from 0 to 1. It seems that there is no threshold. Yet, there is a threshold when $J{'}>J$. The impurity level appears only when the deviation $dev=(J{'}-J)/J{'}$ is greater than $B_{c}$, which is near 0.3 in our calculation.

preprint1993arXiv

The Haldane Energy Gap of A Doped Linear-Chain Heisenberg Antiferromagnet

Using the valence-bond-solid (VBS) approach and the Schwinger boson mean field approximation, we study the dependence of the Haldane gap of a spin-1 linear chain Heisenberg antiferromagnet on impurity doping with different spins. The impurity spins affect the singlet pairing order parameter $Δ$ and the constraint factor $λ$. As a result, the Haldane gap is reduced by a factor $ \sim n_i^{2/3}$, with $n_i$ as the impurity concentration, and eventually collapses at $n_i \sim 1/ξ$ with $ξ$ as the VBS correlation length. This theoretical prediction can be verified by neutron scattering experiments.