Source author record

Jeffrey B. Neaton

Jeffrey B. Neaton 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

39works
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

39 published item(s)

preprint2026arXiv

Discovery of a new weberite-type antiferroelectric: La3NbO7

Antiferroelectrics are antipolar materials which possess an electric field-induced phase transition to a polar, ferroelectric phase and offer significant potential for sensing/actuation and energy-storage applications. Known antiferroelectrics are relatively scarce and mainly based on a limited set of perovskite materials and their alloys (e.g., PbZrO$_3$, AgNbO$_3$, NaNbO$_3$). Here, a new family of lead-free, weberite-type antiferroelectrics, identified through a large-scale, first-principles computational search is introduced. The screening methodology, which connects lattice dynamics to antipolar distortions, predicted that La$_3$NbO$_7$ could exhibit antiferroelectricity. We confirm the prediction through the synthesis and characterization of epitaxial La$_3$NbO$_7$ thin films, which display the signature double hysteresis loops of an antiferroelectric material as well as clear evidence of an antipolar ground state structure from transmission electron microscopy. The antiferroelectricity in La$_3$NbO$_7$ is simpler than most known antiferroelectrics and can be explained by a Kittel-type mechanism involving the movement of niobium atoms in an oxygen octahedron through a single phonon mode which results in a smaller change in the volume during the field-induced phase transition. Additionally, it is found that La$_3$NbO$_7$ combines a high threshold field with a high breakdown field ($\approx$ 6MV/cm) - which opens up opportunities for energy-storage applications. This new weberite-type family of materials offers many opportunities to tune electrical and temperature response especially through substitutions on the rare-earth site. Ultimately, this work demonstrates a successful data-driven theory-to-experiment discovery of an entirely new family of antiferroelectrics and provides a blueprint for the future design of ferroic materials.

preprint2026arXiv

Enhanced Interlayer Coupling and Excitons in Twin-Stacked Two-Dimensional Magnetic CrSBr Bilayers

The degree of electronic coupling between individual layers in few-layer van der Waals heterostructures offers a route to engineer their magnetic, electronic, and optical functionalities. Using state-of-the-art first-principles calculations, we demonstrate that the electronic coupling between two monolayers of CrSBr, an anisotropic two-dimensional magnetic semiconductor, is highly nonlinear and nonmonotonic with respect to their relative twist angle, exhibiting a pronounced maximum at the twin-stacking configuration. The coupling strength scales with both the degree of overlap of Br orbitals adjacent to the van der Waals gap and the cosine of half of the interlayer spin angle. This enhanced interlayer electronic coupling gives rise to excitons delocalized across the two layers with a strong polarization dependence that reflects the details of the interlayer spin alignment. Our results reveal a sensitive interplay between twist angle, magnetism, and excitonic properties in twin-stacked CrSBr bilayers, and they establish twin stacking as an effective route to engineering interlayer coupling and optical response in anisotropic two-dimensional magnets with rectangular lattices.

preprint2022arXiv

An Optimally-Tuned Starting Point for Single-Shot $GW$ Calculations of Solids

The dependence of ab initio many-body perturbation theory within the $GW$ approximation on the eigensystem used in calculating quasiparticle corrections limits this method's predictive power. Here, we investigate the accuracy of the recently developed Wannier-localized optimally tuned screened range-separated hybrid (WOT-SRSH) functional as a generalized Kohn-Sham starting point for single-shot $GW$ ($G_0W_0$) calculations for a range of semiconductors and insulators. Comparison to calculations based on well-established functionals, namely PBE, PBE0, and HSE, as well as to self-consistent $GW$ schemes and to experiment, shows that band gaps computed via $G_0W_0$@WOT-SRSH have a level of precision and accuracy that is comparable to that of more advanced methods such as quasiparticle self-consistent $GW$ (QS$GW$) and eigenvalue self-consistent $GW$ (ev$GW$). We also find that $G_0W_0$@WOT-SRSH improves the description of states deeper in the valence band manifold. Finally, we show that $G_0W_0$@WOT-SRSH significantly reduces the sensitivity of computed band gaps to ambiguities in the underlying WOT-SRSH tuning procedure.

preprint2022arXiv

Density of States Prediction for Materials Discovery via Contrastive Learning from Probabilistic Embeddings

Machine learning for materials discovery has largely focused on predicting an individual scalar rather than multiple related properties, where spectral properties are an important example. Fundamental spectral properties include the phonon density of states (phDOS) and the electronic density of states (eDOS), which individually or collectively are the origins of a breadth of materials observables and functions. Building upon the success of graph attention networks for encoding crystalline materials, we introduce a probabilistic embedding generator specifically tailored to the prediction of spectral properties. Coupled with supervised contrastive learning, our materials-to-spectrum (Mat2Spec) model outperforms state-of-the-art methods for predicting ab initio phDOS and eDOS for crystalline materials. We demonstrate Mat2Spec's ability to identify eDOS gaps below the Fermi energy, validating predictions with ab initio calculations and thereby discovering candidate thermoelectrics and transparent conductors. Mat2Spec is an exemplar framework for predicting spectral properties of materials via strategically incorporated machine learning techniques.

preprint2021arXiv

Chemically-Localized Resonant Excitons in Silver-Pnictogen Halide Double Perovskites

Halide double perovskites with alternating silver and pnictogen cations are an emerging family of photoabsorber materials with robust stability and band gaps in the visible range. However, the nature of optical excitations in these systems is not yet well understood, limiting their utility. Here, we use ab initio many-body perturbation theory within the $GW$ approximation and the Bethe-Salpeter equation approach to calculate the electronic structure and optical excitations of the double perovskite series Cs$_2$AgBX$_6$, with B=Bi$^{3+}$, Sb$^{3+}$, X = Br$^-$, Cl$^-$. We find that these materials exhibit strongly localized resonant excitons with energies from 170 to 434 meV below the direct band gap. In contrast to lead-based perovskites, the Cs$_2$AgBX$_6$ excitons are computed to be non-hydrogenic, with anisotropic effective masses and sensitive to local field effects, a consequence of their chemical heterogeneity. Our calculations demonstrate the limitations of the Wannier-Mott and Elliott models for this class of double perovskites and contribute to a detailed atomistic understanding of their light-matter interactions.

preprint2021arXiv

Imaging gate-tunable Tomonaga-Luttinger liquids in 1H-MoSe$_2$ mirror twin boundaries

One-dimensional electron systems (1DESs) exhibit properties that are fundamentally different from higher-dimensional systems. For example, electron-electron interactions in 1DESs have been predicted to induce Tomonaga-Luttinger liquid behavior. Naturally-occurring grain boundaries in single-layer semiconducting transition metal dichalcogenides provide 1D conducting channels that have been proposed to host Tomonaga-Luttinger liquids, but charge density wave physics has also been suggested to explain their behavior. Clear identification of the electronic ground state of this system has been hampered by an inability to electrostatically gate such boundaries and thereby tune their charge carrier concentration. Here we present a scanning tunneling microscopy/spectroscopy study of gate-tunable mirror twin boundaries (MTBs) in single-layer 1H-MoSe$_2$ devices. Gating here enables STM spectroscopy to be performed for different MTB electron densities, thus allowing precise characterization of electron-electron interaction effects. Visualization of MTB electronic structure under these conditions allows unambiguous identification of collective density wave excitations having two distinct velocities, in quantitative agreement with the spin-charge separation predicted by finite-length Tomonaga-Luttinger-liquid theory.

preprint2020arXiv

Band gap renormalization, carrier mobilities, and the electron-phonon self-energy in crystalline naphthalene

Organic molecular crystals are expected to feature appreciable electron-phonon interactions that influence their electronic properties at zero and finite temperature. In this work, we report first-principles calculations and an analysis of the electron-phonon self-energy in naphthalene crystals. We compute the zero-point renormalization and temperature dependence of the fundamental band gap, and the resulting scattering lifetimes of electronic states near the valence- and conduction-band edges employing density functional theory. Further, our calculated phonon renormalization of the $GW$-corrected quasiparticle band structure predicts a fundamental band gap of 5 eV for naphthalene at room temperature, in good agreement with experiments. From our calculated phonon-induced electron lifetimes, we obtain the temperature-dependent mobilities of electrons and holes in good agreement with experimental measurements at room temperatures. Finally, we show that an approximate energy self-consistent computational scheme for the electron-phonon self-energy leads to the prediction of strong satellite bands in the electronic band structure. We find that a single calculation of the self-energy can reproduce the self-consistent results of the band gap renormalization and electrical mobilities for naphthalene, provided that the on-the-mass-shell approximation is used, i.e., if the self-energy is evaluated at the bare eigenvalues.

preprint2020arXiv

Half-magnetization plateau and the origin of threefold symmetry breaking in an electrically-switchable triangular antiferromagnet

We perform high-field magnetization measurements on the triangular lattice antiferromagnet Fe$_{1/3}$NbS$_2$. We observe a plateau in the magnetization centered at approximately half the saturation magnetization over a wide range of temperature and magnetic field. From density functional theory calculations, we determine a likely set of magnetic exchange constants. Incorporating these constants into a minimal Hamiltonian model of our material, we find that the plateau and of the $Z_3$ symmetry breaking ground state both arise from interplane and intraplane antiferromagnetic interactions acting in competition. These findings are pertinent to the magneto-electric properties of Fe$_{1/3}$NbS$_2$, which allow electrical switching of antiferromagnetic textures at relatively low current densities.

preprint2020arXiv

How Substitutional Point Defects in Two-Dimensional WS$_2$ Induce Charge Localization, Spin-Orbit Splitting, and Strain

Control of impurity concentrations in semiconducting materials is essential to device technology. Because of their intrinsic confinement, the properties of two-dimensional semiconductors such as transition metal dichalcogenides (TMDs) are more sensitive to defects than traditional bulk materials. The technological adoption of TMDs is dependent on the mitigation of deleterious defects and guided incorporation of functional foreign atoms. The first step towards impurity control is the identification of defects and assessment of their electronic properties. Here, we present a comprehensive study of point defects in monolayer tungsten disulfide (WS$_2$) grown by chemical vapor deposition (CVD) using scanning tunneling microscopy/spectroscopy, CO-tip noncontact atomic force microscopy, Kelvin probe force spectroscopy, density functional theory, and tight-binding calculations. We observe four different substitutional defects: chromium (Cr$_{\text{W}}$) and molybdenum (Mo$_{\text{W}}$) at a tungsten site, oxygen at sulfur sites in both bottom and top layers (O$_{\text{S}}$ top/bottom), as well as two negatively charged defects (CDs). Their electronic fingerprints unambiguously corroborate the defect assignment and reveal the presence or absence of in-gap defect states. The important role of charge localization, spin-orbit coupling, and strain for the formation of deep defect states observed at substitutional defects in WS$_2$ as reported here will guide future efforts of targeted defect engineering and doping of TMDs.

preprint2020arXiv

Identifying substitutional oxygen as a prolific point defect in monolayer transition metal dichalcogenides with experiment and theory

Chalcogen vacancies are considered to be the most abundant point defects in two-dimensional (2D) transition-metal dichalcogenide (TMD) semiconductors, and predicted to result in deep in-gap states (IGS). As a result, important features in the optical response of 2D-TMDs have typically been attributed to chalcogen vacancies, with indirect support from Transmission Electron Microscopy (TEM) and Scanning Tunneling Microscopy (STM) images. However, TEM imaging measurements do not provide direct access to the electronic structure of individual defects; and while Scanning Tunneling Spectroscopy (STS) is a direct probe of local electronic structure, the interpretation of the chemical nature of atomically-resolved STM images of point defects in 2D-TMDs can be ambiguous. As a result, the assignment of point defects as vacancies or substitutional atoms of different kinds in 2D-TMDs, and their influence on their electronic properties, has been inconsistent and lacks consensus. Here, we combine low-temperature non-contact atomic force microscopy (nc-AFM), STS, and state-of-the-art ab initio density functional theory (DFT) and GW calculations to determine both the structure and electronic properties of the most abundant individual chalcogen-site defects common to 2D-TMDs. Surprisingly, we observe no IGS for any of the chalcogen defects probed. Our results and analysis strongly suggest that the common chalcogen defects in our 2D-TMDs, prepared and measured in standard environments, are substitutional oxygen rather than vacancies.

preprint2020arXiv

The bright side of defects in MoS$_2$ and WS$_2$ and a generalizable chemical treatment protocol for defect passivation

Structural defects are widely regarded as detrimental to the optoelectronic properties of monolayer transition metal dichalcogenides, leading to concerted efforts to eliminate defects via improved materials growth or post-growth passivation. Here, using steady-state and ultrafast optical spectroscopy, supported by ab initio calculations, we demonstrate that sulfur vacancy defects act as exciton traps. Current chemical treatments do not passivate these sites, leading to decreased mobility and trap-limited photoluminescence. We present a generalizable treatment protocol based on the use of passivating agents such as thiols or sulfides in combination with a Lewis acid to passivate sulfur vacancies in monolayer MoS$_2$ and WS$_2$, increasing photoluminescence up to 275 fold, while maintaining mobilities. Our findings suggest a route for simple and rational defect engineering strategies, where the passivating agent varies the electronic properties, thereby allowing the design of new heterostructures.

preprint2019arXiv

Fermi-crossing Type-II Dirac fermions and topological surface states in NiTe2

Transition-metal dichalcogenides (TMDs) offer an ideal platform to experimentally realize Dirac fermions. However, typically these exotic quasiparticles are located far away from the Fermi level, limiting the contribution of Dirac-like carriers to the transport properties. Here we show that NiTe2 hosts both bulk Type-II Dirac points and topological surface states. The underlying mechanism is shared with other TMDs and based on the generic topological character of the Te p-orbital manifold. However, unique to NiTe2, a significant contribution of Ni d orbital states shifts the energy of the Type-II Dirac point close to the Fermi level. In addition, one of the topological surface states intersects the Fermi energy and exhibits a remarkably large spin splitting of 120 meV. Our results establish NiTe2 as an exciting candidate for next-generation spintronics devices.

preprint2019arXiv

Reproducibility in $G_0W_0$ Calculations for Solids

Ab initio many-body perturbation theory within the $GW$ approximation is a Green's function formalism widely used in the calculation of quasiparticle excitation energies of solids. In what has become an increasingly standard approach, Kohn-Sham eigenenergies, generated from a DFT calculation with a strategically-chosen exchange correlation functional ``starting point'', are used to construct $G$ and $W$, and then perturbatively corrected by the resultant $GW$ self-energy. In practice, there are several ways to construct the $GW$ self-energy, and these can lead to variations in predicted quasiparticle energies. For example, for ZnO and TiO$_2$, reported $GW$ fundamental gaps can vary by more than 1 eV. In this work, we address the convergence and key approximations in contemporary $G_0W_0$ calculations, including frequency-integration schemes and the treatment of the Coulomb divergence in the exact-exchange term. We study several systems,and compare three different $GW$ codes: BerkeleyGW, Abinit and Yambo. We demonstrate, for the first time, that the same quasiparticle energies for systems in the condensed phase can be obtained with different codes, and we provide a comprehensive assessment of implementations of the $GW$ approximation.

preprint2019arXiv

Resonant and Bound States of Charged Defects in Two-Dimensional Semiconductors

A detailed understanding of charged defects in two-dimensional semiconductors is needed for the development of ultrathin electronic devices. Here, we study negatively charged acceptor impurities in monolayer WS$_2$ using a combination of scanning tunnelling spectroscopy and large-scale atomistic electronic structure calculations. We observe several localized defect states of hydrogenic wave function character in the vicinity of the valence band edge. Some of these defect states are bound, while others are resonant. The resonant states result from the multi-valley valence band structure of WS$_2$, whereby localized states originating from the secondary valence band maximum at $Γ$ hybridize with continuum states from the primary valence band maximum at K/K$^{\prime}$. Resonant states have important consequences for electron transport as they can trap mobile carriers for several tens of picoseconds.

preprint2019arXiv

Superlattice-induced ferroelectricity in charge-ordered La$_{1/3}$Sr$_{2/3}$FeO$_{3}$

Charge-order-driven ferroelectrics are an emerging class of functional materials, distinct from conventional ferroelectrics, where electron-dominated switching can occur at high frequency. Despite their promise, only a few systems exhibiting this behavior have been experimentally realized thus far, motivating the need for new materials. Here, we use density functional theory to study the effect of artificial structuring on mixed-valence solid-solution La$_{1/3}$Sr$_{2/3}$FeO$_{3}$ (LSFO), a system well-studied experimentally. Our calculations show that A-site cation (111)-layered LSFO exhibits a ferroelectric charge-ordered phase in which inversion symmetry is broken by changing the registry of the charge order with respect to the superlattice layering. The phase is energetically degenerate with a ground-state centrosymmetric phase, and the computed switching polarization is 39 $μ$C/cm$^{2}$, a significant value arising from electron transfer between Fe ions. Our calculations reveal that artificial structuring of LSFO and other mixed valence oxides with robust charge ordering in the solid solution phase can lead to charge-order-induced ferroelectricity.

preprint2017arXiv

A new class of half-metallic ferromagnets from first principles

Half-metallic ferromagnetism (HMFM) occurs rarely in materials and yet offers great potential for spintronic devices. Recent experiments suggest a class of compounds with the `ThCr$_{2}$Si$_{2}$' (122) structure -- isostructural and containing elements common with Fe pnictide-based superconductors -- can exhibit HMFM. Here we use $ab$ $initio$ density-functional theory calculations to understand the onset of half-metallicity in this family of materials and explain the appearance of ferromagnetism at a quantum critical point. We also predict new candidate materials with HMFM and high Curie temperatures through A-site alloying.

preprint2016arXiv

Effects of quantum confinement on excited state properties of SrTiO$_3$ from ab initio many-body perturbation theory

The Ruddlesden-Popper (RP) homologous series Sr$_{n+1}$Ti$_{n}$O$_{3n+1}$ provides a useful template for the study and control of the effects of dimensionality and quantum confinement on the excited state properties of the complex oxide SrTiO$_3$. We use ab initio many-body perturbation theory within the $GW$ approximation and the Bethe-Salpeter equation approach to calculate quasiparticle energies and absorption spectrum of Sr$_{n+1}$Ti$_{n}$O$_{3n+1}$ for $n=1-5$ and $\infty$. Our computed direct and indirect optical gaps are in excellent agreement with spectroscopic measurements. The calculated optical spectra reproduce the main experimental features and reveal excitonic structure near the gap edge. We find that electron-hole interactions are important across the series, leading to significant exciton binding energies that increase for small $n$ and reach a value of 330~meV for $n=1$, a trend attributed to increased quantum confinement. We find that the lowest-energy singlet exciton of Sr$_2$TiO$_4$ ($n=1$) localizes in the 2D plane defined by the TiO$_2$ layer, and explain the origin of its localization.

preprint2016arXiv

Electric field and strain induced Rashba effect in hybrid halide perovskites

Using first principles density functional theory calculations, we show how Rashba-type energy band splitting in the hybrid organic-inorganic halide perovskites APbX$_3$ (A=CH$_3$NH$_3^+$, CH(NH$_2$)$_2^+$, Cs$^+$ and X=I, Br) can be tuned and enhanced with electric fields and anisotropic strain. In particular, we demonstrate that the magnitude of the Rashba splitting of tetragonal (CH$_3$NH$_3$)PbI$_3$ grows with increasing macroscopic alignment of the organic cations and electric polarization, indicating appreciable tunability with experimentally-feasible applied fields, even at room temperature. Further, we quantify the degree to which this effect can be tuned via chemical substitution at the A and X sites, which alters amplitudes of different polar distortion patterns of the inorganic PbX$_3$ cage that directly impact Rashba splitting. In addition, we predict that polar phases of CsPbI$_3$ and (CH$_3$NH$_3$)PbI$_3$ with $R3c$ symmetry possessing considerable Rashba splitting might be accessible at room temperature via anisotropic strain induced by epitaxy, even in the absence of electric fields.

preprint2016arXiv

Evaluating the GW Approximation with CCSD(T) for Charged Excitations Across the Oligoacenes

Charged excitations of the oligoacene family of molecules, relevant for astrophysics and technological applications, are widely studied and therefore provide an excellent system for benchmarking theoretical methods. In this work, we evaluate the performance of many-body perturbation theory within the GW approximation relative to new high-quality CCSD(T) reference data for charged excitations of the acenes. We compare GW calculations with a number of hybrid density functional theory starting points and with eigenvalue self-consistency. Special focus is given to elucidating the trend of GW-predicted excitations with molecule length increasing from benzene to hexacene. We find that GW calculations with starting points based on an optimally tuned range-separated hybrid (OTRSH) density functional and eigenvalue self-consistency can yield quantitative ionization potentials for the acenes. However, for larger acenes, the predicted electron affinities can deviate considerably from reference values. Our work paves the way for predictive and cost-effective GW calculations of charged excitations of molecules and identifies certain limitations of current GW methods used in practice for larger molecules.

preprint2016arXiv

Interplay of bias-driven charging and the vibrational Stark effect in molecular junctions

We observe large, reversible, bias driven changes in the vibrational energies of PCBM, based on simultaneous transport and surface-enhanced Raman spectroscopy (SERS) measurements on PCBM-gold junctions. A combination of linear and quadratic shifts in vibrational energies with voltage is analyzed and compared with similar measurements involving C60-gold junctions. A theoretical model based on density functional theory (DFT) calculations suggests that both a vibrational Stark effect and bias-induced charging of the junction contribute to the shifts in vibrational energies. In the PCBM case, a linear vibrational Stark effect is observed due to the permanent electric dipole moment of PCBM. The vibrational Stark shifts shown here for PCBM junctions are comparable to or larger than the charging effects that dominate in C60 junctions.

preprint2016arXiv

Observation of charge density wave order in 1D mirror twin boundaries of single-layer MoSe2

Properties of two-dimensional transition metal dichalcogenides are highly sensitive to the presence of defects in the crystal structure. A detailed understanding of defect structure may lead to control of material properties through defect engineering. Here we provide direct evidence for the existence of isolated, one-dimensional charge density waves at mirror twin boundaries in single-layer MoSe2. Our low-temperature scanning tunneling microscopy/spectroscopy measurements reveal a substantial bandgap of 60 - 140 meV opening at the Fermi level in the otherwise one dimensional metallic structure. We find an energy-dependent periodic modulation in the density of states along the mirror twin boundary, with a wavelength of approximately three lattice constants. The modulations in the density of states above and below the Fermi level are spatially out of phase, consistent with charge density wave order. In addition to the electronic characterization, we determine the atomic structure and bonding configuration of the one-dimensional mirror twin boundary by means of high-resolution non-contact atomic force microscopy. Density functional theory calculations reproduce both the gap opening and the modulations of the density of states.

preprint2016arXiv

Structural and excited-state properties of oligoacene crystals from first principles

Molecular crystals are a prototypical class of van der Waals (vdW) bound organic materials with excited state properties relevant for optoelectronics applications. Predicting the structure and excited state properties of molecular crystals presents a challenge for electronic structure theory, as standard approximations to density functional theory (DFT) do not capture long range vdW dispersion interactions and do not yield excited state properties. In this work, we use a combination of DFT including vdW forces) using both non local correlation functionals and pair wise correction methods (together with many body perturbation theory (MBPT) to study the geometry and excited states, respectively, of the entire series of oligoacene crystals, from benzene to hexacene. We find that vdW methods can predict lattice constants within 1 percent of the experimental measurements, on par with the previously reported accuracy of pairwise approximations for the same systems. We further find that excitation energies are sensitive to geometry, but if optimized geometries are used MBPT can yield excited state properties within a few tenths of an eV from experiment. We elucidate trends in MBPT computed charged and neutral excitation energies across the acene series and discuss the role of common approximations used in MBPT.

preprint2015arXiv

A systematic benchmark of the ab initio Bethe-Salpeter equation approach for low-lying optical excitations of small organic molecules

The predictive power of the ab initio Bethe-Salpeter equation (BSE) approach, rigorously based on many-body Green's function theory but incorporating information from density functional theory, has already been demonstrated for the optical gaps and spectra of solid-state systems. Interest in photoactive hybrid organic/inorganic systems has recently increased, and so has the use of the BSE for computing neutral excitations of organic molecules. However, no systematic benchmarks of the BSE for neutral electronic excitations of organic molecules exist. Here, we study the performance of the BSE for the 28 small molecules in Thiel's widely-used time-dependent density functional theory benchmark set [M. Schreiber et al. J. Chem. Phys. 128, 134110 (2008)]. We observe that the BSE produces results that depend critically on the mean-field starting point employed in the perturbative approach. We find that this starting point dependence is mainly introduced through the quasiparticle energies obtained at the intermediate GW step, and that with a judicious choice of starting mean-field, singlet excitation energies obtained from BSE are in excellent quantitative agreement with higher-level wavefunction methods. The quality of the triplet excitations is slightly less satisfactory.

preprint2015arXiv

Solid-state optical absorption from optimally-tuned time-dependent range-separated hybrid density functional theory

We present a framework for obtaining reliable solid-state charge and optical excitations and spectra from optimally-tuned range-separated hybrid density functional theory. The approach, which is fully couched within the formal framework of generalized Kohn-Sham theory, allows for accurate prediction of exciton binding energies. We demonstrate our approach through first principles calculations of one- and two-particle excitations in pentacene, a molecular semiconducting crystal, where our work is in excellent agreement with experiments and prior computations. We further show that with one adjustable parameter, set to produce an accurate bandgap, this method accurately predicts band structures and optical spectra of silicon and lithium flouride, prototypical covalent and ionic solids. Our findings indicate that for a broad range of extended bulk systems, this method may provide a computationally inexpensive alternative to many-body perturbation theory, opening the door to studies of materials of increasing size and complexity.

preprint2014arXiv

Adsorption-Induced Solvent-Based Electrostatic Gating of Charge Transport through Molecular Junctions

Recent experiments have shown that transport properties of molecular-scale devices can be reversibly altered by the surrounding solvent. Here, we use a combination of first-principles calculations and experiment to explain this change in transport properties through a shift in the local electrostatic potential at the junction caused by nearby conducting and solvent molecules chemically bound to the electrodes. This effect is found to alter the conductance of bipyridine - gold junctions by more than 50%. Moreover, we develop a general electrostatic model that quantitatively predicts the relationship between conductance and the binding energies and dipoles of the solvent and conducting molecules. Our work shows that solvent induced effects are a viable route for controlling charge and energy transport at molecular-scale interfaces.

preprint2014arXiv

Effects of self-consistency and plasmon-pole models on GW calculations for closed-shell molecules

We present theoretical calculations of quasiparticle energies in closed-shell molecules using the GW method. We compare three different approaches: a full-frequency $G_0W_0$ (FF-$G_0W_0$) method with density functional theory (DFT-PBE) used as a starting mean field; a full-frequency $GW_0$ (FF-$GW_0$) method where the interacting Green's function is approximated by replacing the DFT energies with self-consistent quasiparticle energies or Hartree-Fock energies; and a $G_0W_0$ method with a Hybertsen-Louie generalized plasmon-pole model (HL GPP-$G_0W_0$). While the latter two methods lead to good agreement with experimental ionization potentials and electron affinities for methane, ozone, and beryllium oxide molecules, FF-$G_0W_0$ results can differ by more than one electron volt from experiment. We trace this failure of the FF-$G_0W_0$ method to the occurrence of incorrect self-energy poles describing shake-up processes in the vicinity of the quasiparticle energies.

preprint2014arXiv

Elastic Properties of Chemical-Vapor-Deposited Monolayer MoS2, WS2, and Their Bilayer Heterostructures

Elastic properties of materials are an important factor in their integration in applications. Chemical vapor deposited (CVD) monolayer semiconductors are proposed as key components in industrial-scale flexible devices and building blocks of 2D van der Waals heterostructures. However, their mechanical and elastic properties have not been fully characterized. Here we report high 2D elastic moduli of CVD monolayer MoS2 and WS2 (~ 170 N/m), which is very close to the value of exfoliated MoS2 monolayers and almost half the value of the strongest material, graphene. The 2D moduli of their bilayer heterostructures are lower than the sum of 2D modulus of each layer, but comparable to the corresponding bilayer homostructure, implying similar interactions between the hetero monolayers as between homo monolayers. These results not only provide deep insight to understanding interlayer interactions in 2D van der Waals structures, but also potentially allow engineering of their elastic properties as desired.

preprint2014arXiv

Energy-dependent resonance broadening in symmetric and asymmetric molecular junctions from an ab initio non-equilibrium Green's function approach

The electronic structure of organic-inorganic interfaces often feature resonances originating from discrete molecular orbitals coupled to continuum lead states. An example are molecular junctions, individual molecules bridging electrodes, where the shape and peak energy of such resonances dictate junction conductance, thermopower, I-V characteristics and related transport properties. In molecular junctions where off-resonance coherent tunneling dominates transport, resonance peaks in the transmission function are often assumed to be Lorentzian functions with an energy-independent broadening parameter $Γ$. Here we define a new energy-dependent resonance broadening function, $Γ(E)$, based on diagonalization of non-Hermitian matrices, which can describe resonances of a more complex, non-Lorentzian nature and can be decomposed into components associated with the left and right lead, respectively. We compute this quantity via an \emph{ab initio} non-equilibrium Green's function approach based on density functional theory for both symmetric and asymmetric molecular junctions, and show that our definition of $Γ(E)$, when combined with Breit-Wigner formula, reproduces the transmission calculated from DFT-NEGF. Through a series of examples, we illustrate how this approach can shed new light on experiments and understanding of junction transport properties in terms of molecular orbitals.

preprint2014arXiv

Numerical integration for ab initio many-electron self energy calculations within the GW approximation

We present a numerical integration scheme for evaluating the convolution of a Green's function with a screened Coulomb potential on the real axis in the GW approximation of the self energy. Our scheme takes the zero broadening limit in Green's function first, replaces the numerator of the integrand with a piecewise polynomial approximation, and performs principal value integration on subintervals analytically. We give the error bound of our numerical integration scheme and show by numerical examples that it is more reliable and accurate than the standard quadrature rules such as the composite trapezoidal rule. We also discuss the benefit of using different self energy expressions to perform the numerical convolution at different frequencies.

preprint2013arXiv

Competing thermodynamic and dynamic factors select molecular assemblies on a gold surface

Controlling the self-assembly of surface-adsorbed molecules into nanostructures requires understanding physical mechanisms that act across multiple length and time scales. By combining scanning tunneling microscopy with hierarchical ab initio and statistical mechanical modeling of 1,4-substituted benzenediamine (BDA) molecules adsorbed on a gold (111) surface, we demonstrate that apparently simple nanostructures are selected by a subtle competition of thermodynamics and dynamics. Of the collection of possible BDA nanostructures mechanically stabilized by hydrogen bonding, the interplay of intermolecular forces, surface modulation, and assembly dynamics select at low temperature a particular subset: low free energy oriented linear chains of monomers, and high free energy branched chains.

preprint2012arXiv

Design of Low Band Gap Double Perovskites from First Principles

Using density functional theory (DFT)-based calculations, we propose a family of metastable, as-yet unmade V5+ and Cr6+ double perovskite compounds with low band gaps spanning much of the visible region of the solar spectrum. Through analysis of a related set of measured optical gaps of d0 ABO3 perovskites and A2B'BO6 double perovskites, an ad hoc procedure is developed to correct DFT and many-body perturbation theory gaps, bringing them into quantitative agreement with experiment for measured compounds, and predicting that V5+ and Cr6+ double perovskites would have gaps ranging from approximately 1.1-2.4 eV, significantly lower than previous materials studied in this class. DFT calculations also establish that these V5+ and Cr6+ compounds are likely able to be synthesized, either in bulk form or as epitaxial thin films. These compounds would comprise a new class of semiconducting double perovskites for potential use in solar energy conversion and other optoelectronic applications.

preprint2012arXiv

Molecular Adsorption on Metal Surfaces with a van der Waals Density Functional

The adsorption of 1,4-benzenediamine (BDA) on the Au(111) surface and azobenzene on the Ag(111) surface is investigated using density functional theory (DFT) with a non-local density functional (vdW-DF) and a semi-local Perdew-Burke-Ernzerhof (PBE) functional. For BDA on Au(111), the inclusion of London dispersion interactions not only dramatically enhances the molecule-substrate binding, resulting in adsorption energies consistent with experimental results, but also significantly alters the BDA binding geometry. For azobenzene on Ag(111), the vdW-DF produces superior adsorption energies compared to those obtained with other dispersion corrected DFT approaches. These results provide evidence for the applicability of the vdW-DF method and serves as a practical benchmark for the investigation of molecules adsorbed on noble metal surfaces.

preprint2012arXiv

Quantitative Current-Voltage Characteristics in Molecular Junctions from First Principles

Using self-energy-corrected density functional theory (DFT) and a coherent scattering-state approach, we explain current-voltage (IV) measurements of four pyridine-Au and amine-Au linked molecular junctions with quantitative accuracy. Parameter-free many-electron self-energy corrections to DFT Kohn-Sham eigenvalues are demonstrated to lead to excellent agreement with experiments at finite bias, improving upon order-of-magnitude errors in currents obtained with standard DFT approaches. We further propose an approximate route for prediction of quantitative IV characteristics for both symmetric and asymmetric molecular junctions based on linear response theory and knowledge of the Stark shifts of junction resonance energies. Our work demonstrates that a quantitative, computationally inexpensive description of coherent transport in molecular junctions is readily achievable, enabling new understanding and control of charge transport properties of molecular-scale interfaces at large bias voltages.

preprint2012arXiv

Quantitative Molecular Orbital Energies within a $G_0W_0$ Approximation

Using many-body perturbation theory within the $G_0W_0$ approximation, we explore routes for computing the ionization potential (IP), electron affinity (EA), and fundamental gap of three gas-phase molecules -- benzene, thiophene, and (1,4) diamino-benzene -- and compare with experiments. We examine the dependence of the IP on the number of unoccupied states used to build the dielectric function and the self energy, as well as the dielectric function plane-wave cutoff. We find that with an effective completion strategy for approximating the unoccupied subspace, and a converged dielectric function kinetic energy cutoff, the computed IPs and EAs are in excellent quantitative agreement with available experiment (within 0.2 eV), indicating that a one-shot $G_0W_0$ approach can be very accurate for calculating addition/removal energies of small organic molecules. Our results indicate that a sufficient dielectric function kinetic energy cutoff may be the limiting step for a wide application of $G_0W_0$ to larger organic systems.

preprint2012arXiv

Quasiparticle and Optical Spectroscopy of Organic Semiconductors Pentacene and PTCDA from First Principles

The broad use of organic semiconductors for optoelectronic applications relies on quantitative understanding and control of their spectroscopic properties. Of paramount importance are the transport gap - the difference between ionization potential and electron affinity - and the exciton binding energy - inferred from the difference between the transport and optical absorption gaps. Transport gaps are commonly established via photoemission and inverse photoemission spectroscopy (PES/IPES). However, PES and IPES are surface-sensitive, average over a dynamic lattice, and are subject to extrinsic effects, leading to significant uncertainty in gaps. Here, we use density functional theory and many-body perturbation theory to calculate the spectroscopic properties of two prototypical organic semiconductors, pentacene and 3,4,9,10-perylene tetracarboxylic dianhydride (PTCDA), quantitatively comparing with measured PES, IPES, and optical absorption spectra. For bulk pentacene and PTCDA, the computed transport gaps are 2.4 and 3.0 eV, and optical gaps are 1.7 and 2.1 eV, respectively. Computed bulk quasiparticle spectra are in excellent agreement with surface-sensitive photoemission measurements over several eV only if the measured gap is reduced by 0.6 eV for pentacene, and 0.6-0.9 eV for PTCDA. We attribute this redshift to several physical effects, including incomplete charge screening at the surface, static and dynamical disorder, and experimental resolution. Optical gaps are in excellent agreement with experiment, with solid-state exciton binding energies of ~0.5 eV for both systems; for pentacene, the exciton is delocalized over several molecules and exhibits significant charge transfer character. Our parameter-free calculations provide new interpretation of spectroscopic properties of organic semiconductors critical to optoelectronics.

preprint2012arXiv

Quasiparticle spectra from a non-empirical optimally-tuned range-separated hybrid density functional

We present a method for obtaining outer valence quasiparticle excitation energies from a DFT-based calculation, with accuracy that is comparable to that of many-body perturbation theory within the GW approximation. The approach uses a range-separated hybrid density functional, with asymptotically exact and short-range fractional Fock exchange. The functional contains two parameters - the range separation and the short-range Fock fraction. Both are determined non-empirically, per system, based on satisfaction of exact physical constraints for the ionization potential and many-electron self-interaction, respectively. The accuracy of the method is demonstrated on four important benchmark organic molecules: perylene, pentacene, 3,4,9,10-perylene-tetracarboxylic-dianydride (PTCDA) and 1,4,5,8-naphthalene-tetracarboxylic dianhydride (NTCDA). We envision that for finite systems the approach could provide an inexpensive alternative to GW, opening the door to the study of presently out of reach large-scale systems.

preprint2011arXiv

Band Gap and Edge Engineering via Ferroic Distortion and Anisotropic Strain: The Case of SrTiO$_{3}$

The effects of ferroic distortion and biaxial strain on the band gap and band edges of SrTiO$_{3}$ (STO) are calculated using density functional theory and many-body perturbation theory. Anisotropic strains are shown to reduce the gap by breaking degeneracies at the band edges. Ferroic distortions are shown to widen the gap by allowing new band edge orbital mixings. Compressive biaxial strains raise band edge energies, while tensile strains lower them. To reduce the STO gap, one must lower the symmetry from cubic while suppressing ferroic distortions. Our calculations indicate that for engineered orientation of the growth direction along [111], the STO gap can be controllably and considerably reduced at room temperature.

preprint2011arXiv

Electronic energy level alignment at metal-molecule interfaces with a GW approach

Using density functional theory and many-body perturbation theory within a GW approximation, we calculate the electronic structure of a metal-molecule interface consisting of benzene diamine (BDA) adsorbed on Au(111). Through direct comparison with photoemission data, we show that a conventional G$_0$W$_0$ approach can underestimate the energy of the adsorbed molecular resonance relative to the Au Fermi level by up to 0.8 eV. The source of this discrepancy is twofold: a 0.7 eV underestimate of the gas phase ionization energy (IE), and a 0.2 eV overestimate of the Au work function. Refinements to self-energy calculations within the GW framework that account for deviations in both the Au work function and BDA gas-phase IE can result in an interfacial electronic level alignment in quantitative agreement with experiment.

preprint2004arXiv

First-principles study of symmetry lowering in relaxed BaTiO3/SrTiO3 superlattices

The crystal structure and local spontaneous polarization of (BaTiO3)m/(SrTiO3)n superlattices is calculated using a first-principles density functional theory method. The in-plane lattice constant is 1% larger than the SrTiO3 substrate to imitate the relaxed superlattice structure and the symmetry is lowered to monoclinic space group Cm which allows polarization to develop along the [110] and [001] directions. The polarization component in the [110] direction is found to develop only in the SrTiO3 layers and falls to zero in the BaTiO3 layers, whereas the polarization in the [001] direction is approximately uniform throughout the superlattice. These findings are consistent with recent experimental data and first-principles results for epitaxially strained BT and ST.