1
|
Rivero Santamaría A, Piris M. Time evolution of natural orbitals in ab initio molecular dynamics. J Chem Phys 2024; 160:071102. [PMID: 38364005 DOI: 10.1063/5.0188491] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/21/2023] [Accepted: 01/28/2024] [Indexed: 02/18/2024] Open
Abstract
This work combines for the first time ab initio molecular dynamics (AIMD) within the Born-Oppenheimer approximation with a global natural orbital functional (GNOF), an approximate functional of the one-particle reduced density matrix. The most prominent feature of GNOF-AIMD is its ability to display the real-time evolution of natural orbitals, providing detailed information on the time-dependent electronic structure of complex systems and processes, including reactive collisions. The quartet ground-state reaction N(4S) + H2(1Σ) → NH(3Σ) + H(2S) is taken as a validation test. Collision energy influences on integral cross sections for different initial rovibrational states of H2 and rotational-state distributions of the NH product are discussed, showing a good agreement with previous high-quality theoretical results.
Collapse
Affiliation(s)
| | - Mario Piris
- Donostia International Physics Center (DIPC), 20018 Donostia, Euskadi, Spain; Euskal Herriko Unibertsitatea (UPV/EHU), PK 1072, 20080 Donostia, Euskadi, Spain; and Basque Foundation for Science (IKERBASQUE), 48009 Bilbao, Euskadi, Spain
| |
Collapse
|
2
|
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
|
3
|
Mazziotti DA. Quantum Many-Body Theory from a Solution of the N-Representability Problem. PHYSICAL REVIEW LETTERS 2023; 130:153001. [PMID: 37115895 DOI: 10.1103/physrevlett.130.153001] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/15/2022] [Revised: 02/13/2023] [Accepted: 03/16/2023] [Indexed: 06/19/2023]
Abstract
Here we present a many-body theory based on a solution of the N-representability problem in which the ground-state two-particle reduced density matrix (2-RDM) is determined directly without the many-particle wave function. We derive an equation that re-expresses physical constraints on higher-order RDMs to generate direct constraints on the 2-RDM, which are required for its derivation from an N-particle density matrix, known as N-representability conditions. The approach produces a complete hierarchy of 2-RDM constraints that do not depend explicitly upon the higher RDMs or the wave function. By using the two-particle part of a unitary decomposition of higher-order constraint matrices, we can solve the energy minimization by semidefinite programming in a form where the low-rank structure of these matrices can be potentially exploited. We illustrate by computing the ground-state electronic energy and properties of the H_{8} ring.
Collapse
Affiliation(s)
- David A Mazziotti
- Department of Chemistry and The James Franck Institute, The University of Chicago, Chicago, Illinois 60637, USA
| |
Collapse
|
4
|
Cioslowski J, Schilling C, Schilling R. 1-Matrix functional for long-range interaction energy of two hydrogen atoms. J Chem Phys 2023; 158:084106. [PMID: 36859076 DOI: 10.1063/5.0139897] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023] Open
Abstract
The leading terms in the large-R asymptotics of the functional of the one-electron reduced density matrix for the ground-state energy of the H2 molecule with the internuclear separation R are derived thanks to the solution of the phase dilemma at the R → ∞ limit. At this limit, the respective natural orbitals (NOs) are given by symmetric and antisymmetric combinations of "half-space" orbitals with the corresponding natural amplitudes having the same amplitudes but opposite signs. Minimization of the resulting explicit functional yields the large-R asymptotics for the occupation numbers of the weakly occupied NOs and the C6 dispersion coefficient. The highly accurate approximates for the radial components of the p-type "half-space" orbitals and the corresponding occupation numbers (that decay like R-6), which are available for the first time thanks to the development of the present formalism, have some unexpected properties.
Collapse
Affiliation(s)
- Jerzy Cioslowski
- Institute of Physics, University of Szczecin, Wielkopolska 15, 70-451 Szczecin, Poland
| | - Christian Schilling
- Department of Physics, Arnold Sommerfeld Center for Theoretical Physics, Ludwig-Maximilians-Universität München, Theresienstrasse 37, 80333 München, Germany
| | - Rolf Schilling
- Institute of Physics, Johannes Gutenberg University, Staudinger Weg 9, 55099 Mainz, Germany
| |
Collapse
|
5
|
Lew-Yee JFH, Del Campo JM, Piris M. Electron Correlation in the Iron(II) Porphyrin by Natural Orbital Functional Approximations. J Chem Theory Comput 2023; 19:211-220. [PMID: 36579972 PMCID: PMC9996833 DOI: 10.1021/acs.jctc.2c01093] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
The relative stability of the singlet, triplet, and quintet spin states of iron(II) porphyrin (FeP) represents a challenging problem for electronic structure methods. While it is currently accepted that the ground state is a triplet, multiconfigurational wave function-based methods predict a quintet, and density functional approximations vary between triplet and quintet states, leading to a prediction that highly depends on the features of the method employed. The recently proposed Global Natural Orbital Functional (GNOF) aims to provide a balanced treatment between static and dynamic correlation, and together with the previous Piris Natural Orbital Functionals (PNOFs), allowed us to explore the importance of each type of correlation in the stability order of the states of FeP with a method that conserves the spin of the system. It is noteworthy that GNOF correlates all electrons in all available orbitals for a given basis set; in the case of the FeP with a double-ζ basis set as used in this work, this means that GNOF can properly correlate 186 electrons in 465 orbitals, significantly increasing the sizes of systems amenable to multiconfigurational treatment. Results show that PNOF5, PNOF7s, and PNOF7 predict the quintet to have a lower energy than the triplet state; however, the addition of dynamic correlation via second-order Møller-Plesset corrections (NOF-MP2) turns the triplet state to be lower than the quintet state, a prediction also reproduced by GNOF that incorporates much more dynamic correlation than its predecessors.
Collapse
Affiliation(s)
- Juan Felipe Huan Lew-Yee
- Departamento de Física y Química Teórica, Facultad de Química, Universidad Nacional Autónoma de México, México CityC.P. 04510, México
| | - Jorge M Del Campo
- Departamento de Física y Química Teórica, Facultad de Química, Universidad Nacional Autónoma de México, México CityC.P. 04510, México
| | - Mario Piris
- Kimika Fakultatea, Euskal Herriko Unibertsitatea (UPV/EHU), P.K. 1072, 20080Donostia, Euskadi, Spain.,Donostia International Physics Center (DIPC), 20018Donostia, Euskadi, Spain.,IKERBASQUE, Basque Foundation for Science, 48013Bilbao, Euskadi, Spain
| |
Collapse
|
6
|
Gedeon J, Schmidt J, Hodgson MJP, Wetherell J, Benavides-Riveros CL, Marques MAL. Machine learning the derivative discontinuity of density-functional theory. MACHINE LEARNING: SCIENCE AND TECHNOLOGY 2022. [DOI: 10.1088/2632-2153/ac3149] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
Abstract
Machine learning is a powerful tool to design accurate, highly non-local, exchange-correlation functionals for density functional theory. So far, most of those machine learned functionals are trained for systems with an integer number of particles. As such, they are unable to reproduce some crucial and fundamental aspects, such as the explicit dependency of the functionals on the particle number or the infamous derivative discontinuity at integer particle numbers. Here we propose a solution to these problems by training a neural network as the universal functional of density-functional theory that (a) depends explicitly on the number of particles with a piece-wise linearity between the integer numbers and (b) reproduces the derivative discontinuity of the exchange-correlation energy. This is achieved by using an ensemble formalism, a training set containing fractional densities, and an explicitly discontinuous formulation.
Collapse
|
7
|
Gibney D, Boyn JN, Mazziotti DA. Toward a Resolution of the Static Correlation Problem in Density Functional Theory from Semidefinite Programming. J Phys Chem Lett 2021; 12:385-391. [PMID: 33356286 DOI: 10.1021/acs.jpclett.0c03371] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
Kohn-Sham density functional theory (DFT) has long struggled with the accurate description of strongly correlated and open shell systems, and improvements have been minor even in the newest hybrid functionals. In this Letter we treat the static correlation in DFT when frontier orbitals are degenerate by the means of using a semidefinite programming (SDP) approach to minimize the system energy as a function of the N-representable, non-idempotent 1-electron reduced density matrix. While showing greatly improved singlet-triplet gaps for local density approximation and generalized gradient approximation (GGA) functionals, the SDP procedure reveals flaws in modern meta and hybrid GGA functionals, which show no major improvements when provided with an accurate electron density.
Collapse
Affiliation(s)
- Daniel Gibney
- The James Franck Institute and The Department of Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| | - Jan-Niklas Boyn
- The James Franck Institute and The Department of Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| | - David A Mazziotti
- The James Franck Institute and The Department of Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
8
|
Cioslowski J. Off-diagonal derivative discontinuities in the reduced density matrices of electronic systems. J Chem Phys 2020; 153:154108. [PMID: 33092376 DOI: 10.1063/5.0023955] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Abstract
An explicit expression relating the magnitude of the fifth-order off-diagonal cusp in the real part of the one-electron reduced density matrix to the "on-top" two-electron density is derived in a rigorous manner from the behavior of the underlying electronic wavefunction at the electron-electron coalescence. The implications of the presence of this cusp upon electronic structure calculations of quantum chemistry and solid-state physics, including the limits imposed upon their accuracy, are elucidated. In particular, the power-law decay of the occupation numbers of the natural orbitals is demonstrated for 1S states of systems composed of arbitrary even numbers of electrons. The practical importance of analogous off-diagonal cusps in many-electron reduced density matrices is briefly discussed.
Collapse
Affiliation(s)
- Jerzy Cioslowski
- Institute of Physics, University of Szczecin, Wielkopolska 15, 70-451 Szczecin, Poland
| |
Collapse
|