1
|
Brown J. Calculating Potential Energy Surfaces with Quantum Computers by Measuring Only the Density Along Adiabatic Transitions. J Chem Theory Comput 2024; 20:3096-3108. [PMID: 38602483 DOI: 10.1021/acs.jctc.3c01177] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/12/2024]
Abstract
We show that chemically accurate potential energy surfaces (PESs) can be generated from quantum computers by measuring only the density along an adiabatic transition between different molecular geometries. In lieu of using phase estimation, the energy is evaluated by performing line-integration using the inverted real-space time-dependent density functional theory Kohn-Sham (KS) potential obtained from the geometry-varying densities of the full wave function. The accuracy of this method depends on the validity of the adiabatic evolution itself and the potential inversion process (which is theoretically exact but can be numerically unstable), whereas the total evolution time is the defining factor for the precision of phase estimation. We examine the method with a one-dimensional system of two electrons for both the ground and first triplet states in first quantization, as well as the ground state of three- and four-electron systems in second quantization. It is shown that few accurate measurements can be utilized to obtain chemical accuracy across the full potential energy curve, with a shorter propagation time than may be required using phase estimation for a similar accuracy. We also show that an accurate potential energy curve can be calculated by making many imprecise density measurements (using a few shots) along the time evolution and smoothing the resulting density evolution. Finally, it is important to note that the method is able to classically provide a check of its own accuracy by comparing the density resulting from a time-independent KS calculation using the inverted potential with the measured density. This can be used to determine whether longer adiabatic evolution times are required to satisfy the adiabatic theorem.
Collapse
Affiliation(s)
- James Brown
- Good Chemistry Company, 200-1285 West Pender Street, Vancouver, British Columbia V6E 4B1, Canada
| |
Collapse
|
2
|
Hait D, Head-Gordon M. When Is a Bond Broken? The Polarizability Perspective. Angew Chem Int Ed Engl 2023; 62:e202312078. [PMID: 37713599 DOI: 10.1002/anie.202312078] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/18/2023] [Revised: 09/15/2023] [Accepted: 09/15/2023] [Indexed: 09/17/2023]
Abstract
The question of when a chemical bond can be said to be broken is of fundamental chemical interest but has not been widely studied. Herein we propose that the maxima of static polarizability along bond dissociation coordinates naturally define cutoff points for bond rupture, as they represent the onset of localization of shared electron density into constituent fragments. Examples of computed polarizability maxima over the course of bond cleavage in main-group and transition metal compounds are provided, across covalent, dative and charge-shift bonds. The behavior along reaction paths is also considered. Overall, the static polarizability is found to be a sensitive reporter of electronic structure reorganization associated with bond stretching, and thus can serve as a metric for describing bond cleavage (or diagnose the absence of a chemical bond).
Collapse
Affiliation(s)
- Diptarka Hait
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, Berkeley, CA 94720, USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA
- Department of Chemistry and PULSE Institute, Stanford University, Stanford, CA 94305, USA
| | - Martin Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, Berkeley, CA 94720, USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA
| |
Collapse
|
3
|
Hellgren M, Baguet L. Strengths and limitations of the adiabatic exact-exchange kernel for total energy calculations. J Chem Phys 2023; 158:2889488. [PMID: 37158324 DOI: 10.1063/5.0146423] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/13/2023] [Accepted: 04/24/2023] [Indexed: 05/10/2023] Open
Abstract
We investigate the adiabatic approximation to the exact-exchange kernel for calculating correlation energies within the adiabatic-connection fluctuation-dissipation framework of time-dependent density functional theory. A numerical study is performed on a set of systems having bonds of different character (H2 and N2 molecules, H-chain, H2-dimer, solid-Ar, and the H2O-dimer). We find that the adiabatic kernel can be sufficient in strongly bound covalent systems, yielding similar bond lengths and binding energies. However, for non-covalent systems, the adiabatic kernel introduces significant errors around equilibrium geometry, systematically overestimating the interaction energy. The origin of this behavior is investigated by studying a model dimer composed of one-dimensional, closed-shell atoms, interacting via soft-Coulomb potentials. The kernel is shown to exhibit a strong frequency dependence at small to intermediate atomic separation that affects both the low-energy spectrum and the exchange-correlation hole obtained from the corresponding diagonal of the two-particle density matrix.
Collapse
Affiliation(s)
- Maria Hellgren
- Sorbonne Université, MNHN, UMR CNRS 7590, IMPMC, 4 place Jussieu, 75005 Paris, France
| | - Lucas Baguet
- CEA, DAM, DIF, F-91297 Arpajon, France
- Université Paris-Saclay, CEA, Laboratoire Matière en Conditions Extrêmes, 91680 Bruyères-le-Châtel, France
| |
Collapse
|
4
|
Giarrusso S, Neugarten R, Baerends EJ, Giesbertz KJH. Secondary Kinetic Peak in the Kohn-Sham Potential and Its Connection to the Response Step. J Chem Theory Comput 2022; 18:4762-4773. [PMID: 35895974 DOI: 10.1021/acs.jctc.2c00332] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We consider a prototypical 1D model Hamiltonian for a stretched heteronuclear molecule and construct individual components of the corresponding KS potential, namely, the kinetic, the N - 1, and the conditional potentials. These components show very special features, such as peaks and steps, in regions where the density is drastically low. Some of these features are quite well-known, whereas others, such as a secondary peak in the kinetic potential or a second bump in the conditional potential, are less or not known at all. We discuss these features building on the analytical model treated in Giarrusso et al. J. Chem. Theory Comput. 2018, 14, 4151. In particular, we provide an explanation for the underlying mechanism which determines the appearance of both peaks in the kinetic potential and elucidate why these peaks delineate the region over which the plateau structure, due to the N - 1 potential, stretches. We assess the validity of the Heitler-London Ansatz at large but finite internuclear distance, showing that, if optimal orbitals are used, this model is an excellent approximation to the exact wave function. Notably, we find that the second natural orbital presents an extra node very far out on the side of the more electronegative atom.
Collapse
Affiliation(s)
- Sara Giarrusso
- Department of Chemistry and Biochemistry, University of California Merced, 5200 North Lake Road, Merced, California 95343, United States
| | - Roeland Neugarten
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling, FEW, Vrije Universiteit, De Boelelaan 1083, 1081HV Amsterdam, The Netherlands
| | - Evert Jan Baerends
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling, FEW, Vrije Universiteit, De Boelelaan 1083, 1081HV Amsterdam, The Netherlands
| | - Klaas J H Giesbertz
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling, FEW, Vrije Universiteit, De Boelelaan 1083, 1081HV Amsterdam, The Netherlands
| |
Collapse
|
5
|
Erhard J, Trushin E, Görling A. Numerically stable inversion approach to construct Kohn-Sham potentials for given electron densities within a Gaussian basis set framework. J Chem Phys 2022; 156:204124. [PMID: 35649824 DOI: 10.1063/5.0087356] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present a Kohn-Sham (KS) inversion approach to construct KS exchange-correlation potentials corresponding to given electron densities. This method is based on an iterative procedure using linear response to update potentials. All involved quantities, i.e., orbitals, potentials, and response functions, are represented by Gaussian basis functions. In contrast to previous KS inversion methods relying on Gaussian basis sets, the method presented here is numerically stable even for standard basis sets from basis set libraries due to a preprocessing of the auxiliary basis used to represent an exchange-correlation charge density that generates the exchange-correlation potential. The new KS inversion method is applied to reference densities of various atoms and molecules obtained by full configuration interaction or CCSD(T) (coupled cluster singles doubles perturbative triples). The considered examples encompass cases known to be difficult, such as stretched hydrogen or lithium hydride molecules or the beryllium isoelectronic series. For the stretched hydrogen molecule, potentials of benchmark quality are obtained by employing large basis sets. For the carbon monoxide molecule, we show that the correlation potential from the random phase approximation (RPA) is in excellent qualitative and quantitative agreement with the correlation potential from the KS inversion of a CCSD(T) reference density. This indicates that RPA correlation potentials, in contrast to those from semi-local density-functionals, resemble the exact correlation potential. Besides providing exchange-correlation potentials for benchmark purposes, the proposed KS inversion method may be used in density-partition-based quantum embedding and in subsystem density-functional methods because it combines numerical stability with computational efficiency.
Collapse
Affiliation(s)
- Jannis Erhard
- Lehrstuhl für Theoretische Chemie, Universität Erlangen-Nürnberg, Egerlandstr. 3, D-91058 Erlangen, Germany
| | - Egor Trushin
- Lehrstuhl für Theoretische Chemie, Universität Erlangen-Nürnberg, Egerlandstr. 3, D-91058 Erlangen, Germany
| | - Andreas Görling
- Lehrstuhl für Theoretische Chemie, Universität Erlangen-Nürnberg, Egerlandstr. 3, D-91058 Erlangen, Germany
| |
Collapse
|
6
|
Baerends EJ. Chemical potential, derivative discontinuity, fractional electrons, jump of the Kohn-Sham potential, atoms as thermodynamic open systems, and other (mis)conceptions of the density functional theory of electrons in molecules. Phys Chem Chem Phys 2022; 24:12745-12766. [PMID: 35593143 DOI: 10.1039/d2cp01585d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Many references exist in the density functional theory (DFT) literature to the chemical potential of the electrons in an atom or a molecule. The origin of this notion has been the identification of the Lagrange multiplier μ = ∂E/∂N in the Euler-Lagrange variational equation for the ground state density as the chemical potential of the electrons. We first discuss why the Lagrange multiplier in this case is an arbitrary constant and therefore cannot be a physical characteristic of an atom or molecule. The switching of the energy derivative ("chemical potential") from -I to -A when the electron number crosses the integer, called integer discontinuity or derivative discontinuity, is not physical but only occurs when the nonphysical noninteger electron systems and the corresponding energy and derivative ∂E/∂N are chosen in a specific discontinuous way. The question is discussed whether in fact the thermodynamical concept of a chemical potential can be defined for the electrons in such few-electron systems as atoms and molecules. The conclusion is that such systems lack important characteristics of thermodynamic systems and do not afford the definition of a chemical potential. They also cannot be considered as analogues of the open systems of thermodynamics that can exchange particles with an environment (a particle bath or other members of a Gibbsian ensemble). Thermodynamical (statistical mechanical) concepts like chemical potential, open systems, grand canonical ensemble etc. are not applicable to a few electron system like an atom or molecule. A number of topics in DFT are critically reviewed in light of these findings: jumps in the Kohn-Sham potential when crossing an integer number of electrons, the band gap problem, the deviation-from-straight-lines error, and the role of ensembles in DFT.
Collapse
|
7
|
Shi Y, Chávez VH, Wasserman A. n2v
: A density‐to‐potential inversion suite. A sandbox for creating, testing, and benchmarking density functional theory inversion methods. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2022. [DOI: 10.1002/wcms.1617] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Affiliation(s)
- Yuming Shi
- Department of Physics and Astronomy Purdue University West Lafayette Indiana USA
| | - Victor H. Chávez
- Department of Chemistry Purdue University West Lafayette Indiana USA
| | - Adam Wasserman
- Department of Physics and Astronomy Purdue University West Lafayette Indiana USA
- Department of Chemistry Purdue University West Lafayette Indiana USA
| |
Collapse
|
8
|
Maitra NT. Double and Charge-Transfer Excitations in Time-Dependent Density Functional Theory. Annu Rev Phys Chem 2021; 73:117-140. [PMID: 34910562 DOI: 10.1146/annurev-physchem-082720-124933] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
Time-dependent density functional theory has emerged as a method of choice for calculations of spectra and response properties in physics, chemistry, and biology, with its system-size scaling enabling computations on systems much larger than otherwise possible. While increasingly complex and interesting systems have been successfully tackled with relatively simple functional approximations, there has also been increasing awareness that these functionals tend to fail for certain classes of approximations. Here I review the fundamental challenges the approximate functionals have in describing double excitations and charge-transfer excitations, which are two of the most common impediments for the theory to be applied in a black-box way. At the same time, I describe the progress made in recent decades in developing functional approximations that give useful predictions for these excitations. Expected final online publication date for the Annual Review of Physical Chemistry, Volume 73 is April 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.
Collapse
Affiliation(s)
- Neepa T Maitra
- Department of Physics, Rutgers University at Newark, Newark, New Jersey, USA;
| |
Collapse
|
9
|
Vieira D. Readdressing molecular dissociation within the Kohn–Sham formalism of density-functional theory: simple models and a different point of view. Mol Phys 2021. [DOI: 10.1080/00268976.2021.2008037] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
Affiliation(s)
- Daniel Vieira
- Departamento de Física, Programa de Pós-Graduação em Física, Universidade do Estado de Santa Catarina, Joinville, SC, Brazil
| |
Collapse
|
10
|
Kraisler E, Hodgson MJP, Gross EKU. From Kohn-Sham to Many-Electron Energies via Step Structures in the Exchange-Correlation Potential. J Chem Theory Comput 2021; 17:1390-1407. [PMID: 33595312 PMCID: PMC8363072 DOI: 10.1021/acs.jctc.0c01093] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
Accurately
describing excited states within Kohn–Sham (KS)
density functional theory (DFT), particularly those which induce ionization
and charge transfer, remains a great challenge. Common exchange-correlation
(xc) approximations are unreliable for excited states owing, in part,
to the absence of a derivative discontinuity in the xc energy (Δ),
which relates a many-electron energy difference to the corresponding
KS energy difference. We demonstrate, analytically and numerically,
how the relationship between KS and many-electron energies leads to
the step structures observed in the exact xc potential in four scenarios:
electron addition, molecular dissociation, excitation of a finite
system, and charge transfer. We further show that steps in the potential
can be obtained also with common xc approximations, as simple as the
LDA, when addressed from the ensemble perspective. The article therefore
highlights how capturing the relationship between KS and many-electron
energies with advanced xc approximations is crucial for accurately
calculating excitations, as well as the ground-state density and energy
of systems which consist of distinct subsystems.
Collapse
Affiliation(s)
- Eli Kraisler
- Fritz Haber Center for Molecular Dynamics and Institute of Chemistry, The Hebrew University of Jerusalem, 9091401 Jerusalem, Israel
| | - M J P Hodgson
- Department of Physics, Durham University, South Road, Durham DH1 3LE, United Kingdom.,Max-Planck-Institut für Mikrostrukturphysik, Weinberg 2, D-06120 Halle, Germany
| | - E K U Gross
- Fritz Haber Center for Molecular Dynamics and Institute of Chemistry, The Hebrew University of Jerusalem, 9091401 Jerusalem, Israel
| |
Collapse
|
11
|
Giarrusso S, Gori-Giorgi P. Exchange-Correlation Energy Densities and Response Potentials: Connection between Two Definitions and Analytical Model for the Strong-Coupling Limit of a Stretched Bond. J Phys Chem A 2020; 124:2473-2482. [PMID: 32118422 PMCID: PMC7104238 DOI: 10.1021/acs.jpca.9b10538] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
![]()
We
analyze in depth two widely used definitions (from the theory
of conditional probability amplitudes and from the adiabatic connection
formalism) of the exchange-correlation energy density and of the response
potential of Kohn–Sham density functional theory. We introduce
a local form of the coupling-constant-dependent Hohenberg–Kohn functional, showing that
the difference between the two definitions is due to a corresponding
local first-order term in the coupling constant, which disappears
globally (when integrated over all space), but not locally. We also
design an analytic representation for the response potential in the
strong-coupling limit of density functional theory for a model single
stretched bond.
Collapse
Affiliation(s)
- Sara Giarrusso
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling, FEW, Vrije Universiteit, De Boelelaan 1083, Amsterdam, 1081HV, The Netherlands
| | - Paola Gori-Giorgi
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling, FEW, Vrije Universiteit, De Boelelaan 1083, Amsterdam, 1081HV, The Netherlands
| |
Collapse
|
12
|
Kraisler E. Asymptotic Behavior of the Exchange‐Correlation Energy Density and the Kohn‐Sham Potential in Density Functional Theory: Exact Results and Strategy for Approximations. Isr J Chem 2020. [DOI: 10.1002/ijch.201900103] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Affiliation(s)
- Eli Kraisler
- Fritz Haber Center for Molecular Dynamics and Institute of Chemistry The Hebrew University of Jerusalem 9091401 Jerusalem Israel
| |
Collapse
|
13
|
Gerolin A, Grossi J, Gori-Giorgi P. Kinetic Correlation Functionals from the Entropic Regularization of the Strictly Correlated Electrons Problem. J Chem Theory Comput 2020; 16:488-498. [PMID: 31855421 PMCID: PMC6964418 DOI: 10.1021/acs.jctc.9b01133] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/13/2019] [Indexed: 11/29/2022]
Abstract
In this work, we study the entropic regularization of the strictly correlated electrons formalism, discussing the implications for density functional theory and establishing a link with earlier works on quantum kinetic energy and classical entropy. We carry out a very preliminary investigation (using simplified models) on the use of the solution of the entropic regularized problem to build approximations for the kinetic correlation functional at large coupling strengths. We also analyze lower and upper bounds to the Hohenberg-Kohn functional using the entropic regularized strictly correlated electrons problem.
Collapse
Affiliation(s)
- Augusto Gerolin
- Department of Theoretical
Chemistry and Amsterdam Center for Multiscale Modeling, FEW, Vrije Universiteit, De Boelelaan 1083, 1081HV Amsterdam, The Netherlands
| | - Juri Grossi
- Department of Theoretical
Chemistry and Amsterdam Center for Multiscale Modeling, FEW, Vrije Universiteit, De Boelelaan 1083, 1081HV Amsterdam, The Netherlands
| | - Paola Gori-Giorgi
- Department of Theoretical
Chemistry and Amsterdam Center for Multiscale Modeling, FEW, Vrije Universiteit, De Boelelaan 1083, 1081HV Amsterdam, The Netherlands
| |
Collapse
|
14
|
Gould T, Vuckovic S. Range-separation and the multiple radii functional approximation inspired by the strongly interacting limit of density functional theory. J Chem Phys 2019; 151:184101. [DOI: 10.1063/1.5125692] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Tim Gould
- Qld Micro- and Nanotechnology Centre, Griffith University, Nathan, Qld 4111, Australia
| | - Stefan Vuckovic
- Department of Chemistry, University of California, Irvine, California 92697, USA
| |
Collapse
|
15
|
Hellgren M, Gould T. Strong Correlation and Charge Localization in Kohn–Sham Theories with Fractional Orbital Occupations. J Chem Theory Comput 2019; 15:4907-4914. [DOI: 10.1021/acs.jctc.9b00477] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Maria Hellgren
- Sorbonne Université, Muséum National d’Histoire Naturelle, UMR CNRS 7590, Institut de Minéralogie, de Physique des Matériaux et de Cosmochimie (IMPMC), 4 place Jussieu, 75005 Paris, France
| | - Tim Gould
- Queensland Micro- and Nanotechnology Centre, Griffith University, Nathan, QLD 4111, Australia
| |
Collapse
|
16
|
Baerends EJ. On derivatives of the energy with respect to total electron number and orbital occupation numbers. A critique of Janak's theorem. Mol Phys 2019. [DOI: 10.1080/00268976.2019.1612955] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
17
|
Fabiano E, Śmiga S, Giarrusso S, Daas KJ, Della Sala F, Grabowski I, Gori-Giorgi P. Investigation of the Exchange-Correlation Potentials of Functionals Based on the Adiabatic Connection Interpolation. J Chem Theory Comput 2019; 15:1006-1015. [PMID: 30620596 DOI: 10.1021/acs.jctc.8b01037] [Citation(s) in RCA: 22] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We have studied the correlation potentials produced by various adiabatic connection models (ACMs) for several atoms and molecules. The results have been compared to accurate reference potentials (coupled cluster and quantum Monte Carlo results) as well as to state-of-the-art ab initio DFT approaches. We have found that all the ACMs yield correlation potentials that exhibit a correct behavior, quite resembling scaled second-order Görling-Levy (GL2) potentials and including most of the physically meaningful features of the accurate reference data. The behavior and contribution of the strong-interaction limit potentials have also been investigated and discussed.
Collapse
Affiliation(s)
- Eduardo Fabiano
- Institute for Microelectronics and Microsystems (CNR-IMM), Via Monteroni , Campus Unisalento, 73100 Lecce , Italy
- Center for Biomolecular Nanotechnologies @UNILE , Istituto Italiano di Tecnologia , Via Barsanti , I-73010 Arnesano , Italy
| | - Szymon Śmiga
- Institute of Physics, Faculty of Physics, Astronomy and Informatics , Nicolaus Copernicus University , Grudziadzka 5 , 87-100 Torun , Poland
| | - Sara Giarrusso
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling, Faculty of Science , Vrije Universiteit , De Boelelaan 1083 , 1081HV Amsterdam , The Netherlands
| | - Kimberly J Daas
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling, Faculty of Science , Vrije Universiteit , De Boelelaan 1083 , 1081HV Amsterdam , The Netherlands
| | - Fabio Della Sala
- Institute for Microelectronics and Microsystems (CNR-IMM), Via Monteroni , Campus Unisalento, 73100 Lecce , Italy
- Center for Biomolecular Nanotechnologies @UNILE , Istituto Italiano di Tecnologia , Via Barsanti , I-73010 Arnesano , Italy
| | - Ireneusz Grabowski
- Institute of Physics, Faculty of Physics, Astronomy and Informatics , Nicolaus Copernicus University , Grudziadzka 5 , 87-100 Torun , Poland
| | - Paola Gori-Giorgi
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling, Faculty of Science , Vrije Universiteit , De Boelelaan 1083 , 1081HV Amsterdam , The Netherlands
| |
Collapse
|
18
|
Gould T, Pittalis S, Toulouse J, Kraisler E, Kronik L. Asymptotic behavior of the Hartree-exchange and correlation potentials in ensemble density functional theory. Phys Chem Chem Phys 2019; 21:19805-19815. [DOI: 10.1039/c9cp03633d] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We show that the Hartree-exchange and correlation potentials of ensemble systems can have unexpected features, including non-vanishing asymptotic constants and non-trivial screening of electrons. These features are demonstrated here on Li, C, and F.
Collapse
Affiliation(s)
- Tim Gould
- Qld Micro- and Nanotechnology Centre
- Griffith University
- Nathan
- Australia
| | | | - Julien Toulouse
- Laboratoire de Chimie Théorique
- Sorbonne Université and CNRS
- F-75005 Paris
- France
| | - Eli Kraisler
- Fritz Haber Center for Molecular Dynamics and Institute of Chemistry
- The Hebrew University of Jerusalem
- 9091401 Jerusalem
- Israel
| | - Leeor Kronik
- Department of Materials and Interfaces
- Weizmann Institute of Science
- Rehovoth 76100
- Israel
| |
Collapse
|
19
|
Kooi DP, Gori-Giorgi P. Local and global interpolations along the adiabatic connection of DFT: a study at different correlation regimes. Theor Chem Acc 2018; 137:166. [PMID: 30464722 PMCID: PMC6223841 DOI: 10.1007/s00214-018-2354-5] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/01/2018] [Accepted: 09/21/2018] [Indexed: 11/16/2022]
Abstract
Interpolating the exchange-correlation energy along the density-fixed adiabatic connection of density functional theory is a promising way to build approximations that are not biased toward the weakly correlated regime. These interpolations can be performed at the global (integrated over all spaces) or at the local level, using energy densities. Many features of the relevant energy densities as well as several different ways to construct these interpolations, including comparisons between global and local variants, are investigated here for the analytically solvable Hooke's atom series, which allows for an exploration of different correlation regimes. We also analyze different ways to define the correlation kinetic energy density, focusing on the peak in the kinetic correlation potential.
Collapse
Affiliation(s)
- Derk P. Kooi
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling, Faculty of Science, Vrije Universiteit, De Boelelaan 1083, 1081 HV Amsterdam, The Netherlands
| | - Paola Gori-Giorgi
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling, Faculty of Science, Vrije Universiteit, De Boelelaan 1083, 1081 HV Amsterdam, The Netherlands
| |
Collapse
|
20
|
Kosugi T, Matsushita YI. One-particle Green's function of interacting two electrons using analytic solutions for a three-body problem: comparison with exact Kohn-Sham system. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2018; 30:435604. [PMID: 30229746 DOI: 10.1088/1361-648x/aae287] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
For a three-electron system with finite-strength interactions confined to a one-dimensional harmonic trap, we solve the Schrödinger equation analytically to obtain the exact solutions, from which we construct explicitly the simultaneous eigenstates of the energy and total spin for the first time. The solutions for the three-electron system allow us to derive analytic expressions for the exact one-particle Green's function (GF) for the corresponding two-electron system. We calculate the GF in frequency domain to examine systematically its behavior depending on the electronic interactions. We also compare the pole structure of non-interacting GF using the exact Kohn-Sham (KS) potential with that of the exact GF to find that the discrepancy of the energy gap between the KS system and the original system is larger for a stronger interaction. We perform numerical examination on the behavior of GFs in real space to demonstrate that the exact and KS GFs can have shapes quite different from each other. Our simple model will help to understand generic characteristics of interacting GFs.
Collapse
Affiliation(s)
- Taichi Kosugi
- Department of Physics, University of Tokyo, Tokyo 113-0033, Japan
| | | |
Collapse
|
21
|
Uniform electron gas limit of an exact expression for the Kohn–Sham exchange-correlation potential. Theor Chem Acc 2018. [DOI: 10.1007/s00214-018-2303-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
|
22
|
Giarrusso S, Vuckovic S, Gori-Giorgi P. Response Potential in the Strong-Interaction Limit of Density Functional Theory: Analysis and Comparison with the Coupling-Constant Average. J Chem Theory Comput 2018; 14:4151-4167. [PMID: 29906106 PMCID: PMC6096453 DOI: 10.1021/acs.jctc.8b00386] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Using the formalism of the conditional amplitude, we study the response part of the exchange-correlation potential in the strong-coupling limit of density functional theory, analyzing its peculiar features and comparing it with the response potential averaged over the coupling constant for small atoms and for the hydrogen molecule. We also use a simple one-dimensional model of a stretched heteronuclear molecule to derive exact properties of the response potential in the strong-coupling limit. The simplicity of the model allows us to unveil relevant features also of the exact Kohn-Sham potential and its different components, namely the appearance of a second peak in the correlation kinetic potential on the side of the most electronegative atom.
Collapse
Affiliation(s)
- Sara Giarrusso
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling, FEW , Vrije Universiteit , De Boelelaan 1083 , 1081HV Amsterdam , The Netherlands
| | - Stefan Vuckovic
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling, FEW , Vrije Universiteit , De Boelelaan 1083 , 1081HV Amsterdam , The Netherlands
| | - Paola Gori-Giorgi
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling, FEW , Vrije Universiteit , De Boelelaan 1083 , 1081HV Amsterdam , The Netherlands
| |
Collapse
|
23
|
Suzuki Y, Lacombe L, Watanabe K, Maitra NT. Exact Time-Dependent Exchange-Correlation Potential in Electron Scattering Processes. PHYSICAL REVIEW LETTERS 2017; 119:263401. [PMID: 29328727 DOI: 10.1103/physrevlett.119.263401] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/29/2017] [Indexed: 06/07/2023]
Abstract
We identify peak and valley structures in the exact exchange-correlation potential of time-dependent density functional theory that are crucial for time-resolved electron scattering in a model one-dimensional system. These structures are completely missed by adiabatic approximations that, consequently, significantly underestimate the scattering probability. A recently proposed nonadiabatic approximation is shown to correctly capture the approach of the electron to the target when the initial Kohn-Sham state is chosen judiciously, and it is more accurate than standard adiabatic functionals but ultimately fails to accurately capture reflection. These results may explain the underestimation of scattering probabilities in some recent studies on molecules and surfaces.
Collapse
Affiliation(s)
- Yasumitsu Suzuki
- Department of Physics, Tokyo University of Science, 1-3 Kagurazaka, Shinjuku-ku, Tokyo 162-8601, Japan
| | - Lionel Lacombe
- Department of Physics and Astronomy, Hunter College and the Graduate Center of the City University of New York, 695 Park Avenue, New York, New York 10065, USA
| | - Kazuyuki Watanabe
- Department of Physics, Tokyo University of Science, 1-3 Kagurazaka, Shinjuku-ku, Tokyo 162-8601, Japan
| | - Neepa T Maitra
- Department of Physics and Astronomy, Hunter College and the Graduate Center of the City University of New York, 695 Park Avenue, New York, New York 10065, USA
| |
Collapse
|
24
|
Hodgson MJP, Kraisler E, Schild A, Gross EKU. How Interatomic Steps in the Exact Kohn-Sham Potential Relate to Derivative Discontinuities of the Energy. J Phys Chem Lett 2017; 8:5974-5980. [PMID: 29179553 DOI: 10.1021/acs.jpclett.7b02615] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
Accurate density functional calculations hinge on reliable approximations to the unknown exchange-correlation (xc) potential. The most popular approximations usually lack features of the exact xc potential that are important for an accurate prediction of the fundamental gap and the distribution of charge in complex systems. Two principal features in this regard are the spatially uniform shift in the potential, as the number of electrons infinitesimally surpasses an integer, and the spatial steps that form, for example, between the atoms of stretched molecules. Although both aforementioned concepts are well known, the exact relationship between them remained unclear. Here we establish this relationship via an analytical derivation. We support our result by numerically solving the many-electron Schrödinger equation to extract the exact Kohn-Sham potential and directly observe its features. Spatial steps in the exact xc potential of a full configuration-interaction (FCI) calculation of a molecule are presented in three dimensions.
Collapse
Affiliation(s)
- M J P Hodgson
- Max-Planck-Institut für Mikrostrukturphysik , Weinberg 2, D-06120 Halle, Germany
| | - Eli Kraisler
- Max-Planck-Institut für Mikrostrukturphysik , Weinberg 2, D-06120 Halle, Germany
| | - Axel Schild
- Laboratorium für Physikalische Chemie, ETH Zürich , 8093 Zürich, Switzerland
| | - E K U Gross
- Max-Planck-Institut für Mikrostrukturphysik , Weinberg 2, D-06120 Halle, Germany
- Fritz Haber Center for Molecular Dynamics, Institute of Chemistry, The Hebrew University of Jerusalem , Jerusalem 91904, Israel
| |
Collapse
|
25
|
Ryabinkin IG, Ospadov E, Staroverov VN. Exact exchange-correlation potentials of singlet two-electron systems. J Chem Phys 2017; 147:164117. [DOI: 10.1063/1.5003825] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Ilya G. Ryabinkin
- Department of Physical and Environmental Sciences, University of Toronto Scarborough, Toronto, Ontario M1C 1A4,
Canada
| | - Egor Ospadov
- Department of Chemistry, The University of Western Ontario, London, Ontario N6A 5B7, Canada
| | - Viktor N. Staroverov
- Department of Chemistry, The University of Western Ontario, London, Ontario N6A 5B7, Canada
| |
Collapse
|
26
|
Maitra NT. Charge transfer in time-dependent density functional theory. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2017; 29:423001. [PMID: 28766507 DOI: 10.1088/1361-648x/aa836e] [Citation(s) in RCA: 58] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
Charge transfer plays a crucial role in many processes of interest in physics, chemistry, and bio-chemistry. In many applications the size of the systems involved calls for time-dependent density functional theory (TDDFT) to be used in their computational modeling, due to its unprecedented balance between accuracy and efficiency. However, although exact in principle, in practise approximations must be made for the exchange-correlation functional in this theory, and the standard functional approximations perform poorly for excitations which have a long-range charge-transfer component. Intense progress has been made in developing more sophisticated functionals for this problem, which we review. We point out an essential difference between the properties of the exchange-correlation kernel needed for an accurate description of charge-transfer between open-shell fragments and between closed-shell fragments. We then turn to charge-transfer dynamics, which, in contrast to the excitation problem, is a highly non-equilibrium, non-perturbative, process involving a transfer of one full electron in space. This turns out to be a much more challenging problem for TDDFT functionals. We describe dynamical step and peak features in the exact functional evolving over time, that are missing in the functionals currently used. The latter underestimate the amount of charge transferred and manifest a spurious shift in the charge transfer resonance position. We discuss some explicit examples.
Collapse
Affiliation(s)
- Neepa T Maitra
- Department of Physics and Astronomy, Hunter College and the Physics Program at the Graduate Center of the City University of New York, 695 Park Avenue, NY 10065, United States of America
| |
Collapse
|
27
|
Kosugi T, Matsushita YI. Quantum Singwi-Tosi-Land-Sjölander approach for interacting inhomogeneous systems under electromagnetic fields: Comparison with exact results. J Chem Phys 2017; 147:114105. [PMID: 28938817 DOI: 10.1063/1.4994720] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
For inhomogeneous interacting electronic systems under a time-dependent electromagnetic perturbation, we derive the linear equation for response functions in a quantum mechanical manner. It is a natural extension of the original semi-classical Singwi-Tosi-Land-Sjölander (STLS) approach for an electron gas. The factorization ansatz for the two-particle distribution is an indispensable ingredient in the STLS approaches for the determination of the response function and the pair correlation function. In this study, we choose an analytically solvable interacting two-electron system as the target for which we examine the validity of the approximation. It is demonstrated that the STLS response function reproduces well the exact one for low-energy excitations. The interaction energy contributed from the STLS response function is also discussed.
Collapse
Affiliation(s)
- Taichi Kosugi
- Department of Applied Physics, The University of Tokyo, Tokyo 113-8656, Japan
| | | |
Collapse
|
28
|
Baerends EJ. From the Kohn–Sham band gap to the fundamental gap in solids. An integer electron approach. Phys Chem Chem Phys 2017; 19:15639-15656. [DOI: 10.1039/c7cp02123b] [Citation(s) in RCA: 40] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
The upshift Δ of the level at the bottom of the conduction band (the LUMO) from the neutral N-electron crystal to the negative N + 1 system, and therefore the fundamental gap εLUMO(N + 1) − εHOMO(N) = I − A, can be calculated simply and cheaply from the response part of vxc.
Collapse
Affiliation(s)
- E. J. Baerends
- Section Theoretical Chemistry
- Vrije Universiteit
- Amsterdam
- The Netherlands
| |
Collapse
|
29
|
Vuckovic S, Irons TJP, Wagner LO, Teale AM, Gori-Giorgi P. Interpolated energy densities, correlation indicators and lower bounds from approximations to the strong coupling limit of DFT. Phys Chem Chem Phys 2017; 19:6169-6183. [DOI: 10.1039/c6cp08704c] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
|
30
|
Komsa DN, Staroverov VN. Elimination of Spurious Fractional Charges in Dissociating Molecules by Correcting the Shape of Approximate Kohn–Sham Potentials. J Chem Theory Comput 2016; 12:5361-5366. [DOI: 10.1021/acs.jctc.6b00798] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Darya N. Komsa
- Department of Chemistry, The University of Western Ontario, London, Ontario N6A 5B7, Canada
| | - Viktor N. Staroverov
- Department of Chemistry, The University of Western Ontario, London, Ontario N6A 5B7, Canada
| |
Collapse
|
31
|
Genova A, Ceresoli D, Pavanello M. Avoiding fractional electrons in subsystem DFT based ab-initio molecular dynamics yields accurate models for liquid water and solvated OH radical. J Chem Phys 2016; 144:234105. [DOI: 10.1063/1.4953363] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Alessandro Genova
- Department of Chemistry, Rutgers University, Newark, New Jersey 07102, USA
| | - Davide Ceresoli
- CNR-ISTM: Institute of Molecular Sciences and Technologies, Milano, Italy
| | - Michele Pavanello
- Department of Chemistry, Rutgers University, Newark, New Jersey 07102, USA
| |
Collapse
|
32
|
Nafziger J, Wasserman A. Fragment-based treatment of delocalization and static correlation errors in density-functional theory. J Chem Phys 2016; 143:234105. [PMID: 26696044 DOI: 10.1063/1.4937771] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Abstract
One of the most important open challenges in modern Kohn-Sham (KS) density-functional theory (DFT) is the correct treatment of systems involving fractional electron charges and spins. Approximate exchange-correlation functionals struggle with such systems, leading to pervasive delocalization and static correlation errors. We demonstrate how these errors, which plague density-functional calculations of bond-stretching processes, can be avoided by employing the alternative framework of partition density-functional theory (PDFT) even using the local density approximation for the fragments. Our method is illustrated with explicit calculations on simple systems exhibiting delocalization and static-correlation errors, stretched H2 (+), H2, He2 (+), Li2 (+), and Li2. In all these cases, our method leads to greatly improved dissociation-energy curves. The effective KS potential corresponding to our self-consistent solutions displays key features around the bond midpoint; these are known to be present in the exact KS potential, but are absent from most approximate KS potentials and are essential for the correct description of electron dynamics.
Collapse
Affiliation(s)
- Jonathan Nafziger
- Department of Chemistry, Purdue University, 560 Oval Dr., West Lafayette, Indiana 47907, USA
| | - Adam Wasserman
- Department of Chemistry, Purdue University, 560 Oval Dr., West Lafayette, Indiana 47907, USA
| |
Collapse
|
33
|
Kohut SV, Polgar AM, Staroverov VN. Origin of the step structure of molecular exchange–correlation potentials. Phys Chem Chem Phys 2016; 18:20938-44. [DOI: 10.1039/c6cp00878j] [Citation(s) in RCA: 29] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The step structure of exact exchange–correlation potentials is linked to the properties of the average local electron energy (ALEE).
Collapse
|
34
|
Cuevas-Saavedra R, Ayers PW, Staroverov VN. Kohn–Sham exchange-correlation potentials from second-order reduced density matrices. J Chem Phys 2015; 143:244116. [DOI: 10.1063/1.4937943] [Citation(s) in RCA: 36] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
| | - Paul W. Ayers
- Department of Chemistry and Chemical Biology, McMaster University, Hamilton, Ontario L8S 4M1, Canada
| | - Viktor N. Staroverov
- Department of Chemistry, The University of Western Ontario, London, Ontario N6A 5B7, Canada
| |
Collapse
|
35
|
Ryabinkin IG, Kohut SV, Staroverov VN. Reduction of Electronic Wave Functions to Kohn-Sham Effective Potentials. PHYSICAL REVIEW LETTERS 2015; 115:083001. [PMID: 26340185 DOI: 10.1103/physrevlett.115.083001] [Citation(s) in RCA: 43] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/09/2015] [Indexed: 06/05/2023]
Abstract
A method for calculating the Kohn-Sham exchange-correlation potential v(XC)(r) from a given electronic wave function is devised and implemented. It requires on input one- and two-electron reduced density matrices and involves construction of the generalized Fock matrix. The method is free from numerical limitations and basis-set artifacts of conventional schemes for constructing v(XC)(r) in which the potential is recovered from a given electron density, and is simpler than various many-body techniques. The chief significance of this development is that it allows one to directly probe the functional derivative of the true exchange-correlation energy functional and to rigorously test and improve various density-functional approximations.
Collapse
Affiliation(s)
- Ilya G Ryabinkin
- Department of Physical and Environmental Sciences, University of Toronto Scarborough, Toronto, Ontario M1C 1A4, Canada
| | - Sviataslau V Kohut
- Department of Chemistry, The University of Western Ontario, London, Ontario N6A 5B7, Canada
| | - Viktor N Staroverov
- Department of Chemistry, The University of Western Ontario, London, Ontario N6A 5B7, Canada
| |
Collapse
|
36
|
Luo K, Fuks JI, Sandoval ED, Elliott P, Maitra NT. Kinetic and interaction components of the exact time-dependent correlation potential. J Chem Phys 2015; 140:18A515. [PMID: 24832323 DOI: 10.1063/1.4867002] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023] Open
Abstract
The exact exchange-correlation (xc) potential of time-dependent density functional theory has been shown to have striking features. For example, step and peak features are generically found when the system is far from its ground-state, and these depend nonlocally on the density in space and time. We analyze the xc potential by decomposing it into kinetic and interaction components and comparing each with their exact-adiabatic counterparts, for a range of dynamical situations in model one-dimensional two-electron systems. We find that often, but not always, the kinetic contribution is largely responsible for these features that are missed by the adiabatic approximation. The adiabatic approximation often makes a smaller error for the interaction component, which we write in two parts, one being the Coulomb potential due to the time-dependent xc hole. Non-adiabatic features of the kinetic component were also larger than those of the interaction component in cases that we studied when there is negligible step structure. In ground-state situations, step and peak structures arise in cases of static correlation, when more than one determinant is essential to describe the interacting state. We investigate the time-dependent natural orbital occupation numbers and find the corresponding relation between these and the dynamical step is more complex than for the ground-state case.
Collapse
Affiliation(s)
- Kai Luo
- Department of Physics and Astronomy, Hunter College and the Graduate Center of the City University of New York, 695 Park Avenue, New York, New York 10065, USA
| | - Johanna I Fuks
- Department of Physics and Astronomy, Hunter College and the Graduate Center of the City University of New York, 695 Park Avenue, New York, New York 10065, USA
| | - Ernesto D Sandoval
- Department of Physics and Astronomy, Hunter College and the Graduate Center of the City University of New York, 695 Park Avenue, New York, New York 10065, USA
| | - Peter Elliott
- Max-Planck-Institut für Mikrostrukturphysik, Weinberg 2, 06120 Halle (Saale), Germany
| | - Neepa T Maitra
- Department of Physics and Astronomy, Hunter College and the Graduate Center of the City University of New York, 695 Park Avenue, New York, New York 10065, USA
| |
Collapse
|
37
|
Vieira D. Strong Correlations in Density-Functional Theory: A Model of Spin-Charge and Spin-Orbital Separations. J Chem Theory Comput 2014; 10:3641-6. [PMID: 26588509 DOI: 10.1021/ct400223d] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
It is known that the separation of electrons into spinons and chargons, the spin-charge separation, plays a decisive role when describing one-dimensional (1D) strongly correlated systems [ Phys. Rev. B 2012 , 86 , 075132 ]. In this paper, within the density-functional theory (DFT) formalism, we extend the investigation by considering a model for the third electron fractionalization: the separation into spinons, chargons and orbitons, the last associated with the electronic orbital degree of freedom. Specifically, we deal with two exact constraints of exchange-correlation (XC) density-functionals: (i) the constancy of the highest occupied (HO) Kohn-Sham (KS) eigenvalues upon fractional electron numbers and (ii) their discontinuities at integers. By means of 1D discrete Hubbard chains and 1D H2 molecules in the continuum, we find that spin-charge separation yields almost constant HO KS eigenvalues, whereas the spin-orbital counterpart can be decisive when describing derivative discontinuities of XC potentials at strong correlations.
Collapse
Affiliation(s)
- Daniel Vieira
- Departamento de Física, Universidade do Estado de Santa Catarina , Joinville, 89219-710 Santa Catarina, Brazil
| |
Collapse
|
38
|
Mori-Sánchez P, Cohen AJ. The derivative discontinuity of the exchange–correlation functional. Phys Chem Chem Phys 2014; 16:14378-87. [DOI: 10.1039/c4cp01170h] [Citation(s) in RCA: 57] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
Abstract
Manifestations of the derivative discontinuity of the energy in density functional theory are demonstrated in simple systems in chemistry and physics.
Collapse
Affiliation(s)
- Paula Mori-Sánchez
- Departamento de Química and Instituto de Física de la Materia Condensada (IFIMAC)
- Universidad Autónoma de Madrid
- Madrid, Spain
| | - Aron J. Cohen
- Department of Chemistry
- University of Cambridge
- Cambridge, UK
| |
Collapse
|
39
|
Fuks JI, Maitra NT. Challenging adiabatic time-dependent density functional theory with a Hubbard dimer: the case of time-resolved long-range charge transfer. Phys Chem Chem Phys 2014; 16:14504-13. [DOI: 10.1039/c4cp00118d] [Citation(s) in RCA: 41] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
Adiabatic TDDFT fails for time-resolved charge-transfer in a simple model, even when linear response spectra is accurate.
Collapse
Affiliation(s)
- Johanna I Fuks
- Department of Physics and Astronomy, Hunter College and the City University of New York, 695 Park Avenue, New York, New York 10065, USA.
| | | |
Collapse
|
40
|
Mosquera MA, Wasserman A. Non-analytic Spin-Density Functionals. Top Curr Chem (Cham) 2014. [DOI: 10.1007/128_2014_619] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
41
|
Armiento R, Kümmel S. Orbital localization, charge transfer, and band gaps in semilocal density-functional theory. PHYSICAL REVIEW LETTERS 2013; 111:036402. [PMID: 23909345 DOI: 10.1103/physrevlett.111.036402] [Citation(s) in RCA: 36] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/12/2012] [Indexed: 06/02/2023]
Abstract
We derive an exchange energy functional of generalized gradient form with a corresponding potential that changes discontinuously at integer particle numbers. The functional is semilocal, yet incorporates key features that are connected to the derivative discontinuity of Kohn-Sham density-functional theory. We validate our construction for several paradigm systems and explain how it addresses central well-known deficiencies of antecedent semilocal methods, i.e., the description of charge transfer, properly localized orbitals, and band gaps. We find, e.g., an improved shell structure for atoms, eigenvalues that more closely correspond to ionization energies, and an improved description of band structure where localized states are lowered in energy.
Collapse
Affiliation(s)
- R Armiento
- Department of Physics, Chemistry and Biology (IFM), Linköping University, SE-58183 Linköping, Sweden
| | | |
Collapse
|
42
|
Fuks JI, Elliott P, Rubio A, Maitra NT. Dynamics of Charge-Transfer Processes with Time-Dependent Density Functional Theory. J Phys Chem Lett 2013; 4:735-739. [PMID: 26281927 DOI: 10.1021/jz302099f] [Citation(s) in RCA: 45] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
We show that whenever an electron transfers between closed-shell molecular fragments, the exact correlation potential of time-dependent density functional theory develops a step and peak structure in the bonding region. This structure has a density dependence that is nonlocal both in space and in time that even the exact adiabatic ground-state exchange-correlation functional fails to capture it. For charge-transfer between open-shell fragments, an initial step and peak vanish as the charge-transfer state is reached. The inability of usual approximations to develop these structures leads to inaccurate charge-transfer dynamics. This is illustrated by the complete lack of Rabi oscillations in the dipole moment under conditions of resonant charge transfer for an exactly solvable model system. The results transcend the model and are applicable to more realistic molecular complexes.
Collapse
Affiliation(s)
- J I Fuks
- †Nano-Bio Spectroscopy Group, Departamento Fı́sica de Materiales, Universidad del Paı́s Vasco, Centro de Fı́sica de Materiales CSIC-UPV/EHU-MPC and DIPC, Av. Tolosa 72, E-20018 San Sebastián, Spain
| | - P Elliott
- ‡Max-Planck-Institut für Mikrostrukturphysik, Weinberg 2, 06120 Halle (Saale), Germany
| | - A Rubio
- †Nano-Bio Spectroscopy Group, Departamento Fı́sica de Materiales, Universidad del Paı́s Vasco, Centro de Fı́sica de Materiales CSIC-UPV/EHU-MPC and DIPC, Av. Tolosa 72, E-20018 San Sebastián, Spain
- §Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, D-14195 Berlin, Germany
| | - N T Maitra
- ⊥Department of Physics and Astronomy, Hunter College and the Graduate Center of the City University of New York, 695 Park Avenue, New York, New York 10065, United States
| |
Collapse
|
43
|
Malet F, Gori-Giorgi P. Strong correlation in Kohn-Sham density functional theory. PHYSICAL REVIEW LETTERS 2012; 109:246402. [PMID: 23368350 DOI: 10.1103/physrevlett.109.246402] [Citation(s) in RCA: 49] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/11/2012] [Indexed: 06/01/2023]
Abstract
We use the exact strong-interaction limit of the Hohenberg-Kohn energy density functional to approximate the exchange-correlation energy of the restricted Kohn-Sham scheme. Our approximation corresponds to a highly nonlocal density functional whose functional derivative can be easily constructed, thus transforming exactly, in a physically transparent way, an important part of the electron-electron interaction into an effective local one-body potential. We test our approach on quasi-one-dimensional systems, showing that it captures essential features of strong correlation that restricted Kohn-Sham calculations using the currently available approximations cannot describe.
Collapse
Affiliation(s)
- Francesc Malet
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling, FEW, Vrije Universiteit, De Boelelaan 1083, 1081HV Amsterdam, The Netherlands
| | | |
Collapse
|
44
|
Hellgren M, Gross EKU. Effect of discontinuities in Kohn-Sham-based chemical reactivity theory. J Chem Phys 2012; 136:114102. [DOI: 10.1063/1.3694103] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
45
|
Hellgren M, Rohr DR, Gross EKU. Correlation potentials for molecular bond dissociation within the self-consistent random phase approximation. J Chem Phys 2012; 136:034106. [DOI: 10.1063/1.3676174] [Citation(s) in RCA: 67] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
46
|
Adiabatic approximation within time-dependent density functional theory using inversion of the ground-state spin-density Kohn–Sham formalism. Chem Phys 2011. [DOI: 10.1016/j.chemphys.2011.02.009] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
47
|
|
48
|
Hesselmann A, Görling A. Correct description of the bond dissociation limit without breaking spin symmetry by a random-phase-approximation correlation functional. PHYSICAL REVIEW LETTERS 2011; 106:093001. [PMID: 21405619 DOI: 10.1103/physrevlett.106.093001] [Citation(s) in RCA: 60] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/24/2010] [Indexed: 05/30/2023]
Abstract
A correlation functional that is termed exact-exchange random phase approximation (EXX-RPA) functional and is obtained with the exact frequency-dependent exchange kernel via the fluctuation-dissipation theorem is shown to correctly describe electron pair bonds in the dissociation limit without the need to resort to symmetry breaking in spin space. Because the functional also yields more accurate electronic energies for molecules in their equilibrium geometry than standard correlation functionals, it combines accuracy at equilibrium bond distances and in dissociation processes with a correct description of spin, something all commonly employed correlation functionals fail to do. The reason why the EXX-RPA correlation functional yields distinctively and qualitatively better results than RPA approaches based on Hartree-Fock and time-dependent Hartree-Fock is explained.
Collapse
Affiliation(s)
- Andreas Hesselmann
- Lehrstuhl für Theoretische Chemie, Universität Erlangen-Nürnberg, Germany
| | | |
Collapse
|
49
|
Helbig N, Tokatly IV, Rubio A. Exact Kohn–Sham potential of strongly correlated finite systems. J Chem Phys 2009; 131:224105. [DOI: 10.1063/1.3271392] [Citation(s) in RCA: 74] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
50
|
Rajam A, Hessler P, Gaun C, Maitra NT. Phase-space explorations in time-dependent density functional theory. ACTA ACUST UNITED AC 2009. [DOI: 10.1016/j.theochem.2009.04.008] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|