1
|
Woźniak AP, Adamowicz L, Pedersen TB, Kvaal S. Gaussians for Electronic and Rovibrational Quantum Dynamics. J Phys Chem A 2024; 128:3659-3671. [PMID: 38687971 PMCID: PMC11089519 DOI: 10.1021/acs.jpca.4c00364] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2024] [Revised: 03/23/2024] [Accepted: 04/08/2024] [Indexed: 05/02/2024]
Abstract
The assumptions underpinning the adiabatic Born-Oppenheimer (BO) approximation are broken for molecules interacting with attosecond laser pulses, which generate complicated coupled electronic-nuclear wave packets that generally will have components of electronic and dissociation continua as well as bound-state contributions. The conceptually most straightforward way to overcome this challenge is to treat the electronic and nuclear degrees of freedom on equal quantum-mechanical footing by not invoking the BO approximation at all. Explicitly correlated Gaussian (ECG) basis functions have proved successful for non-BO calculations of stationary molecular states and energies, reproducing rovibrational absorption spectra with very high accuracy. In this Article, we present a proof-of-principle study of the ability of fully flexible ECGs (FFECGs) to capture the intricate electronic and rovibrational dynamics generated by short, high-intensity laser pulses. By fitting linear combinations of FFECGs to accurate wave function histories obtained on a large real-space grid for a regularized 2D model of the hydrogen atom and for the 2D Morse potential, we demonstrate that FFECGs provide a very compact description of laser-driven electronic and rovibrational dynamics.
Collapse
Affiliation(s)
| | - Ludwik Adamowicz
- Department
of Chemistry and Biochemistry, University
of Arizona, 1306 E University Blvd, Tucson, Arizona 85721-0041, United States
| | - Thomas Bondo Pedersen
- Hylleraas
Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, N-0315 Oslo, Norway
| | - Simen Kvaal
- Hylleraas
Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, N-0315 Oslo, Norway
| |
Collapse
|
2
|
Lang L, Cezar HM, Adamowicz L, Pedersen TB. Quantum Definition of Molecular Structure. J Am Chem Soc 2024; 146:1760-1764. [PMID: 38199236 PMCID: PMC10811664 DOI: 10.1021/jacs.3c11467] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/16/2023] [Revised: 11/21/2023] [Accepted: 11/22/2023] [Indexed: 01/12/2024]
Abstract
Molecular structure, a key concept of chemistry, has remained elusive from the perspective of all-particle quantum mechanics, despite many efforts. Viewing molecular structure as a manifestation of strong statistical correlation between nuclear positions, we propose a practical method based on Markov chain Monte Carlo sampling and unsupervised machine learning. Application to the D3+ molecule unambiguously shows that it possesses an equilateral triangular structure. These results provide a major step forward in our understanding of the molecular structure from fundamental quantum principles.
Collapse
Affiliation(s)
- Lucas Lang
- Hylleraas
Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, 0315 Oslo, Norway
- Technische
Universität Berlin, Institut für Chemie, Theoretische Chemie/Quantenchemie, Sekr. C7, Straße des 17. Juni
135, 10623 Berlin, Germany
| | - Henrique M. Cezar
- Hylleraas
Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, 0315 Oslo, Norway
| | - Ludwik Adamowicz
- Centre
for Advanced Study at the Norwegian Academy of Science and Letters, Drammensveien 78, 0271 Oslo, Norway
- Department
of Chemistry and Biochemistry, University
of Arizona, Tucson, Arizona 85721, United States
| | - Thomas B. Pedersen
- Hylleraas
Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, 0315 Oslo, Norway
- Centre
for Advanced Study at the Norwegian Academy of Science and Letters, Drammensveien 78, 0271 Oslo, Norway
| |
Collapse
|
3
|
Quantum states of a confined hydrogen atom calculated in a basis of explicitly correlated Gaussians. Chem Phys Lett 2022. [DOI: 10.1016/j.cplett.2022.139358] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|
4
|
Adamowicz L, Stanke M, Tellgren E, Helgaker T. A computational quantum-mechanical model of a molecular magnetic trap. J Chem Phys 2018; 149:244112. [PMID: 30599715 DOI: 10.1063/1.5055767] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A model for describing the states of a molecular system trapped in a cavity created by a fast-rotating strong magnetic field is proposed and implemented. All-particle explicitly correlated Gaussian functions with shifted centers are employed in the model to expand the wave functions of the system. Both "internal" states associated with the system's rovibrational and electronic motions and the "external" states associated with translational motion of the center of mass of the system in the cavity are calculated. The states are visualized by density plots. The model is applied to a trapped HD molecule.
Collapse
Affiliation(s)
- Ludwik Adamowicz
- Department of Chemistry and Biochemistry, University of Arizona, Tucson, Arizona 85721, USA and Interdisciplinary Center for Modern Technologies, Nicolaus Copernicus University, ul. Wileńska 4, Toruń, PL 87-100, Poland
| | - Monika Stanke
- Institute of Physics, Faculty of Physics, Astronomy, and Informatics, Nicolaus Copernicus University, ul. Grudzia̧dzka 5, Toruń PL 87-100, Poland
| | - Erik Tellgren
- Hyllerras Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, N-0315 Oslo, Norway
| | - Trygve Helgaker
- Hyllerras Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, N-0315 Oslo, Norway
| |
Collapse
|
5
|
Goli M, Shahbazian S. Hidden aspects of the Structural theory of chemistry: MC-QTAIM analysis reveals "alchemical" transformation from a triatomic to a diatomic structure. Phys Chem Chem Phys 2015; 17:245-55. [PMID: 25388361 DOI: 10.1039/c4cp03722g] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
The Structural theory of chemistry introduces chemical/molecular structure as a combination of relative arrangement and bonding patterns of atoms in a molecule. Nowadays, the structure of atoms in molecules is derived from the topological analysis of the quantum theory of atoms in molecules (QTAIM). In this context, a molecular structure is varied by large geometrical variations and concomitant reorganization of electronic structure that usually take place in chemical reactions or under extreme hydrostatic pressure. In this report, a new mode of structural variation is introduced within the context of the newly proposed multi-component QTAIM (MC-QTAIM) that originates from the mass variation of nuclei. Accordingly, XCN and CNX series of species are introduced where X stands for a quantum particle with a unit of positive charge and a variable mass that is varied in discrete steps between the mass of a proton and a positron. Ab initio non-Born-Oppenheimer (non-BO) calculations are done on both series of species and the resulting non-BO wavefunctions are used for the MC-QTAIM analysis, revealing a triatomic structure for the proton mass and a diatomic structure for the positron mass. In both series of species, a critical mass between that of proton and positron mass is discovered where the transition from triatomic to diatomic structure takes place. This abrupt structural transformation has a topological nature resembling the usual phase transitions in thermodynamics. The discovered mass-induced structural transformation is a hidden aspect of the Structural theory which is revealed only beyond the BO paradigm, when nuclei are treated as quantum waves instead of clamped point charges.
Collapse
Affiliation(s)
- Mohammad Goli
- Faculty of Chemistry, Shahid Beheshti University, G. C., Evin, P.O. Box 19395-4716, Tehran, Iran19839.
| | | |
Collapse
|
6
|
Rodríguez C, Urbina A, Torres F, Cazar D, Ludeña E. Non-Born–Oppenheimer nuclear and electronic densities for a three-particle Hooke–Coulomb model. COMPUT THEOR CHEM 2013. [DOI: 10.1016/j.comptc.2013.05.033] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
7
|
Mátyus E, Reiher M. Molecular structure calculations: A unified quantum mechanical description of electrons and nuclei using explicitly correlated Gaussian functions and the global vector representation. J Chem Phys 2012; 137:024104. [DOI: 10.1063/1.4731696] [Citation(s) in RCA: 53] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
8
|
Goli M, Shahbazian S. The two-component quantum theory of atoms in molecules (TC-QTAIM): foundations. Theor Chem Acc 2012. [DOI: 10.1007/s00214-012-1208-9] [Citation(s) in RCA: 33] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
9
|
Ludeña EV, Echevarría L, Lopez X, Ugalde JM. Non-Born-Oppenheimer electronic and nuclear densities for a Hooke-Calogero three-particle model: non-uniqueness of density-derived molecular structure. J Chem Phys 2012; 136:084103. [PMID: 22380028 DOI: 10.1063/1.3682244] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We consider the calculation of non-Born-Oppenheimer, nBO, one-particle densities for both electrons and nuclei. We show that the nBO one-particle densities evaluated in terms of translationally invariant coordinates are independent of the wavefunction describing the motion of center of mass of the whole system. We show that they depend, however, on an arbitrary reference point from which the positions of the vectors labeling the particles are determined. We examine the effect that this arbitrary choice has on the topology of the one-particle density by selecting the Hooke-Calogero model of a three-body system for which expressions for the one-particle densities can be readily obtained in analytic form. We extend this analysis to the one-particle densities obtained from full Coulomb interaction wavefunctions for three-body systems. We conclude, in view of the fact that there is a close link between the choice of the reference point and the topology of one-particle densities that the molecular structure inferred from the topology of these densities is not unique. We analyze the behavior of one-particle densities for the Hooke-Calogero Born-Oppenheimer, BO, wavefunction and show that topological transitions are also present in this case for a particular mass value of the light particles even though in the BO regime the nuclear masses are infinite. In this vein, we argue that the change in topology caused by variation of the mass ratio between light and heavy particles does not constitute a true indication in the nBO regime of the emergence of molecular structure.
Collapse
Affiliation(s)
- E V Ludeña
- Centro de Química, Instituto Venezolano de Investigaciones Científicas, IVIC, Apartado 21827, Caracas 1020-A, Venezuela.
| | | | | | | |
Collapse
|
10
|
Schmidling D. High accuracy interpolation of diatomic pseudo-potentials with high order Hermite splines. COMPUT THEOR CHEM 2011. [DOI: 10.1016/j.comptc.2010.10.002] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
11
|
Pachucki K, Komasa J. Nonadiabatic corrections to rovibrational levels of H2. J Chem Phys 2009; 130:164113. [DOI: 10.1063/1.3114680] [Citation(s) in RCA: 138] [Impact Index Per Article: 9.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
12
|
Ishimoto T, Tachikawa M, Nagashima U. Simultaneous analytical optimization of variational parameters in Gaussian-type functions with full configuration interaction of multicomponent molecular orbital method by elimination of translational and rotational motions: Application to isotopomers of the hydrogen molecule. J Chem Phys 2008; 128:164118. [DOI: 10.1063/1.2912939] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
13
|
Miyamoto K, Hoshino M, Nakai H. Elimination of Translational and Rotational Motions in Nuclear Orbital Plus Molecular Orbital Theory: Contribution of the First-Order Rovibration Coupling. J Chem Theory Comput 2006; 2:1544-50. [DOI: 10.1021/ct6002065] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Kaito Miyamoto
- Department of Chemistry, School of Science and Engineering, Waseda University, Tokyo 169-8555, Japan
| | - Minoru Hoshino
- Department of Chemistry, School of Science and Engineering, Waseda University, Tokyo 169-8555, Japan
| | - Hiromi Nakai
- Department of Chemistry, School of Science and Engineering, Waseda University, Tokyo 169-8555, Japan
| |
Collapse
|
14
|
Hoshino M, Nakai H. Elimination of translational and rotational motions in nuclear orbital plus molecular orbital theory: Application of Møller-Plesset perturbation theory. J Chem Phys 2006; 124:194110. [PMID: 16729806 DOI: 10.1063/1.2193513] [Citation(s) in RCA: 26] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
The translation- and rotation-free nuclear orbital plus molecular orbital (TRF-NOMO) theory was developed to determine the nonadiabatic nuclear and electronic wave functions. This study presents a formulation of TRF-NOMO second-order Moller-Plesset (MP2) perturbation and Epstein-Nesbet (EN) theory with the use of the TRF Hamiltonian. Numerical assessment of the TRF-NOMO/MP2 and EN is performed for several molecules. We confirm the importance of the elimination of translational and rotational motions in the many-body calculations.
Collapse
Affiliation(s)
- Minoru Hoshino
- Department of Chemistry, School of Science and Engineering, Waseda University, Tokyo 169-8555, Japan
| | | |
Collapse
|