Source author record

Jun-Hyung Cho

Jun-Hyung Cho 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

30works
6topics
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

30 published item(s)

preprint2022arXiv

Emerging Two-Dimensional Magnetism in Nonmagnetic Electrides Hf2X (X = S, Se, Te)

Recent experimental discoveries of two-dimensional (2D) magnets have triggered intense research activities to search for atomically thin magnetic systems. Using first-principles calculations, we predict the emergence of 2D magnetism in the monolayers (MLs), few layers, and surfaces of nonmagnetic layered electrides Hf2X (X = S, Se, Te) consisting of three-atom-thick Hf-X-Hf stacks. It is revealed that each bulk Hf2X hosts a novel quantum state of Dirac nodal lines with a high density of states arising from Hf-5d cationic and interlayer anionic electrons around -0.9 eV below the Fermi level EF. However, for the MLs, few layers, and surfaces of Hf2X, such hybridized states are shifted toward EF to generate van Hove singularities, leading to a Stoner instability. The resulting surface ferromagnetism gives rise to strongly spin-polarized topological surface states at Hf2X(001), demonstrating that anionic electrons, 2D magnetism, and band topology are entangled with each other. Our findings will open new perspectives for the discovery of 2D magnets via exploiting surface effects in nonmagnetic layered electrides.

preprint2021arXiv

Electron-phonon coupling and superconductivity in an alkaline earth hydride CaH$_6$ at high pressures

Recently, an alkaline earth hydride CaH$_6$ having a sodalitelike clathrate structure has been experimentally synthesized at megabar pressures with a maximum $T_c$ of 215 K, comparable to that of a rare earth hydride LaH$_{10}$. Here, based on first-principles calculations, we find that CaH$_6$ exhibits a huge peak in the Eliashberg spectral function $α^{2}F$ around the low-frequency region of H-derived phonon modes, in contrast to LaH$_{10}$ having a widely spreading spectrum of $α^{2}F$ over the whole frequencies of H-derived phonon modes. It is revealed that the huge peak of $α^{2}F$ in CaH$_6$ is associated with an effective electron-phonon coupling (EPC) between low-frequency optical phonons and hybridized H 1$s$ and Ca 3$d$ states near the Fermi energy. As pressure increases, the strengthened H$-$H covalent bonding not only induces a hardening of optical phonon modes but also reduces the electron-phonon matrix elements related to the low-frequency optical modes, thereby leading to a lowering of the EPC constant. It is thus demonstrated that H-derived low-frequency phonon modes play an important role in the pressure-induced variation of $T_c$ in CaH$_6$. Furthermore, unlike the presence of two distinct superconducting gaps in LaH$_{10}$, CaH$_6$ is found to exhibit a single isotropic superconducting gap.

preprint2021arXiv

Formation mechanism of chemically precompressed hydrogen clathrates in metal superhydrides

Recently, the experimental discovery of high-$T_c$ superconductivity in compressed hydrides H$_3$S and LaH$_{10}$ at megabar pressures has triggered searches for various superconducting superhydrides. It was experimentally observed that thorium hydrides, ThH$_{10}$ and ThH$_9$, are stabilized at much lower pressures compared to LaH$_{10}$. Based on first-principles density-functional theory calculations, we reveal that the isolated Th frameworks of ThH$_{10}$ and ThH$_9$ have relatively more excess electrons in interstitial regions than the La framework of LaH$_{10}$. Such interstitial excess electrons easily participate in the formation of anionic H cage surrounding metal atom. The resulting Coulomb attraction between cationic Th atoms and anionic H cages is estimated to be stronger than the corresponding one of LaH$_{10}$, thereby giving rise to larger chemical precompressions in ThH$_{10}$ and ThH$_9$. Such a formation mechanism of H clathrates can also be applied to another experimentally synthesized superhydride CeH$_9$, confirming the experimental evidence that the chemical precompression in CeH$_9$ is larger than that in LaH$_{10}$. Our findings demonstrate that interstitial excess electrons in the isolated metal frameworks of high-pressure superhydrides play an important role in generating the chemical precompression of H clathrates.

preprint2021arXiv

Interlayer Exchange Interaction Driven Topological Phase Transition in Antiferromagnetic Electride Gd$_2$O

Based on first-principles calculations, we discover a two-dimensional layered antiferromagnetic (AFM) electride Gd$_2$O, where anionic excess electrons exist in the interstitial spaces between positively charged cationic layers. It is revealed that each cationic layer composed of three-atom-thick Gd$-$O$-$Gd stacks has in-plane ferromagnetic and out-of-plane AFM superexchange interactions between the localized Gd 4$f$ spins through O 2$p$ orbitals. Furthermore, the interlayer superexchange mediated by the hybridized Gd-5$d$ and interstitial-$s$-like states involves intimate couplings between the spin, lattice, and charge degrees of freedom, thereby inducing simultaneous magnetic, structural, and electronic phase transitions. The resulting ground state with the simple hexagonal lattice hosts massless Dirac fermions protected by nonsymmorphic magnetic symmetry, as well as massive Dirac fermions. We thus demonstrate that the anionic excess electrons in Gd$_2$O play a crucial role in the emergence of magnetic Dirac semimetal states, therefore offering an intriguing interplay between 2D magnetic electrides and topological physics.

preprint2021arXiv

Origin of Charge Density Wave in Layered Kagome Metal CsV$_3$Sb$_5$

Using first-principles calculations, we identify the origin of the observed charge density wave (CDW) formation in a layered kagome metal CsV$_3$Sb$_5$. It is revealed that the structural distortion of kagome lattice forming the trimeric and hexameric V atoms is accompanied by the stabilization of quasimolecular states, which gives rise to the opening of CDW gaps for the V-derived multibands lying around the Fermi level. This Jahn-Teller-like instability having the local lattice distortion and its derived quasimolecular states is a driving force of the CDW order. Specifically, the saddle points of multiple Dirac bands near the Fermi level, located at the $M$ point, are hybridized to disappear along the $k_z$ direction, therefore not supporting the widely accepted Peierls-like electronic instability due to Fermi surface nesting. It is further demonstrated that applied hydrostatic pressure significantly reduces the interlayer spacing to destabilize the quasimolecular states, leading to a disappearance of the CDW phase at a pressure of ${\sim}$2 GPa. The presently proposed underlying mechanism of the CDW order in CsV$_3$Sb$_5$ can also be applicable to other isostructural kagome lattices such as KV$_3$Sb$_5$ and RbV$_3$Sb$_5$.

preprint2020arXiv

Isotope effect in superconducting lanthanum hydride under high compression

Recently, the discovery of room-temperature superconductivity (SC) was experimentally realized in the fcc phase of LaH$_{10}$ under megabar pressure. Specifically, the isotope effect of $T_{\rm c}$ was measured by the replacement of hydrogen (H) with deuterium (D), demonstrating a driving role of phonons in the observed room-temperature SC. Herein, based on the first-principles calculations within the harmonic approximation, we reveal that (i) the identical electron-phonon coupling constants of fcc LaH$_{10}$ and LaD$_{10}$ decrease monotonously with increasing pressure and (ii) the isotope effect of $T_{\rm c}$ is nearly proportional to $M^{-α}$ ($M$: ionic mass) with $α$ ${\approx}$ 0.465, irrespective of pressure. The predicted value of $α$ agrees well with the experimental one ($α=0.46$) measured at around 150 GPa. Thus, our findings provide a theoretical confirmation of the conventional electron-phonon coupling mechanism in a newly discovered room-temperature superconductor of compressed LaH$_{10}$.

preprint2020arXiv

Origin of enhanced chemical precompression in cerium hydride CeH$_{9}$

The rare-earth metal hydrides with clathrate structures have been highly attractive because of their promising high-$T_{\rm c}$ superconductivity at high pressure. Recently, cerium hydride CeH$_9$ composed of Ce-encapsulated clathrate H cages was synthesized at much lower pressures of 80$-$100 GPa, compared to other experimentally synthesized rare-earth hydrides such as LaH$_{10}$ and YH$_6$. Based on density-functional theory calculations, we find that the Ce 5$p$ semicore and 4$f$/5$d$ valence states strongly hybridize with the H 1$s$ state, while a transfer of electrons occurs from Ce to H atoms. Further, we reveal that the delocalized nature of Ce 4$f$ electrons plays an important role in the chemical precompression of clathrate H cages. Our findings not only suggest that the bonding nature between the Ce atoms and H cages is characterized as a mixture of ionic and covalent, but also have important implications for understanding the origin of enhanced chemical precompression that results in the lower pressures required for the synthesis of CeH$_9$.

preprint2020arXiv

Stabilization mechanism of clathrate H cages in a room-temperature superconductor LaH$_{10}$

Lanthanum hydride LaH$_{10}$ with a sodalitelike clathrate structure was experimentally realized to exhibit a room-temperature superconductivity under megabar pressures. Based on first-principles calculations, we reveal that the metal framework of La atoms has the excess electrons at interstitial regions. Such anionic electrons are easily captured to form a stable clathrate structure of H cages. We thus propose that the charge transfer from La to H atoms is mostly driven by the electride property of the La framework. Further, the interaction between La atoms and H cages induces a delocalization of La-5$p$ semicore states to hybridize with H-1$s$ state. Consequently, the bonding nature between La atoms and H cages is characterized as a mixture of ionic and covalent. Our findings demonstrate that anionic and semicore electrons play important roles in stabilizing clathrate H cages in LaH$_{10}$, which can be broadly applicable to other high-pressure rare-earth hydrides with clathrate structures.

preprint2019arXiv

Multiband Nature of the Room-Temperature Superconductivity in Compressed LaH$_{10}$

Recently, the discovery of room-temperature superconductivity (SC) was experimentally realized in the fcc phase of LaH$_{10}$ under megabar pressures. This SC of compressed LaH$_{10}$ has been explained in terms of strong electron-phonon coupling (EPC), but the mechanism of how the large EPC constant and high superconducting transition temperature $T_{\rm c}$ are attained has not yet been clearly identified. Based on the density-functional theory and the Migdal-Eliashberg formalism, we reveal the presence of two nodeless, anisotropic superconducting gaps on the Fermi surface (FS). Here, the small gap is mostly associated with the hybridized states of H $s$ and La $f$ orbitals on the three outer FS sheets, while the large gap arises mainly from the hybridized state of neighboring H $s$ or $p$ orbitals on the one inner FS sheet. Further, we find that the EPC constant of compressed YH$_{10}$ with the same sodalite-like clathrate structure is enhanced due to the two additional FS sheets, leading to a higher $T_{\rm c}$ than LaH$_{10}$. It is thus demonstrated that the multiband pairing of hybridized electronic states is responsible for the large EPC constant and room-temperature SC in compressed hydrides LaH$_{10}$ and YH$_{10}$.

preprint2019arXiv

Two-dimensional topological semimetal states in monolayers Cu$_2$Ge, Fe$_2$Ge, and Fe$_2$Sn

Recent experimental realizations of the topological semimetal states in several monolayer systems are very attractive because of their exotic quantum phenomena and technological applications. Based on first-principles density-functional theory calculations including spin-orbit coupling, we here explore the drastically different two-dimensional (2D) topological semimetal states in three monolayers Cu$_2$Ge, Fe$_2$Ge, and Fe$_2$Sn, which are isostructural with a combination of the honeycomb Cu or Fe lattice and the triangular Ge or Sn lattice. We find that (i) the nonmagnetic (NM) Cu$_{2}$Ge monolayer having a planar geometry exhibits the massive Dirac nodal lines, (ii) the ferromagentic (FM) Fe$_2$Ge monolayer having a buckled geometry exhibits the massive Weyl points, and (iii) the FM Fe$_2$Sn monolayer having a planar geometry and an out-of-plane magnetic easy axis exhibits the massless Weyl nodal lines. It is therefore revealed that mirror symmetry cannot protect the four-fold degenerate Dirac nodal lines in the NM Cu$_{2}$Ge monolayer, but preserves the doubly degenerate Weyl nodal lines in the FM Fe$_{2}$Sn monolayer. Our findings demonstrate that the interplay of crystal symmetry, magnetic easy axis, and band topology is of importance for tailoring various 2D topological states in Cu$_2$Ge, Fe$_2$Ge, and Fe$_2$Sn monlayers.

preprint2016arXiv

Mott- versus Slater-type Insulating Nature of Two-Dimensional Sn Atom Lattice on SiC(0001)

Semiconductor surfaces with narrow surface bands provide unique playgrounds to search for Mott-insulating state. Recently, a combined experimental and theoretical study [Phys. Rev. Lett. 114, 247602 (2015)] of the two-dimensional (2D) Sn atom lattice on a wide-gap SiC(0001) substrate proposed a Mott-type insulator driven by strong on-site Coulomb repulsion U. Our systematic density-functional theory (DFT) study with local, semilocal, and hybrid exchange-correlation functionals shows that the Sn dangling-bond state largely hybridizes with the substrate Si 3p and C 2p states to split into three surface bands due to the crystal field. Such a hybridization gives rise to the stabilization of the antiferromagnetic order via superexchange interactions. The band gap and the density of states predicted by the hybrid DFT calculation agree well with photoemission data. Our findings not only suggest that the Sn/SiC(0001) system can be represented as a Slater-type insulator driven by long-range magnetism, but also have an implication that taking into account long-range interactions beyond the on-site interaction would be of importance for properly describing the insulating nature of Sn/SiC(0001).

preprint2016arXiv

Origin of the Metal-Insulator Transition of Indium Atom Wires on Si(111)

As a prototypical one-dimensional electron system, self-assembled indium (In) nanowires on the Si(111) surface have been believed to drive a metal-insulator transition by a charge-density-wave (CDW) formation due to electron-phonon coupling. Here, our first-principles calculations demonstrate that the structural phase transition from the high-temperature 4x1 phase to the low-temperature 8x2 phase occurs through an exothermic reaction with the consecutive bond-breaking and bond-making processes, giving rise to an energy barrier between the two phases as well as a gap opening. This atomistic picture for the phase transition not only identifies its first-order nature but also solves a long-standing puzzle of the origin of the metal-insulator transition in terms of the x2 periodic lattice reconstruction of In hexagons via bond breakage and new bond formation, not by the Peierls instability-driven CDW formation.

preprint2016arXiv

Spin-orbit coupling effects on the stability of two competing structures in Pb/Si(111) and Pb/Ge(111)

Using first-principles density-functional theory (DFT) calculations, we investigate the 4/3-monolayer structure of Pb on the Si(111) or Ge(111) surface within the two competing structural models termed the H$_3$ and T$_4$ structures. We find that the spin-orbit coupling (SOC) influences the relative stability of the two structures in both the Pb/Si(111) and Pb/Ge(111) systems: i.e., our DFT calculation without including the SOC predicts that the T$_4$ structure is energetically favored over the H$_3$ structure by $ΔE$ = 25 meV for Pb/Si(111) and 22 meV for Pb/Ge(111), but the inclusion of SOC reverses their relative stability as $ΔE$ = $-$12 and $-$7 meV, respectively. Our analysis shows that the SOC-induced switching of the ground state is attributed to a more asymmetric surface charge distribution in the H$_3$ structure, which gives rise to a relatively larger Rashba spin splitting of surface states as well as a relatively larger pseudo-gap opening compared to the T$_4$ structure. By the nudged elastic band calculation, we obtain a sizable energy barrier from the H$_3$ to the T$_4$ structure as ${\sim}$0.59 and ${\sim}$0.27 eV for Pb/Si(111) and Pb/Ge(111), respectively. It is thus likely that the two energetically competing structures can coexist at low temperatures.

preprint2016arXiv

Stacking-sequence-independent band structure and shear exfoliation of two-dimensional electride materials

The electronic band structure of crystals is generally influenced by the periodic arrangement of their constituent atoms. Specifically, the emerging two-dimensional (2D) layered structures have shown different band structures with respect to their stacking configurations. Here, based on first-principles density-functional theory calculations, we demonstrate that the band structure of the recently synthesized 2D Ca$_2$N electride changes little for the stacking sequence as well as the lateral interlayer shift. This intriguing invariance of band structure with respect to geometrical variations can be attributed to a complete screening of [Ca$_2$N]$^{+}$ cationic layers by anionic excess electrons delocalized between the cationic layers. The resulting weak interactions between 2D dressed cationic layers give rise to not only a shallow potential barrier for bilayer sliding but also an electron-doping facilitated shear exfoliation. Our findings open a route for exploration of the peculiar geometry-insensitive electronic properties in 2D electride materials, which will be useful for future thermally stable electronic applications.

preprint2015arXiv

Anomalous doping effect in black phosphorene from first-principles calculations

Using first-principles density functional theory calculations, we investigate the geometries, electronic structures, and thermodynamic stabilities of substitutionally doped phosphorene sheets with group III, IV, V, and VI elements. We find that the electronic properties of phosphorene are drastically modified by the number of valence electrons in dopant atoms. The dopants with even number of valence electrons enable the doped phosphorenes to have a metallic feature, while the dopants with odd number of valence electrons keep the semiconducting feature with a larger band gap than the undoped phosphorene. This even-odd behavior is attributed to the peculiar bonding characteristics of phosphorene and the strong hybridization of sp orbitals between dopants and phosphorene. The calculated formation energies of various substitutional dopants in phosphorene show that such doped systems can be thermodynamically stable.

preprint2015arXiv

Competing Magnetic Orderings and Tunable Topological States in Two-Dimensional Hexagonal Organometallic Lattices

The exploration of topological states is of significant fundamental and practical importance in contemporary condensed matter physics, for which the extension to two-dimensional (2D) organometallic systems is particularly attractive. Using first-principles calculations, we show that a 2D hexagonal triphenyl-lead lattice composed of only main group elements is susceptible to a magnetic instability, characterized by a considerably more stable antiferromagnetic (AFM) insulating state rather than the topologically nontrivial quantum spin Hall state proposed recently. Even though this AFM phase is topologically trivial, it possesses an intricate emergent degree of freedom, defined by the product of spin and valley indices, leading to Berry curvature-induced spin and valley currents under electron or hole doping. Furthermore, such a trivial band insulator can be tuned into a topologically nontrivial matter by the application of an out-of-plane electric field, which destroys the AFM order, favoring instead ferrimagnetic spin ordering and a quantum anomalous Hall state with a non-zero topological invariant. These findings further enrich our understanding of 2D hexagonal organometallic lattices for potential applications in spintronics and valleytronics.

preprint2015arXiv

Contrasting interedge superexchange interactions of graphene nanoribbons embedded in h-BN and graphane

Based on first-principles density-functional theory calculations, we present a comparative study of the elec- tronic structures of ultranarrow zigzag graphene nanoribbons (ZGNRs) embedded in hexagonal boron nitride (BN) sheet and fully hydrogenated graphene (graphane) as a function of their width N (the number of zigzag C chains composing the ZGNRs). We find that ZGNRs/BN have the nonmagnetic ground state except at N = 5 and 6 with a weakly stabilized half-semimetallic state, whereas ZGNRs/graphane with N >= 2 exhibit a strong antiferromagnetic coupling between ferromagnetically ordered edge states on each edge. It is revealed that the disparate magnetic properties of the two classes of ZGNRs are attributed to the contrasting interedge superexchange interactions arising from different interface structures: i.e., the asymmetric interface structure of ZGNRs/BN gives a relatively short-range and weak superexchange interaction between the two inequivalent edge states, while the symmetric interface structure of ZGNRs/graphane gives a long-range, strong interedge superexchange interaction.

preprint2015arXiv

Dilute magnetic semiconductor and half metal behaviors in 3d transition-metal doped black and blue phosphorenes: a first-principles study

We present first-principles density-functional calculations for the structural, electronic, and magnetic properties of substitutional 3d transition metal (TM) impurities in two-dimensional black and blue phosphorenes. We find that the magnetic properties of such substitutional impurities can be understood in terms of a simple model based on the Hund's rule. The TM-doped black phosphorenes with Ti, V, Cr, Mn, Fe and Ni impurities show dilute magnetic semiconductor (DMS) properties while those with Sc and Co impurities show nonmagnetic properties. On the other hand, the TM-doped blue phosphorenes with V, Cr, Mn and Fe impurities show DMS properties, those with Ti and Ni impurities show half-metal properties, whereas Sc and Co doped systems show nonmagnetic properties. We identify two different regimes depending on the occupation of the hybridized electronic states of TM and phosphorous atoms: (i) bonding states are completely empty or filled for Sc- and Co-doped black and blue phosphorenes, leading to non-magnetic; (ii) non-bonding d states are partially occupied for Ti-, V-, Cr-, Mn-, Fe- and Ni-doped black and blue phosphorenes, giving rise to large and localized spin moments. These results provide a new route for the potential applications of dilute magnetic semiconductor and half-metal in spintronic devices by employing black and blue phosphorenes.

preprint2015arXiv

Equivalence of Electronic and Mechanical Stresses in Structural Phase Stabilization: A Case Study of Indium Wires on Si(111)

It was recently proposed that the stress state of a material can also be altered via electron or hole doping, a concept termed electronic stress (ES), which is different from the traditional mechanical stress (MS) due to lattice contraction or expansion. Here we demonstrate the equivalence of ES and MS in structural stabilization, using In wires on Si(111) as a prototypical example. Our systematic density-functional theory calculations reveal that, first, for the same degrees of carrier doping into the In wires, the ES of the high-temperature metallic 4x1 structure is only slightly compressive, while that of the low-temperature insulating 8x2 structure is much larger and highly anisotropic. As a consequence, the intrinsic energy difference between the two phases is significantly reduced towards electronically phase-separated ground states. Our calculations further demonstrate quantitatively that such intriguing phase tunabilities can be achieved equivalently via lattice-contraction induced MS in the absence of charge doping. We also validate the equivalence through our detailed scanning tunneling microscopy experiments. The present findings have important implications in understanding the underlying driving forces involved in various phase transitions of simple and complex systems alike.

preprint2015arXiv

Giant Spin-Orbit-Induced Spin Splitting in Bi Chains on GaAs(110)

The search for one-dimensional electron systems with a giant Rashba-type spin splitting is of importance for the application of spin transport. Here we report, based on a first-principles density-functional-theory calculation, that Bi zigzag chains formed on a heterogeneous GaAs(110) surface have a giant spin splitting of surface states. This giant spin splitting is revealed to originate from spin-orbit coupling (SOC) and electric dipole interaction that are significantly enhanced by (i) the asymmetric surface charge distribution due to the strong SOC-induced hybridization of the Bi px , py , and pz orbitals and (ii) the large out-of-plane and in-plane potential gradients generated by two geometrically and electronically inequivalent Bi atoms bonding to Ga and As atoms. The results demonstrate an important implication of the in-plane and out-of-plane asymmetry of the Bi/GaAs(110) interface system in producing the giant spin splitting with the in-plane and out-of-plane spin components.

preprint2015arXiv

Nature of the Insulating Ground State of the 5d PostPerovskite CaIrO3

The insulating ground state of the 5d transition metal oxide CaIrO3 has been classified as a Mott-type insulator. Based on a systematic density functional theory (DFT) study with local, semilocal, and hybrid exchange-correlation functionals, we reveal that the Ir t2g states exhibit large splittings and one-dimensional electronic states along the c axis due to a tetragonal crystal field. Our hybrid DFT calculation adequately describes the antiferromagnetic (AFM) order along the c direction via a superexchange interaction between Ir4+ spins. Furthermore, the spin-orbit coupling (SOC) hybridizes the t2g states to open an insulating gap. These results indicate that CaIrO3 can be represented as a spin-orbit Slater insulator, driven by the interplay between a long-range AFM order and the SOC. Such a Slater mechanism for the gap formation is also demonstrated by the DFT + dynamical mean field theory calculation, where the metal-insulator transition and the paramagnetic to AFM phase transition are concomitant with each other.

preprint2015arXiv

Origin of Symmetric Dimer Images of Si(001) Observed in Low-Temperature STM

It has been a long-standing puzzle why buckled dimers of the Si(001) surface appeared symmetric below 20 K in scanning tunneling microscopy (STM) experiments. Although such symmetric dimer images were concluded to be due to an artifact induced by STM measurements, its underlying mechanism is still veiled. Here,we demonstrate, based on a first-principles density-functional theory calculation, that the symmetric dimer images are originated from the flip-flop motion of buckled dimers, driven by quantum tunneling (QT). It is revealed that at low temperature the tunneling-induced surface charging with holes reduces the energy barrier for the flipping of buckled dimers, thereby giving rise to a sizable QT-driven frequency of the flip-flop motion. However, such a QT phenomenon becomes marginal in the tunneling-induced surface charging with electrons. Our findings provide an explanation for low-temperature STM data that exhibits apparent symmetric (buckled) dimer structure in the filled-state (empty-state) images.

preprint2014arXiv

Antiferromagnetic Slater Insulator Phase of Na2IrO3

Using a hybrid density-functional theory (DFT) calculation including spin-orbit coupling (SOC), we predict that the zigzag antiferromagnetic (AFM) ground state of the honeycomb layered compound Na2IrO3 opens the observed insulating gap through a long-range magnetic order. We show that the effect of SOC and the correction of self-interaction error inherent in previous local or semilocal DFT calculations play crucial roles in predicting the band gap formation in Na2IrO3. It is revealed that the itinerant AFM order with a strong suppression of the Ir magnetic moment is attributed to a considerable hybridization of the Ir 5d orbitals with the O 2p orbitals. Thus, our results suggest that the insulating phase of Na2IrO3 can be represented as a Slater insulator driven by itinerant magnetism.

preprint2013arXiv

Antiferromagnetic ordering of dangling-bond electrons at the stepped Si(001) surface

Using first-principles density-functional calculations, we explore the possibility of magnetic order at the rebonded $D_B$ step of the Si(001) surface. The rebonded $D_B$ step containing threefold coordinated Si atoms can be treated as a one-dimensional dangling-bond (DB) wire along the step edge. We find that Si atoms composing the step edge are displaced up and down alternatively due to Jahn-Teller-like distortion, but, if Si dimers on the terrace are passivated by H atoms, the antiferromagnetic (AFM) order can be stabilized at the step edge with a suppression of Jahn-Teller-like distortion. We also find that the energy preference of AFM order over Jahn-Teller-like distortion is enhanced in an oscillatory way as the length of DB wires decreases, showing the so-called quantum size effects.

preprint2013arXiv

Driving force of phase transition in Indium nanowires on Si(111)

The precise driving force of the phase transition in indium nanowires on Si(111) has been controversial whether it is driven by a Peierls instability or by a simple energy lowering due to a periodic lattice distortion. The present van der Waals (vdW) corrected hybrid density functional calculation predicts that the low-temperature 8x2 structure whose building blocks are indium hexagons is energetically favored over the room-temperature 4x1 structure. We show that the correction of self-interaction error and the inclusion of vdW interactions play crucial roles in describing the covalent bonding, band-gap opening, and energetics of hexagon structures. The results manifest that the formation of hexagons occurs by a simple energy lowering due to the lattice distortion, not by a charge density wave formation arising from Fermi surface nesting.

preprint2013arXiv

Ferrimagnetic Slater Insulator Phase of the Sn/Ge(111) Surface

We have performed the semilocal and hybrid density-functional theory (DFT) studies of the Sn/Ge(111) surface to identify the origin of the observed insulating ${\sqrt{3}}{\times}{\sqrt{3}}$ phase below ${\sim}$30 K. Contrasting with the semilocal DFT calculation predicting a metallic 3${\times}$3 ground state, the hybrid DFT calculation including van der Waals interactions shows that the insulating ferrimagnetic structure with ${\sqrt{3}}{\times}{\sqrt{3}}$ structural symmetry is energetically favored over the metallic 3${\times}$3 structure. It is revealed that the correction of self-interaction error with a hybrid exchange-correlation functional gives rise to a band-gap opening induced by a ferrimagnetic order. The results manifest that the observed insulating phase is attributed to the Slater mechanism via itinerant magnetic order rather than the hitherto accepted Mott-Hubbard mechanism via electron correlations.

preprint2013arXiv

Fluorine-induced local magnetic moment in graphene: A hybrid DFT study

Recent experimental evidence that fluorinated graphene creates local magnetic moments around F adatoms has not been supported by semilocal density-functional theory (DFT) calculations where the adsorption of an isolated F adatom induces no magnetic moment in graphene. Here, we show that such an incorrect prediction of the nonmagnetic ground state is due to the self-interaction error inherent in semilocal exchange-correlation functionals. The present hybrid DFT calculation for an isolated F adatom on graphene predicts not only a spin-polarized ground state with a spin moment of 1 bohr magneton, but also a long-range spin polarization caused by the bipartite nature of the graphene lattice as well as the induced spin polarization of the graphene states. The results provide support for the experimental observations of local magnetic moments in fluorinated graphene.

preprint2013arXiv

Physisorption of DNA nucleobases on h-BN and graphene: vdW-corrected DFT calculations

We present a comparative study of DNA nucleobases [guanine (G), adenine (A), thymine (T), and cytosine (C)] adsorbed on hexagonal boron nitride (\textit{h}-BN) sheet and graphene, using local, semilocal, and van der Waals (vdW) energy-corrected density-functional theory (DFT) calculations. Intriguingly, despite the very different electronic properties of BN sheet and graphene, we find rather similar binding energies for the various nucleobase molecules when adsorbed on the two types of sheets. The calculated binding energies of the four nucleobases using the local, semilocal, and DFT+vdW schemes are in the range of 0.54 ${\sim}$ 0.75 eV, 0.06 ${\sim}$ 0.15 eV, and 0.93 ${\sim}$ 1.18 eV, respectively. In particular, the DFT+vdW scheme predicts not only a binding energy predominantly determined by vdW interactions between the base molecules and their substrates decreasing in the order of G$>$A$>$T$>$C, but also a very weak hybridization between the molecular levels of the nucleobases and the $π$-states of the BN sheet or graphene. This physisorption of G, A, T, and C on the BN sheet (graphene) induces a small interfacial dipole, giving rise to an energy shift in the work function by 0.11 (0.22), 0.09 (0.15), $-$0.05 (0.01), and 0.06 (0.13) eV, respectively.

preprint2011arXiv

Atomic structure, energetics, and dynamics of topological solitons in Indium chains on Si(111) surfaces

Based on scanning tunneling microscopy and first-principles theoretical studies, we characterize the precise atomic structure of a topological soliton in In chains grown on Si(111) surfaces. Variable-temperature measurements of the soliton population allow us to determine the soliton formation energy to be ~60 meV, smaller than one half of the band gap of ~200 meV. Once created, these solitons have very low mobility, even though the activation energy is only about 20 meV; the sluggish nature is attributed to the exceptionally low attempt frequency for soliton migration. We further demonstrate local electric field-enhanced soliton dynamics.

preprint2011arXiv

Quantum Size Effect and Electronic Stability of Freestanding Metal Atom Wires

Using DFT calculations, we present a thorough study of the quantum size effects on the stability of freestanding metal atom wires. Our systems include Na, Ag, Au, In, Ga and Pb atom wires, i.e. $s$, $sd$, and $sp$ electron prototypes. We found that the total energy always oscillates with the wire length, which clearly indicates the existence of preferred lengths. Increasing the length, the s-system exhibits even-odd oscillations following a $a/x +b/x^2$ decay law in the stability, which can be attributed to electron band filling and quantum confinement along the wire. The $sd$-system exhibits a similar oscillation pattern, even in the presence of $sd$ hybridization. In $sp$-system, the energy oscillations are beyond the simple even-odd nature, likely due to unpaired p orbitals and the corresponding nontrival band filling. Our findings clearly demonstrate that electronic contribution is quite critical to the stability of freestanding wires, and this stability may be important even when wires are deposited on substrates or strained. This study sheds light on the underlying formation mechanism of metal atom wires.