1
|
Kinetic Energy Density Functionals Based on a Generalized Screened Coulomb Potential: Linear Response and Future Perspectives. COMPUTATION 2022. [DOI: 10.3390/computation10020030] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Abstract
We consider kinetic energy functionals that depend, beside the usual semilocal quantities (density, gradient, Laplacian of the density), on a generalized Yukawa potential, that is the screened Coulomb potential of the density raised to some power. These functionals, named Yukawa generalized gradient approximations (yGGA), are potentially efficient real-space semilocal methods that include significant non-local effects and can describe different important exact properties of the kinetic energy. In this work, we focus in particular on the linear response behavior for the homogeneous electron gas (HEG). We show that such functionals are able to reproduce the exact Lindhard function behavior with a very good accuracy, outperforming all other semilocal kinetic functionals. These theoretical advances allow us to perform a detailed analysis of a special class of yGGAs, namely the linear yGGA functionals. Thus, we show how the present approach can generalize the yGGA functionals improving the HEG linear behavior and leading to an extended formula for the kinetic functional. Moreover, testing on several jellium cluster model systems allows highlighting advantages and limitations of the linear yGGA functionals and future perspectives for the development of yGGA kinetic functionals.
Collapse
|
2
|
Jana S, Myneni H, Śmiga S, Constantin LA, Samal P. Benchmark test of a dispersion corrected revised Tao-Mo semilocal functional for thermochemistry, kinetics, and noncovalent interactions of molecules and solids. J Chem Phys 2021; 155:114102. [PMID: 34551544 DOI: 10.1063/5.0060538] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Abstract
In the density functional theory, dispersion corrected semilocal approximations are often used to benchmark weekly interacting finite and extended systems. Here, the focus is on providing a broad overview of the performance of D3 dispersion corrected revised Tao-Mo (revTM) semilocal functionals [A. Patra et al., J. Chem. Phys. 153, 084 117 (2020)] for thermochemistry and kinetics of molecules, molecular crystals, ice polymorphs, metal-organic systems, atom/molecular adsorption on solids, water interacting with nano-materials, binding energies of layered materials, and properties of weekly and strongly bonded solids. We show that the most suitable "optimized power" function for the revTM functional needs a modification to make it suitable for properties related to the diverse nature of finite and extended systems. The present work is an extension of the previously proposed revTM+D3 method with the motivation to design and benchmark the dispersion corrected cost-effective method based on this semilocal approximation. We show that the revised revTM+D3 functional provides various general purpose molecular and solid properties with the closest to experimental findings than its predecessor. The present assessment and benchmarking can be practically useful for performing cost-effective method based simulations of various molecular and solid-state properties.
Collapse
Affiliation(s)
- Subrata Jana
- School of Physical Sciences, National Institute of Science Education and Research, HBNI, Bhubaneswar 752 050, India
| | - Hemanadhan Myneni
- Science Institute and Faculty of Physical Sciences, University of Iceland, VR-III, 107 Reykjavík, Iceland
| | - Szymon Śmiga
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University, Grudziadzka 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 752 050, India
| |
Collapse
|
3
|
Giannone G, Śmiga S, D'Agostino S, Fabiano E, Della Sala F. Plasmon Couplings from Subsystem Time-Dependent Density Functional Theory. J Phys Chem A 2021; 125:7246-7259. [PMID: 34403247 DOI: 10.1021/acs.jpca.1c05384] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Many applications in plasmonics are related to the coupling between metallic nanoparticles (MNPs) or between an emitter and a MNP. The theoretical analysis of such a coupling is thus of fundamental importance to analyze the plasmonic behavior and to design new systems. While classical methods neglect quantum and spill-out effects, time-dependent density functional theory (TD-DFT) considers all of them and with Kohn-Sham orbitals delocalized over the whole system. Thus, within TD-DFT, no definite separation of the subsystems (the single MNP or the emitter) and their couplings is directly available. This important feature is obtained here using the subsystem formulation of TD-DFT, which has been originally developed in the context of weakly interacting organic molecules. In subsystem TD-DFT, interacting MNPs are treated independently, thus allowing us to compute the plasmon couplings directly from the subsystem TD-DFT transition densities. We show that subsystem TD-DFT, as well as a simplified version of it in which kinetic contributions are neglected, can reproduce the reference TD-DFT calculations for gap distances greater than about 6 Å or even smaller in the case of hybrid plasmonic systems (i.e., molecules interacting with MNPs). We also show that the subsystem TD-DFT can be also used as a tool to analyze the impact of charge-transfer effects.
Collapse
Affiliation(s)
- Giulia Giannone
- Center for Biomolecular Nanotechnologies, Istituto Italiano di Tecnologia, Via Barsanti 14, Arnesano (LE) 73010, Italy.,Department of Mathematics and Physics "E. De Giorgi", University of Salento, Via Arnesano, Lecce 73100, Italy
| | - Szymon Śmiga
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University, Grudzia̧dzka 5, Toruń 87-100, Poland
| | - Stefania D'Agostino
- Center for Biomolecular Nanotechnologies, Istituto Italiano di Tecnologia, Via Barsanti 14, Arnesano (LE) 73010, Italy.,Department of Mathematics and Physics "E. De Giorgi", University of Salento, Via Arnesano, Lecce 73100, Italy.,Institute of Nanotechnology, National Research Council (CNR-NANOTEC), c/o Campus Ecotekne, via Monteroni, Lecce 73100, Italy
| | - Eduardo Fabiano
- Center for Biomolecular Nanotechnologies, Istituto Italiano di Tecnologia, Via Barsanti 14, Arnesano (LE) 73010, Italy.,Institute for Microelectronics and Microsystems (CNR-IMM), Via Monteroni, Campus Unisalento, Lecce 73100, Italy
| | - Fabio Della Sala
- Center for Biomolecular Nanotechnologies, Istituto Italiano di Tecnologia, Via Barsanti 14, Arnesano (LE) 73010, Italy.,Institute for Microelectronics and Microsystems (CNR-IMM), Via Monteroni, Campus Unisalento, Lecce 73100, Italy
| |
Collapse
|
4
|
Jana S, Behera SK, Śmiga S, Constantin LA, Samal P. Accurate density functional made more versatile. J Chem Phys 2021; 155:024103. [PMID: 34266258 DOI: 10.1063/5.0051331] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Abstract
We propose a one-electron self-interaction-free correlation energy functional compatible with the order-of-limit problem-free Tao-Mo (TM) semilocal functional (regTM) [J. Tao and Y. Mo, Phys. Rev. Lett. 117, 073001 (2016) and Patra et al., J. Chem. Phys. 153, 184112 (2020)] to be used for general purpose condensed matter physics and quantum chemistry. The assessment of the proposed functional for large classes of condensed matter and chemical systems shows its improvement in most cases compared to the TM functional, e.g., when applied to the relative energy difference of MnO2 polymorphs. In this respect, the present exchange-correction functional, which incorporates the TM technique of the exchange hole model combined with the slowly varying density correction, can achieve broad applicability, being able to solve difficult solid-state problems.
Collapse
Affiliation(s)
- Subrata Jana
- School of Physical Sciences, National Institute of Science Education and Research, HBNI, Bhubaneswar 752050, India
| | - Sushant Kumar Behera
- 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
- 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
|
5
|
Francisco HI, Carmona-Espíndola J, Gázquez JL. Analysis of the kinetic energy functional in the generalized gradient approximation. J Chem Phys 2021; 154:084107. [PMID: 33639771 DOI: 10.1063/5.0040973] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A new density functional for the total kinetic energy in the generalized gradient approximation is developed through an enhancement factor that leads to the correct behavior in the limits when the reduced density gradient tends to 0 and to infinity and by making use of the conjoint conjecture for the interpolation between these two limits, through the incorporation, in the intermediate region of constraints that are associated with the exchange energy functional. The resulting functional leads to a reasonable description of the kinetic energies of atoms and molecules when it is used in combination with Hartree-Fock densities. Additionally, in order to improve the behavior of the kinetic energy density, a new enhancement factor for the Pauli kinetic energy is proposed by incorporating the correct behavior into the limits when the reduced density gradient tends to 0 and to infinity, together with the positivity condition, and imposing through the interpolation function that the sum of its integral over the whole space and the Weiszacker energy must be equal to the value obtained with the enhancement factor developed for the total kinetic energy.
Collapse
Affiliation(s)
- Héctor I Francisco
- Departamento de Química, Universidad Autónoma Metropolitana-Iztapalapa, Av. San Rafael Atlixco 186, México, Ciudad de México 09340, Mexico
| | - Javier Carmona-Espíndola
- Departamento de Química, CONACYT-Universidad Autónoma Metropolitana-Iztapalapa, Av. San Rafael Atlixco 186, Ciudad de México 09340, Mexico
| | - José L Gázquez
- Departamento de Química, Universidad Autónoma Metropolitana-Iztapalapa, Av. San Rafael Atlixco 186, México, Ciudad de México 09340, Mexico
| |
Collapse
|
6
|
Furness JW, Kaplan AD, Ning J, Perdew JP, Sun J. Accurate and Numerically Efficient r 2SCAN Meta-Generalized Gradient Approximation. J Phys Chem Lett 2020; 11:8208-8215. [PMID: 32876454 DOI: 10.1021/acs.jpclett.0c02405] [Citation(s) in RCA: 230] [Impact Index Per Article: 57.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
The recently proposed rSCAN functional [ J. Chem. Phys. 2019 150, 161101] is a regularized form of the SCAN functional [ Phys. Rev. Lett. 2015 115, 036402] that improves SCAN's numerical performance at the expense of breaking constraints known from the exact exchange-correlation functional. We construct a new meta-generalized gradient approximation by restoring exact constraint adherence to rSCAN. The resulting functional maintains rSCAN's numerical performance while restoring the transferable accuracy of SCAN.
Collapse
Affiliation(s)
- James W Furness
- Department of Physics and Engineering Physics, Tulane University, New Orleans, Louisiana 70118, United States
| | - Aaron D Kaplan
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, United States
| | - Jinliang Ning
- Department of Physics and Engineering Physics, Tulane University, New Orleans, Louisiana 70118, United States
| | - John P Perdew
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, United States
- Department of Chemistry, Temple University, Philadelphia, Pennsylvania 19122, United States
| | - Jianwei Sun
- Department of Physics and Engineering Physics, Tulane University, New Orleans, Louisiana 70118, United States
| |
Collapse
|
7
|
Ś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
|
8
|
The Role of the Reduced Laplacian Renormalization in the Kinetic Energy Functional Development. COMPUTATION 2019. [DOI: 10.3390/computation7040065] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
The Laplacian of the electronic density diverges at the nuclear cusp, which complicates the development of Laplacian-level meta-GGA (LLMGGA) kinetic energy functionals for all-electron calculations. Here, we investigate some Laplacian renormalization methods, which avoid this divergence. We developed two different LLMGGA functionals, which improve the kinetic energy or the kinetic potential. We test these KE functionals in the context of Frozen-Density-Embedding (FDE), for a large palette of non-covalently interacting molecular systems. These functionals improve over the present state-of-the-art LLMGGA functionals for the FDE calculations.
Collapse
|
9
|
Witt WC, Jiang K, Carter EA. Upper bound to the gradient-based kinetic energy density of noninteracting electrons in an external potential. J Chem Phys 2019. [DOI: 10.1063/1.5108896] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- William C. Witt
- Department of Mechanical and Aerospace Engineering, Princeton University, Princeton, New Jersey 08544-5263, USA
| | - Kaili Jiang
- Department of Mechanical and Aerospace Engineering, Princeton University, Princeton, New Jersey 08544-5263, USA
| | - Emily A. Carter
- School of Engineering and Applied Science, Princeton University, Princeton, New Jersey 08544-5263, USA
| |
Collapse
|
10
|
Tran F, Kovács P, Kalantari L, Madsen GKH, Blaha P. Orbital-free approximations to the kinetic-energy density in exchange-correlation MGGA functionals: Tests on solids. J Chem Phys 2018; 149:144105. [DOI: 10.1063/1.5048907] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023] Open
Affiliation(s)
- Fabien Tran
- Institute of Materials Chemistry, Vienna University of Technology, Getreidemarkt 9/165-TC, A-1060 Vienna, Austria
| | - Péter Kovács
- Institute of Materials Chemistry, Vienna University of Technology, Getreidemarkt 9/165-TC, A-1060 Vienna, Austria
| | - Leila Kalantari
- Institute of Materials Chemistry, Vienna University of Technology, Getreidemarkt 9/165-TC, A-1060 Vienna, Austria
| | - Georg K. H. Madsen
- Institute of Materials Chemistry, Vienna University of Technology, Getreidemarkt 9/165-TC, A-1060 Vienna, Austria
| | - Peter Blaha
- Institute of Materials Chemistry, Vienna University of Technology, Getreidemarkt 9/165-TC, A-1060 Vienna, Austria
| |
Collapse
|
11
|
Modified Fourth-Order Kinetic Energy Gradient Expansion with Hartree Potential-Dependent Coefficients. J Chem Theory Comput 2017; 13:4228-4239. [DOI: 10.1021/acs.jctc.7b00705] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
12
|
Śmiga S, Fabiano E, Constantin LA, Della Sala F. Laplacian-dependent models of the kinetic energy density: Applications in subsystem density functional theory with meta-generalized gradient approximation functionals. J Chem Phys 2017; 146:064105. [DOI: 10.1063/1.4975092] [Citation(s) in RCA: 39] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022] Open
|