Source author record

Wenhui Duan

Wenhui Duan 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

69works
10topics
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

69 published item(s)

preprint2026arXiv

DeepH-pack: A general-purpose neural network package for deep-learning electronic structure calculations

In computational physics and materials science, first-principles methods, particularly density functional theory, have become central tools for electronic structure prediction and materials design. Recently, rapid advances in artificial intelligence (AI) have begun to reshape the research landscape, giving rise to the emerging field of deep-learning electronic structure calculations. Despite numerous pioneering studies, the field remains in its early stages; existing software implementations are often fragmented, lacking unified frameworks and standardized interfaces required for broad community adoption. Here we present DeepH-pack, a comprehensive and unified software package that integrates first-principles calculations with deep learning. By incorporating fundamental physical principles into neural-network design, such as the nearsightedness principle and the equivariance principle, DeepH-pack achieves robust cross-scale and cross-material generalizability. This allows models trained on small-scale structures to generalize to large-scale and previously unseen materials. The toolkit preserves first-principles accuracy while accelerating electronic structure calculations by several orders of magnitude, establishing an efficient and intelligent computational paradigm for large-scale materials simulation, high-throughput materials database construction, and AI-driven materials discovery.

preprint2026arXiv

Two-dimensional Intrinsic Janus Structures: Design Principle and Anomalous Nonlinear Optics

Two-dimensional Janus structures have garnered rapidly growing attention across multidisciplinary fields. However, despite extensive theoretical and experimental efforts, a principle for designing intrinsic Janus materials remains elusive. Here, we propose a first-principles alloy theory based on cluster expansion, incorporating a strong repulsive interaction of a cation-mediated anion-pair cluster and refined short-range cluster-cluster competitions, to unravel the formation mechanism of intrinsic Janus structures with a distorted 1T phase among numerous competing phases. Our theory not only explains why intrinsic Janus structures are accidentally observed in RhSeCl and BiTeI which are composed of alloyed elements from different groups, but also accurately predicts a wide range of 1T-like intrinsic Janus materials that are ready for synthesis. Intriguingly, as demonstrated in the case of RhSeCl, we reveal that intrinsic Janus materials can exhibit anomalous second-harmonic generation (SHG) with a distinct quantum geometric effect, originating from strong lattice and chemical-potential mirror asymmetry. Furthermore, a novel skin effect unexpectedly emerges in finite-thickness RhSeCl, accompanied by a hidden SHG effect within the bulk region. Our theory paves the way for the ab initio design of intrinsic Janus materials, significantly accelerating progress in Janus science.

preprint2024arXiv

Enhancement of Ising superconductivity in monolayer NbSe$_2$ via surface fluorination

Recently discovered Ising superconductors have garnered considerable interest due to their anomalously large in-plane upper critical fields ($B_{c2}$). However, the requisite strong spin-orbital coupling in the Ising pairing mechanism generally renders these superconductors heavy-element dominant with notably low superconducting transition temperatures ($T_c$). Here, based on the Migdal-Eliashberg theory and the mean-field Bogoliubov-de Gennes Hamiltonian, we demonstrate a significant enhancement of Ising superconductivity in monolayer NbSe$_2$ through surface fluorination, as evidenced by concomitant improvements in $T_c$ and $B_{c2}$. This enhancement arises from three predominant factors. Firstly, fluorine atoms symmetrically and stably adhere to both sides of the monolayer NbSe$_2$, thereby maintaining the out-of-plane mirror symmetry and locking carrier spins out-of-plane. Secondly, fluorination suppresses the charge density wave in monolayer NbSe$_2$ and induces a van Hove singularity in the vicinity of the Fermi level, leading to a marked increase in the number of carriers and, consequently, strengthening the electron-phonon coupling (EPC). Lastly, the appearance of fluorine-related, low-frequency phonon modes further augments the EPC. Our findings suggest a promising avenue to elevate $T_c$ in two-dimensional Ising superconductors without compromising their Ising pairing.

preprint2022arXiv

Berry-curvature-induced linear magnetotransport in magnetic Weyl semimetals

Magnetotransport such as the giant magnetoresistance and Hall effect lies at the heart of fundamental physics and technologies. Recently, some experiments have clearly demonstrated linear magnetotransport (LMT) proportional to magnetic field but the underlying physical mechanism is still unclear. In this work, we show that Berry curvature effect is a new mechanism dominating the LMT. The Berry-curvature-induced LMT widely exists in 66 out of 122 magnetic point groups. For typical magnetic Weyl semimetals Co$_3$Sn$_2$S$_2$ and ferromagnetic MnBi$_2$Te$_4$, Berry curvature induces LMT conductivities reaching orders of $10^4$ and $10^2$ ${\rm Ω^{-1}m^{-1}}$ per tesla, respectively, which are tunable through magnetization canting induced by moderate magnetic fields. We further reveal that Berry-curvature-induced LMT can be detected by Hall effect and especially intrinsic magnetoresistance exceeding $100\%$ per tesla insensitive to the sample quality. Our results agree with recent experiments and uncover the important role of Berry curvature in LMT.

preprint2022arXiv

Coexistence of extended flat band and Kekulé order in Li-intercalated graphene

Doping graphene near the 1/4 filling to shift the extended flat band and van Hove singularity below E$_F$ has been highly desirable. Here we report the experimental observation of an extended flat band below E$_F$ in Li-intercalated graphene. Strong electron-phonon interaction is clearly identified by notable kinks in the band dispersion. Moreover, the evolution of the band structure upon Li intercalation shows that the extended flat band and the Kekulé order emerge simultaneously. Our work provides opportunities for investigating flat band related instabilities and its interplay with the Kekulé order

preprint2022arXiv

Control of phase ordering and elastic properties in phase field crystals through three-point direct correlation

Effects of three-point direct correlation on properties of the phase field crystal (PFC) modeling are examined, for the control of various ordered and disordered phases and their coexistence in both three-dimensional and two-dimensional systems. Such effects are manifested via the corresponding gradient nonlinearity in the PFC free energy functional that is derived from classical density functional theory. Their significant impacts on the stability regimes of ordered phases, phase diagrams, and elastic properties of the system, as compared to those of the original PFC model, are revealed through systematic analyses and simulations. The nontrivial contribution from three-point direct correlation leads to the variation of the critical point of order-disorder transition to which all the phase boundaries in the temperature-density phase diagram converge. It also enables the variation and control of system elastic constants over a substantial range as needed in modeling different types of materials with the same crystalline structure but different elastic properties. The capability of this PFC approach in modeling both solid and soft matter systems is further demonstrated through the effect of three-point correlation on controlling the vapor-liquid-solid coexistence and transitions for body-centered cubic (bcc) phase and on achieving the liquid-stripe or liquid-lamellar phase coexistence. All these provide a valuable and efficient method for the study of structural ordering and evolution in various types of material systems.

preprint2022arXiv

Controllable chirality and band gap of quantum anomalous Hall insulators

Finding guiding principles to optimize properties of quantum anomalous Hall (QAH) insulators is of pivotal importance to fundamental science and applications. Here, we build a first-principles QAH material database of chirality and band gap, explore microscopic mechanisms determining the QAH material properties, and obtain a general physical picture that can comprehensively understand the QAH data. Our results reveal that the usually neglected Coulomb exchange is unexpectedly strong in a large class of QAH materials, which is the key to resolve experimental puzzles. Moreover, we identify simple indicators for property evaluation and suggest material design strategies to control QAH chirality and gap by tuning cooperative or competing contributions via magnetic co-doping, heterostructuring, spin-orbit proximity, etc. The work is valuable to future research of magnetic topological physics and materials.

preprint2022arXiv

Deep-Learning Density Functional Theory Hamiltonian for Efficient ab initio Electronic-Structure Calculation

The marriage of density functional theory (DFT) and deep learning methods has the potential to revolutionize modern computational materials science. Here we develop a deep neural network approach to represent DFT Hamiltonian (DeepH) of crystalline materials, aiming to bypass the computationally demanding self-consistent field iterations of DFT and substantially improve the efficiency of ab initio electronic-structure calculations. A general framework is proposed to deal with the large dimensionality and gauge (or rotation) covariance of DFT Hamiltonian matrix by virtue of locality and is realized by the message passing neural network for deep learning. High accuracy, high efficiency and good transferability of the DeepH method are generally demonstrated for various kinds of material systems and physical properties. The method provides a solution to the accuracy-efficiency dilemma of DFT and opens opportunities to explore large-scale material systems, as evidenced by a promising application to study twisted van der Waals materials.

preprint2022arXiv

Molecular conformer search with low-energy latent space

Identifying low-energy conformers with quantum mechanical accuracy for molecules with many degrees of freedom is challenging. In this work, we use the molecular dihedral angles as features and explore the possibility of performing molecular conformer search in a latent space with a generative model named variational auto-encoder (VAE). We bias the VAE towards low-energy molecular configurations to generate more informative data. In this way, we can effectively build a reliable energy model for the low-energy potential energy surface. After the energy model has been built, we extract local-minimum conformations and refine them with structure optimization. We have tested and benchmarked our low-energy latent-space (LOLS) structure search method on organic molecules with $5-9$ searching dimensions. Our results agree with previous studies.

preprint2021arXiv

Basic formulation and first-principles implementation of nonlinear magneto-optical effects

First-principles calculation of nonlinear magneto-optical effects has become an indispensable tool to reveal the geometric and topological nature of electronic states and to understand light-matter interactions. While intriguingly rich physics could emerge in magnetic materials, further methodological developments are required to deal with time-reversal symmetry breaking, due to the degeneracy and gauge problems caused by symmetry and the low-frequency divergence problem in the existing calculation formalism. Here we present a gauge-covariant and low-frequency convergent formalism for the first-principles computation. Remarkably, this formalism generally works for both non-magnetic and magnetic materials with or without band degeneracy. Reliability and capability of our method are demonstrated by studying example materials (i.e., bilayers of MnBi$_2$Te$_4$ and CrI$_3$) and comparing with published results. Moreover, an importance correction term that ensures gauge covariance of degenerate states is derived, whose influence on physical responses is systematically checked. Our method enables computation of nonlinear magneto-optical effects in magnetic materials and paves the way for exploring rich physics created by the interplay of light and magnetism.

preprint2021arXiv

Symmetry-adapted graph neural networks for constructing molecular dynamics force fields

Molecular dynamics is a powerful simulation tool to explore material properties. Most of the realistic material systems are too large to be simulated with first-principles molecular dynamics. Classical molecular dynamics has lower computational cost but requires accurate force fields to achieve chemical accuracy. In this work, we develop a symmetry-adapted graph neural networks framework, named molecular dynamics graph neural networks (MDGNN), to construct force fields automatically for molecular dynamics simulations for both molecules and crystals. This architecture consistently preserves the translation, rotation and permutation invariance in the simulations. We propose a new feature engineering method including higher order contributions and show that MDGNN accurately reproduces the results of both classical and first-principles molecular dynamics. We also demonstrate that force fields constructed by the model has good transferability. Therefore, MDGNN provides an efficient and promising option for molecular dynamics simulations of large scale systems with high accuracy.

preprint2020arXiv

Electronic states and magnetic response of MnBi2Te4 by scanning tunneling microscopy and spectroscopy

Exotic quantum phenomena have been demonstrated in recently discovered intrinsic magnetic topological insulator MnBi2Te4. At its two-dimensional limit, quantum anomalous Hall (QAH) effect and axion insulator state are observed in odd and even layers of MnBi2Te4, respectively. The measured band structures exhibit intriguing and complex properties. Here we employ low-temperature scanning tunneling microscopy to study its surface states and magnetic response. The quasiparticle interference patterns indicate that the electronic structures on the topmost layer of MnBi2Te4 is different from that of the expected out-of-plane A-type antiferromagnetic phase. The topological surface states may be embedded in deeper layers beneath the topmost surface. Such novel electronic structure presumably related to the modification of crystalline structure during sample cleaving and re-orientation of magnetic moment of Mn atoms near the surface. Mn dopants substituted at the Bi site on the second atomic layer are observed. The ratio of Mn/Bi substitutions is 5%. The electronic structures are fluctuating at atomic scale on the surface, which can affect the magnetism of MnBi2Te4. Our findings shed new lights on the magnetic property of MnBi2Te4 and thus the design of magnetic topological insulators.

preprint2020arXiv

Enhancement of superconductivity in organic-inorganic hybrid topological materials

Inducing or enhancing superconductivity in topological materials is an important route toward topological superconductivity. Reducing the thickness of transition metal dichalcogenides (e.g. WTe2 and MoTe2) has provided an important pathway to engineer superconductivity in topological matters; for instance, emergent superconductivity with Tc=0.82 K was observed in monolayer WTe2 which also hosts intriguing quantum spin Hall effect, although the bulk crystal is nonsuperconducting. However, such monolayer sample is difficult to obtain, unstable in air, and with extremely low Tc, which could pose a grand challenge for practical applications. Here we report an experimentally convenient approach to control the interlayer coupling to achieve tailored topological properties, enhanced superconductivity and good sample stability through organic cation intercalation of the Weyl semimetals MoTe2 and WTe2. The as-formed organic-inorganic hybrid crystals are weak topological insulators with enhanced Tc of 7.0 K for intercalated MoTe2 (0.25 K for pristine crystal) and 2.3 K for intercalated WTe2 (2.8 times compared to monolayer WTe2). Such organic-cationintercalation method can be readily applied to many other layered crystals, providing a new pathway for manipulating their electronic, topological and superconducting properties.

preprint2020arXiv

High-Temperature Quantum Anomalous Hall Insulators in Lithium-Decorated Iron-Based Superconductor Materials

Quantum anomalous Hall (QAH) insulator is the key material to study emergent topological quantum effects, but its ultralow working temperature limits experiments. Here, by first-principles calculations, we find a family of stable two-dimensional (2D) structures generated by lithium decoration of layered iron-based superconductor materials FeX (X = S, Se, Te), and predict room-temperature ferromagnetic semiconductors together with large-gap high-Chern-number QAH insulators in the 2D materials. The extremely robust ferromagnetic order is induced by the electron injection from Li to Fe and stabilized by strong ferromagnetic kinetic exchange in the 2D Fe layer. While in the absence of spin-orbit coupling (SOC), the ferromagnetism polarizes the system into a half Dirac semimetal state protected by mirror symmetry, the SOC effect results in a spontaneous breaking of mirror symmetry and introduces a Dirac mass term, which creates QAH states with sizable gaps (several tens of meV) and multiple chiral edge modes. We also find a 3D QAH insulator phase featured by macroscopic number of chiral conduction channels in bulk LiOH-LiFeX. The findings open new opportunities to realize novel QAH physics and applications at high temperatures.

preprint2020arXiv

Spin-Triplet Excitonic Insulator: The Case of Graphone

While various excitonic insulators have been studied in the literature, due to the perceived too-small spin splitting, spin-triplet excitonic insulator is rare. In two-dimensional systems such as a graphone, however, it is possible, as revealed by first-principles calculations coupled with Bethe-Salpeter equation. The critical temperature, given by an effective Hamiltonian, is 11.5 K. While detecting excitonic insulators is still a daunting challenge, the condensation of triplet excitons will result in spin superfluidity, which can be directly measured by a transport experiment. Nonlocal dielectric screening also leads to an unexpected phenomenon, namely, an indirect-to-direct transition crossover between single-particle band and exciton dispersion in graphone, which offers yet another test by experiment.

preprint2020arXiv

Tunable interlayer magnetism and band topology in van der Waals heterostructures of MnBi2Te4-family materials

Manipulating the interlayer magnetic coupling in van der Waals magnetic materials and heterostructures is the key to tailoring their magnetic and electronic properties for various electronic applications and fundamental studies in condensed matter physics. By utilizing the MnBi2Te4-family compounds and their heterostructures as a model system, we systematically studied the dependence of the sign and strength of interlayer magnetic coupling on constituent elements by using first-principles calculations. It was found that the coupling is a long-range superexchange interaction mediated by the chains of p orbitals between the magnetic atoms of neighboring septuple-layers. The interlayer exchange is always antiferromagnetic in the pure compounds, but can be tuned to ferromagnetic in some combinations of heterostructures, dictated by d orbital occupations. Strong interlayer magnetic coupling can be realized if the medial p electrons are delocalized and the d bands of magnetic atoms are near the Fermi level. The knowledge on the interlayer coupling mechanism enables us to engineer magnetic and topological properties of MnBi2Te4-family materials as well as many other insulating van der Waals magnetic materials and heterostructures.

preprint2019arXiv

Berry Curvature Engineering by Gating Two-Dimensional Antiferromagnets

Recent advances in tuning electronic, magnetic, and topological properties of two-dimensional (2D) magnets have opened a new frontier in the study of quantum physics and promised exciting possibilities for future quantum technologies. In this study, we find that the dual-gate technology can well tune the electronic and topological properties of antiferromagnetic (AFM) even septuple-layer (SL) MnBi$_2$Te$_4$ thin films. Under an out-of-plane electric field that breaks $\mathcal{PT}$ symmetry, the Berry curvature of the thin film could be engineered efficiently, resulting in a huge change of anomalous Hall (AH) signal. Beyond the critical electric field, the double-SL MnBi$_2$Te$_4$ thin film becomes a Chern insulator with a high Chern number of 3. We further demonstrate that such 2D material can be used as an AFM switch via electric-field control of the AH signal. These discoveries inspire the design of low-power memory prototype for future AFM spintronic applications.

preprint2019arXiv

Giant Enhancement of Solid Solubility in Monolayer BNC Alloys by Selective Orbital Coupling

Solid solubility (SS) is one of the most important features of alloys, which is usually difficult to be largely tuned in the entire alloy concentrations by external approaches. Some alloys that were supposed to have promising physical properties could turn out to be much less useful because of their poor SS, e.g., the case for monolayer BNC [(BN)1-x(C2)x] alloys. Until now, an effective approach on significantly enhancing SS of (BN)1-x(C2)x in the entire x is still lacking. In this article, a novel mechanism of selective orbital coupling between high energy wrong-bond states and surface states mediated by the specific substrate has been proposed to stabilize the wrong-bonds and in turn significantly enhance the SS of (BN)1-x(C2)x alloys. Surprisingly, we demonstrate that five ordered alloys, exhibiting variable direct quasi-particle bandgaps from 1.35 to 3.99 eV, can spontaneously be formed at different x when (BN)1-x(C2)x is grown on hcp-phase Cr. Interestingly, the optical transitions around the band edges in these ordered alloys, accompanied by largely tunable exciton binding energies of ~1 eV at different x, are significantly strong due to their unique band structures. Importantly, the disordered (BN)1-x(C2)x alloys, exhibiting fully tunable bandgaps from 0 to ~6 eV in the entire x, can be formed on Cr substrate at the miscibility temperature of ~1200 K, which is greatly reduced compared to that of 4500~5600 K in free-standing form or on other substrates. Our discovery not only may resolve the long-standing SS problem of BNC alloys, but also could significantly extend the applications of BNC alloys for various optoelectronic applications.

preprint2019arXiv

Type-II Ising Pairing in Few-Layer Stanene

Spin-orbit coupling has proven indispensable in realizing topological materials and more recently Ising pairing in two-dimensional superconductors. This pairing mechanism relies on inversion symmetry breaking and sustains anomalously large in-plane polarizing magnetic fields whose upper limit is expected to diverge at low temperatures, although experimental demonstration of this has remained elusive due to the required fields. In this work, the recently discovered superconductor few-layer stanene, i.e. epitaxially strained $α$-Sn, is shown to exhibit a new type of Ising pairing between carriers residing in bands with different orbital indices near the $Γ$-point. The bands are split as a result of spin-orbit locking without the participation of inversion symmetry breaking. The in-plane upper critical field is strongly enhanced at ultra-low temperature and reveals the sought for upturn.

preprint2016arXiv

Electronic properties of SnTe-class topological crystalline insulator materials

The rise of topological insulators in recent years has broken new ground both in the conceptual cognition of condensed matter physics and the promising revolution of the electronic devices. It also stimulates the explorations of more topological states of matter. Topological crystalline insulator is a new topological phase, which combines the electronic topology and crystal symmetry together. In this article, we review the recent progress in the studies of SnTe-class topological crystalline insulator materials. Starting from the topological identifications in the aspects of the bulk topology, surface states calculations and experimental observations, we present the electronic properties of topological crystalline insulators under various perturbations, including native defect, chemical doping, strain, and thickness-dependent confinement effects, and then discuss their unique quantum transport properties, such as valley-selective filtering and helicity-resolved functionalities for Dirac fermions. The rich properties and high tunability make SnTe-class materials promising candidates for novel quantum devices.

preprint2016arXiv

Energy gaps of atomically precise armchair graphene nanoribbons

Graphene nanoribbons (GNRs) are one-dimensional (1D) structures that exhibit a rich variety of electronic properties1-17. Therefore, they are predicted to be the building blocks in next-generation nanoelectronic devices. Theoretically, it has been demonstrated that armchair GNRs can be divided into three families, i.e., Na = 3p, Na = 3p + 1, and Na = 3p + 2 (here Na is the number of dimer lines across the ribbon width and p is an integer), according to their electronic structures, and the energy gaps for the three families are quite different even with the same p1,3-6. However, a systematic experimental verification of this fundamental prediction is still lacking, owing to very limited atomic-level control of the width of the armchair GNRs investigated7,9,10,13,17. Here, we studied electronic structures of the armchair GNRs with atomically well-defined widths ranging from Na = 6 to Na = 26 by using scanning tunnelling microscope (STM). Our result demonstrated explicitly that all the studied armchair GNRs exhibit semiconducting gaps due to quantum confinement and, more importantly, the observed gaps as a function of Na are well grouped into the three categories, as predicted by density-functional theory calculations3. Such a result indicated that we can tune the electronic properties of the armchair GNRs dramatically by simply adding or cutting one carbon dimer line along the ribbon width.

preprint2016arXiv

Experimental observation of topological Fermi arcs in type-II Weyl semimetal MoTe2

Weyl semimetal is a new quantum state of matter [1-12] hosting the condensed matter physics counterpart of relativisticWeyl fermion [13] originally introduced in high energy physics. The Weyl semimetal realized in the TaAs class features multiple Fermi arcs arising from topological surface states [10, 11, 14-16] and exhibits novel quantum phenomena, e.g., chiral anomaly induced negative mag-netoresistance [17-19] and possibly emergent supersymmetry [20]. Recently it was proposed theoretically that a new type (type-II) of Weyl fermion [21], which does not have counterpart in high energy physics due to the breaking of Lorentz invariance, can emerge as topologically-protected touching between electron and hole pockets. Here, we report direct spectroscopic evidence of topological Fermi arcs in the predicted type-II Weyl semimetal MoTe2 [22-24]. The topological surface states are confirmed by directly observing the surface states using bulk-and surface-sensitive angle-resolved photoemission spectroscopy (ARPES), and the quasi-particle interference (QPI) pattern between the two putative Fermi arcs in scanning tunneling microscopy (STM). Our work establishes MoTe2 as the first experimental realization of type-II Weyl semimetal, and opens up new opportunities for probing novel phenomena such as exotic magneto-transport [21] in type-II Weyl semimetals.

preprint2016arXiv

Heavy Dirac fermions in a graphene/topological insulator hetero-junction

The low energy physics of both graphene and surface states of three-dimensional topological insulators is described by gapless Dirac fermions with linear dispersion. In this work, we predict the emergence of a "heavy" Dirac fermion in a graphene/topological insulator hetero-junction, where the linear term almost vanishes and the corresponding energy dispersion becomes highly non-linear. By combining {\it ab initio} calculations and an effective low-energy model, we show explicitly how strong hybridization between Dirac fermions in graphene and the surface states of topological insulators can reduce the Fermi velocity of Dirac fermions. Due to the negligible linear term, interaction effects will be greatly enhanced and can drive "heavy" Dirac fermion states into the half quantum Hall state with non-zero Hall conductance.

preprint2016arXiv

Stable Dirac semi-metal in the allotrope of IV elements

Three dimensional topological Dirac semi-metals represent a novel state of quantum matter with exotic electronic properties, in which a pair of Dirac points with the linear dispersion along all momentum directions exist in the bulk. Herein, by using the first principles calculations, we discover a new metastable allotrope of Ge and Sn in the staggered layered dumbbell structure, named as germancite and stancite, to be Dirac semi-metals with a pair of Dirac points on its rotation axis. On the surface parallel to the rotation axis, a pair of topologically non-trivial Fermi arcs are observed and a Lifshitz transition is found by tuning the Fermi level. Furthermore, the quantum thin film of germancite is found to be an intrinsic quantum spin Hall insulator. These discoveries suggest novel physical properties and future applications of the new metastable allotrope of Ge and Sn.

preprint2016arXiv

Topological nodal-line semimetals in alkaline-earth stannides, germanides and silicides

Based on first-principles calculations and an effective Hamiltonian analysis, we systematically investigate the electronic and topological properties of alkaline-earth compounds $AX_2$ ($A$=Ca, Sr, Ba; $X$=Si, Ge, Sn). Taking BaSn$_2$ as an example, we find that when spin-orbit coupling is ignored, these materials are three-dimensional topological nodal-line semimetals characterized by a snake-like nodal loop in three-dimensional momentum space. Drumhead-like surface states emerge either inside or outside the loop circle on the (001) surface depending on surface termination, while complicated double-drumhead-like surface states appear on the (010) surface. When spin-orbit coupling is included, the nodal line is gapped and the system becomes a topological insulator with Z$_2$ topological invariants (1;001). Since spin-orbit coupling effects are weak in light elements, the nodal-line semimetal phase is expected to be achievable in some alkaline-earth germanides and silicides.

preprint2015arXiv

Band engineering of Dirac surface states in topological insulators-based van der Waals heterostructures

The existence of gapless Dirac surface band of a three dimensional (3D) topological insulator (TI) is guaranteed by the non-trivial topological character of the bulk band, yet the surface band dispersion is mainly determined by the environment near the surface. In this Letter, through in-situ angle-resolved photoemission spectroscopy (ARPES) and the first-principles calculation on 3D TI-based van der Waals heterostructures, we demonstrate that one can engineer the surface band structures of 3D TIs by surface modifications without destroying their topological non-trivial property. The result provides an accessible method to independently control the surface and bulk electronic structures of 3D TIs, and sheds lights in designing artificial topological materials for electronic and spintronic purposes.

preprint2015arXiv

Converting normal insulators into topological insulators via tuning orbital levels

Tuning the spin-orbit coupling strength via foreign element doping and/or modifying bonding strength via strain engineering are the major routes to convert normal insulators to topological insulators. We here propose an alternative strategy to realize topological phase transition by tuning the orbital level. Following this strategy, our first-principles calculations demonstrate that a topological phase transition in some cubic perovskite-type compounds CsGeBr$_3$ and CsSnBr$_3$ could be facilitated by carbon substitutional doping. Such unique topological phase transition predominantly results from the lower orbital energy of the carbon dopant, which can pull down the conduction bands and even induce band inversion. Beyond conventional approaches, our finding of tuning the orbital level may greatly expand the range of topologically nontrivial materials.

preprint2015arXiv

Creating and probing wide-bandgap nanoribbon-like structures in a continuous metallic graphene sheet

The light-like dispersion of graphene monolayer results in many novel electronic properties in it1, however, this gapless feature also limits the applications of graphene monolayer in digital electronics2. A rare working solution to generate a moderate bandgap in graphene monolayer is to cut it into one-dimensional (1D) nanometre-wide ribbons3-13. Here we show that a wide bandgap can be created in a unique 1D strained structure, i.e., graphene-nanoribbon-like (GNR-like) structure, of a continuous graphene sheet via strong interaction between graphene and the metal substrate, instead of cutting graphene monolayer. The GNR-like structures with width of only a few nanometers are observed in a continuous graphene sheet grown on Rh foil by using thermal strain engineering. Spatially-resolved scanning tunnelling spectroscopy revealed bandgap opening of a few hundreds meV in the GNR-like structure in an otherwise continuous metallic graphene sheet, directly demonstrating the realization of a metallic-semiconducting-metallic junction entirely in a graphene monolayer. We also show that it is possible to tailor the structure and electronic properties of the GNR-like structure by using scanning tunnelling microscope.

preprint2015arXiv

Effects of ferroelectric polarization on surface phase diagram: an evolutionary algorithm study of the BaTiO$_{3}$(001) surface

We have constructed the surface phase diagram of the BaTiO$_{3}$(001) surface by employing an evolutionary algorithm for surface structure prediction, where the ferroelectric polarization is included as a degree of freedom. Among over 1000 candidate structures explored, a surface reconstruction of (2$\times$1)-TiO is discovered to be thermodynamically stable and have the \emph{p2mm} plane group symmetry as observed experimentally. We find that the influence of ferroelectric polarization on the surface free energy can be either negligibly small or sizably large (over 1 eV per ($2 \times 1$) supercell), depending strongly on the surface structure and resulting in a significant distinction of surface phase diagram with varying ferroelectric polarization. It is therefore feasible to control the surface stability by applying an external electric field. Our results may have important implications in understanding the surface reconstruction of ferroelectric materials and tuning surface properties.

preprint2015arXiv

Emergence of a Chern-insulating state from a semi-Dirac dispersion

A Chern insulator (quantum anomalous Hall insulator) phase is demonstrated to exist in a typical semi-Dirac system, the TiO2/VO2 heterostructure. By combining first-principles calculations with Wannier-based tight-binding model, we calculate the Berry curvature distribution, finding a Chern number of -2 for the valence bands, and demonstrate the existence of gapless chiral edge states, ensuring quantization of the Hall conductivity to 2e^2/h. A new semi-Dirac model, where each semi-Dirac cone is formed by merging three conventional Dirac points, is proposed to reveal how the nontrivial topology with finite Chern number is compatible with a semi-Dirac electronic spectrum.

preprint2015arXiv

Graphene/g-C2N bilayer:gap opening, enhanced visible light response and electrical field tuning band structure

Opening up a band gap of the graphene and finding a suitable substrate are two challenges for constituting the nano-electronic equipment. A new two-dimensional layered crystal g-C2N (Nat. Commun. 2015, 6, 1--7) with novel electronic and optical properties can be effectively synthesized via a wet-chemical reaction. And g-C2N can be used as a suitable substrate to open the band gap of graphene as much as 0.239 eV, which is large enough for the band gap opening at room temperature. The physics behind the band gap opening is that g-C2N substrate can produce the inhomogeneous electrostatic potential over the graphene layer. The imposition of external electrical field can tune the band gap of the hybrid of graphene/g-C2N effectively from the semiconductors to the metal. The hybrid graphene/g-C2N displays an enhanced optical activity compared with the pure g-C2N monolayer.

preprint2015arXiv

High-$T_c$ $d^{0}$ ferromagnetism in a doped Mott insulator: the case of hydrogenated epitaxial graphene on SiC(0001)

We show that the $d^{0}$ ferromagnetism with high Curie temperature ($T_c$) can be achieved in the electron doped hydrogenated epitaxial graphene on some certain SiC substrates through first-principles calculations. The pristine systems are found to be a Mott insulator independent of SiC polytypes (2$H$, 4$H$ or 6$H$) which, however, play a significant role in the modulation of magnetic interaction. Carrier doping enhances the ferromagnetic coupling due to the double exchange mechanism and thus realizes the phase transition from antiferromagnetism to ferromagnetism. A $T_c$ of around 400 K is predicted on the 2$H$-SiC. We employ a non-degenerate Hubbard model to demonstrate how the SiC affects the interfacial magnetism in intra-atomic Coulomb repulsion and intersite hopping interactions.

preprint2015arXiv

Structural phase transition and electronic structure evolution in Ir1-xPtxTe2 studied by scanning tunneling microscopy

The IrTe2 transition metal dichalcogenide undergoes a series of structural and electronic phase transitions when doped with Pt. The nature of each phase and the mechanism of the phase transitions have attracted much attention. In this paper, we report scanning tunneling microscopy and spectroscopy studies of Pt doped IrTe2 with varied Pt contents. In pure IrTe2, we find that the ground state has a 1/6 superstructure, and the electronic structure is inconsistent with Fermi surface nesting induced charge density wave order. Upon Pt doping, the crystal structure changes to a 1/5 superstructure and then to a quasi-periodic hexagonal phase. First principles calculations show that the superstructures and electronic structures are determined by the global chemical strain and local impurity states that can be tuned systematically by Pt doping.

preprint2015arXiv

Theory of the Dirac Half Metal and Quantum Anomalous Hall Effect in Mn Intercalated Epitaxial Graphene

The prospect of a Dirac half metal, a material which is characterized by a bandstructure with a gap in one spin channel but a Dirac cone in the other, is of both fundamental interest and a natural candidate for use in spin-polarized current applications. However, while the possibility of such a material has been reported based on model calculations[H. Ishizuka and Y. Motome, Phys. Rev. Lett. 109, 237207 (2012)], it remains unclear what material system might realize such an exotic state. Using first-principles calculations, we show that the experimentally accessible Mn intercalated epitaxial graphene on SiC(0001) transits to a Dirac half metal when the coverage is > 1/3 monolayer. This transition results from an orbital-selective breaking of quasi-2D inversion symmetry, leading to symmetry breaking in a single spin channel which is robust against randomness in the distribution of Mn intercalates. Furthermore, the inclusion of spin-orbit interaction naturally drives the system into the quantum anomalous Hall (QAH) state. Our results thus not only demonstrate the practicality of realizing the Dirac half metal beyond a toy model but also open up a new avenue to the realization of the QAH effect.

preprint2014arXiv

A general group theoretical method to unfold band structures and its application

We present a general method to unfold energy bands of supercell calculations to primitive Brillouin zone using group theoretical techniques, where an isomorphic factor group is introduced to connect the primitive translation group with the supercell translation group via a direct product. Originating from the translation group symmetry, our method gives an uniform description of unfolding approaches based on various basis sets, and therefore, should be easy to implement in both tight-binding model and existing ab initio code packages using different basis sets. This makes the method applicable to a variety of problems involving the use of supercells, such as defects, disorder, and interfacial reconstructions. As a realistic example, we calculate electronic properties of an monolayer FeSe on SrTiO$_3$ in checkerboard and collinear antiferromagnetic spin configurations, illustrating the potential of our method.

preprint2014arXiv

Chemical potential dependent gap-opening at the Dirac surface states of Bi2Se3 induced by aggregated substitutional Cr atoms

With angle-resolved photoemission spectroscopy, gap-opening is resolved at up to room temperature in the Dirac surface states of molecular beam epitaxy grown Cr-doped Bi2Se3 topological insulator films, which however show no long-range ferromagnetic order down to 1.5 K. The gap size is found decreasing with increasing electron doping level. Scanning tunneling microscopy and first principles calculations demonstrate that substitutional Cr atoms aggregate into superparamagnetic multimers in Bi2Se3 matrix, which contribute to the observed chemical potential dependent gap-opening in the Dirac surface states without long-range ferromagnetic order.

preprint2014arXiv

Experimental observation of Dirac-like surface states and topological phase transition in Pb$_{1-x}$Sn$_x$Te(111) films

The surface of a topological crystalline insulator (TCI) carries an even number of Dirac cones protected by crystalline symmetry. We epitaxially grew high quality Pb$_{1-x}$Sn$_x$Te(111) films and investigated the TCI phase by in-situ angle-resolved photoemission spectroscopy. Pb$_{1-x}$Sn$_x$Te(111) films undergo a topological phase transition from trivial insulator to TCI via increasing the Sn/Pb ratio, accompanied by a crossover from n-type to p-type doping. In addition, a hybridization gap is opened in the surface states when the thickness of film is reduced to the two-dimensional limit. The work demonstrates an approach to manipulating the topological properties of TCI, which is of importance for future fundamental research and applications based on TCI.

preprint2014arXiv

Functionalized Germanene as a Prototype of Large-Gap Two-Dimensional Topological Insulators

We propose new two-dimensional (2D) topological insulators (TIs) in functionalized germanenes (GeX, X=H, F, Cl, Br or I) using first-principles calculations. We find GeI is a 2D TI with a bulk gap of about 0.3 eV, while GeH, GeF, GeCl and GeBr can be transformed into TIs with sizeable gaps under achievable tensile strains. A unique mechanism is revealed to be responsible for large topologically-nontrivial gap obtained: owing to the functionalization, the $σ$ orbitals with stronger spin-orbit coupling (SOC) dominate the states around the Fermi level, instead of original $π$ orbitals with weaker SOC; thereinto, the coupling of the $p_{xy}$ orbitals of Ge and heavy halogens in forming the $σ$ orbitals also plays a key role in the further enlargement of the gaps in halogenated germanenes. Our results suggest a realistic possibility for the utilization of topological effects at room temperature.

preprint2014arXiv

Prediction of a Stable Post-Post-Perovskite Structure from First Principles

A novel stable crystallographic structure is discovered in a variety of ABO3, ABF3 and A2O3 compounds (including materials of geological relevance, prototypes of multiferroics, exhibiting strong spin-orbit effects, etc...), via the use of first principles. This novel structure appears under hydrostatic pressure, and is the first "post-post-perovskite" phase to be found. It provides a successful solution to experimental puzzles in important systems, and is characterized by one-dimensional chains linked by group of two via edge-sharing oxygen/fluorine octahedra. Such unprecedented organization automatically results in anisotropic elastic properties and new magnetic arrangements. Depending on the system of choice, this post-post-perovskite structure also possesses electronic band gaps ranging from zero to ~ 10 eV being direct or indirect in nature, which emphasizes its "universality" and its potential to have striking, e.g., electrical or transport phenomena.

preprint2014arXiv

Quantum Unfolding: A program for unfolding electronic energy bands of materials

We present Quantum Unfolding, a Fortran90 program for unfolding first-principles electronic energy bands. It unfolds energy bands accurately by handling the Fourier components of Bloch wavefunctions, which are reconstructed from Wannier functions from Wannier90. Due to the wide application of Wannier90 package and the possibility of focusing only on the most important energy bands, the present code works very conveniently.

preprint2014arXiv

Replying "the comment on Interlayer interactions in graphites"

For the moment, there is no exact description of van der Waals (vdW) interactions. ACFD-RPA \cite{Gould1} is expected to better describe vdW bonding, but it is not exact. The PBE/DFT-D2 method is less satisfactory, however, its results are in good agreement with experimental data. Although our fitting technique may weaken (not neglect) the vdW interactions and produce interlayer potentials with weakened vdW, the obtained interlayer potentials reproduce energetics of graphite near the equilibrium interlayer distance very well, as shown in Ref. \cite{Chen}. If having inputs which fully include vdW interactions and having better fitting functions, we believe that interlayer potentials can also fully include vdW interactions in graphite system.

preprint2014arXiv

Stable two-dimensional dumbbell stanene: a quantum spin Hall insulator

We predict from first-principles calculations a novel structure of stanene with dumbbell units (DB), and show that it is a two-dimensional topological insulator with inverted band gap which can be tuned by compressive strain. Furthermore, we propose that the boron nitride sheet and reconstructed ($2\times2$) InSb(111) surfaces are ideal substrates for the experimental realization of DB stanene, maintaining its non-trivial topology. Combined with standard semiconductor technologies, such as magnetic doping and electrical gating, the quantum anomalous Hall effect, Chern half metallicity and topological superconductivity can be realized in DB stanene on those substrates. These properties make the two-dimensional supported stanene a good platform for the study of new quantum spin Hall insulator as well as other exotic quantum states of matter.

preprint2014arXiv

Thermal and Thermoelectric Properties of Graphene

The subject of thermal transport at the mesoscopic scale and in low-dimensional systems is interesting for both fundamental research and practical applications. As the first example of truly two-dimensional materials, graphene has exceptionally high thermal conductivity, and thus provides an ideal platform for the research. Here we review recent studies on thermal and thermoelectric properties of graphene, with an emphasis on experimental progresses. A general physical picture based on the Landauer transport formalism is introduced to understand underlying mechanisms. We show that the superior thermal conductivity of graphene is contributed not only by large ballistic thermal conductance but also by very long phonon mean free path (MFP). The long phonon MFP, explained by the low-dimensional nature and high sample purity of graphene, results in important isotope effects and size effects on thermal conduction. In terms of various scattering mechanisms in graphene, several approaches are suggested to control thermal conductivity. Among them, introducing rough boundaries and weakly-coupled interfaces are promising ways to suppress thermal conduction effectively. We also discuss the Seebeck effect of graphene. Graphene itself might not be a good thermoelectric material. However, the concepts developed by graphene research might be applied to improve thermoelectric performance of other materials.

preprint2014arXiv

Tunable Thermal Conduction in Graphane Nanoribbons

Graphane and graphene are both two-dimensional materials but of different bonding configurations, which can result in distinct thermal conduction properties. We simulate thermal conduction in graphane nanoribbons (GANRs) using the nonequilibrium Green's function method. It is found that GANRs have lower ballistic thermal conductance and stronger thermal conductance anisotropy than the graphene counterparts. Furthermore, hydrogen vacancies of GANRs considerably suppress thermal conduction, accompanied by enhanced thermal conductance anisotropy. The tunable thermal conduction, realized by controlling the width, edge shape and hydrogen vacancy concentration of GANRs, could be useful for thermal management and thermoelectric applications.

preprint2013arXiv

Large-gap quantum spin Hall insulators in tin films

The search of large-gap quantum spin Hall (QSH) insulators and effective approaches to tune QSH states is important for both fundamental and practical interests. Based on first-principles calculations we find two-dimensional tin films are QSH insulators with sizable bulk gaps of 0.3 eV, sufficiently large for practical applications at room temperature. These QSH states can be effectively tuned by chemical functionalization and by external strain. The mechanism for the QSH effect in this system is band inversion at the Γpoint, similar to the case of HgTe quantum well. With surface doping of magnetic elements, the quantum anomalous Hall effect could also be realized.

preprint2013arXiv

On the electronic structure of silicene on Ag substrate: strong hybridization effects

The electronic structure of the recently synthesised (3x3) reconstructed silicene on (4x4) Ag(111) is investigated by first-principles calculations. New states emerge due to the strong hybridization between silicene and Ag. Analyzing the nature and composition of these hybridized states, we show that i) it is possible to clearly distinguish them from states coming from the Dirac cone of free-standing silicene or from the sp-bands of bulk Ag and ii) assign their contribution to the description of the linearly dispersing band observed in photoemission. Furthermore, we show that silicene atoms contribute to the Fermi level, which leads to similar STM patterns as observed below or above the Fermi level. Our findings are crucial for the proper interpretation of experimental observations.

preprint2013arXiv

Spin-filtered Edge States with an Electrically Tunable Gap in a Two-Dimensional Topological Crystalline Insulator

Three-dimensional topological crystalline insulators were recently predicted and observed in the SnTe class of IV-VI semiconductors, which host metallic surface states protected by crystal symmetries. In this work, we study thin films of these materials and expose their potential for device applications. We demonstrate that thin films of SnTe and Pb(1-x)Sn(x)Se(Te) grown along the (001) direction are topologically nontrivial in a wide range of film thickness and carry conducting spin-filtered edge states that are protected by the (001) mirror symmetry via a topological invariant. Application of an electric field perpendicular to the film will break the mirror symmetry and generate a band gap in these edge states. This functionality motivates us to propose a novel topological transistor device, in which charge and spin transport are maximally entangled and simultaneously controlled by an electric field. The high on/off operation speed and coupling of spin and charge in such a device may lead to electronic and spintronic applications for topological crystalline insulators.

preprint2013arXiv

Superconducting Graphene: the conspiracy of doping and strain

Graphene has exhibited a wealth of fascinating properties, but is also known not to be a superconductor. Remarkably, we show that graphene can be made a conventional Bardeen-Cooper-Schrieffer superconductor by the combined effect of charge doping and tensile strain. While the effect of doping is obvious to enlarge Fermi surface, the effect of strain is profound to greatly increase the electron-phonon coupling. At the experimental accessible doping (4E+14cm-2) and strain (~16%) levels, the superconducting critical temperature Tc reaches as high as ~30 K, the highest for a single-element material above the liquid hydrogen temperature. This significantly makes graphene a commercially viable superconductor.

preprint2013arXiv

Surface States of Topological Crystalline Insulators in IV-VI Semiconductors

Topological crystalline insulators (TCI) are new topological phases of matter protected by crystal symmetry of solids. Recently, the first realization of TCI has been predicted and observed in IV-VI semiconductor SnTe and related alloys Pb_{1-x}Sn_{x}(Te, Se). By combining k.p theory and band structure calculation, we present a unified approach to study topological surface states on various crystal surfaces of TCI in IV-VI semiconductors. We explicitly derive k.p Hamiltonian for topological surface states from electronic structure of the bulk, thereby providing a microscopic understanding of bulk-boundary correspondence in TCI. Depending on the surface orientation, we find two types of surface states with qualitatively different properties. In particular, we predict that (111) surface states consist of four Dirac cones centered at time-reversal-invariant momenta Γ and M, while (110) surface states consist of Dirac cones at non-time-reversal-invariant momenta, similar to (001). Moreover, both (001) and (110) surface states exhibit a Lifshitz transition as a function of Fermi energy, which is accompanied by a Van-Hove singularity in density of states arising from saddle points in the band structure.

preprint2013arXiv

Weak Topological Insulators in PbTe/SnTe Superlattices

It is desirable to realize topological phases in artificial structures by engineering electronic band structures. In this paper, we investigate $(PbTe)_m(SnTe)_{2n-m}$ superlattices along [001] direction and find a robust weak topological insulator phase for a large variety of layer numbers m and 2n-m. We confirm this topologically non-trivial phase by calculating Z2 topological invariants and topological surface states based on the first-principles calculations. We show that the folding of Brillouin zone due to the superlattice structure plays an essential role in inducing topologically non-trivial phases in this system. This mechanism can be generalized to other systems in which band inversion occurs at multiple momenta, and gives us a brand-new way to engineer topological materials in artificial structures.

preprint2013arXiv

Weak topological insulators induced by the inter-layer coupling: A first-principles study of stacked Bi2TeI

Based on first-principles calculations, we predict Bi2TeI, a stoichiometric compound synthesized, to be a weak topological insulator (TI) in layered subvalent bismuth telluroiodides. Within a bulk energy gap of 80 meV, two Dirac-cone-like topological surface states exist on the side surface perpendicular to BiTeI layer plane. These Dirac cones are relatively isotropic due to the strong inter-layer coupling, distinguished from those of previously reported weak TI candidates. Moreover, with chemically stable cladding layers, the BiTeI-Bi2-BiTeI sandwiched structure is a robust quantum spin Hall system, which can be obtained by simply cleaving the bulk Bi2TeI.

preprint2012arXiv

Dirac Fermion in Strongly-Bound Graphene Systems

It is highly desirable to integrate graphene into existing semiconductor technology, where the combined system is thermodynamically stable yet maintain a Dirac cone at the Fermi level. Firstprinciples calculations reveal that a certain transition metal (TM) intercalated graphene/SiC(0001), such as the strongly-bound graphene/intercalated-Mn/SiC, could be such a system. Different from free-standing graphene, the hybridization between graphene and Mn/SiC leads to the formation of a dispersive Dirac cone of primarily TM d characters. The corresponding Dirac spectrum is still isotropic, and the transport behavior is nearly identical to that of free-standing graphene for a bias as large as 0.6 V, except that the Fermi velocity is half that of graphene. A simple model Hamiltonian is developed to qualitatively account for the physics of the transfer of the Dirac cone from a dispersive system (e.g., graphene) to an originally non-dispersive system (e.g., TM).

preprint2012arXiv

Electronic Strengthening of Graphene by Charge Doping

Graphene is known as the strongest 2D material in nature, yet we show that moderate charge doping of either electrons or holes can further enhance its ideal strength by up to ~17%, based on first principles calculations. This unusual electronic enhancement, versus conventional structural enhancement, of material's strength is achieved by an intriguing physical mechanism of charge doping counteracting on strain induced enhancement of Kohn anomaly, which leads to an overall stiffening of zone boundary K1 phonon mode whose softening under strain is responsible for graphene failure. Electrons and holes work in the same way due to the high electron-hole symmetry around the Dirac point of graphene, while over doping may weaken the graphene by softening other phonon modes. Our findings uncover another fascinating property of graphene with broad implications in graphene-based electromechanical devices.

preprint2012arXiv

Topo-electronic transitions in Sb(111) nanofilm: the interplay between quantum confinement and surface effect

When the dimension of a solid structure is reduced, there will be two emerging effects, quantum confinement and surface effect, which dominate at nanoscale. Based on first-principles calculations, we demonstrate that due to an intriguing interplay between these two dominating effects, the topological and electronic (topo-electronic) properties of Sb (111) nanofilms undergo a series of transitions as a function of the reducing film thickness: transforming from a topological semimetal to a topological insulator at 7.8 nm (22 bilayer), then to a quantum spin hall (QSH) phase at 2.7 nm (8 bilayer), and finally to a normal (topological trivial) semiconductor at 1.0 nm (3 bilayer). Our theoretical findings for the first time identify the existence of the QSH in the Sb (111) nanofilms within a narrow range of thickness and suggest that the Sb (111) nanofilms provide an ideal test bed for experimental study of topo-electronic phase transitions.

preprint2012arXiv

Topological Crystalline Insulators in the SnTe Material Class

Topological crystalline insulators are new states of matter in which the topological nature of electronic structures arises from crystal symmetries. Here we predict the first material realization of topological crystalline insulator in the semiconductor SnTe, by identifying its nonzero topological index. We predict that as a manifestation of this nontrivial topology, SnTe has metallic surface states with an even number of Dirac cones on high-symmetry crystal surfaces such as {001}, {110} and {111}. These surface states form a new type of high-mobility chiral electron gas, which is robust against disorder and topologically protected by reflection symmetry of the crystal with respect to {110} mirror plane. Breaking this mirror symmetry via elastic strain engineering or applying an in-plane magnetic field can open up a continuously tunable band gap on the surface, which may lead to wide-ranging applications in thermoelectrics, infrared detection, and tunable electronics. Closely related semiconductors PbTe and PbSe also become topological crystalline insulators after band inversion by pressure, strain and alloying.

preprint2011arXiv

Activated dissociation of O2 on Pb(111) surfaces by Pb adatoms

We investigate the dissociation of O2 on Pb(111) surface using first-principles calculations. It is found that in a practical high-vacuum environment, the adsorption of molecular O2 takes place on clean Pb surfaces only at low temperatures such as 100 K, but the O2 easily desorbs at (elevated) room temperatures. It is further found that the Pb adatoms enhance the molecular adsorption and activate the adsorbed O2 to dissociate during subsequent room-temperature annealing. Our theory explains the observation of a two-step oxidation process on the Pb surfaces by the unique role of Pb adatoms.

preprint2011arXiv

Klein-tunneling-enhanced directional coupler for Dirac electron wave in graphene

Using the coupled-mode theory in guided-wave optics and electronics, we explore a directional coupling structure composed of two parallel waveguides electrostatically induced by the split-gate technique in bulk graphene. Our results show that Klein tunneling can greatly enhance the coupling strength of the structure. By adjusting a gate voltage, the probability density of Dirac electron wave function initially in one waveguide can be completely transferred into the other waveguide within several hundred nanometers. Our findings could not only lead to functional coherent coupling devices for quantum-based electronic signal processing and on-chip device integration in graphene, but also shrink the size of the devices to facilitate the fabrication of graphene-based large-scale integrated logic circuits.

preprint2010arXiv

Ab initio study of beryllium-decorated fullerenes for hydrogen storage

We have found that a beryllium (Be) atom on nanostructured materials with H2 molecules generates a Kubas-like dihydrogen complex [H. Lee et al. arXiv:1002.2247v1 (2010)]. Here, we investigate the feasibility of Be-decorated fullerenes for hydrogen storage using ab initio calculations. We find that the aggregation of Be atoms on pristine fullerenes is energetically preferred, resulting in the dissociation of the dihydrogen. In contrast, for boron (B)-doped fullerenes, Be atoms prefer to be individually attached to B sites of the fullerenes, and a maximum of one H2 molecule binds to each Be atom in a form of dihydrogen with a binding energy of ~0.3 eV. Our results show that individual dispersed Be-decorated B-doped fullerenes can serve as a room-temperature hydrogen storage medium.

preprint2010arXiv

Beryllium-dihydrogen complexes on nanostructures

Using the pseudopotential density functional method, we find that a Be atom on a nanostructure with H2 molecules forms a Be-dihydrogen complex through the hybridization of the Be s or p orbits with the H2 sigma orbits and the binding energy of the H2 molecules is in the range of ~0.3 - 0.8 eV/H2. We also study Be-dihydrogen complexes on various nanostructures and demonstrate the feasibility of the application of the complexes to a hydrogen storage medium that operates near room temperature and ambient pressure.

preprint2010arXiv

Edge Stability of BN sheets and Its Application for Designing Hybrid BNC Structures

First-principles investigations on the edge energies and edge stresses of single-layer hexagonal boron-nitride (BN) are presented. The armchair edges of BN nanoribbons (BNNRs) are more stable in energy than zigzag ones. Armchair BNNRs are under compressive edge stress while zigzag BNNRs are under tensile edge stress. The intrinsic spin-polarization and edge saturation play important roles in modulating the edge stability of BNNRs. The edge energy difference between BN and graphene could be used to guide the design of the specific hybrid BNC structures: in an armchair BNC nanoribbon (BNCNR), BN domains are expected to grow outside of C domains, while the opposite occurs in a zigzag BNCNR. More importantly, armchair BNCNRs can reproduce unique electronic properties of armchair graphene nanoribbons (GNRs), which are expected to be robust against edge functionalization or disorder. Within certain C/BN ratio, zigzag BNCNRs may exhibit intrinsic half-metallicity without any external constraints. These unexpected electronic properties of BNCNRs may offer unique opportunities to develop nanoscale electronics and spintronics beyond individual graphene and BN. Generally, these principles for designing BNC could be extended to other hybrid nanostructures.

preprint2010arXiv

Multiple Localized States and Magnetic Orderings in Partially Open Zigzag Carbon Nanotube Superlattices: An Ab Initio Study

Using first-principles calculations, we examine the electronic and magnetic properties of partially open zigzag carbon nanotube (CNT) superlattices. It is found that depending on their opening degree, these superlattices can exhibit multiple localized states around the Fermi energy. More importantly, some electronic states confined in some parts of the structure even have special magnetic orderings. We demonstrate that, as a proof of principle, some partially open zigzag CNT superlattices are by themselves giant (100%) magnetoresistive devices. Furthermore, the localized(and spin-polarized) states as well as the band gaps of the superlattices could be further modulated by external electric fields perpendicular to the tube axis, and a bias voltage along the tube axis may be used to control the conductance of two spin states. We believe that these results will open the way to the production of novel nanoscale electronic and spintronic devices.

preprint2010arXiv

Quantum Manifestations of Graphene Edge Stress and Edge Instability: A First-Principles Study

We have performed first-principles calculations of graphene edge stresses, which display two interesting quantum manifestations absent from the classical interpretation: the armchair edge stress oscillates with a nanoribbon width, and the zigzag edge stress is noticeably reduced by spin polarization. Such quantum stress effects in turn manifest in mechanical edge twisting and warping instability, showing features not captured by empirical potentials or continuum theory. Edge adsorption of H and Stone-Wales reconstruction are shown to provide alternative mechanisms in relieving the edge compression and hence to stabilize the planar edge structure.

preprint2010arXiv

Releasing H2 molecules with a partial pressure difference without the use of temperature

Using the pseudopotential density functional method as well as equilibrium thermodynamic functions, we explore the process of releasing H2 molecules adsorbed on a transition metal atom caused by the hydrogen-ammonia partial pressure difference. The H2 molecules bind to a transition metal atom at H2 pressure-NH3 pressure-temperature 50 atm-10-9 atm-25 °C, and they are released at 3 atm-10-6 atm-25 °C. This process involves the same mechanism responsible for carbon monoxide poisoning of hemoglobin with the O2-CO partial pressure difference. We show that our findings can be applicable to an approach to induce hydrogen desorption on nanostructured hydrogen storage materials without the need for increasing temperature.

preprint2010arXiv

Role of Symmetry in the Transport Properties of Graphene Nanoribbons under Bias

The intrinsic transport properties of zigzag graphene nanoribbons (ZGNRs) are investigated using first principles calculations. It is found that although all ZGNRs have similar metallic band structure, they show distinctly different transport behaviors under bias voltages, depending on whether they are mirror symmetric with respect to the midplane between two edges. Asymmetric ZGNRs behave as conventional conductors with linear current-voltage dependence, while symmetric ZGNRs exhibit unexpected very small currents with the presence of a conductance gap around the Fermi level. This difference is revealed to arise from different coupling between the conducting subbands around the Fermi level, which is dependent on the symmetry of the systems.

preprint2010arXiv

Spontaneous edge-defect formation and defect-induced conductance suppression in graphene nanoribbons

We present a first-principles study of the migration and recombination of edge defects (carbon adatom and/or vacancy) and their influence on electrical conductance in zigzag graphene nanoribbons (ZGNRs). It is found that at room temperature, the adatom is quite mobile while the vacancy is almost immobile along the edge of ZGNRs. The recombination of an adatom-vacancy pair leads to a pentagon-heptagon ring defect structure having a lower energy than the perfect edge, implying that such an edge-defect can be formed spontaneously. This edge defect can suppresses the conductance of ZGNRs drastically, which provides some useful hints for understanding the observed semiconducting behavior of the fabricated narrow GNRs.

preprint2010arXiv

The Half-Metallicity of Zigzag Graphene Nanoribbons with Asymmetric Edge Terminations

The spin-polarized electronic structure and half-metallicity of zigzag graphene nanoribbons (ZGNRs) with asymmetric edge terminations are investigated by using first principles calculations. It is found that compared with symmetric hydrogen-terminated counterparts, such ZGNRs maintain a spin-polarized ground state with the anti-ferromagnetic configuration at opposite edges, but their energy bands are no longer spin degenerate. In particular, the energy gap of one spin orientation decreases remarkably. Consequently, the ground state of such ZGNRs is very close to half-metallic state, and thus a smaller critical electric field is required for the systems to achieve the half-metallic state. Moreover, two kinds of studied ZGNRs present massless Dirac-fermion band structure when they behave like half-metals.

preprint2010arXiv

Towards Graphene Nanoribbon-based Electronics

The successful fabrication of single layer graphene has greatly stimulated the progress of the research on graphene. In this article, focusing on the basic electronic and transport properties of graphene nanoribbons (GNRs), we review the recent progress of experimental fabrication of GNRs, and the theoretical and experimental investigations of physical properties and device applications of GNRs. We also briefly discuss the research efforts on the spin polarization of GNRs in relation to the edge states.

preprint2009arXiv

Intrinsic anisotropy of thermal conductance in graphene nanoribbons

Thermal conductance of graphene nanoribbons (GNRs) with the width varying from 0.5 to 35 nm is systematically investigated using nonequilibrium Green's function method. Anisotropic thermal conductance is observed with the room temperature thermal conductance of zigzag GNRs up to ~ 30% larger than that of armchair GNRs. At room temperature, the anisotropy is found to disappear until the width is larger than 100 nm. This intrinsic anisotropy originate from different boundary condition at ribbon edges, and can be used to tune thermal conductance, which have important implications for the applications of GNRs in nanoelectronics and thermoelectricity.

preprint2003arXiv

Metal-to-semiconductor transition in squashed armchair carbon nanotubes

We investigate electronic transport properties of the squashed armchair carbon nanotubes, using tight-binding molecular dynamics and Green's function method. We demonstrate a metal-to-semiconductor transistion while squashing the nanotubes and a general mechanism for such transistion. It is the distinction of the two sublattices in the nanotube that opens an energy gap near the Fermi energy. We show that the transition has to be achieved by a combined effect of breaking of mirror symmetry and bond formation between the flattened faces in the squashed nanotubes.