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
|
Vuckovic S, Bahmann H. Nonlocal Functionals Inspired by the Strongly Interacting Limit of DFT: Exact Constraints and Implementation. J Chem Theory Comput 2023; 19:6172-6184. [PMID: 37611177 DOI: 10.1021/acs.jctc.3c00437] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/25/2023]
Abstract
Capturing strong correlation effects remains a key challenge for the development of improved exchange-correlation (XC) functionals in density functional theory. The recently proposed multiple radii functional (MRF) [J. Phys. Chem. Lett. 2017, 8, 2799; J. Chem. Theory Comput. 2019, 15, 3580] was designed to capture strong correlation effects seamlessly, as its mathematical structure draws from that of the exact XC functional in the limit of infinite correlations. The MRF functional provides a framework for building approximations along the density-fixed adiabatic connection, delivers accurate XC energy densities in the standard DFT gauge (same as that of the exact exchange energy density), and is free of one-electron self-interaction errors. To facilitate the development of XC functionals based on the MRF, we examine the behavior of the MRF functional when applied to uniform and scaled densities and consider how it can be made exact for the uniform electron gas. These theoretical insights are then used to build improved forms for the fluctuation function, an object that defines XC energy densities within the MRF framework. We also show how the MRF fluctuation function for physical correlation can be easily readjusted to accurately capture the XC functional in the limit of infinite correlations, demonstrating the versatility of MRF for building approximations for different correlation regimes. We describe the implementation of MRF using densities expanded on Gaussian basis sets, which improves the efficiency of previous grid-based MRF implementations.
Collapse
Affiliation(s)
- Stefan Vuckovic
- Department of Chemistry, University of Fribourg, 1700 Fribourg, Switzerland
| | - Hilke Bahmann
- Physical and Theoretical Chemistry, University of Wuppertal, Gaußstr. 20, 42119 Wuppertal, Germany
| |
Collapse
|
3
|
Giarrusso S, Pribram-Jones A. Møller-Plesset and Density-Fixed Adiabatic Connections for a Model Diatomic System at Different Correlation Regimes. J Chem Theory Comput 2023; 19:5835-5850. [PMID: 37642270 DOI: 10.1021/acs.jctc.3c00529] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/31/2023]
Abstract
In recent years, adiabatic connection (AC) interpolations developed within density functional theory (DFT) have been found to provide good performances in the calculation of interaction energies when used with Hartree-Fock (HF) ingredients. The physical and mathematical reasons for such unanticipated performance have been clarified, to some extent, by studying the strong-interaction limit of the Møller-Plesset (MP) AC. In this work, we calculate both the MP and the DFT AC integrand for the asymmetric Hubbard dimer, which allows for a systematic investigation of different correlation regimes by varying two simple parameters in the Hamiltonian: the external potential, Δv, and the interaction strength, U. Notably, we find that, while the DFT AC integrand appears to be convex in the full parameter space, the MP integrand may change curvature twice. Furthermore, we discuss different aspects of the second-order expansion of the correlation energy in each AC, and we demonstrate why the derivative of the λ-dependent density in the MP AC at λ = 0 (i.e., at the HF density) is zero in the model. Concerning the strong-interaction limit of both ACs in the Hubbard dimer setting, we show that the asymptotic value of the MP AC, W∞HF, is lower than (or equal to) its DFT analogue, W∞KS, if the two are compared at a given density, just like in real space. However, we also show that this is not always the case if the two quantities are compared at a given external potential.
Collapse
Affiliation(s)
- Sara Giarrusso
- Department of Chemistry and Biochemistry, University of California Merced, 5200 North Lake Rd., Merced, California 95343, United States
| | - Aurora Pribram-Jones
- Department of Chemistry and Biochemistry, University of California Merced, 5200 North Lake Rd., Merced, California 95343, United States
| |
Collapse
|
4
|
Cuierrier E, Roy PO, Ernzerhof M. The factorization ansatz for non-local approximations to the exchange-correlation hole. J Chem Phys 2022; 156:184110. [PMID: 35568557 DOI: 10.1063/5.0077287] [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
Among the various types of approximations to the exchange-correlation energy (EXC), the completely non-local approach is one of the lesser explored approximation schemes. It has not yet reached the predictive power of the widely used generalized gradient approximations, meta-generalized gradient approximations, hybrids, etc. In non-local functionals pursued here, the electron density at every point in space is employed to express the exchange-correlation energy per particle ϵXC(r) at a given position r. Here, we use the non-local, spherical-averaged density ρ(r,u)=∫dΩu4πρ(r+u) as a starting point to construct approximate exchange-correlation holes through the factorization ansatz ρXC(r, u) = f(r, u)ρ(r, u). We present upper and lower bounds to the exchange energy per particle ϵX(r) in terms of ρ(r, u). The factor f(r, u) is then designed to satisfy various conditions that represent important exchange and correlation effects. We assess the resulting approximations and find that the complex, oscillatory structure of ρ(r, u) makes the construction of a corresponding f(r, u) very challenging. This conclusion, identifying the main issue of the non-local approximation, is supported by a detailed analysis of the resulting exchange-correlation holes.
Collapse
Affiliation(s)
- Etienne Cuierrier
- Département de Chimie, Université de Montréal, C.P. 6128 Succursale A, Montréal, Québec H3C 3J7, Canada
| | - Pierre-Olivier Roy
- Département de Chimie, Université de Montréal, C.P. 6128 Succursale A, Montréal, Québec H3C 3J7, Canada
| | - Matthias Ernzerhof
- Département de Chimie, Université de Montréal, C.P. 6128 Succursale A, Montréal, Québec H3C 3J7, Canada
| |
Collapse
|
5
|
Via-Nadal M, Rodríguez Mayorga MA, Ramos Cordoba E, Matito E. Natural Range Separation of the Coulomb Hole. J Chem Phys 2022; 156:184106. [DOI: 10.1063/5.0085284] [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
A natural range separation of the Coulomb hole into two components, one of them being predominant at long interelectronic separations (hcI ) and the other at short distances (hcII ), is exhaustively analyzed throughout various examples that put forward the most relevant features of this approach and how they can be used to develop efficient ways to capture electron correlation. We show that hcI, which only depends on the first-order reduced density matrix, can be used to identify molecules with a predominant nondynamic correlation regime and differentiate between two types of nondynamic correlation, types A and B. Through the asymptotic properties of the hole components, we explain how hcI can retrieve the long-range part of electron correlation. We perform an exhaustive analysis of the hydrogen molecule in a minimal basis set, dissecting the hole contributions into spin components. We also analyze the simplest molecule presenting a dispersion interaction and how hcII helps identify it. The study of several atoms in different spin states reveals that the Coulomb hole components distinguish correlation regimes that are not apparent from the entire hole. The results of this work hold out the promise to aid in developing new electronic structure methods that efficiently capture electron correlation.
Collapse
Affiliation(s)
| | | | - Eloy Ramos Cordoba
- Theoretical Chemistry Group, Donostia International Physics Center, Spain
| | - Eduard Matito
- Donostia International Physics Center, Donostia International Physics Center, Spain
| |
Collapse
|
6
|
Daas T, Fabiano E, Della Sala F, Gori-Giorgi P, Vuckovic S. Noncovalent Interactions from Models for the Møller-Plesset Adiabatic Connection. J Phys Chem Lett 2021; 12:4867-4875. [PMID: 34003655 PMCID: PMC8280728 DOI: 10.1021/acs.jpclett.1c01157] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/10/2021] [Accepted: 05/13/2021] [Indexed: 05/08/2023]
Abstract
Given the omnipresence of noncovalent interactions (NCIs), their accurate simulations are of crucial importance across various scientific disciplines. Here we construct accurate models for the description of NCIs by an interpolation along the Møller-Plesset adiabatic connection (MP AC). Our interpolation approximates the correlation energy, by recovering MP2 at small coupling strengths and the correct large-coupling strength expansion of the MP AC, recently shown to be a functional of the Hartree-Fock density. Our models are size consistent for fragments with nondegenerate ground states, have the same cost as double hybrids, and require no dispersion corrections to capture NCIs accurately. These interpolations greatly reduce large MP2 errors for typical π-stacking complexes (e.g., benzene-pyridine dimers) and for the L7 data set. They are also competitive with state-of-the-art dispersion enhanced functionals and can even significantly outperform them for a variety of data sets, such as CT7 and L7.
Collapse
Affiliation(s)
- Timothy
J. Daas
- 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), Via Monteroni, Campus Unisalento, 73100 Lecce, Italy
- Center
for Biomolecular Nanotechnologies, Istituto
Italiano di Tecnologia, Via Barsanti 14, 73010 Arnesano (LE), Italy
| | - 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
| | - 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
| | - Stefan Vuckovic
- 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
- Physical
and Theoretical Chemistry, University of
Saarland, 66123 Saarbrücken, Germany
- Department
of Chemistry, University of California, Irvine, California 92697, United States
| |
Collapse
|
7
|
Jana S, Śmiga S, Constantin LA, Samal P. Generalizing Double-Hybrid Density Functionals: Impact of Higher-Order Perturbation Terms. J Chem Theory Comput 2020; 16:7413-7430. [PMID: 33205659 PMCID: PMC7735712 DOI: 10.1021/acs.jctc.0c00823] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/06/2020] [Indexed: 11/28/2022]
Abstract
Connections between the Görling-Levy (GL) perturbation theory and the parameters of double-hybrid (DH) density functional are established via adiabatic connection formalism. Moreover, we present a more general DH density functional theory, where the higher-order perturbation terms beyond the second-order GL2 one, such as GL3 and GL4, also contribute. It is shown that a class of DH functionals including previously proposed ones can be formed using the present construction. Based on the proposed formalism, we assess the performance of higher-order DH and long-range corrected DH formed on the Perdew-Burke-Ernzerhof (PBE) semilocal functional and second-order GL2 correlation energy. The underlying construction of DH functionals based on the generalized many-body perturbation approaches is physically appealing in terms of the development of the non-local forms using more accurate and sophisticated semilocal functionals.
Collapse
Affiliation(s)
- Subrata Jana
- School
of Physical Sciences, National Institute
of Science Education and Research, HBNI, Bhubaneswar 752050, India
| | - Szymon Śmiga
- Institute
of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University, Grudziadzka 5, 87-100 Toruń, Poland
| | - Lucian A. Constantin
- Consiglio
Nazionale delle Ricerche CNR-NANO, Istituto
di Nanoscienze, 41125 Modena, Italy
| | - Prasanjit Samal
- School
of Physical Sciences, National Institute
of Science Education and Research, HBNI, Bhubaneswar 752050, India
| |
Collapse
|
8
|
Daas TJ, Grossi J, Vuckovic S, Musslimani ZH, Kooi DP, Seidl M, Giesbertz KJH, Gori-Giorgi P. Large coupling-strength expansion of the Møller–Plesset adiabatic connection: From paradigmatic cases to variational expressions for the leading terms. J Chem Phys 2020; 153:214112. [DOI: 10.1063/5.0029084] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Timothy J. Daas
- Department of Chemistry and Pharmaceutical Sciences, Amsterdam Institute of Molecular and Life Sciences (AIMMS), Faculty of Science, Vrije Universiteit, De Boelelaan 1083, 1081HV Amsterdam, The Netherlands
| | - Juri Grossi
- Department of Chemistry and Pharmaceutical Sciences, Amsterdam Institute of Molecular and Life Sciences (AIMMS), Faculty of Science, Vrije Universiteit, De Boelelaan 1083, 1081HV Amsterdam, The Netherlands
| | - Stefan Vuckovic
- Department of Chemistry, University of California, Irvine, California 92697, USA
| | - Ziad H. Musslimani
- Department of Chemistry and Pharmaceutical Sciences, Amsterdam Institute of Molecular and Life Sciences (AIMMS), Faculty of Science, Vrije Universiteit, De Boelelaan 1083, 1081HV Amsterdam, The Netherlands
| | - Derk P. Kooi
- Department of Chemistry and Pharmaceutical Sciences, Amsterdam Institute of Molecular and Life Sciences (AIMMS), Faculty of Science, Vrije Universiteit, De Boelelaan 1083, 1081HV Amsterdam, The Netherlands
| | - Michael Seidl
- Department of Chemistry and Pharmaceutical Sciences, Amsterdam Institute of Molecular and Life Sciences (AIMMS), Faculty of Science, Vrije Universiteit, De Boelelaan 1083, 1081HV Amsterdam, The Netherlands
| | - Klaas J. H. Giesbertz
- Department of Chemistry and Pharmaceutical Sciences, Amsterdam Institute of Molecular and Life Sciences (AIMMS), Faculty of Science, Vrije Universiteit, De Boelelaan 1083, 1081HV Amsterdam, The Netherlands
| | - Paola Gori-Giorgi
- Department of Chemistry and Pharmaceutical Sciences, Amsterdam Institute of Molecular and Life Sciences (AIMMS), Faculty of Science, Vrije Universiteit, De Boelelaan 1083, 1081HV Amsterdam, The Netherlands
| |
Collapse
|
9
|
Alipour M, Khorrami M. Pauli energy and information-theoretic approach for evaluating dynamic and nondynamic electron correlation. Theor Chem Acc 2020. [DOI: 10.1007/s00214-020-02689-7] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
10
|
Ś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
|
11
|
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
|
12
|
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
|
13
|
Singh R, Patra B, Patra A, Harbola MK, Samal P. Adiabatic connection in density functional theory in two-dimensions: A semi-analytic wavefunction based study for two-electron atomic systems. J Chem Phys 2019; 151:204104. [PMID: 31779332 DOI: 10.1063/1.5120853] [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
This work focuses on studying the adiabatic-connection in density functional theory in two dimensions. It employs a recently developed accurate form of wavefunction for two-electron systems. The explicit semianalytic form of the wavefunction makes it possible to calculate ground state wavefunctions, energies, densities, and the resulting properties for the scaled Coulomb interaction between the electrons at fixed density accurately. The results so obtained for the correlation energies are then used as the reference values for studying the performance of two-dimensional correlation energy functionals.
Collapse
Affiliation(s)
- Rabeet Singh
- School of Physical Sciences, National Institute of Science Education and Research, HBNI, Bhubaneswar 752050, India
| | - Bikash Patra
- School of Physical Sciences, National Institute of Science Education and Research, HBNI, Bhubaneswar 752050, India
| | - Abhilash Patra
- School of Physical Sciences, National Institute of Science Education and Research, HBNI, Bhubaneswar 752050, India
| | - Manoj K Harbola
- Department of Physics, Indian Institute of Technology Kanpur, Kanpur 208016, India
| | - Prasanjit Samal
- School of Physical Sciences, National Institute of Science Education and Research, HBNI, Bhubaneswar 752050, India
| |
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
|
Vuckovic S, Song S, Kozlowski J, Sim E, Burke K. Density Functional Analysis: The Theory of Density-Corrected DFT. J Chem Theory Comput 2019; 15:6636-6646. [DOI: 10.1021/acs.jctc.9b00826] [Citation(s) in RCA: 37] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Stefan Vuckovic
- Departments of Chemistry and of Physics, University of California, Irvine, California 92697, United States
| | - Suhwan Song
- Department of Chemistry, Yonsei University, 50 Yonsei-ro Seodaemun-gu, Seoul 03722, Korea
| | - John Kozlowski
- Departments of Chemistry and of Physics, University of California, Irvine, California 92697, United States
| | - Eunji Sim
- Department of Chemistry, Yonsei University, 50 Yonsei-ro Seodaemun-gu, Seoul 03722, Korea
| | - Kieron Burke
- Departments of Chemistry and of Physics, University of California, Irvine, California 92697, United States
| |
Collapse
|
16
|
Via-Nadal M, Rodríguez-Mayorga M, Ramos-Cordoba E, Matito E. Singling Out Dynamic and Nondynamic Correlation. J Phys Chem Lett 2019; 10:4032-4037. [PMID: 31276421 DOI: 10.1021/acs.jpclett.9b01376] [Citation(s) in RCA: 22] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
The correlation part of the pair density is separated into two components, one of them being predominant at short electronic ranges and the other at long ranges. The analysis of the intracular part of these components permits to classify molecular systems according to the prevailing correlation: dynamic or nondynamic. The study of the long-range asymptotics reveals the key component of the pair density that is responsible for the description of London dispersion forces and a universal decay with the interelectronic distance. The natural range-separation, the identification of the dispersion forces, and the kind of predominant correlation type that arise from this analysis are expected to be important assets in the development of new electronic structure methods in wave function, density, and reduced density-matrix functional theories.
Collapse
Affiliation(s)
- Mireia Via-Nadal
- Donostia International Physics Center (DIPC) , 20018 Donostia , Euskadi , Spain
- Kimika Fakultatea , Euskal Herriko Unibertsitatea (UPV/EHU) , Donostia , Euskadi , Spain
| | - Mauricio Rodríguez-Mayorga
- Donostia International Physics Center (DIPC) , 20018 Donostia , Euskadi , Spain
- Kimika Fakultatea , Euskal Herriko Unibertsitatea (UPV/EHU) , Donostia , Euskadi , Spain
- Institut de Química Computacional i Catàlisi (IQCC) and Departament de Química , Universitat de Girona , C/Maria Aurèlia Capmany, 69 , 17003 Girona , Catalonia , Spain
| | - Eloy Ramos-Cordoba
- Donostia International Physics Center (DIPC) , 20018 Donostia , Euskadi , Spain
- Kimika Fakultatea , Euskal Herriko Unibertsitatea (UPV/EHU) , Donostia , Euskadi , Spain
| | - Eduard Matito
- Donostia International Physics Center (DIPC) , 20018 Donostia , Euskadi , Spain
- IKERBASQUE , Basque Foundation for Science , 48013 Bilbao , Euskadi , Spain
| |
Collapse
|
17
|
Vuckovic S. Density Functionals from the Multiple-Radii Approach: Analysis and Recovery of the Kinetic Correlation Energy. J Chem Theory Comput 2019; 15:3580-3590. [DOI: 10.1021/acs.jctc.9b00129] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Stefan Vuckovic
- Department of Chemistry, University of California, Irvine, California 92697, United States
| |
Collapse
|
18
|
Seidl M, Giarrusso S, Vuckovic S, Fabiano E, Gori-Giorgi P. Communication: Strong-interaction limit of an adiabatic connection in Hartree-Fock theory. J Chem Phys 2018; 149:241101. [DOI: 10.1063/1.5078565] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Michael Seidl
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling, Faculty of Science, Vrije Universiteit, De Boelelaan 1083, 1081HV Amsterdam, The Netherlands
| | - Sara Giarrusso
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling, Faculty of Science, Vrije Universiteit, De Boelelaan 1083, 1081HV Amsterdam, The Netherlands
| | - Stefan Vuckovic
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling, Faculty of Science, Vrije Universiteit, De Boelelaan 1083, 1081HV Amsterdam, The Netherlands
- Department of Chemistry, University of California, Irvine, California 92697, USA
| | - 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
| | - 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
|
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
|
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
|
21
|
Vuckovic S, Gori-Giorgi P, Della Sala F, Fabiano E. Restoring Size Consistency of Approximate Functionals Constructed from the Adiabatic Connection. J Phys Chem Lett 2018; 9:3137-3142. [PMID: 29787273 PMCID: PMC5994725 DOI: 10.1021/acs.jpclett.8b01054] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/05/2018] [Accepted: 05/22/2018] [Indexed: 05/24/2023]
Abstract
Approximate exchange-correlation functionals built by modeling in a nonlinear way the adiabatic connection (AC) integrand of density functional theory have many attractive features, being virtually parameter-free and satisfying different exact properties, but they also have a fundamental flaw: they violate the size-consistency condition, crucial to evaluate interaction energies of molecular systems. We show that size consistency in the AC-based functionals can be restored in a very simple way at no extra computational cost. Results on a large set of benchmark molecular interaction energies show that functionals based on the interaction strength interpolation approximations are significantly more accurate than second-order perturbation theory.
Collapse
Affiliation(s)
- 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
| | - 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
| | - 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
| |
Collapse
|
22
|
Giarrusso S, Gori-Giorgi P, Della Sala F, Fabiano E. Assessment of interaction-strength interpolation formulas for gold and silver clusters. J Chem Phys 2018; 148:134106. [DOI: 10.1063/1.5022669] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022] Open
Affiliation(s)
- Sara Giarrusso
- 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
| | - 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
| | - 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
| |
Collapse
|
23
|
Vuckovic S, Levy M, Gori-Giorgi P. Augmented potential, energy densities, and virial relations in the weak- and strong-interaction limits of DFT. J Chem Phys 2017; 147:214107. [DOI: 10.1063/1.4997311] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022] Open
Affiliation(s)
- Stefan Vuckovic
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling, FEW, Vrije Universiteit, De Boelelaan 1083, 1081HV Amsterdam, The Netherlands
| | - Mel Levy
- Department of Chemistry and Quantum Theory Group, Tulane University, New Orleans, Louisiana 70118, USA
| | - Paola Gori-Giorgi
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling, FEW, Vrije Universiteit, De Boelelaan 1083, 1081HV Amsterdam, The Netherlands
| |
Collapse
|
24
|
Grossi J, Kooi DP, Giesbertz KJH, Seidl M, Cohen AJ, Mori-Sánchez P, Gori-Giorgi P. Fermionic Statistics in the Strongly Correlated Limit of Density Functional Theory. J Chem Theory Comput 2017; 13:6089-6100. [PMID: 29111724 PMCID: PMC5729548 DOI: 10.1021/acs.jctc.7b00998] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Exact pieces of information on the adiabatic connection integrand, Wλ[ρ], which allows evaluation of the exchange-correlation energy of Kohn-Sham density functional theory, can be extracted from the leading terms in the strong coupling limit (λ → ∞, where λ is the strength of the electron-electron interaction). In this work, we first compare the theoretical prediction for the two leading terms in the strong coupling limit with data obtained via numerical implementation of the exact Levy functional in the simple case of two electrons confined in one dimension, confirming the asymptotic exactness of these two terms. We then carry out a first study on the incorporation of the Fermionic statistics at large coupling λ, both numerical and theoretical, confirming that spin effects enter at orders ∼e-√λ.
Collapse
Affiliation(s)
- Juri Grossi
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling, FEW, Vrije Universiteit , De Boelelaan 1083, 1081HV Amsterdam, The Netherlands
| | - Derk P Kooi
- 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
| | - Michael Seidl
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling, FEW, Vrije Universiteit , De Boelelaan 1083, 1081HV Amsterdam, The Netherlands
| | - Aron J Cohen
- Max-Planck Institute for Solid State Research , Heisenbergstrasse 1, 70569 Stuttgart, Germany
| | - Paula Mori-Sánchez
- Departamento de Quimíca and Instituto de Física de la Materia Condensada (IFIMAC), Universidad Autónoma de Madrid , 28049 Madrid, Spain
| | - Paola Gori-Giorgi
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling, FEW, Vrije Universiteit , De Boelelaan 1083, 1081HV Amsterdam, The Netherlands
| |
Collapse
|
25
|
Vuckovic S, Gori-Giorgi P. Simple Fully Nonlocal Density Functionals for Electronic Repulsion Energy. J Phys Chem Lett 2017; 8:2799-2805. [PMID: 28581751 PMCID: PMC5502414 DOI: 10.1021/acs.jpclett.7b01113] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/07/2017] [Accepted: 06/05/2017] [Indexed: 05/24/2023]
Abstract
From a simplified version of the mathematical structure of the strong coupling limit of the exact exchange-correlation functional, we construct an approximation for the electronic repulsion energy at physical coupling strength, which is fully nonlocal. This functional is self-interaction free and yields energy densities within the definition of the electrostatic potential of the exchange-correlation hole that are locally accurate and have the correct asymptotic behavior. The model is able to capture strong correlation effects that arise from chemical bond dissociation, without relying on error cancellation. These features, which are usually missed by standard density functional theory (DFT) functionals, are captured by the highly nonlocal structure, which goes beyond the "Jacob's ladder" framework for functional construction, by using integrals of the density as the key ingredient. Possible routes for obtaining the full exchange-correlation functional by recovering the missing kinetic component of the correlation energy are also implemented and discussed.
Collapse
|