1
|
Conte R, Aieta C, Cazzaniga M, Ceotto M. A Perspective on the Investigation of Spectroscopy and Kinetics of Complex Molecular Systems with Semiclassical Approaches. J Phys Chem Lett 2024; 15:7566-7576. [PMID: 39024505 DOI: 10.1021/acs.jpclett.4c01338] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/20/2024]
Abstract
In this Perspective we show that semiclassical methods provide a rigorous hierarchical way to study the vibrational spectroscopy and kinetics of complex molecular systems. The time averaged approach to spectroscopy and the semiclassical transition state theory for kinetics, which have been first adopted and then further developed in our group, provide accurate quantum results on rigorous physical grounds and can be applied even when dealing with a large number of degrees of freedom. In spectroscopy, the multiple coherent, divide-and-conquer, and adiabatically switched semiclassical approaches have practically permitted overcoming issues related to the convergence of results. In this Perspective we demonstrate the possibility of studying the semiclassical vibrational spectroscopy of a molecule adsorbed on an anatase (101) surface, a system made of 51 atoms. In kinetics, the semiclassical transition state theory is able to account for anharmonicity and the coupling between the reactive and bound modes. Our group has developed this technique for practical applications involving the study of phenomena like kinetic isotope effect, heavy atom tunneling, and elusive conformer lifetimes. Here, we show that our multidimensional anharmonic quantum approach is able to tackle on-the-fly the thermal kinetic rate constant of a 135 degree-of-freedom system. Overall, semiclassical methods open up the possibility to describe at the quantum mechanical level systems characterized by hundreds of degrees of freedom leading to the accurate spectroscopic and kinetic description of biomolecules and complex molecular systems.
Collapse
Affiliation(s)
- Riccardo Conte
- Dipartimento di Chimica, Università degli Studi di Milano, via Golgi 19, 20133 Milano, Italy
| | - Chiara Aieta
- Dipartimento di Chimica, Università degli Studi di Milano, via Golgi 19, 20133 Milano, Italy
| | - Marco Cazzaniga
- Dipartimento di Chimica, Università degli Studi di Milano, via Golgi 19, 20133 Milano, Italy
| | - Michele Ceotto
- Dipartimento di Chimica, Università degli Studi di Milano, via Golgi 19, 20133 Milano, Italy
| |
Collapse
|
2
|
Singh K, Lee KH, Peláez D, Bande A. Accelerating wavepacket propagation with machine learning. J Comput Chem 2024. [PMID: 39031712 DOI: 10.1002/jcc.27443] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/15/2023] [Revised: 05/13/2024] [Accepted: 05/16/2024] [Indexed: 07/22/2024]
Abstract
In this work, we discuss the use of a recently introduced machine learning (ML) technique known as Fourier neural operators (FNO) as an efficient alternative to the traditional solution of the time-dependent Schrödinger equation (TDSE). FNOs are ML models which are employed in the approximated solution of partial differential equations. For a wavepacket propagating in an anharmonic potential and for a tunneling system, we show that the FNO approach can accurately and faithfully model wavepacket propagation via the density. Additionally, we demonstrate that FNOs can be a suitable replacement for traditional TDSE solvers in cases where the results of the quantum dynamical simulation are required repeatedly such as in the case of parameter optimization problems (e.g., control). The speed-up from the FNO method allows for its combination with the Markov-chain Monte Carlo approach in applications that involve solving inverse problems such as optimal and coherent laser control of the outcome of dynamical processes.
Collapse
Affiliation(s)
- Kanishka Singh
- Theory of Electron Dynamics and Spectroscopy, Helmholtz-Zentrum Berlin für Materialien und Energie GmbH, Berlin, Germany
- Institute of Chemistry and Biochemistry, Freie Universität Berlin, Berlin, Germany
| | - Ka Hei Lee
- Theory of Electron Dynamics and Spectroscopy, Helmholtz-Zentrum Berlin für Materialien und Energie GmbH, Berlin, Germany
- Fachbereich Physik, Freie Universität Berlin, Berlin, Germany
| | - Daniel Peláez
- CNRS, Institut des Sciences Moléculaires d'Orsay, Université Paris-Saclay, Orsay, France
| | - Annika Bande
- Theory of Electron Dynamics and Spectroscopy, Helmholtz-Zentrum Berlin für Materialien und Energie GmbH, Berlin, Germany
- Institute of Inorganic Chemistry, Leibniz University Hannover, Hannover, Germany
- Cluster of Excellence PhoenixD, Leibniz University Hannover, Hannover, Germany
| |
Collapse
|
3
|
Lanzi C, Aieta C, Ceotto M, Conte R. A time averaged semiclassical approach to IR spectroscopy. J Chem Phys 2024; 160:214107. [PMID: 38828809 DOI: 10.1063/5.0214037] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/15/2024] [Accepted: 05/14/2024] [Indexed: 06/05/2024] Open
Abstract
We propose a new semiclassical approach to the calculation of molecular IR spectra. The method employs the time averaging technique of Kaledin and Miller upon symmetrization of the quantum dipole-dipole autocorrelation function. Spectra at high and low temperatures are investigated. In the first case, we are able to point out the possible presence of hot bands in the molecular absorption line shape. In the second case, we are able to reproduce accurate IR spectra as demonstrated by a calculation of the IR spectrum of the water molecule, which is within 4% of the exact intensity. Our time averaged IR spectra can be directly compared to time averaged semiclassical power spectra as shown in an application to the CO2 molecule, which points out the differences between IR and power spectra and demonstrates that our new approach can identify active IR transitions correctly. Overall, the method features excellent accuracy in calculating absorption intensities and provides estimates for the frequencies of vibrations in agreement with the corresponding power spectra. In perspective, this work opens up the possibility to interface the new method with the semiclassical techniques developed for power spectra, such as the divide-and-conquer one, to get accurate IR spectra of complex and high-dimensional molecular systems.
Collapse
Affiliation(s)
- Cecilia Lanzi
- Dipartimento di Chimica, Università degli Studi di Milano, via Golgi 19, 20133 Milano, Italy
| | - Chiara Aieta
- Dipartimento di Chimica, Università degli Studi di Milano, via Golgi 19, 20133 Milano, Italy
| | - Michele Ceotto
- Dipartimento di Chimica, Università degli Studi di Milano, via Golgi 19, 20133 Milano, Italy
| | - Riccardo Conte
- Dipartimento di Chimica, Università degli Studi di Milano, via Golgi 19, 20133 Milano, Italy
| |
Collapse
|
4
|
Moscato D, Mandelli G, Bondanza M, Lipparini F, Conte R, Mennucci B, Ceotto M. Unraveling Water Solvation Effects with Quantum Mechanics/Molecular Mechanics Semiclassical Vibrational Spectroscopy: The Case of Thymidine. J Am Chem Soc 2024; 146:8179-8188. [PMID: 38470354 DOI: 10.1021/jacs.3c12700] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/13/2024]
Abstract
We introduce a quantum mechanics/molecular mechanics semiclassical method for studying the solvation process of molecules in water at the nuclear quantum mechanical level with atomistic detail. We employ it in vibrational spectroscopy calculations because this is a tool that is very sensitive to the molecular environment. Specifically, we look at the vibrational spectroscopy of thymidine in liquid water. We find that the C═O frequency red shift and the C═C frequency blue shift, experienced by thymidyne upon solvation, are mainly due to reciprocal polarization effects, that the molecule and the water solvent exert on each other, and nuclear zero-point energy effects. In general, this work provides an accurate and practical tool to study quantum vibrational spectroscopy in solution and condensed phase, incorporating high-level and computationally affordable descriptions of both electronic and nuclear problems.
Collapse
Affiliation(s)
- Davide Moscato
- Dipartimento di Chimica, Università degli Studi di Milano, Via Golgi, 19, 20133 Milano, Italy
| | - Giacomo Mandelli
- Dipartimento di Chimica, Università degli Studi di Milano, Via Golgi, 19, 20133 Milano, Italy
| | - Mattia Bondanza
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via Giuseppe Moruzzi, 13, 56124 Pisa, Italy
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via Giuseppe Moruzzi, 13, 56124 Pisa, Italy
| | - Riccardo Conte
- Dipartimento di Chimica, Università degli Studi di Milano, Via Golgi, 19, 20133 Milano, Italy
| | - Benedetta Mennucci
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via Giuseppe Moruzzi, 13, 56124 Pisa, Italy
| | - Michele Ceotto
- Dipartimento di Chimica, Università degli Studi di Milano, Via Golgi, 19, 20133 Milano, Italy
| |
Collapse
|
5
|
Mino L, Cazzaniga M, Moriggi F, Ceotto M. Elucidating NO x Surface Chemistry at the Anatase (101) Surface in TiO 2 Nanoparticles. THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2023; 127:437-449. [PMID: 36660096 PMCID: PMC9841571 DOI: 10.1021/acs.jpcc.2c07489] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 10/25/2022] [Revised: 12/12/2022] [Indexed: 06/17/2023]
Abstract
Understanding NO x chemistry at titania nanoparticle surfaces is important for photocatalytic environmental remediation processes. We focus on this problem and put forward an experimental-computational approach based on vibrational spectroscopy grounds. Temperature-dependent IR experiments of NO x adsorption on shape-engineered nanoparticle (101) anatase surfaces are paired with power spectra obtained from Born-Oppenheimer trajectories. Then, the harmonic versus anharmonic vibrational frequencies of several adsorption scenarios are directly compared with the IR experiments. We conclude that molecules are adsorbed mainly by the N-end side and both the intermolecular interactions between adsorbed molecules and (NO)2 dimer formation are responsible for the main NO adsorption spectroscopic features. We also investigate the spectroscopy and the mechanism of formation on defective anatase surfaces of the long-lived greenhouse gas N2O.
Collapse
Affiliation(s)
- Lorenzo Mino
- Department
of Chemistry and NIS Centre, University
of Torino, Via Giuria
7, I-10125 Torino, Italy
| | - Marco Cazzaniga
- Dipartimento
di Chimica, Università Degli Studi
di Milano, Via Golgi 19, I-20133 Milano, Italy
| | - Francesco Moriggi
- Dipartimento
di Chimica, Università Degli Studi
di Milano, Via Golgi 19, I-20133 Milano, Italy
| | - Michele Ceotto
- Dipartimento
di Chimica, Università Degli Studi
di Milano, Via Golgi 19, I-20133 Milano, Italy
| |
Collapse
|
6
|
Conte R, Nandi A, Qu C, Yu Q, Houston PL, Bowman JM. Semiclassical and VSCF/VCI Calculations of the Vibrational Energies of trans- and gauche-Ethanol Using a CCSD(T) Potential Energy Surface. J Phys Chem A 2022; 126:7709-7718. [PMID: 36240438 DOI: 10.1021/acs.jpca.2c06322] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
A recent full-dimensional Δ-Machine learning potential energy surface (PES) for ethanol is employed in semiclassical and vibrational self-consistent field (VSCF) and virtual-state configuration interaction (VCI) calculations, using MULTIMODE, to determine the anharmonic vibrational frequencies of vibration for both the trans and gauche conformers of ethanol. Both semiclassical and VSCF/VCI energies agree well with the experimental data. We find significant mixing between the VSCF basis states due to Fermi resonances between bending and stretching modes. The same effects are also accurately described by the full-dimensional semiclassical calculations. These are the first high-level anharmonic calculations using a PES, in particular a "gold-standard" CCSD(T) one.
Collapse
Affiliation(s)
- Riccardo Conte
- Dipartimento di Chimica, Università degli Studi di Milano, via Golgi 19, 20133 Milano, Italy
| | - Apurba Nandi
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - Chen Qu
- Independent Researcher, Toronto, Ontario M9B0E3, Canada
| | - Qi Yu
- Department of Chemistry Yale University, New Haven, Connecticut 06520, United States
| | - Paul L Houston
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States.,Department of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332, United States
| | - Joel M Bowman
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| |
Collapse
|
7
|
Galimberti DR. Vibrational Circular Dichroism from DFT Molecular Dynamics: The AWV Method. J Chem Theory Comput 2022; 18:6217-6230. [PMID: 36112978 PMCID: PMC9558311 DOI: 10.1021/acs.jctc.2c00736] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
The paper illustrates the Activity Weighted Velocities
(AWV) methodology
to compute Vibrational Circular Dichroism (VCD) anharmonic spectra
from Density Functional Theory (DFT) molecular dynamics. AWV calculates
the spectra by the Fourier Transform of the time correlation functions
of velocities, weighted by specific observables: the Atomic Polar
Tensors (APTs) and the Atomic Axial Tensors (AATs). Indeed, AWV shows
to correctly reproduce the experimental spectra for systems in the
gas and liquid phases, both in the case of weakly and strongly interacting
systems. The comparison with the experimental spectra is striking
especially in the fingerprint region, as demonstrated by the three
benchmark systems discussed: (1S)-Fenchone in the
gas phase, (S)-(−)-Propylene oxide in the
liquid phase, and (R)-(−)-2-butanol in the
liquid phase. The time evolution of APTs and AATs can be adequately
described by a linear combination of the tensors of a small set of
appropriate reference structures, strongly reducing the computational
cost without compromising accuracy. Additionally, AWV allows the partition
of the spectral signal in its molecular components without any expensive
postprocessing and any localization of the charge density or the wave
function.
Collapse
Affiliation(s)
- Daria Ruth Galimberti
- Institute for Molecules and Materials, Radboud University, Heyendaalseweg 135, 6525 AJ Nijmegen, The Netherlands
| |
Collapse
|
8
|
Nandi A, Conte R, Qu C, Houston PL, Yu Q, Bowman JM. Quantum Calculations on a New CCSD(T) Machine-Learned Potential Energy Surface Reveal the Leaky Nature of Gas-Phase Trans and Gauche Ethanol Conformers. J Chem Theory Comput 2022; 18:5527-5538. [PMID: 35951990 PMCID: PMC9476654 DOI: 10.1021/acs.jctc.2c00760] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
![]()
Ethanol is a molecule of fundamental interest in combustion,
astrochemistry,
and condensed phase as a solvent. It is characterized by two methyl
rotors and trans (anti) and gauche conformers, which are known to be very close in energy.
Here we show that based on rigorous quantum calculations of the vibrational
zero-point state, using a new ab initio potential
energy surface (PES), the ground state resembles the trans conformer, but substantial delocalization to the gauche conformer is present. This explains experimental issues about identification
and isolation of the two conformers. This “leak” effect
is partially quenched when deuterating the OH group, which further
demonstrates the need for a quantum mechanical approach. Diffusion
Monte Carlo and full-dimensional semiclassical dynamics calculations
are employed. The new PES is obtained by means of a Δ-machine
learning approach starting from a pre-existing low level density functional
theory surface. This surface is brought to the CCSD(T) level of theory
using a relatively small number of ab initio CCSD(T)
energies. Agreement between the corrected PES and direct ab
initio results for standard tests is excellent. One- and
two-dimensional discrete variable representation calculations focusing
on the trans–gauche torsional
motion are also reported, in reasonable agreement with experiment.
Collapse
Affiliation(s)
- Apurba Nandi
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - Riccardo Conte
- Dipartimento di Chimica, Università Degli Studi di Milano, via Golgi 19, 20133 Milano, Italy
| | - Chen Qu
- Independent Researcher, Toronto 66777, Canada
| | - Paul L Houston
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States.,Department of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332, United States
| | - Qi Yu
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | - Joel M Bowman
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| |
Collapse
|
9
|
Cazzaniga M, Micciarelli M, Gabas F, Finocchi F, Ceotto M. Quantum Anharmonic Calculations of Vibrational Spectra for Water Adsorbed on Titania Anatase(101) Surface: Dissociative versus Molecular Adsorption. THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2022; 126:12060-12073. [PMID: 35928238 PMCID: PMC9340806 DOI: 10.1021/acs.jpcc.2c02137] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 03/29/2022] [Revised: 06/10/2022] [Indexed: 06/15/2023]
Abstract
The interaction of water molecules and hydroxyl groups with titanium dioxide (TiO2) surfaces is ubiquitous and very important in anatase nanoparticle photocatalytic processes. Infrared spectroscopy, assisted by ab initio calculations of vibrational frequencies, can be a powerful tool to elucidate the mechanisms behind water adsorption. However, a straightforward comparison between measurements and calculations remains a challenging task because of the complexity of the physical phenomena occurring on nanoparticle surfaces. Consequently, severe computational approximations, such as harmonic vibrational ones, are usually employed. In the present work we partially address this complexity issue by overcoming some of the standard approximations used in theoretical simulations and employ the Divide and Conquer Semiclassical Initial Value Representation (DC-SCIVR) molecular dynamics. This method allows to perform simulations of vibrational spectra of large dimensional systems accounting not only for anharmonicities, but also for nuclear quantum effects. We apply this computational method to water and deuterated water adsorbed on the ideal TiO2 anatase(101) surface, contemplating both the molecular and the dissociated adsorption processes. The results highlight not only the presence of an anharmonic shift of the frequencies in agreement with the experiments, but also complex quantum mechanical spectral signatures induced by the coupling of molecular vibrational modes with the surface ones, which are different in the hydrogenated case from the deuterated one. These couplings are further analyzed by exploiting the mode subdivision performed during the divide and conquer procedure.
Collapse
Affiliation(s)
- Marco Cazzaniga
- Dipartimento
di Chimica, Universitá degli Studi
di Milano, via Golgi 19, 20133 Milano, Italy
| | - Marco Micciarelli
- Dipartimento
di Chimica, Universitá degli Studi
di Milano, via Golgi 19, 20133 Milano, Italy
| | - Fabio Gabas
- Dipartimento
di Chimica, Universitá degli Studi
di Milano, via Golgi 19, 20133 Milano, Italy
| | - Fabio Finocchi
- Sorbonne
Université, CNRS, Institut des NanoSciences
de Paris (INSP), 4 Place
Jussieu, Paris F- 75005, France
| | - Michele Ceotto
- Dipartimento
di Chimica, Universitá degli Studi
di Milano, via Golgi 19, 20133 Milano, Italy
| |
Collapse
|
10
|
Bowman JM, Qu C, Conte R, Nandi A, Houston PL, Yu Q. The MD17 datasets from the perspective of datasets for gas-phase “small” molecule potentials. J Chem Phys 2022; 156:240901. [DOI: 10.1063/5.0089200] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
There has been great progress in developing methods for machine-learned potential energy surfaces. There have also been important assessments of these methods by comparing so-called learning curves on datasets of electronic energies and forces, notably the MD17 database. The dataset for each molecule in this database generally consists of tens of thousands of energies and forces obtained from DFT direct dynamics at 500 K. We contrast the datasets from this database for three “small” molecules, ethanol, malonaldehyde, and glycine, with datasets we have generated with specific targets for the potential energy surfaces (PESs) in mind: a rigorous calculation of the zero-point energy and wavefunction, the tunneling splitting in malonaldehyde, and, in the case of glycine, a description of all eight low-lying conformers. We found that the MD17 datasets are too limited for these targets. We also examine recent datasets for several PESs that describe small-molecule but complex chemical reactions. Finally, we introduce a new database, “QM-22,” which contains datasets of molecules ranging from 4 to 15 atoms that extend to high energies and a large span of configurations.
Collapse
Affiliation(s)
- Joel M. Bowman
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| | - Chen Qu
- Independent Researcher, Toronto, Canada
| | - Riccardo Conte
- Dipartimento di Chimica, Università Degli Studi di Milano, via Golgi 19, 20133 Milano, Italy
| | - Apurba Nandi
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| | - Paul L. Houston
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
- Department of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
| | - Qi Yu
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, USA
| |
Collapse
|
11
|
Botti G, Aieta C, Conte R. The complex vibrational spectrum of proline explained through the adiabatically switched semiclassical initial value representation. J Chem Phys 2022; 156:164303. [DOI: 10.1063/5.0089720] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Proline, a 17-atom amino acid with a closed-ring side chain, has a complex potential energy surface characterized by several minima. Its IR experimental spectrum, reported in the literature, is of difficult and controversial assignment. In particular, the experimental signal at 3559 cm−1 associated with the OH stretch is interesting because it is inconsistent with the global minimum, trans-proline conformer. This suggests the possibility that multiple conformers may contribute to the IR spectrum. The same conclusion is obtained by investigating the splitting of the CO stretch at 1766 and 1789 cm−1 and other, more complex spectroscopic features involving CH stretches and COH/CNH bendings. In this work, we perform full-dimensional, on-the-fly adiabatically switched semiclassical initial value representation simulations employing the ab initio dft-d3-B3LYP level of theory with aug-cc-pVDZ basis set. We reconstruct the experimental spectrum of proline in its main features by studying the vibrational features of trans-proline and cis1-proline, and provide a new assignment for the OH stretch of trans-proline.
Collapse
Affiliation(s)
| | - Chiara Aieta
- Dipartimento di Chimica, Università degli Studi di Milano, Italy
| | - Riccardo Conte
- Dipartimento di Chimica, Università degli Studi di Milano, Italy
| |
Collapse
|
12
|
Gabas F, Conte R, Ceotto M. Quantum Vibrational Spectroscopy of Explicitly Solvated Thymidine in Semiclassical Approximation. J Phys Chem Lett 2022; 13:1350-1355. [PMID: 35109652 PMCID: PMC8842300 DOI: 10.1021/acs.jpclett.1c04087] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/15/2021] [Accepted: 01/28/2022] [Indexed: 06/14/2023]
Abstract
In this paper, we demonstrate the possibility to perform spectroscopy simulations of solvated biological species taking into consideration quantum effects and explicit solvation. We achieve this goal by interfacing our recently developed divide-and-conquer approach for semiclassical initial value representation molecular dynamics with the polarizable AMOEBABIO18 force field. The method is applied to the study of solvation of the thymidine nucleoside in two different polar solvents, water and N,N-dimethylformamide. Such systems are made of up to 2476 atoms. Experimental evidence concerning the different behavior of thymidine in the two solvents is well reproduced by our study, even though quantitative estimates are hampered by the limited accuracy of the classical force field employed. Overall, this study shows that semiclassically approximate quantum dynamical studies of explicitly solvated biological systems are both computationally affordable and insightful.
Collapse
|
13
|
Product Vibrational State Distributions of the F CH 3OH Reaction on a Full-Dimensional Accurate Potential Energy Surface. CHINESE J CHEM PHYS 2022. [DOI: 10.1063/1674-0068/cjcp2111252] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022]
|
14
|
Botti G, Ceotto M, Conte R. On-the-fly adiabatically switched semiclassical initial value representation molecular dynamics for vibrational spectroscopy of biomolecules. J Chem Phys 2021; 155:234102. [PMID: 34937370 DOI: 10.1063/5.0075220] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Semiclassical (SC) vibrational spectroscopy is a technique capable of reproducing quantum effects (such as zero-point energies, quantum resonances, and anharmonic overtones) from classical dynamics runs even in the case of very large dimensional systems. In a previous study [Conte et al. J. Chem. Phys. 151, 214107 (2019)], a preliminary sampling based on adiabatic switching has been shown to be able to improve the precision and accuracy of semiclassical results for challenging model potentials and small molecular systems. In this paper, we investigate the possibility to extend the technique to larger (bio)molecular systems whose dynamics must be integrated by means of ab initio "on-the-fly" calculations. After some preliminary tests on small molecules, we obtain the vibrational frequencies of glycine improving on pre-existing SC calculations. Finally, the new approach is applied to 17-atom proline, an amino acid characterized by a strong intramolecular hydrogen bond.
Collapse
Affiliation(s)
- Giacomo Botti
- Dipartimento di Chimica, Università degli Studi di Milano, via Golgi 19, 20133 Milano, Italy
| | - Michele Ceotto
- Dipartimento di Chimica, Università degli Studi di Milano, via Golgi 19, 20133 Milano, Italy
| | - Riccardo Conte
- Dipartimento di Chimica, Università degli Studi di Milano, via Golgi 19, 20133 Milano, Italy
| |
Collapse
|
15
|
Gandolfi M, Ceotto M. Unsupervised Machine Learning Neural Gas Algorithm for Accurate Evaluations of the Hessian Matrix in Molecular Dynamics. J Chem Theory Comput 2021; 17:6733-6746. [PMID: 34705463 PMCID: PMC8582248 DOI: 10.1021/acs.jctc.1c00707] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2021] [Indexed: 11/29/2022]
Abstract
The Hessian matrix of the potential energy of molecular systems is employed not only in geometry optimizations or high-order molecular dynamics integrators but also in many other molecular procedures, such as instantaneous normal mode analysis, force field construction, instanton calculations, and semiclassical initial value representation molecular dynamics, to name a few. Here, we present an algorithm for the calculation of the approximated Hessian in molecular dynamics. The algorithm belongs to the family of unsupervised machine learning methods, and it is based on the neural gas idea, where neurons are molecular configurations whose Hessians are adopted for groups of molecular dynamics configurations with similar geometries. The method is tested on several molecular systems of different dimensionalities both in terms of accuracy and computational time versus calculating the Hessian matrix at each time-step, that is, without any approximation, and other Hessian approximation schemes. Finally, the method is applied to the on-the-fly, full-dimensional simulation of a small synthetic peptide (the 46 atom N-acetyl-l-phenylalaninyl-l-methionine amide) at the level of DFT-B3LYP-D/6-31G* theory, from which the semiclassical vibrational power spectrum is calculated.
Collapse
Affiliation(s)
- Michele Gandolfi
- Dipartimento di Chimica, Università
degli Studi di Milano, via Golgi 19, 20133 Milano, Italy
| | - Michele Ceotto
- Dipartimento di Chimica, Università
degli Studi di Milano, via Golgi 19, 20133 Milano, Italy
| |
Collapse
|
16
|
Rognoni A, Conte R, Ceotto M. Caldeira-Leggett model vs ab initio potential: A vibrational spectroscopy test of water solvation. J Chem Phys 2021; 154:094106. [PMID: 33685187 DOI: 10.1063/5.0040494] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present a semiclassically approximate quantum treatment of solvation with the purpose of investigating the accuracy of the Caldeira-Leggett model. We do that by simulating the vibrational features of water solvation by means of two different approaches. One is entirely based on the adoption of an accurate ab initio potential to describe water clusters of increasing dimensionality. The other one consists of a model made of a central water molecule coupled to a high-dimensional Caldeira-Leggett harmonic bath. We demonstrate the role of quantum effects in the detection of water solvation and show that the computationally cheap approach based on the Caldeira-Leggett bath is only partially effective. The main conclusion of the study is that quantum methods associated with high-level potential energy surfaces are necessary to correctly study solvation features, while simplified models, even if attractive owing to their reduced computational cost, can provide some useful insights but are not able to come up with a comprehensive description of the solvation phenomenon.
Collapse
Affiliation(s)
- Alessandro Rognoni
- Dipartimento di Chimica, Universitá degli Studi di Milano, via Golgi 19, 20133 Milano, Italy
| | - Riccardo Conte
- Dipartimento di Chimica, Universitá degli Studi di Milano, via Golgi 19, 20133 Milano, Italy
| | - Michele Ceotto
- Dipartimento di Chimica, Universitá degli Studi di Milano, via Golgi 19, 20133 Milano, Italy
| |
Collapse
|
17
|
Conte R, Houston PL, Qu C, Li J, Bowman JM. Full-dimensional, ab initio potential energy surface for glycine with characterization of stationary points and zero-point energy calculations by means of diffusion Monte Carlo and semiclassical dynamics. J Chem Phys 2020; 153:244301. [DOI: 10.1063/5.0037175] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023] Open
Affiliation(s)
- Riccardo Conte
- Dipartimento di Chimica, Università degli Studi di Milano, via Golgi 19, 20133 Milano, Italy
| | - Paul L. Houston
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA and Department of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
| | - Chen Qu
- Department of Chemistry and Biochemistry, University of Maryland, College Park, Maryland 20742, USA
| | - Jeffrey Li
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| | - Joel M. Bowman
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| |
Collapse
|
18
|
Gandolfi M, Rognoni A, Aieta C, Conte R, Ceotto M. Machine learning for vibrational spectroscopy via divide-and-conquer semiclassical initial value representation molecular dynamics with application to N-methylacetamide. J Chem Phys 2020; 153:204104. [DOI: 10.1063/5.0031892] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Michele Gandolfi
- Dipartimento di Chimica, Università degli Studi di Milano, Via Golgi 19, 20133 Milano, Italy,
| | - Alessandro Rognoni
- Dipartimento di Chimica, Università degli Studi di Milano, Via Golgi 19, 20133 Milano, Italy,
| | - Chiara Aieta
- Dipartimento di Chimica, Università degli Studi di Milano, Via Golgi 19, 20133 Milano, Italy,
| | - Riccardo Conte
- Dipartimento di Chimica, Università degli Studi di Milano, Via Golgi 19, 20133 Milano, Italy,
| | - Michele Ceotto
- Dipartimento di Chimica, Università degli Studi di Milano, Via Golgi 19, 20133 Milano, Italy,
| |
Collapse
|
19
|
Li J, Zhao B, Xie D, Guo H. Advances and New Challenges to Bimolecular Reaction Dynamics Theory. J Phys Chem Lett 2020; 11:8844-8860. [PMID: 32970441 DOI: 10.1021/acs.jpclett.0c02501] [Citation(s) in RCA: 32] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
Dynamics of bimolecular reactions in the gas phase are of foundational importance in combustion, atmospheric chemistry, interstellar chemistry, and plasma chemistry. These collision-induced chemical transformations are a sensitive probe of the underlying potential energy surface(s). Despite tremendous progress in past decades, our understanding is still not complete. In this Perspective, we survey the recent advances in theoretical characterization of bimolecular reaction dynamics, stimulated by new experimental observations, and identify key new challenges.
Collapse
Affiliation(s)
- Jun Li
- School of Chemistry and Chemical Engineering & Chongqing Key Laboratory of Theoretical and Computational Chemistry, Chongqing University, Chongqing 401331, China
| | - Bin Zhao
- Theoretische Chemie, Fakultät für Chemie, Universität Bielefeld, Universitätsstraße 25, D-33615 Bielefeld, Germany
| | - Daiqian Xie
- Institute of Theoretical and Computational Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| |
Collapse
|
20
|
Gabas F, Conte R, Ceotto M. Semiclassical Vibrational Spectroscopy of Biological Molecules Using Force Fields. J Chem Theory Comput 2020; 16:3476-3485. [PMID: 32374992 PMCID: PMC7901649 DOI: 10.1021/acs.jctc.0c00127] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
Semiclassical spectroscopy is a practical way to get an accurately approximate quantum
description of spectral features starting from ab initio molecular
dynamics simulations. The computational bottleneck for the method is represented by the
cost of ab initio potential, gradient, and Hessian matrix estimates.
This drawback is particularly severe for biological systems due to their unique
complexity and large dimensionality. The main goal of this manuscript is to demonstrate
that quantum dynamics and spectroscopy, at the level of semiclassical approximation, are
doable even for sizable biological systems. To this end, we investigate the possibility
of performing semiclassical spectroscopy simulations when ab initio
calculations are replaced by computationally cheaper force field evaluations. Both
polarizable (AMOEBABIO18) and nonpolarizable (AMBER14SB) force fields are tested.
Calculations of some particular vibrational frequencies of four nucleosides, i.e.,
uridine, thymidine, deoxyguanosine, and adenosine, show that ab initio
simulations are accurate and widely applicable. Conversely, simulations based on
AMBER14SB are limited to harmonic approximations, but those relying on AMOEBABIO18 yield
acceptable semiclassical values if the investigated conformation has been included in
the force field parametrization. The main conclusion is that AMOEBABIO18 may provide a
viable route to assist semiclassical spectroscopy in the study of large biological
molecules for which an ab initio approach is not computationally
affordable.
Collapse
Affiliation(s)
- Fabio Gabas
- Dipartimento di Chimica, Università degli Studi di Milano, via Golgi 19, 20133 Milano, Italy
| | - Riccardo Conte
- Dipartimento di Chimica, Università degli Studi di Milano, via Golgi 19, 20133 Milano, Italy
| | - Michele Ceotto
- Dipartimento di Chimica, Università degli Studi di Milano, via Golgi 19, 20133 Milano, Italy
| |
Collapse
|
21
|
Cazzaniga M, Micciarelli M, Moriggi F, Mahmoud A, Gabas F, Ceotto M. Anharmonic calculations of vibrational spectra for molecular adsorbates: A divide-and-conquer semiclassical molecular dynamics approach. J Chem Phys 2020; 152:104104. [PMID: 32171221 DOI: 10.1063/1.5142682] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The vibrational spectroscopy of adsorbates is becoming an important investigation tool for catalysis and material science. This paper presents a semiclassical molecular dynamics method able to reproduce the vibrational energy levels of systems composed by molecules adsorbed on solid surfaces. Specifically, we extend our divide-and-conquer semiclassical method for power spectra calculations to gas-surface systems and interface it with plane-wave electronic structure codes. The Born-Oppenheimer classical dynamics underlying the semiclassical calculation is full dimensional, and our method includes not only the motion of the adsorbate but also those of the surface and the bulk. The vibrational spectroscopic peaks related to the adsorbate are accounted together with the most coupled phonon modes to obtain spectra amenable to physical interpretations. We apply the method to the adsorption of CO, NO, and H2O on the anatase-TiO2 (101) surface. We compare our semiclassical results with the single-point harmonic estimates and the classical power spectra obtained from the same trajectory employed in the semiclassical calculation. We find that CO and NO anharmonic effects of fundamental vibrations are similarly reproduced by the classical and semiclassical dynamics and that H2O adsorption is fully and properly described in its overtone and combination band relevant components only by the semiclassical approach.
Collapse
Affiliation(s)
- Marco Cazzaniga
- Dipartimento di Chimica, Università degli Studi di Milano, via Golgi 19, 20133 Milano, Italy
| | - Marco Micciarelli
- Dipartimento di Chimica, Università degli Studi di Milano, via Golgi 19, 20133 Milano, Italy
| | - Francesco Moriggi
- Dipartimento di Chimica, Università degli Studi di Milano, via Golgi 19, 20133 Milano, Italy
| | - Agnes Mahmoud
- Dipartimento di Chimica, Università degli Studi di Milano, via Golgi 19, 20133 Milano, Italy
| | - Fabio Gabas
- Dipartimento di Chimica, Università degli Studi di Milano, via Golgi 19, 20133 Milano, Italy
| | - Michele Ceotto
- Dipartimento di Chimica, Università degli Studi di Milano, via Golgi 19, 20133 Milano, Italy
| |
Collapse
|