BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Efthimios Kaxiras

Efthimios Kaxiras appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.mtrl-scicond-mat.mes-hallcond-mat.str-elcond-mat.supr-conphysics.comp-phcond-mat.softmtrl-thBiological Physicscond-matcond-mat.otherphysics.chem-phPopulations and Evolutioncond-mat.dis-nnMachine Learningphysics.flu-dynphysics.opticsquant-ph

Catalog footprint

What is connected

52works
17topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

52 published item(s)

preprint2022arXiv

Hydrogen Atoms on Zigzag Graphene Nanoribbons: Chemistry and Magnetism Meet at the Edge

Although the unconventional $π$-magnetism at the zigzag edges of graphene holds promise for a wide array of applications, whether and to what degree it plays a role in their chemistry remains poorly understood. Here, we investigate the addition of a hydrogen atom $-$ the simplest yet the most experimentally relevant adsorbate $-$ to zigzag graphene nanoribbons (ZGNRs). We show that the $π$-magnetism governs the chemistry of ZGNRs, giving rise to a site-dependent reactivity of the carbon atoms and driving the hydrogenation process to the nanoribbon edges. Conversely, the chemisorbed hydrogen atom governs the $π$-magnetism of ZGNRs, acting as a spin-$\frac{1}{2}$ paramagnetic center in the otherwise antiferromagnetic ground state and spin-polarizing the charge carriers at the band extrema. Our findings establish a comprehensive picture of the peculiar interplay between chemistry and magnetism that emerges at the zigzag edges of graphene.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Non-Abelian topological defects and strain mapping in 2D moiré materials

We present a general method to analyze the topological nature of the domain boundary connectivity that appeared in relaxed moiré superlattice patterns at the interface of 2-dimensional (2D) van der Waals (vdW) materials. At large enough moiré lengths, all moiré systems relax into commensurated 2D domains separated by networks of dislocation lines. The nodes of the 2D dislocation line network can be considered as vortex-like topological defects. We find that a simple analogy to common topological systems with an $S^1$ order parameter, such as a superconductor or planar ferromagnet, cannot correctly capture the topological nature of these defects. For example, in twisted bilayer graphene, the order parameter space for the relaxed moiré system is homotopy equivalent to a punctured torus. Here, the nodes of the 2D dislocation network can be characterized as elements of the fundamental group of the punctured torus, the free group on two generators, endowing these network nodes with non-Abelian properties. Extending this analysis to consider moiré patterns generated from any relative strain, we find that antivortices occur in the presence of anisotropic heterostrain, such as shear or anisotropic expansion, while arrays of vortices appear under twist or isotropic expansion between vdW materials. Experimentally, utilizing the dark field imaging capability of transmission electron microscopy (TEM), we demonstrate the existence of vortex and antivortex pair formation in a moiré system, caused by competition between different types of heterostrains in the vdW interfaces. We also present a methodology for mapping the underlying heterostrain of a moiré structure from experimental TEM data, which provides a quantitative relation between the various components of heterostrain and vortex-antivortex density in moiré systems.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Outstanding thermoelectric performance predicted for out-of-plane p-doped GeSe

The record-breaking thermoelectric performance of tin selenide (SnSe) has motivated the investigation of analogue compounds with the same structure. A promising candidate that emerged recently is germanium selenide (GeSe). Here, using extensive first-principles calculations of the hole-phonon and hole-impurity scattering, we investigate the thermoelectric transport properties of the orthorhombic phase of p-doped GeSe. We predict outstanding thermoelectric performance for GeSe over a broad range of temperatures due to its high Seebeck coefficients, extremely low Lorenz numbers, ultralow total thermal conductivity, and relatively large band gap. In particular, the out-of-plane direction in GeSe presents equivalent or even higher performance than SnSe for temperatures above 500 K. By extending the analysis to 900 K, we obtained an ultrahigh value for the thermoelectric figure of merit (zT = 3.2) at the optimal hole density of 4x10^19 cm^-3. Our work provides strong motivation for continued experimental work focusing on improving the GeSe doping efficiency in order to achieve this optimal hole density.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Twistronics of Janus transition metal dichalcogenide bilayers

Twisted multilayers of two-dimensional (2D) materials are an increasingly important platform for investigating quantum phases of matter, and in particular, strongly correlated electrons. The moiré pattern introduced by the relative twist between layers creates effective potentials of long-wavelength, leading to electron localization. However, in contrast to the abundance of 2D materials, few twisted heterostructures have been studied until now. Here we develop a first-principle continuum theory to study the electronic bands introduced by moire patterns of twisted Janus transition metal dichalcogenides (TMD) homo- and hetero-bilayers. The model includes lattice relaxation, stacking-dependent effective mass, and Rashba spin-orbit coupling. We then perform a high-throughput generation and characterization of DFT-extracted continuum models for more than a hundred possible combinations of materials and stackings. Our model predicts that the moiré physics and emergent symmetries depend on chemical composition, vertical layer orientation, and twist angle, so that the minibands wavefunctions can form triangular, honeycomb, and Kagome networks. Rashba spin-orbit effects, peculiar of these systems, can dominate the moiré bandwidth at small angles. Our work enables the detailed investigation of Janus twisted heterostructures, allowing the discovery and control of novel electronic phenomena.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Ferromagnetic helical nodal line and Kane-Mele spin-orbit coupling in kagome metal Fe3Sn2

The two-dimensional kagome lattice hosts Dirac fermions at its Brillouin zone corners K and K', analogous to the honeycomb lattice. In the density functional theory electronic structure of ferromagnetic kagome metal Fe$_3$Sn$_2$, without spin-orbit coupling we identify two energetically split helical nodal lines winding along $z$ in the vicinity of K and K' resulting from the trigonal stacking of the kagome layers. We find that hopping across A-A stacking introduces a layer splitting in energy while that across A-B stacking controls the momentum space amplitude of the helical nodal lines. The effect of spin-orbit coupling is found to resemble that of a Kane-Mele term, where the nodal lines can either be fully gapped to quasi-two-dimensional massive Dirac fermions, or remain gapless at discrete Weyl points depending on the ferromagnetic moment orientation. Aside from numerically establishing Fe$_3$Sn$_2$ as a model Dirac kagome metal, our results provide insights into materials design of topological phases from the lattice point of view, where paradigmatic low dimensional lattice models often find realizations in crystalline materials with three-dimensional stacking.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Imprinting Tunable $π$-Magnetism in Graphene Nanoribbons via Edge Extensions

Magnetic carbon nanostructures are currently under scrutiny for a wide spectrum of applications. Here, we theoretically investigate armchair graphene nanoribbons patterned with asymmetric edge extensions consisting of laterally fused naphtho groups, as recently fabricated via on-surface synthesis. We show that an individual edge extension acts as a spin-$\frac{1}{2}$ center and develops a sizable spin-polarization of the conductance around the band edges. The Heisenberg exchange coupling between a pair of edge extensions is dictated by the position of the second naphtho group in the carbon backbone, thus enabling ferromagnetic, antiferromagnetic, or non-magnetic states. The periodic arrangement of edge extensions yields full spin-polarization at the band extrema, and the accompanying ferromagnetic ground state can be driven into non-magnetic or antiferromagnetic phases through external stimuli. Overall, our work reveals precise tunability of the $π$-magnetism in graphene nanoribbons induced by naphtho groups, thereby establishing these one-dimensional architectures as suitable platforms for logic spintronics.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Modeling the effect of the vaccination campaign on the Covid-19 pandemic

Population-wide vaccination is critical for containing the SARS-CoV-2 (Covid-19) pandemic when combined with restrictive and prevention measures. In this study, we introduce SAIVR, a mathematical model able to forecast the Covid-19 epidemic evolution during the vaccination campaign. SAIVR extends the widely used Susceptible-Infectious-Removed (SIR) model by considering the Asymptomatic (A) and Vaccinated (V) compartments. The model contains several parameters and initial conditions that are estimated by employing a semi-supervised machine learning procedure. After training an unsupervised neural network to solve the SAIVR differential equations, a supervised framework then estimates the optimal conditions and parameters that best fit recent infectious curves of 27 countries. Instructed by these results, we performed an extensive study on the temporal evolution of the pandemic under varying values of roll-out daily rates, vaccine efficacy, and a broad range of societal vaccine hesitancy/denial levels. The concept of herd immunity is questioned by studying future scenarios which involve different vaccination efforts and more infectious Covid-19 variants.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Twistons in a Sea of Magic

Magic angle twisted trilayer graphene (TTG) has recently emerged as a new platform to engineer strongly correlated flat bands. Here, we reveal the structural and electronic properties of TTG using low temperature scanning tunneling microscopy at twist angles for which superconductivity has been observed. Real trilayer samples deviate from their idealized structure due to a strong reconstruction of the moiré lattice, which locks layers into near-magic angle, mirror symmetric domains comparable in size to the superconducting coherence length. The price for this magic relaxation is the introduction of an array of localized twist angle faults, termed twistons. These novel, gate-tunable moiré defects offer a natural explanation for the superconducting dome observed in transport and provide an avenue to probe superconducting pairing mechanisms through disorder tuning.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Boosting the efficiency of ab initio electron-phonon coupling calculations through dual interpolation

The coupling between electrons and phonons in solids plays a central role in describing many phenomena, including superconductivity and thermoelecric transport. Calculations of this coupling are exceedingly demanding as they necessitate integrations over both the electron and phonon momenta, both of which span the Brillouin zone of the crystal, independently. We present here an ab initio method for efficiently calculating electron-phonon mediated transport properties by dramatically accelerating the computation of the double integrals with a dual interpolation technique that combines maximally localized Wannier functions with symmetry-adapted plane waves. The performance gain in relation to the current state-of-the-art Wannier-Fourier interpolation is approximately 2n_s \times M, where n_s is the number of crystal symmetry operations and M, a number in the range 5 - 60, governs the expansion in star functions. We demonstrate with several examples how our method performs some ab initio calculations involving electron-phonon interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Duality between atomic configurations and Bloch states in twistronic materials

The relative orientation (twist) of successive layers of stacked two-dimensional (2D) materials creates variations in the interlayer atomic registry. The variations often form a super lattice, called a moiré pattern, which can alter electronic properties. In this work we introduce a classification of the single-particle electronic structures that can occur in twisted stacks of 2D layers by characterizing them as "moiré molecules" or "moiré crystals". The molecules generate localized electronic states and moiré flat bands, while the crystals are sometimes unconventional and produce electronic banding in the configuration basis. The underpinning of this classification is the duality between interlayer configuration and monolayer Bloch momentum in moiré Hamiltonians. We apply this understanding to diagrams of local electron density in untwisted geometries to produce intuitive and quantitative predictions of twistronic properties. We provide a conceptual introduction to this framework through a one-dimensional model, and then apply it to 2D twisted bilayers of the semi-metal graphene, and of MoS$_2$, a representative material of the transition metal dichalcogenide (TMDC) family of semiconductors. This level of thorough understanding of twistronic phenomena is vital in the search for new material platforms for localized moiré electrons.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Effects of Structural Distortions on the Electronic Structure of T-type Transition Metal Dichalcogenides

Single-layer transition metal dichalcogenides (TMDCs) can adopt two distinct structures corresponding to different coordination of the metal atoms. TMDCs adopting the T-type structure exhibit a rich and diverse set of phenomena, including charge density waves (CDW) in a $\sqrt{13}\times\sqrt{13}$ supercell pattern in TaS$_2$ and TaSe$_2$, and a possible excitonic insulating phase in TiSe$_2$. These properties make the T-TMDCs desirable components of layered heterostructure devices. In order to predict the emergent properties of combinations of different layered materials, one needs simple and accurate models for the constituent layers which can take into account potential effects of lattice mismatch, relaxation, strain, and structural distortion. Previous studies have developed ab initio tight-binding Hamiltonians for H-type TMDCs [arXiv:1709.07510]. Here we extend this work to include T-type TMDCs. We demonstrate the capabilities of our model using three example systems: a 1-dimensional sinusoidal ripple, the 2$\times$2 CDW in TiSe$_2$, and the $\sqrt{13}\times\sqrt{13}$ CDW in TaS$_2$. Using the technique of band unfolding we compare the electronic structure of the distorted crystals to the pristine band structure and find excellent agreement with direct DFT calculations, provided the magnitude of the distortions remains in the linear regime.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Electronic structure calculations of twisted multi-layer graphene superlattices

Quantum confinement endows two-dimensional (2D) layered materials with exceptional physics and novel properties compared to their bulk counterparts. Although certain two- and few-layer configurations of graphene have been realized and studied, a systematic investigation of the properties of arbitrarily layered graphene assemblies is still lacking. We introduce theoretical concepts and methods for the processing of materials information, and as a case study, apply them to investigate the electronic structure of multi-layer graphene-based assemblies in a high-throughput fashion. We provide a critical discussion of patterns and trends in tight binding band structures and we identify specific layered assemblies using low-dispersion electronic bands as indicators of potentially interesting physics like strongly correlated behavior. A combination of data-driven models for visualization and prediction is used to intelligently explore the materials space. This work more generally aims to increase confidence in the combined use of physics-based and data-driven modeling for the systematic refinement of knowledge about 2D layered materials, with implications for the development of novel quantum devices.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

How carbon vacancies can affect the properties of group IV color centers in diamond: A study of thermodynamics and kinetics

Recently there has been much interest in using Group IV elements from the Periodic Table to fabricate and study X$V$ color centers in diamond where X = Si, Ge, Sn, or Pb and $V$ is a carbon vacancy. These Group IV color centers have a number of interesting spin and optical properties which could potentially make them better candidates than N$V^-$ centers for important applications in quantum computing and quantum information processing. Unfortunately, the very same ion implantation process that is required to create these X$V$ color centers in diamond necessarily also produces many carbon vacancies ($V_{\rm C}$) which can form complexes with these color centers ($V_{\rm C}-$X$V$) that can dramatically affect the properties of the isolated X$V$ color centers. The main focus of this work is to use density-functional theory (DFT) to study the thermodynamics and kinetics of the formation of these $V_{\rm C}-$X$V$ complexes and to suggest experimental ways to impede this process such as varying the Fermi level of the host diamond material through chemical doping or applying an external electrical bias. We also include a discussion of how the simple presence of many $V_{\rm C}$ can negatively impact the spin coherence times ($T_2$) of Group IV color centers through the presence of acoustic phonons.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

LAN -- A materials notation for 2D layered assemblies

Two-dimensional (2D) layered materials offer intriguing possibilities for novel physics and applications. Before any attempt at exploring the materials space in a systematic fashion, or combining insights from theory, computation and experiment, a formal description of information about an assembly of arbitrary composition is required. Here, we introduce a domain-generic notation that is used to describe the space of 2D layered materials from monolayers to twisted assemblies of arbitrary composition, existent or not-yet-fabricated. The notation corresponds to a theoretical materials concept of stepwise assembly of layered structures using a sequence of rotation, vertical stacking, and other operations on individual 2D layers. Its scope is demonstrated with a number of example structures using common single-layer materials as building blocks. This work overall aims to contribute to the systematic codification, capture and transfer of materials knowledge in the area of 2D layered materials.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Moiré metrology of energy landscapes in van der Waals heterostructures

The emerging field of twistronics, which harnesses the twist angle between two-dimensional materials, represents a promising route for the design of quantum materials, as the twist-angle-induced superlattices offer means to control topology and strong correlations. At the small twist limit, and particularly under strain, as atomic relaxation prevails, the emergent moiré superlattice encodes elusive insights into the local interlayer interaction. Here we introduce moiré metrology as a combined experiment-theory framework to probe the stacking energy landscape of bilayer structures at the 0.1 meV/atom scale, outperforming the gold-standard of quantum chemistry. Through studying the shapes of moiré domains with numerous nano-imaging techniques, and correlating with multi-scale modelling, we assess and refine first-principle models for the interlayer interaction. We document the prowess of moiré metrology for three representative twisted systems: bilayer graphene, double bilayer graphene and H-stacked $MoSe_2/WSe_2$. Moiré metrology establishes sought after experimental benchmarks for interlayer interaction, thus enabling accurate modelling of twisted multilayers.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Moiré of Moiré: Modeling Mechanical Relaxation in Incommensurate Trilayer van der Waals Heterostructures

The incommensurate stacking of multi-layered two-dimensional materials is a challenging problem from a theoretical perspective and an intriguing avenue for manipulating their physical properties. Here we present a multi-scale model to obtain the mechanical relaxation pattern of twisted trilayer van der Waals (vdW) heterostructures with two independent twist angles, a generally incommensurate system without a supercell description. We adopt the configuration space as a natural description of such incommensurate layered materials, based on the local environment of atomic positions, bypassing the need for commensurate approximations. To obtain the relaxation pattern, we perform energy minimization with respect to the relaxation displacement vectors. We use a continuum model in combination with the Generalized Stacking Fault energy to describe the interlayer coupling, obtained from first-principles calculations based on Density Functional Theory. We show that the relaxation patterns of twisted trilayer graphene and $\mathrm{WSe_2}$ are "moiré of moiré", as a result of the incommensurate coupling two bilayer moiré patterns. We also show that, in contrast to the symmetry-preserving in-plane relaxation in twisted bilayers, trilayer relaxation can break the two-fold rotational symmetry about the xy-plane when the two twist angles are equal.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The first 100 days: modeling the evolution of the COVID-19 pandemic

A simple analytical model for modeling the evolution of the 2020 COVID-19 pandemic is presented. The model is based on the numerical solution of the widely used Susceptible-Infectious-Removed (SIR) populations model for describing epidemics. We consider an expanded version of the original Kermack-McKendrick model, which includes a decaying value of the parameter $β$ (the effective contact rate) due to externally imposed conditions, to which we refer as the forced-SIR (FSIR) model. We introduce an approximate analytical solution to the differential equations that represent the FSIR model which gives very reasonable fits to real data for a number of countries over a period of 100 days (from the first onset of exponential increase, in China). The proposed model contains 3 adjustable parameters which are obtained by fitting actual data (up to April 28, 2020). We analyze these results to infer the physical meaning of the parameters involved. We use the model to make predictions about the total expected number of infections in each country as well as the date when the number of infections will have reached 99% of this total. We also compare key findings of the model with recently reported results on the high contagiousness and rapid spread of the disease.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Topological flat bands in frustrated kagome lattice CoSn

Electronic flat bands in momentum space, arising from strong localization of electrons in real space, are an ideal stage to realize strong correlation phenomena. In certain lattices with built-in geometrical frustration, electronic confinement and flat bands can naturally arise from the destructive interference of electronic hopping pathways. Such lattice-borne flat bands are often endowed with nontrivial topology if combined with spin-orbit coupling, while their experimental realization in condensed matter system has been elusive so far. Here, we report the direct observation of topological flat bands in the vicinity of the Fermi level in frustrated kagome system CoSn, using angle-resolved photoemission spectroscopy and band structure calculations. The flat band manifests itself as a dispersionless electronic excitation along the G-M high symmetry direction, with an order of magnitude lower bandwidth (below 150 meV) compared to the Dirac bands originating from the same orbitals. The frustration-driven nature of the flat band is directly confirmed by the real-space chiral d-orbital texture of the corresponding effective Wannier wave functions. Spin-orbit coupling opens a large gap of 80 meV at the quadratic band touching point between the Dirac and flat bands, endowing a nonzero Z2 topological invariant to the flat band in the two-dimensional Brillouin zone. Our observation of lattice-driven topological flat band opens a promising route to engineer novel emergent phases of matter at the crossroad between strong correlation physics and electronic topology.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Unifying framework for strong and fragile liquids via machine learning: a study of liquid silica

The fragility of a glassforming liquid characterizes how rapidly its relaxation dynamics slow down with cooling. The viscosity of strong liquids follows an Arrhenius law with a temperature-independent barrier height to rearrangements responsible for relaxation, whereas fragile liquids experience a much faster increase in their dynamics, suggesting a barrier height that increases with decreasing temperature. Strong glassformers are typically network glasses, while fragile glassformers are typically molecular or hard-sphere-like. As a result of these differences at the microscopic level, strong and fragile glassformers are usually treated separately from a theoretical point of view. Silica is the archetypal strong glassformer at low temperatures, but also exhibits a mysterious strong-to-fragile crossover at higher temperatures. Here we show that softness, a structure-based machine learned parameter that has previously been applied to fragile glassformers provides a useful description of model liquid silica in the strong and fragile regimes, and through the strong-to-fragile crossover. Just as for fragile glassformers, the relationship between softness and dynamics is invariant and Arrhenius in all regimes, but the average softness changes with temperature. The strong-to-fragile crossover in silica is not due to a sudden, qualitative change in structure, but can be explained by a simple Arrhenius form with a continuously and linearly changing local structure. Our results unify the study of liquid silica under a single simple conceptual picture.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Coexistence of ultraheavy and ultrarelativistic Dirac quasiparticles in sandwiched trilayer graphene

Electrons in quantum materials exhibiting coexistence of dispersionless (flat) bands piercing dispersive (steep) bands can give rise to strongly correlated phenomena, and are associated with unconventional superconductivity. It is known that in twisted trilayer graphene steep Dirac cones can coexist with band flattening, but the phenomenon is not stable under layer misalignments. Here we show that such a twisted sandwiched graphene (TSWG) -- a three-layer van der Waals heterostructure with a twisted middle layer -- can have very stable flat bands coexisting with Dirac cones near the Fermi energy when twisted to 1.5 degrees. These flat bands require a specific high-symmetry stacking order, and our atomistic calculations predict that TSWG always relaxes to it. Additionally, with external fields, we can control the relative energy offset between the Dirac cone vertex and the flat bands. Our work establishes twisted sandwiched graphene as a new platform for research into strongly interacting phases, and topological transport beyond Dirac and Weyl semimetals.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Creating Weyl nodes and controlling their energy by magnetization rotation

As they do not rely on the presence of any crystal symmetry, Weyl nodes are robust topological features of an electronic structure that can occur at any momentum and energy. Acting as sinks and sources of Berry curvature, Weyl nodes have been predicted to strongly affect the transverse electronic response, like in the anomalous Hall or Nernst effects. However, to observe large anomalous effects the Weyl nodes need to be close to or at the Fermi-level, which implies the band structure must be tuned by an external parameter, e.g. chemical doping or pressure. Here we show that in a ferromagnetic metal tuning of the Weyl node energy and momentum can be achieved by rotation of the magnetization. Taking Co$_3$Sn$_2$S$_2$ as an example, we use electronic structure calculations based on density-functional theory to show that not only new Weyl fermions can be created by canting the magnetization away from the easy axis, but also that the Weyl nodes can be driven exactly to the Fermi surface. We also show that the dynamics in energy and momentum of the Weyl nodes strongly affect the calculated anomalous Hall and Nernst conductivities.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Effects of Lithium Intercalation in Twisted Bilayer Graphene

We investigate the effects of lithium intercalation in twisted bilayers of graphene, using first-principles electronic structure calculations. To model this system we employ commensurate supercells that correspond to twist angles of 7.34$^\circ$ and 2.45$^\circ$. From the energetics of lithium absorption we demonstrate that for low Li concentration the intercalants cluster in the AA regions with double the density of a uniform distribution. The charge donated by the Li atoms to the graphene layers results in modifications to the band structure that can be qualitatively captured using a continuum model with modified interlayer couplings in a region of parameter space that has yet to be explored either experimentally or theoretically. Thus, the combination of intercalation and twisted layers simultaneously provides the means for spatial control over material properties and an additional knob with which to tune moiré physics in twisted bilayers of graphene, with potential applications ranging from energy storage and conversion to quantum information.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Imaginary-time time-dependent density functional theory and its application for robust convergence of electronic states

Reliable and robust convergence to the electronic ground state within density functional theory (DFT) Kohn-Sham (KS) calculations remains a thorny issue in many systems of interest. In such cases, charge sloshing can delay or completely hinder the convergence. Here, we use an approach based on transforming the time-dependent DFT equations to imaginary time, followed by imaginary-time evolution, as a reliable alternative to the self-consistent field (SCF) procedure for determining the KS ground state. We discuss the theoretical and technical aspects of this approach and show that the KS ground state should be expected to be the long-imaginary-time output of the evolution, independent of the exchange-correlation functional or the level of theory used to simulate the system. By maintaining self-consistency between the single-particle wavefunctions and the electronic density throughout the determination of the stationary state, our method avoids the typical difficulties encountered in SCF. To demonstrate dependability of our approach, we apply it to selected systems which struggle to converge with SCF schemes. In addition, through the van Leeuwen theorem, we affirm the physical meaningfulness of imaginary time TDDFT, justifying its use in certain topics of statistical mechanics such as in computing imaginary time path integrals.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Bounds on Nanoscale Nematicity in Single-Layer FeSe/SrTiO$_3$

We use scanning tunneling microscopy (STM) and quasiparticle interference (QPI) imaging to investigate the low-energy orbital texture of single-layer FeSe/SrTiO$_3$. We develop a $T$-matrix model of multi-orbital QPI to disentangle scattering intensities from Fe $3d_{xz}$ and $3d_{yz}$ bands, enabling the use of STM as a nanoscale detection tool of nematicity. By sampling multiple spatial regions of a single-layer FeSe/SrTiO$_3$ film, we quantitatively exclude static $xz/yz$ orbital ordering with domain size larger than $δr^2$ = 20 nm $\times$ 20 nm, $xz/yz$ Fermi wave vector difference larger than $δk$ = 0.014 $π$, and energy splitting larger than $δE$ = 3.5 meV. The lack of detectable ordering pinned around defects places qualitative constraints on models of fluctuating nematicity.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Dumbbell Defects in FeSe Films: A Scanning Tunneling Microscopy and First-Principles Investigation

The properties of iron-based superconductors (Fe-SCs) can be varied dramatically with the introduction of dopants and atomic defects. As a pressing example, FeSe, parent phase of the highest-$T_c$ Fe-SC, exhibits prevalent defects with atomic-scale "dumbbell" signatures as imaged by scanning tunneling microscopy (STM). These defects spoil superconductivity when their concentration exceeds 2.5%. Resolving their chemical identity is prerequisite to applications such as nanoscale patterning of superconducting/nonsuperconducting regions in FeSe, as well as fundamental questions such as the mechanism of superconductivity and the path by which the defects destroy it. We use STM and density functional theory to characterize and identify the dumbbell defects. In contrast to previous speculations about Se adsorbates or substitutions, we find that an Fe-site vacancy is the most energetically favorable defect in Se-rich conditions, and reproduces our observed STM signature. Our calculations shed light more generally on the nature of Se capping, the removal of Fe vacancies via annealing, and their ordering into a $\sqrt{5}$$\times$$\sqrt{5}$ superstructure in FeSe and related alkali-doped compounds.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Electric Field Tuning of Band Offsets in Transition Metal Dichalcogenides

We use first-principles calculations to investigate the band structure evolution of W$X_2$/Mo$X_2$ ($X$ = S, Se) heterobilayers under a perpendicular electric field. We characterize the extent to which the type-II band alignment in these compounds can be tuned or inverted electrostatically. Our results demonstrate two effects of stacking configuration. First, different stackings produce different net dipole moments, resulting in band offset variations that are larger than 0.1 eV. Second, based on symmetry constraints that depend on stacking, a perpendicular electric field may hybridize W$X_2$ and Mo$X_2$ bands that cross at the Brillouin Zone corner $K$. Our results suggest that external electric fields can be used to tune the physics of intralayer and interlayer excitons in heterobilayers of transition metal dichalcogenides.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Electronic Structure Theory of Weakly Interacting Bilayers

We derive electronic structure models for weakly interacting bilayers such as graphene-graphene and graphene-hexagonal boron nitride, based on density functional theory calculations followed by Wannier transformation of electronic states. These transferable interlayer coupling models can be applied to investigate the physics of bilayers with arbitrary translations and twists. The functional form, in addition to the dependence on the distance, includes the angular dependence that results from higher angular momentum components in the Wannier $p_z$ orbitals. We demonstrate the capabilities of the method by applying it to a rotated graphene bilayer, which produces the analytically predicted renormalization of the Fermi velocity, van Hove singularities in the density of states, and Moiré pattern of the electronic localization at small twist angles. We further extend the theory to obtain the effective couplings by integrating out neighboring layers. This approach is instrumental for the design of van der Walls heterostructures with desirable electronic features and transport properties and for the derivation of low-energy theories for graphene stacks, including proximity effects from other layers.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Epsilon-Near-Zero behavior from plasmonic Dirac point: theory and realization using two-dimensional materials

The electromagnetic response of a two-dimensional metal embedded in a periodic array of a dielectric host can give rise to a plasmonic Dirac point that emulates Epsilon-Near-Zero (ENZ) behavior. This theoretical result is extremely sensitive to tructural features like periodicity of the dielectric medium and thickness imperfections. We propose that such a device can actually be realized by using graphene as the 2D metal and materials like the layered semiconducting transition-metal dichalcogenides or hexagonal boron nitride as the dielectric host. We propose a systematic approach, in terms of design characteristics, for constructing metamaterials with linear, elliptical and hyperbolic dispersion relations which produce ENZ behavior, normal or negative diffraction.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Mapping Reactive Flow Patterns in Monolithic Nanoporous Catalysts

The development of high-efficiency porous catalyst membranes critically depends on our understanding of where the majority of the chemical conversions occur within the porous structure. This requires mapping of chemical reactions and mass transport inside the complex nano-scale architecture of porous catalyst membranes which is a multiscale problem in both the temporal and spatial domain. To address this problem, we developed a multi-scale mass transport computational framework based on the Lattice Boltzmann Method (LBM) that allows us to account for catalytic reactions at the gas-solid interface by introducing a new boundary condition. In good agreement with experiments, the simulations reveal that most catalytic reactions occur near the gas-flow facing side of the catalyst membrane if chemical reactions are fast compared to mass transport within the porous catalyst membrane.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Stability of metallic edges and Fermi-level pinning in transition-metal dichalcogenide nanoribbons

Nanoribbons of MoS$_2$ present a unique electronic structure that consists of a semiconducting bulk bounded by metallic edges; same holds for other Transition-Metal Dichalcogenides (TMDs) (Mo-,W-,S$_{2}$,Se$_{2}$). We perform first-principles calculations for TMD nanoribbons with reconstructed zig-zag metal terminated edges that contain chalcogen adatoms. All nanorobbons have possitive edge energies when the chemical potential of chalcogens is close to the energy of solids, and negative edge energies for high chemical potential. The reconstruction with two chalcogen adatoms is expected to be the most stable one. In all nanoribbons, a metallic phase is found near their edges, with the Fermi level of this metallic phase being lower than the Fermi level of the 2D material.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

\textit{Ab-initio} Tight-Binding Hamiltonian for Transition Metal Dichalcogenides

We present an accurate \textit{ab-initio} tight-binding hamiltonian for the transition-metal dichalcogenides, MoS$_2$, MoSe$_2$, WS$_2$, WSe$_2$, with a minimal basis (the \textit{d} orbitals for the metal atoms and \textit{p} orbitals for the chalcogen atoms) based on a transformation of the Kohn-Sham density function theory (DFT) hamiltonian to a basis of maximally localized Wannier functions (MLWF). The truncated tight-binding hamiltonian (TBH), with only on-site, first and partial second neighbor interactions, including spin-orbit coupling, provides a simple physical picture and the symmetry of the main band-structure features. Interlayer interactions between adjacent layers are modeled by transferable hopping terms between the chalcogen \textit{p} orbitals. The full-range tight-binding hamiltonian (FTBH) can be reduced to hybrid-orbital k $\cdot$ p effective hamiltonians near the band extrema that captures important low-energy excitations. These \textit{ab-initio} hamiltonians can serve as the starting point for applications to interacting many-body physics including optical transitions and Berry curvature of bands, of which we give some examples.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A structural approach to relaxation in glassy liquids

When a liquid freezes, a change in the local atomic structure marks the transition to the crystal. When a liquid is cooled to form a glass, however, no noticeable structural change marks the glass transition. Indeed, characteristic features of glassy dynamics that appear below an onset temperature, T_0, are qualitatively captured by mean field theory, which assumes uniform local structure at all temperatures. Even studies of more realistic systems have found only weak correlations between structure and dynamics. This raises the question: is structure important to glassy dynamics in three dimensions? Here, we answer this question affirmatively by using machine learning methods to identify a new field, that we call softness, which characterizes local structure and is strongly correlated with rearrangement dynamics. We find that the onset of glassy dynamics at T_0 is marked by the onset of correlations between softness (i.e. structure) and dynamics. Moreover, we use softness to construct a simple model of slow glassy relaxation that is in excellent agreement with our simulation results, showing that a theory of the evolution of softness in time would constitute a theory of glassy dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Ab-initio Studies of (Li$_{0.8}$Fe$_{0.2}$)OHFeSe Superconductors: Revealing the Dual Roles of Fe$_{0.2}$ in Structural Stability and Charge Transfer

The recently discovered (Li$_{0.8}$Fe$_{0.2}$)OHFeSe superconductor provides a new platform for exploiting the microscopic mechanisms of high-$T_c$ superconductivity in FeSe-derived systems. Using density functional theory calculations, we first show that substitution of Li by Fe not only significantly strengthens the attraction between the (Li$_{0.8}$Fe$_{0.2}$)OH spacing layers and the FeSe superconducting layers along the \emph{c} axis, but also minimizes the lattice mismatch between the two in the \emph{ab} plane, both favorable for stabilizing the overall structure. Next we explore the electron injection into FeSe from the spacing layers, and unambiguously identify the Fe$_{0.2}$ components to be the dominant atomic origin of the dramatically enhanced interlayer charge transfer. We further reveal that the system strongly favors collinear antiferromagnetic ordering in the FeSe layers, but the spacing layers can be either antiferromagnetic or ferromagnetic depending on the Fe$_{0.2}$ spatial distribution. Based on these understandings, we also predict (Li$_{0.8}$Co$_{0.2}$)OHFeSe to be structurally stable with even larger electron injection and potentially higher $T_c$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

High-Temperature Quantum Anomalous Hall Effect in n-p Codoped Topological Insulators

The quantum anomalous Hall effect (QAHE) is a fundamental quantum transport phenomenon that manifests as a quantized transverse conductance in response to a longitudinally applied electric field in the absence of an external magnetic field, and promises to have immense application potentials in future dissipation-less quantum electronics. Here we present a novel kinetic pathway to realize the QAHE at high temperatures by $n$-$p$ codoping of three-dimensional topological insulators. We provide proof-of-principle numerical demonstration of this approach using vanadium-iodine (V-I) codoped Sb$_2$Te$_3$ and demonstrate that, strikingly, even at low concentrations of $\sim$2\% V and $\sim$1\% I, the system exhibits a quantized Hall conductance, the tell-tale hallmark of QAHE, at temperatures of at least $\sim$ 50 Kelvin, which is three orders of magnitude higher than the typical temperatures at which it has been realized so far. The proposed approach is conceptually general and may shed new light in experimental realization of high-temperature QAHE.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Revealing the Empty-State Electronic Structure of Single-Unit-Cell FeSe/SrTiO$_{3}$

We use scanning tunneling spectroscopy to investigate the filled and empty electronic states of superconducting single-unit-cell FeSe deposited on SrTiO$_3$(001). We map the momentum-space band structure by combining quasiparticle interference imaging with decay length spectroscopy. In addition to quantifying the filled-state bands, we discover a $Γ$-centered electron pocket 75 meV above the Fermi energy. Our density functional theory calculations show the orbital nature of empty states at $Γ$ and suggest that the Se height is a key tuning parameter of their energies, with broad implications for electronic properties.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Theory of Graphene Raman Spectroscopy

Raman spectroscopy plays a key role in studies of graphene and related carbon systems. Graphene is perhaps the most promising material of recent times for many novel applications, including electronics. In this paper, the traditional and well established Kramers-Heisenberg-Dirac (KHD) Raman scattering theory (1925-1927) is extended to crystalline graphene for the first time. It demands different phonon production mechanisms and phonon energies than does the popular "double resonance" Raman scattering model. The latter has never been compared to KHD. Within KHD, phonons are produced instantly along with electrons and holes, in what we term an electron-hole-phonon triplet, which does not suffer Pauli blocking. A new mechanism for double phonon production we name "transition sliding" explains the brightness of the 2D mode and other overtones, as a result of linear (Dirac cone) electron dispersion. Direct evidence for sliding resides in hole doping experiments performed in 2011 \cite{chenCrommie}. Whole ranges of electronic transitions are permitted and may even constructively interfere for the same laser energy and phonon q, explaining the dispersion, bandwidth, and strength of many two phonon Raman bands. Graphene's entire Raman spectrum, including dispersive and fixed bands, missing bands not forbidden by symmetries, weak bands, overtone bands, Stokes anti-Stokes anomalies, individual bandwidths, trends with doping, and D-2D band spacing anomalies emerge naturally and directly in KHD theory.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Identifying structural flow defects in disordered solids using machine learning methods

We use machine learning methods on local structure to identify flow defects - or regions susceptible to rearrangement - in jammed and glassy systems. We apply this method successfully to two disparate systems: a two dimensional experimental realization of a granular pillar under compression, and a Lennard-Jones glass in both two and three dimensions above and below its glass transition temperature. We also identify characteristics of flow defects that differentiate them from the rest of the sample. Our results show it is possible to discern subtle structural features responsible for heterogeneous dynamics observed across a broad range of disordered materials.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Topological Proximity Effects in Graphene Nanoribbon Heterostructures

Topological insulators (TI) are bulk insulators that possess robust chiral conducting states along their interfaces with normal insulators. A tremendous research effort has recently been devoted to TI-based heterostructures, in which conventional proximity effects give rise to many exotic physical phenomena. Here we establish the potential existence of "topological proximity effects" at the interface of a topological graphene nanoribbon (GNR) and a normal GNR. Specifically, we show that the location of the topological edge states exhibits versatile tunability as a function of the interface orientation, as well as the strengths of the interface coupling and spin-orbit coupling in the normal GNR. For zigzag and bearded GNRs, the topological edge state can be tuned to be either at the interface or outer edge of the normal ribbon. For armchair GNR, the potential location of the topological edge state can be further enriched to be at the edge of or within the normal ribbon, at the interface, or diving into the topological GNR. We also discuss potential experimental realization of the predicted topological proximity effects, which may pave the way for integrating the salient functionality of TI and graphene in future device applications.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Approaching the Intrinsic Bandgap in Suspended High-Mobility Graphene Nanoribbons

We report electrical transport measurements on a suspended ultra-low-disorder graphene nanoribbon(GNR) with nearly atomically smooth edges that reveal a high mobility exceeding 3000 cm2 V-1 s-1 and an intrinsic band gap. The experimentally derived bandgap is in quantitative agreement with the results of our electronic-structure calculations on chiral GNRs with comparable width taking into account the electron-electron interactions, indicating that the origin of the bandgap in non-armchair GNRs is partially due to the magnetic zigzag edges.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Direct correlation of crystal structure and optical properties in wurtzite/zinc-blende GaAs nanowire heterostructures

A novel method for the direct correlation at the nanoscale of structural and optical properties of single GaAs nanowires is reported. Nanowires consisting of 100% wurtzite and nanowires presenting zinc-blende/wurtzite polytypism are investigated by photoluminescence spectroscopy and transmission electron microscopy. The photoluminescence of wurtzite GaAs is consistent with a band gap of 1.5 eV. In the polytypic nanowires, it is shown that the regions that are predominantly composed of either zinc-blende or wurtzite phase show photoluminescence emission close to the bulk GaAs band gap, while regions composed of a nonperiodic superlattice of wurtzite and zinc-blende phases exhibit a redshift of the photoluminescence spectra as low as 1.455 eV. The dimensions of the quantum heterostructures are correlated with the light emission, allowing us to determine the band alignment between these two crystalline phases. Our first-principles electronic structure calculations within density functional theory, employing a hybrid-exchange functional, predict band offsets and effective masses in good agreement with experimental results.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

On Graphene Hydrate

Using first-principles calculations, we show that the formation of carbohydrate directly from carbon and water is energetically favored when graphene membrane is subjected to aqueous environment with difference in chemical potential across the two sides. The resultant carbohydrate is two-dimensional, where the hydrogen atoms are exclusively attached on one side of graphene while the hydroxyl groups on the other side form a herringbone reconstruction that optimizes hydrogen bonding. We show that graphene undergos semi-metal-insulator transition upon hydration which is readily detectable from the significant shift in the vibration spectrum. The hydrate form of graphene suggests new applications for graphene in electronics, either deposited on a substrate or in solution.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Properties of nitrogen-vacancy centers in diamond: group theoretic approach

We present a procedure that makes use of group theory to analyze and predict the main properties of the negatively charged nitrogen-vacancy (NV) center in diamond. We focus on the relatively low temperatures limit where both the spin-spin and spin-orbit effects are important to consider. We demonstrate that group theory may be used to clarify several aspects of the NV structure, such as ordering of the singlets in the ($e^2$) electronic configuration, the spin-spin and the spin-orbit interactions in the ($ae$) electronic configuration. We also discuss how the optical selection rules and the response of the center to electric field can be used for spin-photon entanglement schemes. Our general formalism is applicable to a broad class of local defects in solids. The present results have important implications for applications in quantum information science and nanomagnetometry.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Numerical simulation of conformational variability in biopolymer translocation through wide nanopores

Numerical results on the translocation of long biopolymers through mid-sized and wide pores are presented. The simulations are based on a novel methodology which couples molecular motion to a mesoscopic fluid solvent. Thousands of events of long polymers (up to 8000 monomers) are monitored as they pass through nanopores. Comparison between the different pore sizes shows that wide pores can host a larger number of multiple biopolymer segments, as compared to smaller pores. The simulations provide clear evidence of folding quantization in the translocation process as the biopolymers undertake multi-folded configurations, characterized by a well-defined integer number of folds. Accordingly, the translocation time is no longer represented by a single-exponent power law dependence on the length, as it is the case for single-file translocation through narrow pores. The folding quantization increases with the biopolymer length, while the rate of translocated beads at each time step is linearly correlated to the number of resident beads in the pore. Finally, analysis of the statistics over the translocation work unravels the importance of the hydrodynamic interactions in the process.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Theory of Spin-Conserving Excitation of the $N-V^-$ Center in Diamond

The negatively charged nitrogen-vacancy defect ($N-V^-$ center) in diamond is an important atomic-scale structure that can be used as a qubit in quantum computing and as a marker in biomedical applications. Its usefulness relies on the ability to optically excite electrons between well-defined gap states, which requires clear and detailed understanding of the relevant states and excitation processes. Here we show that by using hybrid density-functional-theory calculations in a large supercell we can reproduce the zero-phonon line and the Stokes and anti-Stokes shifts, yielding a complete picture of the spin-conserving excitation of this defect.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Topological Frustration in Graphene Nanoflakes: Magnetic Order and Spin Logic Devices

Magnetic order in graphene-related structures can arise from size effects or from topological frustration. We introduce a rigorous classification scheme for the types of finite graphene structures (nano-flakes) which lead to large net spin or to antiferromagnetic coupling between groups of electron spins. Based on this scheme, we propose specific examples of structures that can serve as the fundamental (NOR and NAND) logic gates for the design of high-density ultra-fast spintronic devices. We demonstrate, using ab initio electronic structure calculations, that these gates can in principle operate at room temperature with very low and correctable error rates.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Multiscale modeling of biopolymer translocation through a nanopore

We employ a multiscale approach to model the translocation of biopolymers through nanometer size pores. Our computational scheme combines microscopic Langevin molecular dynamics (MD) with a mesoscopic lattice Boltzmann (LB) method for the solvent dynamics, explicitly taking into account the interactions of the molecule with the surrounding fluid. Both dynamical and statistical aspects of the translocation process were investigated, by simulating polymers of various initial configurations and lengths. For a representative molecule size, we explore the effects of important parameters that enter in the simulation, paying particular attention to the strength of the molecule-solvent coupling and of the external electric field which drives the translocation process. Finally, we explore the connection between the generic polymers modeled in the simulation and DNA, for which interesting recent experimental results are available.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

Electronic structure of overstretched DNA

Minuscule molecular forces can transform DNA into a structure that is elongated by more than half its original length. We demonstrate that this pronounced conformational transition is of relevance to ongoing experimental and theoretical efforts to characterize the conducting properties of DNA wires. We present quantum mechanical calculations for acidic, dry, poly(CG).poly(CG) DNA which has undergone elongation of up to 90 % relative to its natural length, along with a method for visualizing the effects of stretching on the electronic eigenstates. We find that overstretching leads to a drastic drop of the hopping matrix elements between localized occupied electronic states suggesting a dramatic decrease in the conductivity through holes.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

Environment Dependent Interatomic Potential for Bulk Silicon

We use recent theoretical advances to develop a new functional form for interatomic forces in bulk silicon. The theoretical results underlying the model include a novel analysis of elastic properties for the diamond and graphitic structures and inversions of ab initio cohesive energy curves. The interaction model includes two-body and three-body terms which depend on the local atomic environment through an effective coordination number. This formulation is able to capture successfully: (i) the energetics and elastic properties of the ground state diamond lattice; (ii) the covalent re-hybridization of undercoordinated atoms; (iii) and a smooth transition to metallic bonding for overcoordinated atoms. Because the essential features of chemical bonding in the bulk are built into the functional form, this model promises to be useful for describing interatomic forces in silicon bulk phases and defects. Although this functional form is remarkably realistic by usual standards, it contains a small number of fitting parameters and requires computational effort comparable to the most efficient existing models. In a companion paper, a complete parameterization of the model is given, and excellent performance for condensed phases and bulk defects is demonstrated.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

Interatomic Potential for Silicon Defects and Disordered Phases

We develop an empirical potential for silicon which represents a considerable improvement over existing models in describing local bonding for bulk defects and disordered phases. The model consists of two- and three-body interactions with theoretically motivated functional forms that capture chemical and physical trends as explained in a companion paper. The numerical parameters in the functional form are obtained by fitting to a set of ab initio results from quantum mechanical calculations based on density functional theory in the local density approximation, which include various bulk phases and defect structures. We test the potential by applying it to the relaxation of point defects, core properties of partial dislocations and the structure of disordered phases, none of which are included in the fitting procedure. For dislocations, our model makes predictions in excellent agreement with ab initio and tight-binding calculations. It is the only potential known to describe both the 30$^\circ$- and 90$^\circ$-partial dislocations in the glide set {111}. The structural and thermodynamic properties of the liquid and amorphous phases are also in good agreement with experimental and ab initio results. Our potential is the first capable of simulating a quench directly from the liquid to the amorphous phase, and the resulting amorphous structure is more realistic than with existing empirical preparation methods. These advances in transferability come with no extra computational cost, since force evaluation with our model is faster than with the popular potential of Stillinger-Weber, thus allowing reliable atomistic simulations of very large atomic systems.

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

Atomistic Aspects of Diffusion and Growth on the Si and Ge (111) Surfaces

The stability of interfaces and the mechanisms of thin film growth on semiconductors are issues of central importance in electronic devices. These issues can only be understood through detailed study of the relevant microscopic processes. Experimental studies are able to provide detailed, atomic scale information for model systems. Theoretical analysis of experimental results is essential in explaining certain surprising observations and in providing guidance for optimizing conditions and methods of growth. We review recent theoretical work on the diffusion of adatoms, the structure of adsorbate monolayers, and their implications for growth on the Si and Ge (111) surfaces. The theoretical analysis consists of first-principles calculations of the total-energy and entropy factors for stable, metastable and saddle-point configurations. These calculations are supplemented by Monte Carlo simulations of simple models that afford direct contact with experimental observations.

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

Theoretical Study of Cubic Structures Based on Fullerene Carbon Clusters: C$_{28}$C and (C$_{28})_{2}$

We study a new hypothetical form of solid carbon \csc, with a unit cell which is composed of the \cs \ fullerene cluster and an additional single carbon atom arranged in the zincblende structure. Using {\it ab initio} calculations, we show that this new form of solid carbon has lower energy than hyperdiamond, the recently proposed form composed of \cs \ units in the diamond structure. To understand the bonding character of of these cluster-based solids, we analyze the electronic structure of \csc \ and of hyperdiamond and compare them to the electronic states of crystalline cubic diamond.

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

Use of the Generalized Gradient Approximation in Pseudopotential Calculations of Solids

We present a study of the equilibrium properties of $sp$-bonded solids within the pseudopotential approach, employing recently proposed generalized gradient approximation (GGA) exchange correlation functionals. We analyze the effects of the gradient corrections on the behavior of the pseudopotentials and discuss possible approaches for constructing pseudopotentials self-consistently in the context of gradient corrected functionals. The calculated equilibrium properties of solids using the GGA functionals are compared to the ones obtained through the local density approximation (LDA) and to experimental data. A significant improvement over the LDA results is achieved with the use of the GGA functionals for cohesive energies. For the lattice constant, the same accuracy as in LDA can be obtained when the nonlinear coupling between core and valence electrons introduced by the exchange correlation functionals is properly taken into account. However, GGA functionals give bulk moduli that are too small compared to experiment.

0 citations0 reviewsNo code linkedOpen access