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
|
Huan Lew-Yee JF, Piris M, Del Campo JM. Outstanding improvement in removing the delocalization error by global natural orbital functional. J Chem Phys 2023; 158:084110. [PMID: 36859086 DOI: 10.1063/5.0137378] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/11/2023] Open
Abstract
This work assesses the performance of the recently proposed global natural orbital functional (GNOF) against the charge delocalization error. GNOF provides a good balance between static and dynamic electronic correlations leading to accurate total energies while preserving spin, even for systems with a highly multi-configurational character. Several analyses were applied to the functional, namely, (i) how the charge is distributed in super-systems of two fragments, (ii) the stability of ionization potentials while increasing the system size, and (iii) potential energy curves of a neutral and charged diatomic system. GNOF was found to practically eliminate the charge delocalization error in many of the studied systems or greatly improve the results obtained previously with PNOF7.
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, Mexico City C.P. 04510, Mexico
| | - 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
| | - Jorge M Del Campo
- Departamento de Física y Química Teórica, Facultad de Química, Universidad Nacional Autónoma de México, Mexico City C.P. 04510, Mexico
| |
Collapse
|
3
|
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
|
4
|
Mercero JM, Grande-Aztatzi R, Ugalde JM, Piris M. Natural orbital functional theory studies of all-metal aromaticity: The Al 3−anion. ADVANCES IN QUANTUM CHEMISTRY 2023. [DOI: 10.1016/bs.aiq.2023.02.006] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 03/28/2023]
|
5
|
Lemke Y, Kussmann J, Ochsenfeld C. Efficient Integral-Direct Methods for Self-Consistent Reduced Density Matrix Functional Theory Calculations on Central and Graphics Processing Units. J Chem Theory Comput 2022; 18:4229-4244. [DOI: 10.1021/acs.jctc.2c00231] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Y. Lemke
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstr. 5−13, D-81377 Munich, Germany
| | - J. Kussmann
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstr. 5−13, D-81377 Munich, Germany
| | - C. Ochsenfeld
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstr. 5−13, D-81377 Munich, Germany
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, D-70569 Stuttgart, Germany
| |
Collapse
|
6
|
Quintero-Monsebaiz R, Perea-Ramírez LI, Piris M, Vela A. Spectroscopic properties of open shell diatomic molecules using Piris natural orbital functionals. Phys Chem Chem Phys 2021; 23:2953-2963. [PMID: 33480904 DOI: 10.1039/d0cp05430e] [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/21/2022]
Abstract
Spectroscopic properties such as equilibrium distances, vibrational constants, rotational constants, dissociation energies, and excitation energies are calculated for nine heteronuclear diatomic molecules (PH, NF, NH, NO, CS, AlF, ClF, BeO and CF) using an interactive pair model (PNOF7s), that has been generalized for spin multiplet states, and its second order perturbation variant, NOF-MP2, which was also generalized for multiplets. The results obtained are compared with Complete Active Space (CASSCF) and Complete Active Space Perturbation Theory (CASPT2). It is shown that the potential energy curves provided by the PNOF functional for open shell diatomic molecules are in acceptable agreement with those from CASSCF and CASPT2. The spectroscopic constants depending at most on the second derivative of the potential energy are in good agreement with experiment, while those requiring the evaluation of the third and fourth derivatives show larger deviations from experiment and from those predicted by CASPT2. Thus, it is shown that the PNOF functional extension to multiplets is an alternative approach in predicting spectroscopic constants of molecules where static correlation plays an important role, like the open shell heteronuclear diatomic molecules studied in this work.
Collapse
Affiliation(s)
- Raúl Quintero-Monsebaiz
- Centro de Investigación y de Estudios Avanzados del Instituto Politécnico Nacional, Av Instituto Politécnico Nacional 2508, San Pedro Zacatenco, Gustavo A. Madero, 07360 Ciudad de México, Mexico.
| | - Luis Ignacio Perea-Ramírez
- Centro de Investigación y de Estudios Avanzados del Instituto Politécnico Nacional, Av Instituto Politécnico Nacional 2508, San Pedro Zacatenco, Gustavo A. Madero, 07360 Ciudad de México, Mexico.
| | - Mario Piris
- Donostia International Physics Center (DIPC), 20018 Donostia, Euskadi, Spain and Euskal Herriko Unibertsitatea (UPV/EHU), PK 1072, 20080 Donostia, Euskadi, Spain and Basque Foundation for Science (IKERBASQUE), Plaza Euskadi 5, 48009 Bilbao, Euskadi, Spain.
| | - Alberto Vela
- Centro de Investigación y de Estudios Avanzados del Instituto Politécnico Nacional, Av Instituto Politécnico Nacional 2508, San Pedro Zacatenco, Gustavo A. Madero, 07360 Ciudad de México, Mexico.
| |
Collapse
|