-
Electronic structure and optical properties of halide double perovskites from a Wannier-localized optimally-tuned screened range-separated hybrid functional
Authors:
Francisca Sagredo,
Stephen E. Gant,
Guy Ohad,
Jonah B. Haber,
Marina R. Filip,
Leeor Kronik,
Jeffrey B. Neaton
Abstract:
Halide double perovskites are a chemically-diverse and growing class of compound semiconductors that are promising for optoelectronic applications. However, the prediction of their fundamental gaps and optical properties with density functional theory (DFT) and {\it ab initio} many-body perturbation theory has been a significant challenge. Recently, a nonempirical Wannier-localized optimally-tuned…
▽ More
Halide double perovskites are a chemically-diverse and growing class of compound semiconductors that are promising for optoelectronic applications. However, the prediction of their fundamental gaps and optical properties with density functional theory (DFT) and {\it ab initio} many-body perturbation theory has been a significant challenge. Recently, a nonempirical Wannier-localized optimally-tuned screened range-separated hybrid (WOT-SRSH) functional has been shown to accurately produce the fundamental band gaps of a wide set of semiconductors and insulators, including lead halide perovskites. Here we apply the WOT-SRSH functional to five halide double perovskites, and compare the results with those obtained from other known functionals and previous $GW$ calculations. We also use the approach as a starting point for $GW$ calculations and we compute the band structures and optical absorption spectrum for Cs\textsubscript{2}Ag{Bi}Br\textsubscript{6}, using both time-dependent DFT and the $GW$-Bethe-Salpeter equation approach. We show that the WOT-SRSH functional leads to accurate fundamental and optical band gaps, as well as optical absorption spectra, consistent with spectroscopic measurements, thereby establishing WOT-SRSH as a viable method for the accurate prediction of optoelectronic properties of halide double perovskites.
△ Less
Submitted 7 August, 2024;
originally announced August 2024.
-
Polarization Saturation in Multi-layered Interfacial Ferroelectrics
Authors:
Wei Cao,
Swarup Deb,
Maayan Vizner Stern,
Noam Raab,
Michael Urbakh,
Oded Hod,
Leeor Kronik,
Moshe Ben Shalom
Abstract:
Van der Waals (vdW) polytypes of broken inversion and mirror symmetries were recently shown to exhibit switchable electric polarization even at the ultimate two-layer thin limit. Their out-of-plane polarization was found to accumulate in a ladder-like fashion with each successive layer, offering 2D building blocks for the bottom-up construction of 3D ferroelectrics. Here, we demonstrate experiment…
▽ More
Van der Waals (vdW) polytypes of broken inversion and mirror symmetries were recently shown to exhibit switchable electric polarization even at the ultimate two-layer thin limit. Their out-of-plane polarization was found to accumulate in a ladder-like fashion with each successive layer, offering 2D building blocks for the bottom-up construction of 3D ferroelectrics. Here, we demonstrate experimentally that beyond a critical stack thickness, the accumulated polarization in rhombohedral polytypes of molybdenum disulfide (r-MoS2) saturates. The underlying saturation mechanism, deciphered via DFT and self-consistent Poisson-Schrödinger calculations, point to a purely electronic redistribution involving: (1) polarization-induced bandgap closure that allows for cross-stack charge transfer and the emergence of free surface charge; (2) reduction of the polarization saturation value, as well as the critical thickness at which it is obtained, by the presence of free carriers. The resilience of polar layered structures to atomic surface reconstruction, which is essentially unavoidable in polar 3D crystals, potentially allows for the design of new devices with mobile surface charges. Our findings, which are of general nature, should be accounted for when designing switching and/or conductive devices based on ferroelectric layered materials.
△ Less
Submitted 29 July, 2024;
originally announced July 2024.
-
State-specific density functionals for excited states from ensembles
Authors:
Tim Gould,
Stephen G Dale,
Leeor Kronik,
Stefano Pittalis
Abstract:
We present a first principles strategy for developing state-specific density functional approximations for excited states. We first clarify why approaches based on conventional ground state approximations miss density-driven correlations, by considering excited state physics through the lens of ensemble density functional theory. To solve this issue we gain insights on density driven correlations…
▽ More
We present a first principles strategy for developing state-specific density functional approximations for excited states. We first clarify why approaches based on conventional ground state approximations miss density-driven correlations, by considering excited state physics through the lens of ensemble density functional theory. To solve this issue we gain insights on density driven correlations by exploiting the recently understood low-density limit of electrons in excited states. The theory developments are then combined to produce a proof-of-concept excited state approximation that resolves urgent paradigmatic failures (double excitations, charge transfer excitations, piecewise linearity) of existing state-of-art density-functional approaches, directly from differences of self-consistent field calculations; i.e., $Δ$SCF. In light of its observed impressive performance, we conclude that the approach represents a major step toward unified and accurate modelling of neutral and charged excitations.
△ Less
Submitted 14 October, 2024; v1 submitted 26 June, 2024;
originally announced June 2024.
-
Electronic and Optical Excitations in van der Waals Materials from a Non-Empirical Wannier-Localized Optimally-Tuned Screened Range-Separated Hybrid Functional
Authors:
María Camarasa-Gómez,
Stephen E. Gant,
Guy Ohad,
Jeffrey B. Neaton,
Ashwin Ramasubramanian,
Leeor Kronik
Abstract:
Accurate prediction of electronic and optical excitations in van der Waals (vdW) materials is a long-standing challenge for density functional theory. The recently proposed Wannier-localized optimally-tuned screened range-separated hybrid (WOT-SRSH) functional has proven successful in non-empirical determination of electronic band gaps and optical absorption spectra for various covalent and ionic…
▽ More
Accurate prediction of electronic and optical excitations in van der Waals (vdW) materials is a long-standing challenge for density functional theory. The recently proposed Wannier-localized optimally-tuned screened range-separated hybrid (WOT-SRSH) functional has proven successful in non-empirical determination of electronic band gaps and optical absorption spectra for various covalent and ionic crystals. However, for vdW materials the tuning of the material- and structure-dependent functional parameters has, until now, only been attained semi-empirically. Here, we present a non-empirical WOT-SRSH approach applicable to vdW materials, with the optimal functional parameters transferable between monolayer and bulk. We apply this methodology to prototypical vdW materials: black phosphorus, molybdenum disulfide, and hexagonal boron nitride (in the latter case including zero-point renormalization). We show that the WOT-SRSH approach consistently achieves accuracy levels comparable to experiments and ab initio many-body perturbation theory (MBPT) calculations for band structures and optical absorption spectra, both on its own and as an optimal starting point for MBPT calculations.
△ Less
Submitted 1 May, 2024;
originally announced May 2024.
-
Non-empirical prediction of the length-dependent ionization potential in molecular chains
Authors:
Guy Ohad,
Michal Hartstein,
Tim Gould,
Jeffrey B. Neaton,
Leeor Kronik
Abstract:
The ionization potential of molecular chains is well-known to be a tunable nano-scale property that exhibits clear quantum confinement effects. State-of-the-art methods can accurately predict the ionization potential in the small molecule limit and in the solid-state limit, but for intermediate, nano-sized systems prediction of the evolution of the electronic structure between the two limits is mo…
▽ More
The ionization potential of molecular chains is well-known to be a tunable nano-scale property that exhibits clear quantum confinement effects. State-of-the-art methods can accurately predict the ionization potential in the small molecule limit and in the solid-state limit, but for intermediate, nano-sized systems prediction of the evolution of the electronic structure between the two limits is more difficult. Recently, optimal tuning of range-separated hybrid functionals has emerged as a highly accurate method for predicting ionization potentials. This was first achieved for molecules using the ionization potential theorem (IPT) and more recently extended to solid-state systems, based on an \textit{ansatz} that generalizes the IPT to the removal of charge from a localized Wannier function. Here, we study one-dimensional molecular chains of increasing size, from the monomer limit to the infinite polymer limit using this approach. By comparing our results with other localization-based methods and where available with experiment, we demonstrate that Wannier-localization-based optimal tuning is highly accurate in predicting ionization potentials for any chain length, including the nano-scale regime.
△ Less
Submitted 27 March, 2024;
originally announced March 2024.
-
Optical absorption spectra of metal oxides from time-dependent density functional theory and many-body perturbation theory based on optimally-tuned hybrid functionals
Authors:
Guy Ohad,
Stephen E. Gant,
Dahvyd Wing,
Jonah B. Haber,
María Camarasa-Gómez,
Francisca Sagredo,
Marina R. Filip,
Jeffrey B. Neaton,
Leeor Kronik
Abstract:
Using both time-dependent density functional theory (TDDFT) and the ``single-shot" $GW$ plus Bethe-Salpeter equation ($GW$-BSE) approach, we compute optical band gaps and optical absorption spectra from first principles for eight common binary and ternary closed-shell metal oxides (MgO, Al$_2$O$_3$, CaO, TiO$_2$, Cu$_2$O, ZnO, BaSnO$_3$, and BiVO$_4$), based on the non-empirical Wannier-localized…
▽ More
Using both time-dependent density functional theory (TDDFT) and the ``single-shot" $GW$ plus Bethe-Salpeter equation ($GW$-BSE) approach, we compute optical band gaps and optical absorption spectra from first principles for eight common binary and ternary closed-shell metal oxides (MgO, Al$_2$O$_3$, CaO, TiO$_2$, Cu$_2$O, ZnO, BaSnO$_3$, and BiVO$_4$), based on the non-empirical Wannier-localized optimally-tuned screened range-separated hybrid functional. Overall, we find excellent agreement between our TDDFT and $GW$-BSE results and experiment, with a mean absolute error less than 0.4 eV, including for Cu$_2$O and ZnO, traditionally considered to be challenging for both methods.
△ Less
Submitted 5 September, 2023;
originally announced September 2023.
-
Anisotropic Interlayer Force Field for Group-VI Transition Metal Dichalcogenides
Authors:
Wenwu Jiang,
Reut Sofer,
Xiang Gao,
Alexandre Tkatchenko,
Leeor Kronik,
Wengen Ouyang,
Michael Urbakh,
Oded Hod
Abstract:
An anisotropic interlayer force field that describes the interlayer interactions in homogeneous and heterogeneous interfaces of group-VI transition metal dichalcogenides (MX2 where M = Mo, W and X = S, Se) is presented. The force field is benchmarked against density functional theory calculations for bilayer systems within the Heyd-Scuseria-Ernzerhof hybrid density functional approximation, augmen…
▽ More
An anisotropic interlayer force field that describes the interlayer interactions in homogeneous and heterogeneous interfaces of group-VI transition metal dichalcogenides (MX2 where M = Mo, W and X = S, Se) is presented. The force field is benchmarked against density functional theory calculations for bilayer systems within the Heyd-Scuseria-Ernzerhof hybrid density functional approximation, augmented by a nonlocal many-body dispersion treatment of long-range correlation. The parametrization yields good agreement with reference calculations of binding energy curves and sliding potential energy surfaces. It is found to be transferable to TMD junctions outside the training set that contain the same atom types. Calculated bulk moduli agree with most previous dispersion corrected DFT predictions, which underestimate available experimental values. Calculated phonon spectra of the various junctions under consideration demonstrate the importance of appropriately treating the anisotropic nature of layered interfaces. Considering our previous parameterization for MoS2, the interlayer potential enables accurate and efficient large-scale simulations of the dynamical, tribological, and thermal transport properties of a large set of homogeneous and heterogeneous TMD interfaces.
△ Less
Submitted 23 July, 2023; v1 submitted 20 July, 2023;
originally announced July 2023.
-
Effects of Quantum and Dielectric Confinement on the Emission of Cs-Pb-Br Composites
Authors:
Sebastián Caicedo-Dávila,
Pietro Caprioglio,
Frederike Lehmann,
Sergiu Levcenco,
Martin Stolterfoht,
Dieter Neher,
Leeor Kronik,
Daniel Abou-Ras
Abstract:
The halide perovskite CsPbBr$_3$ belongs to the Cs-Pb-Br material system, which features two additional thermodynamically stable ternary phases, Cs$_4$PbBr$_6$ and CsPb$_2$Br$_5$. The coexistence of these phases and their reportedly similar photoluminescence have resulted in a debate on the nature of the emission in these systems. Here, we combine optical and microscopic characterization with an e…
▽ More
The halide perovskite CsPbBr$_3$ belongs to the Cs-Pb-Br material system, which features two additional thermodynamically stable ternary phases, Cs$_4$PbBr$_6$ and CsPb$_2$Br$_5$. The coexistence of these phases and their reportedly similar photoluminescence have resulted in a debate on the nature of the emission in these systems. Here, we combine optical and microscopic characterization with an effective mass, correlated electron-hole model of excitons in confined systems, to investigate the emission properties of the ternary phases in the Cs-Pb-Br system. We find that all Cs-Pb-Br phases exhibit green emission and the non-perovskite phases exhibit photoluminescence quantum yields orders of magnitude larger than CsPbBr$_3$. In particular, we measure blue- and red-shifted emission for the Cs- and Pb-rich phases, respectively, stemming from embedded CsPbBr$_3$ nanocrystals. Our model reveals that the difference in emission shift is caused by the combined effects of nanocrystal size and different band mismatch. Furthermore, we demonstrate the importance of including the dielectric mismatch in the calculation of the emission energy for Cs-Pb-Br composites. Our results explain the reportedly limited blue shift in CsPbBr$_3$@Cs$_4$PbBr$_6$ composites and rationalize some of its differences with CsPb$_2$Br$_5$.
△ Less
Submitted 21 June, 2023; v1 submitted 14 June, 2023;
originally announced June 2023.
-
Transferable screened range-separated hybrid functionals for electronic and optical properties of van der Waals materials
Authors:
María Camarasa-Gómez,
Ashwin Ramasubramaniam,
Jeffrey B. Neaton,
Leeor Kronik
Abstract:
The accurate description of electronic properties and optical absorption spectra is a long-standing challenge for density functional theory. Recently, the introduction of screened range-separated hybrid (SRSH) functionals for solid-state materials has allowed for the calculation of fundamental band gaps and optical absorption spectra that are in very good agreement with many-body perturbation theo…
▽ More
The accurate description of electronic properties and optical absorption spectra is a long-standing challenge for density functional theory. Recently, the introduction of screened range-separated hybrid (SRSH) functionals for solid-state materials has allowed for the calculation of fundamental band gaps and optical absorption spectra that are in very good agreement with many-body perturbation theory. However, since solid-state SRSH functionals are typically tuned to reproduce the properties of bulk phases, their transferability to low-dimensional structures, which experience substantially different screening than in the bulk, remains an open question. In this work, we explore the transferability of SRSH functionals to several prototypical van der Waals materials, including transition-metal sulfides and selenides, indium selenide, black phosphorus, and hexagonal boron nitride. Considering the bulk and a monolayer of these materials as limiting cases, we show that the parameters of the SRSH functional can be determined systematically, using only the band-edge quasiparticle energies of these extremal structural phases as fitting targets. The resulting SRSH functionals can describe both electronic bandstructures and optical absorption spectra with accuracy comparable to more demanding ab initio many-body perturbation theory (GW and Bethe-Salpeter equation) approaches. Selected examples also demonstrate that the SRSH parameters, obtained from the bulk and monolayer reference structures, display good accuracy for bandstructures and optical spectra of bilayers, indicating a degree of transferability that is independent of the fitting procedure.
△ Less
Submitted 22 May, 2023;
originally announced May 2023.
-
Spontaneous Electric Polarization in Graphene Polytypes
Authors:
Simon Salleh Atri,
Wei Cao,
Bar Alon,
Nirmal Roy,
Maayan Vizner Stern,
Vladimir Falko,
Moshe Goldstein,
Leeor Kronik,
Michael Urbakh,
Oded Hod,
Moshe Ben Shalom
Abstract:
A crystalline solid is a periodic sequence of identical cells, each containing one or more atoms. If the constituting unit cell is not centrosymmetric, charge may distribute unevenly between the atoms, resulting in internal electric polarization. This effect serves as the basis for numerous ferroelectric, piezoelectric, and pyroelectric phenomena. In nearly all polar materials, including multilaye…
▽ More
A crystalline solid is a periodic sequence of identical cells, each containing one or more atoms. If the constituting unit cell is not centrosymmetric, charge may distribute unevenly between the atoms, resulting in internal electric polarization. This effect serves as the basis for numerous ferroelectric, piezoelectric, and pyroelectric phenomena. In nearly all polar materials, including multilayered van der Waals stacks that were recently found to exhibit interfacial polarization, inversion symmetry is broken by having two or more atomic species within the unit cell. Here, we show that even elemental crystals, consisting of one type of atom, and composed of non-polar centrosymmetric layers, exhibit electric polarization if arranged in an appropriate three-dimensional architecture. This concept is demonstrated here for inversion and mirror asymmetric mixed-stacking tetra-layer polytypes of non-polar graphene sheets. Furthermore, we find that the room temperature out-of-plane electric polarization increases with external electrostatic doping, rather than decreases owing to screening. Using first-principles calculations, as well as tight-binding modeling, we unveil the origin of polytype-induced polarization and its dependence on doping. Extension of this idea to graphene multilayers suggests that solely by lateral shifts of constituent monolayers one can obtain multiple meta-stable interlayer stacking sequences that may allow for even larger electrical polarization.
△ Less
Submitted 18 May, 2023;
originally announced May 2023.
-
Two-dimensional carrier gas at a polar interface without surface band gap states: A first-principles perspective
Authors:
Federico Brivio,
Andrew Rappe,
Leeor Kronik,
Dan Ritter
Abstract:
We present first-principles calculations of the interface between GaN and strained AlN, using a slab model in which polarization is compensated via surface fractional-charge pseudo-hydrogen atoms. We show that an interface two-dimensional carrier electron or hole gas emerges naturally in response to different compensating surface charges, but that this need not involve in-gap surface states.
We present first-principles calculations of the interface between GaN and strained AlN, using a slab model in which polarization is compensated via surface fractional-charge pseudo-hydrogen atoms. We show that an interface two-dimensional carrier electron or hole gas emerges naturally in response to different compensating surface charges, but that this need not involve in-gap surface states.
△ Less
Submitted 7 March, 2023;
originally announced March 2023.
-
Range-Separated Hybrid Functional Pseudopotentials
Authors:
Yang Yang,
Georgia Prokopiou,
Tian Qiu,
Aaron M. Schankler,
Andrew M. Rappe,
Leeor Kronik,
Robert A. DiStasio Jr
Abstract:
Consistency between the exchange-correlation (xc) functional used during pseudopotential construction and planewave-based electronic structure calculations is important for an accurate and reliable description of the structure and properties of condensed-phase systems. In this work, we present a general scheme for constructing pseudopotentials with range-separated hybrid (RSH) xc functionals based…
▽ More
Consistency between the exchange-correlation (xc) functional used during pseudopotential construction and planewave-based electronic structure calculations is important for an accurate and reliable description of the structure and properties of condensed-phase systems. In this work, we present a general scheme for constructing pseudopotentials with range-separated hybrid (RSH) xc functionals based on the solution of the all-electron radial integro-differential equation for a spherical atomic configuration. As proof-of-principle, we demonstrate pseudopotential construction with the PBE, PBE0, HSE, and sRSH (based on LC-$ω$PBE0) xc functionals for a select set of atoms, and then investigate the importance of pseudopotential consistency when computing band gaps, equilibrium lattice parameters, bulk moduli, and atomization energies of several solid-state systems. In doing so, we find that pseudopotential consistency errors (PSCE) tend to be systematic and can be as large as $1.4\%$ when computing these properties.
△ Less
Submitted 22 February, 2023;
originally announced February 2023.
-
Uncovering Multiple Intrinsic Chiral Phases in (PEA)2PbI4 halide Perovskites
Authors:
Shahar Zuri,
Arthur Shapiro,
Leeor Kronik,
Efrat Lifshitz
Abstract:
Two-dimensional (2D) halide perovskites offer a unique platform to investigate the ground-state of materials possessing significant anharmonicity. In contrast to three-dimensional perovskites, their 2D counterparts offer substantially fewer degrees of freedom, resulting in multiple well-defined crystal structures. In this work, we thoroughly investigate the anharmonic ground-state of the benchmark…
▽ More
Two-dimensional (2D) halide perovskites offer a unique platform to investigate the ground-state of materials possessing significant anharmonicity. In contrast to three-dimensional perovskites, their 2D counterparts offer substantially fewer degrees of freedom, resulting in multiple well-defined crystal structures. In this work, we thoroughly investigate the anharmonic ground-state of the benchmark (PEA)2PbI4 compound, using complementary information from low-temperature x-ray diffraction (XRD) and photoluminescence spectroscopy, supported by density functional theory (DFT) calculations. We extrapolate four crystallographic configurations from the low-temperature XRD. These configurations imply that the ground-state has an intrinsic disorder stemming from two coexisting chiral sub-lattices, each with a bi-oriented organic spacer molecule. We further show evidence that these chiral structures form unevenly populated ground-states, portraying uneven anharmonicity, where the state population may be tuned by surface effects. Our results uncover a disordered ground-state that may induce intrinsic grain boundaries, which cannot be ignored in practical applications.
△ Less
Submitted 23 November, 2022;
originally announced November 2022.
-
Nuclear cusps and singularities in the non-additive kinetic potential bi-functional from analytical inversion
Authors:
Mojdeh Banafsheh,
Tomasz A. Wesolowski,
Tim Gould,
Leeor Kronik,
David A. Strubbe
Abstract:
The non-additive kinetic potential $v^{\text{NAD}}$ is a key quantity in density-functional theory (DFT) embedding methods, such as frozen density embedding theory and partition DFT. $v^{\text{NAD}}$ is a bi-functional of electron densities $ρ_{\rm B}$ and $ρ_{\rm tot} = ρ_{\rm A} + ρ_{\rm B}$. It can be evaluated using approximate kinetic-energy functionals, but accurate approximations are challe…
▽ More
The non-additive kinetic potential $v^{\text{NAD}}$ is a key quantity in density-functional theory (DFT) embedding methods, such as frozen density embedding theory and partition DFT. $v^{\text{NAD}}$ is a bi-functional of electron densities $ρ_{\rm B}$ and $ρ_{\rm tot} = ρ_{\rm A} + ρ_{\rm B}$. It can be evaluated using approximate kinetic-energy functionals, but accurate approximations are challenging. The behavior of $v^{\text{NAD}}$ in the vicinity of the nuclei has long been questioned, and singularities were seen in some approximate calculations. In this article, the existence of singularities in $v^{\text{NAD}}$ is analyzed analytically for various choices of $ρ_{\rm B}$ and $ρ_{\rm tot}$, using the nuclear cusp conditions for the density and Kohn-Sham potential. It is shown that no singularities arise from smoothly partitioned ground-state Kohn-Sham densities. We confirm this result by numerical calculations on diatomic test systems HeHe, HeLi$^+$, and H$_2$, using analytical inversion to obtain a numerically exact $v^{\rm NAD}$ for the local density approximation. We examine features of $v^{\rm NAD}$ which can be used for development and testing of approximations to $v^{\rm NAD}[ρ_{\rm B},ρ_{\rm tot}]$ and kinetic-energy functionals.
△ Less
Submitted 10 October, 2022; v1 submitted 8 July, 2022;
originally announced July 2022.
-
Cumulative Polarization Coexisting with Conductivity at Interfacial Ferroelectrics
Authors:
Swarup Deb,
Wei Cao,
Noam Raab,
Kenji Watanabe,
Takashi Taniguchi,
Moshe Goldstein,
Leeor Kronik,
Michael Urbakh,
Oded Hod,
Moshe Ben Shalom
Abstract:
Ferroelectricity in atomically thin bilayer structures has been recently predicted1 and measured[2-4] in two-dimensional (2D) materials with hexagonal non-centrosymmetric unit-cells. Interestingly, the crystal symmetry translates lateral shifts between parallel 2D layers to a change of sign in their out-of-plane electric polarization, a mechanism referred to as "Slide-Tronics"[4]. These observatio…
▽ More
Ferroelectricity in atomically thin bilayer structures has been recently predicted1 and measured[2-4] in two-dimensional (2D) materials with hexagonal non-centrosymmetric unit-cells. Interestingly, the crystal symmetry translates lateral shifts between parallel 2D layers to a change of sign in their out-of-plane electric polarization, a mechanism referred to as "Slide-Tronics"[4]. These observations, however, have been restricted to switching between only two polarization states under low charge carrier densities[5-12], strongly limiting the practical application of the revealed phenomena[13]. To overcome these issues, one needs to explore the nature of the polarization that arises in multi-layered van der Waals (vdW) stacks, how it is governed by intra- and inter-layer charge redistribution, and to which extent it survives the introduction of mobile charge carriers, all of which are presently unknown14. To explore these questions, we conduct surface potential measurements of parallel WSe2 and MoS2 multi-layers with aligned and anti-aligned configurations of the polar interfaces. We find evenly spaced, nearly decoupled potential steps, indicating highly confined interfacial electric fields, which provide means to design multi-state "ladder ferroelectrics". Furthermore, we find that the internal polarization remains significant upon electrostatic doping of a mobile charge carrier density as high as 1013 cm-2, with substantial in-plane conductivity. Using first-principles calculations based on density functional theory (DFT), we trace the extra charge redistribution in real and momentum space and identify an eventual doping-induced depolarization mechanism.
△ Less
Submitted 24 June, 2022;
originally announced June 2022.
-
An Optimally-Tuned Starting Point for Single-Shot $GW$ Calculations of Solids
Authors:
Stephen E. Gant,
Jonah B. Haber,
Marina R. Filip,
Francisca Sagredo,
Dahvyd Wing,
Guy Ohad,
Leeor Kronik,
Jeffrey B. Neaton
Abstract:
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…
▽ More
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.
△ Less
Submitted 24 May, 2022; v1 submitted 1 February, 2022;
originally announced February 2022.
-
Theory of Chirality Induced Spin Selectivity: Progress and Challenges
Authors:
Ferdinand Evers,
Amnon Aharony,
Nir Bar-Gill,
Ora Entin-Wohlman,
Per Hedegård,
Oded Hod,
Pavel Jelinek,
Grzegorz Kamieniarz,
Mikhail Lemeshko,
Karen Michaeli,
Vladimiro Mujica,
Ron Naaman,
Yossi Paltiel,
Sivan Refaely-Abramson,
Oren Tal,
Jos Thijssen,
Michael Thoss,
Jan M. van Ruitenbeek,
Latha Venkataraman,
David H. Waldeck,
Binghai Yan,
Leeor Kronik
Abstract:
We provide a critical overview of the theory of the chirality-induced spin selectivity (CISS) effect, i.e., phenomena in which the chirality of molecular species imparts significant spin selectivity to various electron processes. Based on discussions in a recently held workshop, and further work published since, we review the status of CISS effects - in electron transmission, electron transport, a…
▽ More
We provide a critical overview of the theory of the chirality-induced spin selectivity (CISS) effect, i.e., phenomena in which the chirality of molecular species imparts significant spin selectivity to various electron processes. Based on discussions in a recently held workshop, and further work published since, we review the status of CISS effects - in electron transmission, electron transport, and chemical reactions. For each, we provide a detailed discussion of the state-of-the-art in theoretical understanding and identify remaining challenges and research opportunities.
△ Less
Submitted 23 August, 2021;
originally announced August 2021.
-
Band gaps of crystalline solids from Wannier-localization based optimal tuning of a screened range-separated hybrid functional
Authors:
Dahvyd Wing,
Guy Ohad,
Jonah B. Haber,
Marina R. Filip,
Stephen E. Gant,
Jeffrey B. Neaton,
Leeor Kronik
Abstract:
Accurate prediction of fundamental band gaps of crystalline solid state systems entirely within density functional theory is a long standing challenge. Here, we present a simple and inexpensive method that achieves this by means of non-empirical optimal tuning of the parameters of a screened range-separated hybrid functional. The tuning involves the enforcement of an ansatz that generalizes the io…
▽ More
Accurate prediction of fundamental band gaps of crystalline solid state systems entirely within density functional theory is a long standing challenge. Here, we present a simple and inexpensive method that achieves this by means of non-empirical optimal tuning of the parameters of a screened range-separated hybrid functional. The tuning involves the enforcement of an ansatz that generalizes the ionization potential theorem to the removal of an electron in an occupied state described by a localized Wannier function in a modestly sized supercell calculation. The method is benchmarked against experiment for a set of systems ranging from narrow band gap semiconductors to large band gap insulators, spanning a range of fundamental band gaps from 0.2 to 14.2 eV and is found to yield quantitative accuracy across the board, with a mean absolute error of $\sim$0.1 eV and a maximal error of $\sim$0.2 eV.
△ Less
Submitted 11 October, 2021; v1 submitted 6 December, 2020;
originally announced December 2020.
-
The pursuit of stability in halide perovskites: the monovalent cation and the key for surface and bulk self-repair
Authors:
D R Ceratti,
A V Cohen,
R Tenne,
Y Rakita,
L Snarski,
L Cremonesi,
I Goldian,
I Kaplan-Ashiri,
T Bendikov,
V Kalchenko,
M Elbaum,
M A C Potenza,
L Kronik,
G Hodes,
D Cahen
Abstract:
We find significant differences between degradation and healing at the surface or in the bulk for each of the different APbBr3 single crystals (A=CH3NH3+, methylammonium (MA); HC(NH2)2+, formamidinium (FA); and cesium, Cs+). Using 1- and 2-photon microscopy and photobleaching we conclude that kinetics dominate the surface, and thermodynamics the bulk stability. Fluorescence-lifetime imaging micros…
▽ More
We find significant differences between degradation and healing at the surface or in the bulk for each of the different APbBr3 single crystals (A=CH3NH3+, methylammonium (MA); HC(NH2)2+, formamidinium (FA); and cesium, Cs+). Using 1- and 2-photon microscopy and photobleaching we conclude that kinetics dominate the surface, and thermodynamics the bulk stability. Fluorescence-lifetime imaging microscopy, as well as results from several other methods, relate the (damaged) state of the halide perovskite (HaP) after photobleaching to its modified optical and electronic properties. The A cation type strongly influences both the kinetics and the thermodynamics of recovery and degradation: FA heals best the bulk material with faster self-healing; Cs+ protects the surface best, being the least volatile of the A cations and possibly through O-passivation; MA passivates defects via methylamine from photo-dissociation, which binds to Pb2+. DFT simulations not only provide insight into the latter conclusion, but also show the importance and stability of the Br3- defect. These results rationalize the use of mixed A-cation materials for optimizing both solar cell stability and overall performance of HaP-based devices, and provide a basis for designing new HaP variants.
△ Less
Submitted 30 September, 2020;
originally announced September 2020.
-
Mechanical and Tribological Properties of Layered Materials Under High Pressure: Assessing the Importance of Many-Body Dispersion Effects
Authors:
Wengen Ouyang,
Ido Azuri,
Davide Mandelli,
Alexandre Tkatchenko,
Leeor Kronik,
Michael Urbakh,
Oded Hod
Abstract:
The importance of many-body dispersion effects in layered materials subjected to high external loads is evaluated. State-of-the-art many-body dispersion density functional theory calculations performed for graphite, hexagonal boron nitride, and their hetero-structures were used to fit the parameters of a classical registry-dependent interlayer potential. Using the latter, we performed extensive eq…
▽ More
The importance of many-body dispersion effects in layered materials subjected to high external loads is evaluated. State-of-the-art many-body dispersion density functional theory calculations performed for graphite, hexagonal boron nitride, and their hetero-structures were used to fit the parameters of a classical registry-dependent interlayer potential. Using the latter, we performed extensive equilibrium molecular dynamics simulations and studied the mechanical response of homogeneous and heterogeneous bulk models under hydrostatic pressures up to 30 GPa. Comparison with experimental data demonstrates that the reliability of the many-body dispersion model extends deep into the sub-equilibrium regime. Friction simulations demonstrate the importance of many-body dispersion effects for the accurate description of the tribological properties of layered materials interfaces under high pressure.
△ Less
Submitted 28 December, 2019; v1 submitted 12 December, 2019;
originally announced December 2019.
-
Breakdown of the static picture of defect energetics in halide perovskites: the case of the Br vacancy in CsPbBr3
Authors:
Ayala V. Cohen,
David A. Egger,
Andrew M. Rappe,
Leeor Kronik
Abstract:
We consider the Br vacancy in CsPbBr3 as a prototype for the impact of structural dynamics on defect energetics in halide perovskites (HaPs). Using first-principles molecular dynamics based on density functional theory, we find that the static picture of defect energetics breaks down; the energy of the Br vacancy level is found to be intrinsically dynamic, oscillating by as much as 1 eV on the ps…
▽ More
We consider the Br vacancy in CsPbBr3 as a prototype for the impact of structural dynamics on defect energetics in halide perovskites (HaPs). Using first-principles molecular dynamics based on density functional theory, we find that the static picture of defect energetics breaks down; the energy of the Br vacancy level is found to be intrinsically dynamic, oscillating by as much as 1 eV on the ps time scale at room temperature. These significant energy fluctuations are correlated with the distance between the neighboring Pb atoms across the vacancy and with the electrostatic potential at these Pb atomic sites. We expect this unusually strong coupling of structural dynamics and defect energetics to bear important implications for both experimental and theoretical analysis of defect characteristics in HaPs. It may also hold significant ramifications for carrier transport and defect tolerance in this class of photovoltaic materials.
△ Less
Submitted 27 June, 2019; v1 submitted 10 October, 2018;
originally announced October 2018.
-
Prediction of electronic couplings for molecular charge transfer using optimally-tuned range-separated hybrid functionals
Authors:
Debashree Manna,
Jochen Blumberger,
Jan M. L. Martin,
Leeor Kronik
Abstract:
Electronic coupling matrix elements are important to the theoretical description of electron transfer processes. However, they are notoriously difficult to obtain accurately from time- dependent density functional theory (TDDFT). Here, we use the HAB11 benchmark dataset of coupling matrix elements to assess whether TDDFT using optimally-tuned range-separated hybrid functionals, already known to be…
▽ More
Electronic coupling matrix elements are important to the theoretical description of electron transfer processes. However, they are notoriously difficult to obtain accurately from time- dependent density functional theory (TDDFT). Here, we use the HAB11 benchmark dataset of coupling matrix elements to assess whether TDDFT using optimally-tuned range-separated hybrid functionals, already known to be successful for the description of charge transfer excitation energies, also allows for an improved accuracy in the prediction of coupling matrix elements. We find that this approach outperforms all previous TDDFT calculations, based on semi-local, hybrid, or non-tuned range-separated hybrid functionals, with a remaining average deviation as low as ~12%. We discuss potential sources for the remaining error.
△ Less
Submitted 5 June, 2018;
originally announced June 2018.
-
Charge transfer excitations from exact and approximate ensemble Kohn-Sham theory
Authors:
Tim Gould,
Leeor Kronik,
Stefano Pittalis
Abstract:
By studying the lowest excitations of an exactly solvable one-dimensional molecular model, we show that components of Kohn-Sham ensembles can be used to describe charge transfers. Furthermore, we compute the approximate excitation energies obtained by using thee exact ensemble densities in the recently formulated ensemble Hartree-exchange theory [Gould and Pittalis, Phys. Rev. Lett. 119, 243001 (2…
▽ More
By studying the lowest excitations of an exactly solvable one-dimensional molecular model, we show that components of Kohn-Sham ensembles can be used to describe charge transfers. Furthermore, we compute the approximate excitation energies obtained by using thee exact ensemble densities in the recently formulated ensemble Hartree-exchange theory [Gould and Pittalis, Phys. Rev. Lett. 119, 243001 (2017)]. Remarkably, our results show that triplet excitations are accurately reproduced across a dissociation curve in all cases tested, even in systems where ground state energies are poor due to strong static correlations. Singlet excitations exhibit larger deviations from exact results but are still reproduced semi-quantitatively.
△ Less
Submitted 19 January, 2018;
originally announced January 2018.
-
The Effect of Ionic Composition on Acoustic Phonon Speeds in Hybrid Perovskites from Brillouin Spectroscopy and Density Functional Theory
Authors:
Irina V. Kabakova,
Ido Azuri,
Zhuoying Chen,
Pabitra K. Nayak,
Henry J. Snaith,
Leeor Kronik,
Carl Paterson,
Artem A. Bakulin,
David A. Egger
Abstract:
Hybrid organic-inorganic perovskites (HOIPs) have recently emerged as highly promising solution-processable materials for photovoltaic (PV) and other optoelectronic devices. HOIPs represent a broad family of materials with properties highly tuneable by the ions that make up the perovskite structure as well as their multiple combinations. Interestingly, recent high-efficiency PV devices using HOIPs…
▽ More
Hybrid organic-inorganic perovskites (HOIPs) have recently emerged as highly promising solution-processable materials for photovoltaic (PV) and other optoelectronic devices. HOIPs represent a broad family of materials with properties highly tuneable by the ions that make up the perovskite structure as well as their multiple combinations. Interestingly, recent high-efficiency PV devices using HOIPs with substantially improved long-term stability have used combinations of different ionic compositions. The structural dynamics of these systems are unique for semiconducting materials and are currently argued to be central to HOIPs stability and charge-transport properties. Here, we studied the impact of ionic composition on phonon speeds of HOIPs from Brillouin spectroscopy experiments and density functional theory calculations for FAPbBr$_3$, MAPbBr$_3$, MAPbCl$_3$, and the mixed halide MAPbBr$_{1.25}$Cl$_{1.75}$. Our results show that the acoustic phonon speeds can be strongly modified by ionic composition, which we explain by analysing the lead-halide sublattice in detail. The vibrational properties of HOIPs are therefore tuneable by using targeted ionic compositions in the perovskite structure. This tuning can be rationalized with non-trivial effects, for example, considering the influence of the shape and dipole moment of organic cations. This has an important implication to further improvements in the stability and charge-transport properties of these systems.
△ Less
Submitted 10 January, 2018;
originally announced January 2018.
-
Energy Level Alignment at Molecule-Metal Interfaces from an Optimally-Tuned Range-Separated Hybrid Functional
Authors:
Zhenfei Liu,
David A. Egger,
Sivan Refaely-Abramson,
Leeor Kronik,
Jeffrey B. Neaton
Abstract:
The alignment of the frontier orbital energies of an adsorbed molecule with the substrate Fermi level at metal-organic interfaces is a fundamental observable of significant practical importance in nanoscience and beyond. Typical density functional theory calculations, especially those using local and semi-local functionals, often underestimate level alignment leading to inaccurate electronic struc…
▽ More
The alignment of the frontier orbital energies of an adsorbed molecule with the substrate Fermi level at metal-organic interfaces is a fundamental observable of significant practical importance in nanoscience and beyond. Typical density functional theory calculations, especially those using local and semi-local functionals, often underestimate level alignment leading to inaccurate electronic structure and charge transport properties. In this work, we develop a new fully self-consistent predictive scheme to accurately compute level alignment at certain classes of complex heterogeneous molecule-metal interfaces based on optimally-tuned range-separated hybrid functionals. Starting from a highly accurate description of the gas-phase electronic structure, our method by construction captures important nonlocal surface polarization effects via tuning of the long-range screened exchange in a range-separated hybrid in a non-empirical and system-specific manner. We implement this functional in a plane-wave code and apply it to several physisorbed and chemisorbed molecule-metal interface systems. Our results are in quantitative agreement with experiments, both the level alignment and work function changes. Our approach constitutes a new practical scheme for accurate and efficient calculations of the electronic structure of molecule-metal interfaces.
△ Less
Submitted 6 March, 2017;
originally announced March 2017.
-
Optical Phonons in Methylammonium Lead Halide Perovskites and Implications for Charge Transport
Authors:
Michael Sendner,
Pabitra K. Nayak,
David A. Egger,
Sebastian Beck,
Christian Müller,
Bernd Epding,
Wolfgang Kowalsky,
Leeor Kronik,
Henry J. Snaith,
Annemarie Pucci,
Robert Lovrinčić
Abstract:
Lead-halide perovskites are promising materials for opto-electronic applications. Recent reports indicated that their mechanical and electronic properties are strongly affected by the lattice vibrations. Herein we report far-infrared spectroscopy measurements of CH$_{3}$NH$_{3}$Pb(I/Br/Cl)$_{3}$ thin films and single crystals at room temperature and a detailed quantitative analysis of the spectra.…
▽ More
Lead-halide perovskites are promising materials for opto-electronic applications. Recent reports indicated that their mechanical and electronic properties are strongly affected by the lattice vibrations. Herein we report far-infrared spectroscopy measurements of CH$_{3}$NH$_{3}$Pb(I/Br/Cl)$_{3}$ thin films and single crystals at room temperature and a detailed quantitative analysis of the spectra. We find strong broadening and anharmonicity of the lattice vibrations for all three halide perovskites, which indicates dynamic disorder of the lead-halide cage at room temperature. We determine the frequencies of the transversal and longitudinal optical phonons, and use them to calculate, via appropriate models, the static dielectric constants, polaron masses, electron-phonon coupling constants, and upper limits for the phonon-scattering limited charge carrier mobilities. Within the limitations of the model used, we can place an upper limit of 200$\,$cm$^{2}$V$^{-1}$s$^{-1}$ for the room temperature charge carrier mobility in MAPbI$_{3}$ single crystals. Our findings are important for the basic understanding of charge transport processes and mechanical properties in metal halide perovskites.
△ Less
Submitted 11 October, 2016; v1 submitted 28 July, 2016;
originally announced July 2016.
-
Interplay of bias-driven charging and the vibrational Stark effect in molecular junctions
Authors:
Yajing Li,
Pavlo Zolotavin,
Peter Doak,
Leeor Kronik,
Jeffrey B. Neaton,
Douglas Natelson
Abstract:
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 densi…
▽ More
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.
△ Less
Submitted 10 May, 2016;
originally announced May 2016.
-
Local polar fluctuations in lead halide perovskite crystals
Authors:
Omer Yaffe,
Yinsheng Guo,
Liang Z. Tan,
David A. Egger,
Trevor Hull,
Constantinos C. Stoumpos,
Fan Zheng,
Tony F. Heinz,
Leeor Kronik,
Mercouri G. Kanatzidis,
Jonathan S Owen,
Andrew M. Rappe,
Marcos A. Pimenta,
Louis E. Brus
Abstract:
Hybrid lead-halide perovskites have emerged as an excellent class of photovoltaic materials. Recent reports suggest that the organic molecular cation is responsible for local polar fluctuations that inhibit carrier recombination. We combine low frequency Raman scattering with first-principles molecular dynamics (MD) to study the fundamental nature of these local polar fluctuations. Our observation…
▽ More
Hybrid lead-halide perovskites have emerged as an excellent class of photovoltaic materials. Recent reports suggest that the organic molecular cation is responsible for local polar fluctuations that inhibit carrier recombination. We combine low frequency Raman scattering with first-principles molecular dynamics (MD) to study the fundamental nature of these local polar fluctuations. Our observations of a strong central peak in both hybrid (CH$_3$NH$_3$PbBr$_3$) and all-inorganic (CsPbBr$_3$) lead-halide perovskites show that anharmonic, local polar fluctuations are intrinsic to the general lead-halide perovskite structure, and not unique to the dipolar organic cation. MD simulations show that head-to-head Cs motion coupled to Br face expansion, on a few hundred femtosecond time scale, drives the local polar fluctuations in CsPbBr$_3$.
△ Less
Submitted 10 February, 2017; v1 submitted 27 April, 2016;
originally announced April 2016.
-
Structural and excited-state properties of oligoacene crystals from first principles
Authors:
Tonatiuh Rangel,
Kristian Berland,
Sahar Sharifzadeh,
Florian Brown-Altvater,
Kyuho Lee,
Per Hyldgaard,
Leeor Kronik,
Jeffrey B. Neaton
Abstract:
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 int…
▽ More
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.
△ Less
Submitted 31 March, 2016;
originally announced April 2016.
-
Are Mobilities in Hybrid Organic-Inorganic Halide Perovskites Actually 'High'?
Authors:
Thomas M. Brenner,
David A. Egger,
Andrew M. Rappe,
Leeor Kronik,
Gary Hodes,
David Cahen
Abstract:
We present an experimental and theoretical viewpoint on the electronic carrier mobilities of typical hybrid organic-inorganic perovskites (HOIPs). While these mobilities are often quoted as high, a review of them shows that although otherwise the semiconducting properties of HOIPs are impressively good, mobilities of HOIPs used in most solar cells are actually not that high. This is especially app…
▽ More
We present an experimental and theoretical viewpoint on the electronic carrier mobilities of typical hybrid organic-inorganic perovskites (HOIPs). While these mobilities are often quoted as high, a review of them shows that although otherwise the semiconducting properties of HOIPs are impressively good, mobilities of HOIPs used in most solar cells are actually not that high. This is especially apparent if they are compared to those of inorganic semiconductors used in other high efficiency solar cells. We critically examine possible causes and focus on electron-lattice coupling mechanisms that are active at room temperature, and can lead to carrier scattering. From this, we propose scattering due to acoustic phonons or polarons as possible causes, but also point out the difficulties with each of these in view of additional experimental and theoretical findings in the literature. Further research in this direction will contribute to making HOIP solar cells even more efficient than they already are.
△ Less
Submitted 29 October, 2015;
originally announced October 2015.
-
Probing the Orbital Origin of Conductance Oscillations in Atomic Chains
Authors:
Ran Vardimon,
Tamar Yelin,
Marina Klionsky,
Soumyajit Sarkar,
Ariel Biller,
Leeor Kronik,
Oren Tal
Abstract:
We investigate periodical oscillations in the conductance of suspended Au and Pt atomic chains during elongation under mechanical stress. Analysis of conductance and shot noise measurements reveals that the oscillations are mainly related to variations in a specific conduction channel as the chain undergoes transitions between zigzag and linear atomic configurations. The calculated local electroni…
▽ More
We investigate periodical oscillations in the conductance of suspended Au and Pt atomic chains during elongation under mechanical stress. Analysis of conductance and shot noise measurements reveals that the oscillations are mainly related to variations in a specific conduction channel as the chain undergoes transitions between zigzag and linear atomic configurations. The calculated local electronic structure shows that the oscillations originate from varying degrees of hybridization between the atomic orbitals along the chain as a function of the zigzag angle. These variations are highly dependent on the directionality and symmetry of the relevant orbitals, in agreement with the order-of-magnitude difference between the Pt and Au oscillation amplitudes observed in experiment. Our results demonstrate that the sensitivity of conductance to structural variations can be controlled by designing atomic-scale conductors in view of the directional interactions between atomic orbitals.
△ Less
Submitted 25 September, 2015;
originally announced September 2015.
-
One-electron self-interaction and the asymptotics of the Kohn-Sham potential: an impaired relation
Authors:
Tobias Schmidt,
Eli Kraisler,
Leeor Kronik,
Stephan Kümmel
Abstract:
One-electron self-interaction and an incorrect asymptotic behavior of the Kohn-Sham exchange-correlation potential are among the most prominent limitations of many present-day density functionals. However, a one-electron self-interaction-free energy does not necessarily lead to the correct long-range potential. This is here shown explicitly for local hybrid functionals. Furthermore, carefully stud…
▽ More
One-electron self-interaction and an incorrect asymptotic behavior of the Kohn-Sham exchange-correlation potential are among the most prominent limitations of many present-day density functionals. However, a one-electron self-interaction-free energy does not necessarily lead to the correct long-range potential. This is here shown explicitly for local hybrid functionals. Furthermore, carefully studying the ratio of the von Weizsäcker kinetic energy density to the (positive) Kohn-Sham kinetic energy density, $τ_\mathrm{W}/τ$, reveals that this ratio, which frequently serves as an iso-orbital indicator and is used to eliminate one-electron self-interaction effects in meta-generalized-gradient approximations and local hybrid functionals, can fail to approach its expected value in the vicinity of orbital nodal planes. This perspective article suggests that the nature and consequences of one-electron self-interaction and some of the strategies for its correction need to be reconsidered.
△ Less
Submitted 2 September, 2015;
originally announced September 2015.
-
Effect of ensemble generalization on the highest-occupied Kohn-Sham eigenvalue
Authors:
Eli Kraisler,
Tobias Schmidt,
Stephan Kümmel,
Leeor Kronik
Abstract:
There are several approximations to the exchange-correlation functional in density-functional theory that accurately predict total energy-related properties of many-electron systems, such as binding energies, bond lengths, and crystal structures. Other approximations are designed to describe potential-related processes, such as charge transfer and photoemission. However, the development of a funct…
▽ More
There are several approximations to the exchange-correlation functional in density-functional theory that accurately predict total energy-related properties of many-electron systems, such as binding energies, bond lengths, and crystal structures. Other approximations are designed to describe potential-related processes, such as charge transfer and photoemission. However, the development of a functional which can serve the two purposes simultaneously is a long-standing challenge. Trying to address it, we employ in the current work the ensemble generalization procedure proposed in Phys. Rev. Lett. 110, 126403 (2013). Focusing on the prediction of the ionization potential via the highest occupied Kohn-Sham eigenvalue, we examine a variety of exchange-correlation approximations: the local spin-density approximation, semi-local generalized gradient approximations, and global and local hybrid functionals. Results for a test set of 26 diatomic molecules and single atoms are presented. We find that the aforementioned ensemble generalization systematically improves the prediction of the ionization potential, for various systems and exchange-correlation functionals, without compromising the accuracy of total energy-related properties. We specifically examine hybrid functionals. These depend on a parameter controlling the ratio of semi-local to non-local functional components. The ionization potential obtained with ensemble-generalized functionals is found to depend only weakly on the parameter value, contrary to common experience with non-generalized hybrids, thus eliminating one aspect of the so-called `parameter dilemma' of hybrid functionals.
△ Less
Submitted 31 August, 2015;
originally announced August 2015.
-
Elimination of the asymptotic fractional dissociation problem in Kohn-Sham density functional theory using the ensemble-generalization approach
Authors:
Eli Kraisler,
Leeor Kronik
Abstract:
Many approximations within density-functional theory spuriously predict that a many-electron system can dissociate into fractionally charged fragments. Here, we revisit the case of dissociated diatomic molecules, known to exhibit this problem when studied within standard approaches, including the local spin-density approximation (LSDA). By employing our recently proposed [E. Kraisler and L. Kronik…
▽ More
Many approximations within density-functional theory spuriously predict that a many-electron system can dissociate into fractionally charged fragments. Here, we revisit the case of dissociated diatomic molecules, known to exhibit this problem when studied within standard approaches, including the local spin-density approximation (LSDA). By employing our recently proposed [E. Kraisler and L. Kronik, Phys. Rev. Lett. 110, 126403 (2013)] ensemble-generalization we find that asymptotic fractional dissociation is eliminated in all systems examined, even if the underlying exchange-correlation (xc) is still the LSDA. Furthermore, as a result of the ensemble generalization procedure, the Kohn-Sham potential develops a spatial step between the dissociated atoms, reflecting the emergence of the derivative discontinuity in the xc energy functional. This step, predicted in the past for the exact Kohn-Sham potential and observed in some of its more advanced approximate forms, is a desired feature that prevents any fractional charge transfer between the system's fragments. It is usually believed that simple xc approximations such as the LSDA cannot develop this step. Our findings show, however, that ensemble generalization to fractional electron densities automatically introduces the desired step even to the most simple approximate xc functionals and correctly predicts asymptotic integer dissociation.
△ Less
Submitted 30 August, 2015;
originally announced August 2015.
-
Fundamental gaps with approximate density functionals: the derivative discontinuity revealed from ensemble considerations
Authors:
Eli Kraisler,
Leeor Kronik
Abstract:
The fundamental gap is a central quantity in the electronic structure of matter. Unfortunately, the fundamental gap is not generally equal to the Kohn-Sham gap of density functional theory (DFT), even in principle. The two gaps differ precisely by the derivative discontinuity, namely, an abrupt change in slope of the exchange-correlation (xc) energy as a function of electron number, expected acros…
▽ More
The fundamental gap is a central quantity in the electronic structure of matter. Unfortunately, the fundamental gap is not generally equal to the Kohn-Sham gap of density functional theory (DFT), even in principle. The two gaps differ precisely by the derivative discontinuity, namely, an abrupt change in slope of the exchange-correlation (xc) energy as a function of electron number, expected across an integer-electron point. Popular approximate functionals are thought to be devoid of a derivative discontinuity, strongly compromising their performance for prediction of spectroscopic properties. Here we show that, in fact, all exchange-correlation functionals possess a derivative discontinuity, which arises naturally from the application of ensemble considerations within DFT, without any empiricism. This derivative discontinuity can be expressed in closed form using only quantities obtained in the course of a standard DFT calculation of the neutral system. For small, finite systems, addition of this derivative discontinuity indeed results in a greatly improved prediction for the fundamental gap, even when based on the most simple approximate exchange-correlation density functional - the local density approximation (LDA). For solids, the same scheme is exact in principle, but when applied to LDA it results in a vanishing derivative discontinuity correction. This failure is shown to be directly related to the failure of LDA in predicting fundamental gaps from total energy differences in extended systems.
△ Less
Submitted 30 August, 2015;
originally announced August 2015.
-
Experimental and theoretical electronic structure of quinacridone
Authors:
Daniel Lüftner,
Sivan Refaely-Abramson,
Michael Pachler,
Roland Resel,
Michael G. Ramsey,
Leeor Kronik,
Peter Puschnig
Abstract:
The energy positions of frontier orbitals in organic electronic materials are often studied experimentally by (inverse) photoemission spectroscopy and theoretically within density functional theory. However, standard exchange-correlation functionals often result in too small fundamental gaps, may lead to wrong orbital energy ordering, and do not capture polarization-induced gap renormalization. He…
▽ More
The energy positions of frontier orbitals in organic electronic materials are often studied experimentally by (inverse) photoemission spectroscopy and theoretically within density functional theory. However, standard exchange-correlation functionals often result in too small fundamental gaps, may lead to wrong orbital energy ordering, and do not capture polarization-induced gap renormalization. Here, we examine these issues and a strategy for overcoming them by studying the gas phase and bulk electronic structure of the organic molecule quinacridone (5Q), a promising material with many interesting properties for organic devices. Experimentally, we perform angle-resolved photoemission spectroscopy (ARUPS) on thin films of the crystalline $β$-phase of 5Q. Theoretically, we employ an optimally-tuned range-separated hybrid functional (OT-RSH) within density functional theory. For the gas phase molecule, our OT-RSH result for the ionization potential (IP) represents a substantial improvement over the semi-local PBE and the PBE0 hybrid functional results, producing an IP in quantitative agreement with experiment. For the bulk crystal, we take into account the correct screening in the bulk, using the recently developed OT-SRSH approach, while retaining the optimally-tuned parameters for the range-separation and the short-range Fock exchange. This leads to a band gap narrowing due to polarization effects and results in a valence band spectrum in excellent agreement with experimental ARUPS data, with respect to both peak positions and heights. Finally, full-frequency $G_0W_0$ results based on a hybrid functional starting point are shown to agree with the OT-SRSH approach, improving substantially on the PBE-starting point.
△ Less
Submitted 19 August, 2015;
originally announced August 2015.
-
Atomically wired molecular junctions: Connecting a single organic molecule by chains of metal atoms
Authors:
Tamar Yelin,
Ran Vardimon,
Natalia Kuritz,
Richard Korytár,
Alexei Bagrets,
Ferdinand Evers,
Leeor Kronik,
Oren Tal
Abstract:
Using a break junction technique, we find a clear signature for the formation of conducting hybrid junctions composed of a single organic molecule (benzene, naphthalene or anthracene) connected to chains of platinum atoms. The hybrid junctions exhibit metallic-like conductance (~0.1-1G0), which is rather insensitive to further elongation by additional atoms. At low bias voltage the hybrid junction…
▽ More
Using a break junction technique, we find a clear signature for the formation of conducting hybrid junctions composed of a single organic molecule (benzene, naphthalene or anthracene) connected to chains of platinum atoms. The hybrid junctions exhibit metallic-like conductance (~0.1-1G0), which is rather insensitive to further elongation by additional atoms. At low bias voltage the hybrid junctions can be elongated significantly beyond the length of the bare atomic chains. Ab initio calculations reveal that benzene based hybrid junctions have a significant binding energy and high structural flexibility that may contribute to the survival of the hybrid junction during the elongation process. The fabrication of hybrid junctions opens the way for combining the different properties of atomic chains and organic molecules to realize a new class of atomic scale interfaces.
△ Less
Submitted 22 June, 2015;
originally announced June 2015.
-
Solid-state optical absorption from optimally-tuned time-dependent range-separated hybrid density functional theory
Authors:
Sivan Refaely-Abramson,
Manish Jain,
Sahar Sharifzadeh,
Jeffrey B. Neaton,
Leeor Kronik
Abstract:
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 o…
▽ More
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.
△ Less
Submitted 14 August, 2015; v1 submitted 7 May, 2015;
originally announced May 2015.
-
A multi-scale approach to the electronic structure of doped semiconductor surfaces
Authors:
Ofer Sinai,
Oliver T. Hofmann,
Patrick Rinke,
Matthias Scheffler,
Georg Heimel,
Leeor Kronik
Abstract:
The inclusion of the global effects of semiconductor doping poses a unique challenge for first-principles simulations, because the typically low concentration of dopants renders an explicit treatment intractable. In systems which do not display long-range band bending, a satisfactory remedy is offered by the use of "pseudoatoms", with a fractional nuclear charge matching the bulk doping concentrat…
▽ More
The inclusion of the global effects of semiconductor doping poses a unique challenge for first-principles simulations, because the typically low concentration of dopants renders an explicit treatment intractable. In systems which do not display long-range band bending, a satisfactory remedy is offered by the use of "pseudoatoms", with a fractional nuclear charge matching the bulk doping concentration. However, this alone is not always sufficient for charged surfaces, where the width of the space-charge region (SCR) often exceeds realistic supercell dimensions. One generalization of the pseudoatom approach which overcomes this difficulty relies on the introduction of an artificially high doping level within a slab calculation, in conjunction with a multi-scale electrostatic energy correction. Here, we present an alternative technique that naturally extends the pseudoatom approach while bypassing the need for calculations with an unrealistically high doping level. It is based on the introduction of excess charge, mimicking free charge carriers from the SCR, along with a fixed sheet of counter-charge mimicking the SCR-related field. Self-consistency is obtained by imposing charge conservation and Fermi level equilibration between the bulk, treated semi-classically, and the electronic states of the slab/surface, which are treated quantum-mechanically. The method, which we call CREST - the Charge-Reservoir Electrostatic Sheet Technique - can be used with standard electronic structure codes. We validate the approach using a simple tight-binding model, which allows for comparison of its results with calculations encompassing the full SCR explicitly. We then employ it with density functional theory, where it is used to obtain insights into the electronic structures of the "clean-cleaved" Si(111) surface and its buckled (2x1) reconstruction, at various doping densities.
△ Less
Submitted 11 November, 2014;
originally announced November 2014.
-
Deviations from piecewise linearity in the solid-state limit with approximate density functionals
Authors:
Vojtěch Vlček,
Helen R. Eisenberg,
Gerd Steinle-Neumann,
Leeor Kronik,
Roi Baer
Abstract:
In exact density functional theory (DFT) the total ground-state energy is a series of linear segments between integer electron points, a condition known as "piecewise linearity". Deviation from this condition is indicative of poor predictive capabilities for electronic structure, in particular of ionization energies, fundamental gaps, and charge transfer. In this article, we take a new look at the…
▽ More
In exact density functional theory (DFT) the total ground-state energy is a series of linear segments between integer electron points, a condition known as "piecewise linearity". Deviation from this condition is indicative of poor predictive capabilities for electronic structure, in particular of ionization energies, fundamental gaps, and charge transfer. In this article, we take a new look at the deviation from linearity (i.e., curvature) in the solid-state limit by considering two different ways of approaching it: a large finite system of increasing size and a crystal represented by an increasingly large reference cell with periodic boundary conditions. We show that the curvature approaches vanishing values in both limits, even for functionals which yield poor predictions of electronic structure, and therefore can not be used as a diagnostic or constructive tool in solids. We find that the approach towards zero curvature is different in each of the two limits, owing to the presence of a compensating background charge in the periodic case. Based on these findings, we present a new criterion for functional construction and evaluation, derived from the size-dependence of the curvature, along with a practical method for evaluating this criterion. For large finite systems we further show that the curvature is dominated by the self-interaction of the highest occupied eigenstate. These findings are illustrated by computational studies of various solids, semiconductor nanocrystals, and long alkane chains.
△ Less
Submitted 16 September, 2014;
originally announced September 2014.
-
A self-interaction-free local hybrid functional: Accurate binding energies vis-à-vis accurate ionization potentials from Kohn-Sham eigenvalues
Authors:
Tobias Schmidt,
Eli Kraisler,
Adi Makmal,
Leeor Kronik,
Stephan Kümmel
Abstract:
We present and test a new approximation for the exchange-correlation (xc) energy of Kohn-Sham density functional theory. It combines exact exchange with a compatible non-local correlation functional. The functional is by construction free of one-electron self-interaction, respects constraints derived from uniform coordinate scaling, and has the correct asymptotic behavior of the xc energy density.…
▽ More
We present and test a new approximation for the exchange-correlation (xc) energy of Kohn-Sham density functional theory. It combines exact exchange with a compatible non-local correlation functional. The functional is by construction free of one-electron self-interaction, respects constraints derived from uniform coordinate scaling, and has the correct asymptotic behavior of the xc energy density. It contains one parameter that is not determined ab initio. We investigate whether it is possible to construct a functional that yields accurate binding energies and affords other advantages, specifically Kohn-Sham eigenvalues that reliably reflect ionization potentials. Tests for a set of atoms and small molecules show that within our local-hybrid form accurate binding energies can be achieved by proper optimization of the free parameter in our functional, along with an improvement in dissociation energy curves and in Kohn-Sham eigenvalues. However, the correspondence of the latter to experimental ionization potentials is not yet satisfactory, and if we choose to optimize their prediction, a rather different value of the functional's parameter is obtained. We put this finding in a larger context by discussing similar observations for other functionals and possible directions for further functional development that our findings suggest.
△ Less
Submitted 22 May, 2014;
originally announced May 2014.
-
Voltage tuning of vibrational mode energies in single-molecule junctions
Authors:
Y. Li,
P. Doak,
L. Kronik,
J. B. Neaton,
D. Natelson
Abstract:
Vibrational modes of molecules are fundamental properties determined by intramolecular bonding, atomic masses, and molecular geometry, and often serve as important channels for dissipation in nanoscale processes. Although single-molecule junctions have been employed to manipulate electronic structure and related functional properties of molecules, electrical control of vibrational mode energies ha…
▽ More
Vibrational modes of molecules are fundamental properties determined by intramolecular bonding, atomic masses, and molecular geometry, and often serve as important channels for dissipation in nanoscale processes. Although single-molecule junctions have been employed to manipulate electronic structure and related functional properties of molecules, electrical control of vibrational mode energies has remained elusive. Here we use simultaneous transport and surface-enhanced Raman spectroscopy measurements to demonstrate large, reversible, voltage-driven shifts of vibrational mode energies of C60 molecules in gold junctions. C60 mode energies are found to vary approximately quadratically with bias, but in a manner inconsistent with a simple vibrational Stark effect. Our theoretical model suggests instead that the mode shifts are a signature of bias-driven addition of electronic charge to the molecule. These results imply that voltage-controlled tuning of vibrational modes is a general phenomenon at metal-molecule interfaces and is a means of achieving significant shifts in vibrational energies relative to a pure Stark effect.
△ Less
Submitted 4 March, 2014;
originally announced March 2014.
-
A State Representation Approach for Atomistic Time-Dependent Transport Calculations in Molecular Junctions
Authors:
Tamar Zelovich,
Leeor Kronik,
Oded Hod
Abstract:
We propose a new method for simulating electron dynamics in open quantum systems out of equilibrium, using a finite atomistic model. The proposed method is motivated by the intuitive and practical nature of the driven Liouville von-Neumann equation approach of Sánchez et al. [J. Chem. Phys., 124, 214708 (2006)]. A key ingredient of our approach is a transformation of the Hamiltonian matrix from an…
▽ More
We propose a new method for simulating electron dynamics in open quantum systems out of equilibrium, using a finite atomistic model. The proposed method is motivated by the intuitive and practical nature of the driven Liouville von-Neumann equation approach of Sánchez et al. [J. Chem. Phys., 124, 214708 (2006)]. A key ingredient of our approach is a transformation of the Hamiltonian matrix from an atomistic to a state representation of the molecular junction. This allows us to uniquely define the bias voltage across the system while maintaining a proper thermal electronic distribution within the finite lead models. Furthermore, it allows us to investigate complex molecular junctions, including non-linear setups and multi-lead configurations. A heuristic derivation of our working equation leads to explicit expressions for the damping and driving terms, which serve as appropriate electron sources and sinks that effectively "open" the finite model system. Although the method does not forbid it, in practice we find neither violation of Pauli's exclusion principles nor deviation from density matrix positivity throughout our numerical simulations of various tight-binding model systems. We believe that the new approach offers a practical and physically sound route for performing atomistic time-dependent transport calculations in realistic molecular junctions.
△ Less
Submitted 20 February, 2014;
originally announced February 2014.
-
A real-space method for highly parallelizable electronic transport calculations
Authors:
Baruch Feldman,
Tamar Seideman,
Oded Hod,
Leeor Kronik
Abstract:
We present a real-space method for first-principles nano-scale electronic transport calculations. We use the non-equilibrium Green's function method with density functional theory and implement absorbing boundary conditions (ABCs, also known as complex absorbing potentials, or CAPs) to represent the effects of the semi-infinite leads. In real space, the Kohn-Sham Hamiltonian matrix is highly spars…
▽ More
We present a real-space method for first-principles nano-scale electronic transport calculations. We use the non-equilibrium Green's function method with density functional theory and implement absorbing boundary conditions (ABCs, also known as complex absorbing potentials, or CAPs) to represent the effects of the semi-infinite leads. In real space, the Kohn-Sham Hamiltonian matrix is highly sparse. As a result, the transport problem parallelizes naturally and can scale favorably with system size, enabling the computation of conductance in relatively large molecular junction models. Our use of ABCs circumvents the demanding task of explicitly calculating the leads' self-energies from surface Green's functions, and is expected to be more accurate than the use of the jellium approximation. In addition, we take advantage of the sparsity in real space to solve efficiently for the Green's function over the entire energy range relevant to low-bias transport. We illustrate the advantages of our method with calculations on several challenging test systems and find good agreement with reference calculation results.
△ Less
Submitted 2 July, 2014; v1 submitted 4 January, 2014;
originally announced January 2014.
-
Inter-layer Potential for Hexagonal Boron Nitride
Authors:
Itai Leven,
Ido Azuri,
Leeor Kronik,
Oded Hod
Abstract:
A new interlayer force-field for layered hexagonal boron nitride (h-BN) based structures is presented. The force-field contains three terms representing the interlayer attraction due to dispersive interactions, repulsion due to anisotropic overlaps of electron clouds, and monopolar electrostatic interactions. With appropriate parameterization, the potential is able to simultaneously capture well t…
▽ More
A new interlayer force-field for layered hexagonal boron nitride (h-BN) based structures is presented. The force-field contains three terms representing the interlayer attraction due to dispersive interactions, repulsion due to anisotropic overlaps of electron clouds, and monopolar electrostatic interactions. With appropriate parameterization, the potential is able to simultaneously capture well the binding and lateral sliding energies of planar h-BN based dimer systems as well as the interlayer telescoping and rotation of double walled boron-nitride nanotubes of different crystallographic orientations. The new potential thus enables the accurate and efficient modeling and simulation of large-scale h-BN based layered structures.
△ Less
Submitted 10 October, 2013;
originally announced October 2013.
-
Piecewise Linearity of Approximate Density Functionals Revisited: Implications for Frontier Orbital Energies
Authors:
Eli Kraisler,
Leeor Kronik
Abstract:
In the exact Kohn-Sham density-functional theory (DFT), the total energy versus the number of electrons is a series of linear segments between integer points. However, commonly used approximate density functionals produce total energies that do not exhibit this piecewise-linear behavior. As a result, the ionization potential theorem, equating the highest occupied eigenvalue with the ionization pot…
▽ More
In the exact Kohn-Sham density-functional theory (DFT), the total energy versus the number of electrons is a series of linear segments between integer points. However, commonly used approximate density functionals produce total energies that do not exhibit this piecewise-linear behavior. As a result, the ionization potential theorem, equating the highest occupied eigenvalue with the ionization potential, is grossly disobeyed. Here, we show that, contrary to conventional wisdom, most of the required piecewise-linearity of an arbitrary approximate density functional can be restored by careful consideration of the ensemble generalization of DFT. Furthermore, the resulting formulation introduces the desired derivative discontinuity to any approximate exchange-correlation functional, even one that is explicitly density-dependent. This opens the door to calculations of the ionization potential and electron affinity even without explicit electron removal or addition. All these advances are achieved while neither introducing empiricism nor changing the underlying functional form. The power of the approach is demonstrated on benchmark systems using the local density approximation as an illustrative example.
△ Less
Submitted 2 April, 2013; v1 submitted 26 November, 2012;
originally announced November 2012.
-
Curvature and frontier orbital energies in density functional theory
Authors:
Tamar Stein,
Jochen Autschbach,
Niranjan Govind,
Leeor Kronik,
Roi Baer
Abstract:
Perdew et al. [Phys. Rev. Lett 49, 1691 (1982)] discovered and proved two different properties of exact Kohn-Sham density functional theory (DFT): (i) The exact total energy versus particle number is a series of linear segments between integer electron points; (ii) Across an integer number of electrons, the exchange-correlation potential may "jump" by a constant, known as the derivative discontinu…
▽ More
Perdew et al. [Phys. Rev. Lett 49, 1691 (1982)] discovered and proved two different properties of exact Kohn-Sham density functional theory (DFT): (i) The exact total energy versus particle number is a series of linear segments between integer electron points; (ii) Across an integer number of electrons, the exchange-correlation potential may "jump" by a constant, known as the derivative discontinuity (DD). Here, we show analytically that in both the original and the generalized Kohn-Sham formulation of DFT, the two are in fact two sides of the same coin. Absence of a derivative discontinuity necessitates deviation from piecewise linearity, and the latter can be used to correct for the former, thereby restoring the physical meaning of the orbital energies. Using selected small molecules, we show that this results in a simple correction scheme for any underlying functional, including semi-local and hybrid functionals as well as Hartree-Fock theory, suggesting a practical correction for the infamous gap problem of DFT. Moreover, we show that optimally-tuned range-separated hybrid functionals can inherently minimize both DD and curvature, thus requiring no correction, and show that this can be used as a sound theoretical basis for novel tuning strategies.
△ Less
Submitted 7 August, 2012;
originally announced August 2012.
-
Quasiparticle spectra from a non-empirical optimally-tuned range-separated hybrid density functional
Authors:
Sivan Refaely-Abramson,
Sahar Sharifzadeh,
Niranjan Govind,
Jochen Autschbach,
Jeffrey B. Neaton,
Roi Baer,
Leeor Kronik
Abstract:
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 sep…
▽ More
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.
△ Less
Submitted 6 September, 2012; v1 submitted 11 March, 2012;
originally announced March 2012.
-
Quasiparticle and Optical Spectroscopy of Organic Semiconductors Pentacene and PTCDA from First Principles
Authors:
Sahar Sharifzadeh,
Ariel Biller,
Leeor Kronik,
Jeffrey B. Neaton
Abstract:
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…
▽ More
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.
△ Less
Submitted 1 February, 2012; v1 submitted 21 October, 2011;
originally announced October 2011.
-
Charge transport across metal/molecular (alkyl) monolayer-Si junctions is dominated by the LUMO level
Authors:
Omer Yaffe,
Yabing Qi,
Lior Segev,
Luc Scheres,
Sreenivasa Reddy Puniredd,
Tal Ely,
Hossam Haick,
Han Zuilhof,
Leeor Kronik,
Antoine Kahn,
Ayelet Vilan,
David Cahen
Abstract:
We compare the charge transport characteristics of heavy doped p- and n-Si-alkyl chain/Hg junctions. Photoelectron spectroscopy (UPS, IPES and XPS) results for the molecule-Si band alignment at equilibrium show the Fermi level to LUMO energy difference to be much smaller than the corresponding Fermi level to HOMO one. This result supports the conclusion we reach, based on negative differential res…
▽ More
We compare the charge transport characteristics of heavy doped p- and n-Si-alkyl chain/Hg junctions. Photoelectron spectroscopy (UPS, IPES and XPS) results for the molecule-Si band alignment at equilibrium show the Fermi level to LUMO energy difference to be much smaller than the corresponding Fermi level to HOMO one. This result supports the conclusion we reach, based on negative differential resistance in an analogous semiconductor-inorganic insulator/metal junction, that for both p- and n-type junctions the energy difference between the Fermi level and LUMO, i.e., electron tunneling, controls charge transport. The Fermi level-LUMO energy difference, experimentally determined by IPES, agrees with the non-resonant tunneling barrier height deduced from the exponential length-attenuation of the current.
△ Less
Submitted 26 October, 2011; v1 submitted 20 June, 2011;
originally announced June 2011.