Automatic computation of global intermolecular potential energy surfaces for non-covalently bound systems

Ramón Lorenzo Panadés-Barrueta (1) , Emilio Martinez-Nunez (2) , Daniel Peláez (1)

1 PCMT - Physico-Chimie Moléculaire Théorique
2 Departamento de Química Física, Facultade de Química

We aim at simulating full quantum mechanically (nuclei and electrons) the processes of adsorption and photoreactivity of NO2 adsorbed on soot (modeled as large Polycyclic Aromatic Hydrocarbons, PAHs) in atmospheric conditions. A detailed description of these processes is necessary to understand the differential day-nighttime behavior of the production of HONO [1, 2], which is a precursor of the hydroxyl radical (OH) [3]. In particular, the specific mechanism of the soot-mediated interconversion between NO2 and HONO is to date not fully understood. Due to its relevance, we have chosen the NO2-Pyrene system [1]. The first stage in this study has consisted in the determination of the stable configurations (transition states and minima) of the NO2-Pyrene system. To this end, we have used the van der Waals Transition State Search Using Chemical Dynamics (vdW-TSSCDS) method [4], the non-covalent version of the TSSCDS algorithm [5]. Starting from a single input geometry, vdW-TSSCDS permits the characterization of the topography of a intermolecular Potential Energy Surface (PES) in a fully automated fashion. This topographical information will be used to obtain a global description (fit) of the interaction potential, necessary for the dynamical elucidation of spectroscopic properties and reactivity of the adsorbed species. For this, we have developed the Specific Reaction Parameter Multigrid POTFIT (SRP-MGPF) algorithm and software (SRPTucker) [6]. This method computes chemically accurate (intermolecular) PESs through reparametrization of semiempirical methods which are subsequently tensor decomposed using MGPF [7]. This software has been interfaced with the Heidelberg version of the Multiconfiguration Time-Dependent Hartree (MCTDH) software [8]. References [1] C. Guan, X. Li, W. Zhang and Z. Huang, J. Phys. Chem. A, 121, 482-492 (2017). [2] M. E. Monge et al , PNAS, 107, 6605 (2010) [3] A. M. Holloway and R. P. Wayne, Atmos. Chem. (2010) [4] S. Kopec, E. Martínez-Núñez, J. Soto, D. Peláez, Int. J. Quantum Chem., (accepted) [5] E. Martínez-Núñez, Phys. Chem. Chem. Phys. ,17.22 (2015) [6] R.L.Panadés-Barrueta, E. Martínez-Núñez, D. Peláez, Front. Chem., (submitted) [7] D. Peláez, H.-D. Meyer, J. Chem. Phys. 138, 014108 (2013) [8] M. H. Beck, A. Jäckle, G. A. Worth, and H.-D. Meyer, Phys. Rep., 324 (2000)

Domaines: Chimie théorique et/ou physique; Chimie-Physique [physics.chem-ph]

Ramón Lorenzo Panadés-Barrueta, Emilio Martinez-Nunez, Daniel Peláez. Automatic computation of global intermolecular potential energy surfaces for non-covalently bound systems. 10th International Meeting on Atomic and Molecular Physics and Chemistry 2019, Jun 2019, Madrid, Spain. , 2019. ⟨hal-02157098⟩