1
|
Li X, Qian Y, Chen J. Electric Polarization from a Many-Body Neural Network Ansatz. PHYSICAL REVIEW LETTERS 2024; 132:176401. [PMID: 38728714 DOI: 10.1103/physrevlett.132.176401] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/15/2023] [Revised: 12/01/2023] [Accepted: 03/22/2024] [Indexed: 05/12/2024]
Abstract
Ab initio calculation of dielectric response with high-accuracy electronic structure methods is a long-standing problem, for which mean-field approaches are widely used and electron correlations are mostly treated via approximated functionals. Here we employ a neural network wave function ansatz combined with quantum Monte Carlo method to incorporate correlations into polarization calculations. On a variety of systems, including isolated atoms, one-dimensional chains, two-dimensional slabs, and three-dimensional cubes, the calculated results outperform conventional density functional theory and are consistent with the most accurate calculations and experimental data. Furthermore, we have studied the out-of-plane dielectric constant of bilayer graphene using our method and reestablished its thickness dependence. Overall, this approach provides a powerful tool to accurately describe electron correlation in the modern theory of polarization.
Collapse
Affiliation(s)
- Xiang Li
- ByteDance Research, Zhonghang Plaza, No. 43, North 3rd Ring West Road, Haidian District, Beijing, People's Republic of China
| | - Yubing Qian
- ByteDance Research, Zhonghang Plaza, No. 43, North 3rd Ring West Road, Haidian District, Beijing, People's Republic of China
- School of Physics, Peking University, Beijing 100871, People's Republic of China
| | - Ji Chen
- School of Physics, Peking University, Beijing 100871, People's Republic of China
- Interdisciplinary Institute of Light-Element Quantum Materials, Frontiers Science Center for Nano-Optoelectronics, Peking University, Beijing 100871, People's Republic of China
| |
Collapse
|
2
|
Kitamura N, Kurashige Y. Required nearest-neighbor Coulomb interactions for a charge-ordered phase transition in (TMTTF)2MF6 with inversion symmetry breaking in crystal. Chem Phys Lett 2022. [DOI: 10.1016/j.cplett.2021.139254] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
3
|
Kehrer J, Richter R, Foerster JM, Schelter I, Kümmel S. Self-interaction correction, electrostatic, and structural influences on time-dependent density functional theory excitations of bacteriochlorophylls from the light-harvesting complex 2. J Chem Phys 2020; 153:144114. [PMID: 33086803 DOI: 10.1063/5.0014938] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
First-principles calculations offer the chance to obtain a microscopic understanding of light-harvesting processes. Time-dependent density functional theory can have the computational efficiency to allow for such calculations. However, the (semi-)local exchange-correlation approximations that are computationally most efficient fail to describe charge-transfer excitations reliably. We here investigate whether the inexpensive average density self-interaction correction (ADSIC) remedies the problem. For the systems that we study, ADSIC is even more prone to the charge-transfer problem than the local density approximation. We further explore the recently reported finding that the electrostatic potential associated with the chromophores' protein environment in the light-harvesting complex 2 beneficially shifts spurious excitations. We find a great sensitivity on the chromophores' atomistic structure in this problem. Geometries obtained from classical molecular dynamics are more strongly affected by the spurious charge-transfer problem than the ones obtained from crystallography or density functional theory. For crystal structure geometries and density-functional theory optimized ones, our calculations confirm that the electrostatic potential shifts the spurious excitations out of the energetic range that is most relevant for electronic coupling.
Collapse
Affiliation(s)
- Juliana Kehrer
- Theoretical Physics IV, University of Bayreuth, D-95440 Bayreuth, Germany
| | - Rian Richter
- Theoretical Physics IV, University of Bayreuth, D-95440 Bayreuth, Germany
| | | | - Ingo Schelter
- Theoretical Physics IV, University of Bayreuth, D-95440 Bayreuth, Germany
| | - Stephan Kümmel
- Theoretical Physics IV, University of Bayreuth, D-95440 Bayreuth, Germany
| |
Collapse
|
4
|
Hofmann F, Kümmel S. Molecular excitations from meta-generalized gradient approximations in the Kohn-Sham scheme. J Chem Phys 2020; 153:114106. [PMID: 32962375 DOI: 10.1063/5.0023657] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023] Open
Abstract
Meta-Generalized Gradient Approximations (meta-GGAs) can, in principle, include spatial and temporal nonlocality in time-dependent density functional theory at a much lower computational cost than functionals that use exact exchange. We here test whether a meta-GGA that has recently been developed with a focus on capturing nonlocal response properties and the particle number discontinuity can realize such features in practice. To this end, we extended the frequency-dependent Sternheimer formalism to the meta-GGA case. Using the Krieger-Li-Iafrate (KLI) approximation, we calculate the optical response for the selected paradigm molecular systems and compare the meta-GGA Kohn-Sham response to the one found with exact exchange and conventional (semi-)local functionals. We find that the new meta-GGA captures important properties of the nonlocal exchange response. The KLI approximation, however, emerges as a limiting factor in the evaluation of charge-transfer excitations.
Collapse
Affiliation(s)
- Fabian Hofmann
- Theoretical Physics IV, University of Bayreuth, D-95440 Bayreuth, Germany
| | - Stephan Kümmel
- Theoretical Physics IV, University of Bayreuth, D-95440 Bayreuth, Germany
| |
Collapse
|
5
|
Samanta PK, Blunt NS, Booth GH. Response Formalism within Full Configuration Interaction Quantum Monte Carlo: Static Properties and Electrical Response. J Chem Theory Comput 2018; 14:3532-3546. [DOI: 10.1021/acs.jctc.8b00454] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Pradipta Kumar Samanta
- Institut für Theoretische Chemie, Universität Stuttgart, D-70569 Stuttgart, Germany
- Max Planck Institute for Solid State Research, D-70569 Stuttgart, Germany
| | - Nick S. Blunt
- University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - George H. Booth
- Department of Physics, King’s College London, Strand, London WC2R 2LS, United Kingdom
| |
Collapse
|
6
|
Zhao L, Neuscamman E. A Blocked Linear Method for Optimizing Large Parameter Sets in Variational Monte Carlo. J Chem Theory Comput 2017; 13:2604-2611. [DOI: 10.1021/acs.jctc.7b00119] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Luning Zhao
- Department
of Chemistry, University of California, Berkeley, California 94720, United States
| | - Eric Neuscamman
- Department
of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical
Science Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
7
|
Vargas J, Springborg M, Kirtman B. Electronic responses of long chains to electrostatic fields: Hartree-Fock vs. density-functional theory: a model study. J Chem Phys 2014; 140:054117. [PMID: 24511932 DOI: 10.1063/1.4864038] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023] Open
Abstract
The response to an electrostatic field is determined through simple model calculations, within both the restricted Hartree-Fock and density functional theory methods, for long, finite as well as infinite, periodic chains. The permanent dipole moment, μ0, the polarizability, α, and the hyperpolarizabilities β and γ, calculated using a finite-field approach, are extensively analyzed. Our simple model allows for treatment of large systems and for separation of the properties into atomic and unit-cell contributions. That part of the response properties attributable to the terminations of the finite system change into delocalized current contributions in the corresponding infinite periodic system. Special emphasis is placed on analyzing the reasons behind the dramatic overestimation of the response properties found with density functional theory methods presently in common use.
Collapse
Affiliation(s)
- Jorge Vargas
- Physical and Theoretical Chemistry, University of Saarland, 66123 Saarbrücken, Germany
| | - Michael Springborg
- Physical and Theoretical Chemistry, University of Saarland, 66123 Saarbrücken, Germany
| | - Bernard Kirtman
- Department of Chemistry and Biochemistry, University of California, Santa Barbara, California 93106, USA
| |
Collapse
|
8
|
Wouters S, Limacher PA, Van Neck D, Ayers PW. Longitudinal static optical properties of hydrogen chains: Finite field extrapolations of matrix product state calculations. J Chem Phys 2012; 136:134110. [DOI: 10.1063/1.3700087] [Citation(s) in RCA: 51] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
9
|
Umari P, Marzari N. Linear and nonlinear susceptibilities from diffusion quantum Monte Carlo: Application to periodic hydrogen chains. J Chem Phys 2009; 131:094104. [DOI: 10.1063/1.3213567] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
10
|
Murillo MS, Dharma-wardana MWC. Temperature relaxation in hot dense hydrogen. PHYSICAL REVIEW LETTERS 2008; 100:205005. [PMID: 18518546 DOI: 10.1103/physrevlett.100.205005] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/07/2007] [Indexed: 05/26/2023]
Abstract
Temperature equilibration of hydrogen is studied for conditions relevant to inertial confinement fusion. New molecular-dynamics simulations and results from quantum many-body theory are compared with Landau-Spitzer predictions for temperatures T with 50<T<4000 eV and densities with Wigner-Seitz radii rs=1.0 and 0.5. The relaxation is slower than the Landau-Spitzer result, even for T in the kilo-electron-volt range, but converge to agreement in the high-T limit.
Collapse
Affiliation(s)
- Michael S Murillo
- Physics Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA.
| | | |
Collapse
|
11
|
Körzdörfer T, Mundt M, Kümmel S. Electrical response of molecular systems: the power of self-interaction corrected kohn-sham theory. PHYSICAL REVIEW LETTERS 2008; 100:133004. [PMID: 18517945 DOI: 10.1103/physrevlett.100.133004] [Citation(s) in RCA: 37] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/31/2007] [Indexed: 05/26/2023]
Abstract
The accurate prediction of electronic response properties of extended molecular systems has been a challenge for conventional, explicit density functionals. We demonstrate that a self-interaction correction (SIC) implemented rigorously within Kohn-Sham theory via the optimized effective potential (OEP) yields polarizabilities close to the ones from highly accurate wave-function-based calculations and exceeding the quality of exact-change OEP. The orbital structure obtained with the OEP-SIC functional and approximations to it are discussed.
Collapse
Affiliation(s)
- T Körzdörfer
- Physics Institute, University of Bayreuth, D-95440 Bayreuth, Germany
| | | | | |
Collapse
|