1
|
Constantin LA, Jana S, Śmiga S, Della Sala F. Adiabatic connection interaction strength interpolation method made accurate for the uniform electron gas. J Chem Phys 2023; 159:244111. [PMID: 38149733 DOI: 10.1063/5.0178800] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/28/2023] [Accepted: 11/27/2023] [Indexed: 12/28/2023] Open
Abstract
The adiabatic connection interaction strength interpolation (ISI)-like method provides a high-level expression for the correlation energy, being, in principle, exact not only in the weak-interaction limit, where it recovers the second-order Görling-Levy perturbation term, but also in the strong-interaction limit that is described by the strictly correlated electron approach. In this work, we construct a genISI functional made accurate for the uniform electron gas, a solid-state physics paradigm that is a very difficult test for ISI-like correlation functionals. We assess the genISI functional for various jellium spheres with the number of electrons Z ≤ 912 and for the non-relativistic noble atoms with Z ≤ 290. For the jellium clusters, the genISI is remarkably accurate, while for the noble atoms, it shows a good performance, similar to other ISI-like methods. Then, the genISI functional can open the path using the ISI-like method in solid-state calculations.
Collapse
Affiliation(s)
- Lucian A Constantin
- Institute for Microelectronics and Microsystems (CNR-IMM), Via Monteroni, Campus Unisalento, 73100 Lecce, Italy
| | - Subrata Jana
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, Rehovoth 76100, Israel
| | - Szymon Śmiga
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University in Toruń, ul. Grudziądzka 5, 87-100 Toruń, Poland
| | - Fabio Della Sala
- Institute for Microelectronics and Microsystems (CNR-IMM), Via Monteroni, Campus Unisalento, 73100 Lecce, Italy
- Center for Biomolecular Nanotechnologies, Istituto Italiano di Tecnologia, Via Barsanti 14, 73010 Arnesano (LE), Italy
| |
Collapse
|
2
|
Jana S, Śmiga S, Constantin LA, Samal P. Semilocal Meta-GGA Exchange-Correlation Approximation from Adiabatic Connection Formalism: Extent and Limitations. J Phys Chem A 2023; 127:8685-8697. [PMID: 37811903 PMCID: PMC10591512 DOI: 10.1021/acs.jpca.3c03976] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/12/2023] [Revised: 08/24/2023] [Indexed: 10/10/2023]
Abstract
The incorporation of a strong-interaction regime within the approximate semilocal exchange-correlation functionals still remains a very challenging task for density functional theory. One of the promising attempts in this direction is the recently proposed adiabatic connection semilocal correlation (ACSC) approach [Constantin, L. A.; Phys. Rev. B 2019, 99, 085117] allowing one to construct the correlation energy functionals by interpolation of the high and low-density limits for the given semilocal approximation. The current study extends the ACSC method to the meta-generalized gradient approximations (meta-GGA) level of theory, providing some new insights in this context. As an example, we construct the correlation energy functional on the basis of the high- and low-density limits of the Tao-Perdew-Staroverov-Scuseria (TPSS) functional. Arose in this way, the TPSS-ACSC functional is one-electron self-interaction free and accurate for the strictly correlated and quasi-two-dimensional regimes. Based on simple examples, we show the advantages and disadvantages of ACSC semilocal functionals and provide some new guidelines for future developments in this context.
Collapse
Affiliation(s)
- Subrata Jana
- Department
of Chemistry & Biochemistry, The Ohio
State University, Columbus, Ohio 43210, United States
| | - Szymon Śmiga
- Institute
of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University in Toruń, ul. Grudzikadzka 5, 87-100 Toruń, Poland
| | - Lucian A. Constantin
- Istituto
di Nanoscienze, Consiglio Nazionale delle
Ricerche CNR-NANO, 41125 Modena, Italy
| | - Prasanjit Samal
- School
of Physical Sciences, National Institute of Science Education and
Research, HBNI, Bhubaneswar 752050, India
| |
Collapse
|
3
|
Cioslowski J, Englert BG, Trappe MI, Hue JH. Contactium: A strongly correlated model system. J Chem Phys 2023; 158:2890212. [PMID: 37171198 DOI: 10.1063/5.0150523] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2023] [Accepted: 04/21/2023] [Indexed: 05/13/2023] Open
Abstract
At the limit of an infinite confinement strength ω, the ground state of a system that comprises two fermions or bosons in harmonic confinement interacting through the Fermi-Huang pseudopotential remains strongly correlated. A detailed analysis of the one-particle description of this "contactium" reveals several peculiarities that are not encountered in conventional model systems (such as the two-electron harmonium atom, ballium, and spherium) involving Coulombic interparticle interactions. First of all, none of the natural orbitals (NOs) {ψn(ω;r)} of the contactium is unoccupied, which implies nonzero collective occupancies for all the angular momenta. Second, the NOs and their non-ascendingly ordered occupation numbers {νn} turn out to be related to the eigenfunctions and eigenvalues of a zero-energy Schrödinger equation with an attractive Gaussian potential. This observation enables the derivation of their properties, such as the n-4/3 asymptotic decay of νn at the n→∞ limit (which differs from that of n-8/3 in the Coulombic systems), the independence of the confinement energy vn=⟨ψn(ω;r)|12ω2r2|ψn(ω;r)⟩ of n, and the n-2/3 asymptotic decay of the respective contribution νntn to the kinetic energy. Upon suitable scaling, the weakly occupied NOs of the contactium turn out to be virtually identical to those of the two-electron harmonium atom at the ω → ∞ limit, despite the entirely different interparticle interactions in these systems.
Collapse
Affiliation(s)
- Jerzy Cioslowski
- Institute of Physics, University of Szczecin, Wielkopolska 15, 70-451 Szczecin, Poland
| | - Berthold-Georg Englert
- Centre for Quantum Technologies, National University of Singapore, 3 Science Drive 2, Singapore 117543, Singapore
- Department of Physics, National University of Singapore, 2 Science Drive 3, Singapore 117542, Singapore
- MajuLab, CNRS-UCA-SU-NUS-NTU International, Joint Research Unit, Singapore
| | - Martin-Isbjörn Trappe
- Centre for Quantum Technologies, National University of Singapore, 3 Science Drive 2, Singapore 117543, Singapore
| | - Jun Hao Hue
- Centre for Quantum Technologies, National University of Singapore, 3 Science Drive 2, Singapore 117543, Singapore
| |
Collapse
|
4
|
Śmiga S, Della Sala F, Gori-Giorgi P, Fabiano E. Self-Consistent Implementation of Kohn-Sham Adiabatic Connection Models with Improved Treatment of the Strong-Interaction Limit. J Chem Theory Comput 2022; 18:5936-5947. [PMID: 36094908 PMCID: PMC9558377 DOI: 10.1021/acs.jctc.2c00352] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
Adiabatic connection
models (ACMs), which interpolate between the
limits of weak and strong interaction, are powerful tools to build
accurate exchange–correlation functionals. If the exact weak-interaction
expansion from the second-order perturbation theory is included, a
self-consistent implementation of these functionals is challenging
and still absent in the literature. In this work, we fill this gap
by presenting a fully self-consistent-field (SCF) implementation of
some popular ACM functionals. While using second-order perturbation
theory at weak interactions, we have also introduced new generalized
gradient approximations (GGAs), beyond the usual point-charge-plus-continuum
model, for the first two leading terms at strong interactions, which
are crucial to ensure robustness and reliability. We then assess the
SCF–ACM functionals for molecular systems and for prototypical
strong-correlation problems. We find that they perform well for both
the total energy and the electronic density and that the impact of
SCF orbitals is directly connected to the accuracy of the ACM functional
form. For the H2 dissociation, the SCF–ACM functionals
yield significant improvements with respect to standard functionals
also thanks to the use of the new GGAs for the strong-coupling functionals.
Collapse
Affiliation(s)
- Szymon Śmiga
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University in Toruń, ul. Grudziądzka 5, 87-100 Toruń, Poland
| | - Fabio Della Sala
- Institute for Microelectronics and Microsystems (CNR-IMM), Campus Unisalento, Lecce, Via Monteroni 73100, Italy.,Center for Biomolecular Nanotechnologies, Istituto Italiano di Tecnologia, Via Barsanti 14, Arnesano, Lecce 73010, Italy
| | - Paola Gori-Giorgi
- Department of Chemistry & Pharmaceutical Sciences and Amsterdam Institute of Molecular and Life Sciences (AIMMS), Faculty of Science, Vrije Universiteit, De Boelelaan 1083, 1081HV Amsterdam, The Netherlands
| | - Eduardo Fabiano
- Institute for Microelectronics and Microsystems (CNR-IMM), Campus Unisalento, Lecce, Via Monteroni 73100, Italy.,Center for Biomolecular Nanotechnologies, Istituto Italiano di Tecnologia, Via Barsanti 14, Arnesano, Lecce 73010, Italy
| |
Collapse
|
5
|
Everhart LM, Derteano JA, Bates JE. Tension between predicting accurate ground state correlation energies and excitation energies from adiabatic approximations in TDDFT. J Chem Phys 2022; 156:084116. [DOI: 10.1063/5.0080382] [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
The connection between the adiabatic excitation energy of time-dependent density functional theory and the ground state correlation energy from the adiabatic connection fluctuation–dissipation theorem (ACFDT) is explored in the limiting case of one excited state. An exact expression is derived for any adiabatic Hartree-exchange–correlation kernel that connects the excitation energy and the potential contribution to correlation. The resulting formula is applied to the asymmetric Hubbard dimer, a system where this limit is exact. Results from a hierarchy of approximations to the kernel, including the random phase approximation (RPA) with and without exchange and the adiabatically exact (AE) approximation, are compared to the exact ones. At full coupling, the numerical results indicate a tension between predicting an accurate excitation energy and an accurate potential contribution to correlation. The AE approximation is capable of making accurate predictions of both quantities, but only in parts of the parameter space that classify as weakly correlated, while RPA tends to be unable to accurately predict these properties simultaneously anywhere. For a strongly correlated dimer, the AE approximation greatly overestimates the excitation energy yet continues to yield an accurate ground state correlation energy due to its accurate prediction of the adiabatic connection integrand. If similar trends hold for real systems, the development of correlation kernels will be important for applications of the ACFDT in systems with large potential contributions to correlation.
Collapse
Affiliation(s)
- Lucas M. Everhart
- Department of Chemistry and Fermentation Sciences, Appalachian State University, Boone, North Carolina 28607, USA
| | - Julio A. Derteano
- Department of Chemistry and Fermentation Sciences, Appalachian State University, Boone, North Carolina 28607, USA
| | - Jefferson E. Bates
- Department of Chemistry and Fermentation Sciences, Appalachian State University, Boone, North Carolina 28607, USA
| |
Collapse
|
6
|
Śmiga S, Constantin LA. Modified Interaction-Strength Interpolation Method as an Important Step toward Self-Consistent Calculations. J Chem Theory Comput 2020; 16:4983-4992. [PMID: 32559078 PMCID: PMC7588043 DOI: 10.1021/acs.jctc.0c00328] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
The modified point charge plus continuum (mPC) model [ConstantinL. A.; Phys. Rev. B2019, 99, 085117] solves
the important failures of the original counterpart, namely, the divergences
when the reduced gradient of the density is large, such as in the
tail of the density and in quasi-dimensional density regimes. The
mPC allows us to define a modified interaction-strength interpolation
(mISI) method inheriting these good features, which are important
steps toward the full self-consistent treatment. Here, we provide
an assessment of mISI for molecular systems (i.e.,
considering thermochemistry properties, correlation energies, vertical
ionization potentials, and several noncovalent interactions), harmonium
atoms, and functional derivatives in the strong-interaction limit.
For all our tests, mISI provides a systematic improvement over the
original ISI method. Semilocal approximations of the second-order
Görling–Levy (GL2) perturbation theory are also considered
in the mISI method, showing considerable worsening of the results.
Possible further development of mISI is briefly discussed.
Collapse
Affiliation(s)
- Szymon Śmiga
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University, 87-100 Toruń, Poland
| | - Lucian A Constantin
- Consiglio Nazionale delle Ricerche CNR-NANO, Istituto di Nanoscienze, 41125 Modena, Italy
| |
Collapse
|
7
|
Vuckovic S, Fabiano E, Gori-Giorgi P, Burke K. MAP: An MP2 Accuracy Predictor for Weak Interactions from Adiabatic Connection Theory. J Chem Theory Comput 2020; 16:4141-4149. [PMID: 32379454 DOI: 10.1021/acs.jctc.0c00049] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Second-order Møller-Plesset perturbation theory (MP2) approximates the exact Hartree-Fock (HF) adiabatic connection (AC) curve by a straight line. Thus, by using the deviation of the exact curve from the linear behavior, we construct an indicator for the accuracy of MP2. We then use an interpolation along the HF AC to transform the exact form of our indicator into a highly practical MP2 accuracy predictor (MAP) that comes at a negligible additional computational cost. We show that this indicator is already applicable to systems that dissociate into fragments with a nondegenerate ground state, and we illustrate its usefulness by applying it to the S22 and S66 datasets.
Collapse
Affiliation(s)
- Stefan Vuckovic
- Department of Chemistry, University of California, Irvine, California 92697, United States
| | - Eduardo Fabiano
- Institute for Microelectronics and Microsystems (CNR-IMM), Via Monteroni, Campus Unisalento, Lecce 73100, Italy
| | - Paola Gori-Giorgi
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling, FEW, Vrije Universiteit, De Boelelaan 1083, Amsterdam 1081HV, The Netherlands
| | - Kieron Burke
- Department of Chemistry, University of California, Irvine, California 92697, United States
| |
Collapse
|
8
|
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
|