1
|
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
|
2
|
Di Liberto G, Pacchioni G. Modeling Single-Atom Catalysis. ADVANCED MATERIALS (DEERFIELD BEACH, FLA.) 2023; 35:e2307150. [PMID: 37749881 DOI: 10.1002/adma.202307150] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/19/2023] [Revised: 08/17/2023] [Indexed: 09/27/2023]
Abstract
Electronic structure calculations represent an essential complement of experiments to characterize single-atom catalysts (SACs), consisting of isolated metal atoms stabilized on a support, but also to predict new catalysts. However, simulating SACs with quantum chemistry approaches is not as simple as often assumed. In this work, the essential factors that characterize a reliable simulation of SACs activity are examined. The Perspective focuses on the importance of precise atomistic characterization of the active site, since even small changes in the metal atom's surroundings can result in large changes in reactivity. The dynamical behavior and stability of SACs under working conditions, as well as the importance of adopting appropriate methods to solve the Schrödinger equation for a quantitative evaluation of reaction energies are addressed. The Perspective also focuses on the relevance of the model adopted. For electrocatalysis this must include the effects of the solvent, the presence of electrolytes, the pH, and the external potential. Finally, it is discussed how the similarities between SACs and coordination compounds may result in reaction intermediates that usually are not observed on metal electrodes. When these aspects are not adequately considered, the predictive power of electronic structure calculations is quite limited.
Collapse
Affiliation(s)
- Giovanni Di Liberto
- Dipartimento di Scienza dei Materiali, Università degli studi di Milano Bicocca, Via R. Cozzi 55, Milano, 20125, Italy
| | - Gianfranco Pacchioni
- Dipartimento di Scienza dei Materiali, Università degli studi di Milano Bicocca, Via R. Cozzi 55, Milano, 20125, Italy
| |
Collapse
|
3
|
Botti G, Ceotto M, Conte R. Investigating the Spectroscopy of the Gas Phase Guanine-Cytosine Pair: Keto versus Enol Configurations. J Phys Chem Lett 2023; 14:8940-8947. [PMID: 37768143 PMCID: PMC10577776 DOI: 10.1021/acs.jpclett.3c02073] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/26/2023] [Accepted: 09/26/2023] [Indexed: 09/29/2023]
Abstract
We report on a vibrational study of the guanine-cytosine dimer tautomers using state-of-the-art quasiclassical trajectory and semiclassical vibrational spectroscopy. The latter includes possible quantum mechanical effects. Through an accurate comparison to the experimental spectra, we are able to shine a light on the hydrogen bond network of one of the main subunits of DNA and put the experimental assignment on a solid footing. Our calculations corroborate the experimental conclusion that the global minimum Watson-and-Crick structure is not detected in the spectra, and there is no evidence of tunnel-effect-based double proton hopping. Our accurate assignment of the spectral features may also serve as a basis for the development of precise force fields to study the guanine-cytosine dimer.
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
|
4
|
Gandolfi M, Ceotto M. Molecular Dynamics of Artificially Pair-Decoupled Systems: An Accurate Tool for Investigating the Importance of Intramolecular Couplings. J Chem Theory Comput 2023; 19:6093-6108. [PMID: 37698951 PMCID: PMC10536992 DOI: 10.1021/acs.jctc.3c00553] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/25/2023] [Indexed: 09/14/2023]
Abstract
We propose a numerical technique to accurately simulate the vibrations of organic molecules in the gas phase, when pairs of atoms (or, in general, groups of degrees of freedom) are artificially decoupled, so that their motion is instantaneously decorrelated. The numerical technique we have developed is a symplectic integration algorithm that never requires computation of the force but requires estimates of the Hessian matrix. The theory we present to support our technique postulates a pair-decoupling Hamiltonian function, which parametrically depends on a decoupling coefficient α ∈ [0, 1]. The closer α is to 0, the more decoupled the selected atoms. We test the correctness of our numerical method on small molecular systems, and we apply it to study the vibrational spectroscopic features of salicylic acid at the Density Functional Theory ab initio level on a fitted potential. Our pair-decoupled simulations of salicylic acid show that decoupling hydrogen-bonded atoms do not significantly influence the frequencies of stretching modes, but enhance enormously the out-of-plane wagging and twisting motions of the hydroxyl and carboxyl groups to the point that the carboxyl and hydroxyl groups may overcome high potential energy barriers and change the salicylic acid conformation after a short simulation time. In addition, we found that the acidity of salicylic acid is more influenced by the dynamical couplings of the proton of the carboxylic group with the carbon ring than with the hydroxyl group.
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
|
5
|
Barbiero D, Bertaina G, Ceotto M, Conte R. Anharmonic Assignment of the Water Octamer Spectrum in the OH Stretch Region. J Phys Chem A 2023; 127:6213-6221. [PMID: 37477983 PMCID: PMC10405218 DOI: 10.1021/acs.jpca.3c02902] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/03/2023] [Revised: 07/03/2023] [Indexed: 07/23/2023]
Abstract
We interface the quasi-classical trajectory approach with an ab initio potential energy surface for water to assign the vibrational spectroscopical features of the OH stretch region of the water octamer cluster, which is considered to be a precursor of ice. An attempt by Li et al. to assign their recent reference experiment involved lower-level calculations based on an ad hoc scaled harmonic approach. Differently from the conclusions of this previous assignment, which invoked the contribution of 5 conformers and a solvated form of the water heptamer in the spectrum, we find out that the spectroscopic features can be related to the 4 conformers of the octamer lying lower in energy.
Collapse
Affiliation(s)
- Davide Barbiero
- Dipartimento
di Chimica, Università degli Studi
di Milano, via Golgi 19, 20133 Milano, Italy
| | - Gianluca Bertaina
- Istituto
Nazionale di Ricerca Metrologica, Strada delle Cacce 91, I-10135 Torino, 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
|
6
|
Moscato D, Gabas F, Conte R, Ceotto M. Vibrational spectroscopy simulation of solvation effects on a G-quadruplex. J Biomol Struct Dyn 2023; 41:14248-14258. [PMID: 36856120 DOI: 10.1080/07391102.2023.2180435] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/26/2022] [Accepted: 02/07/2023] [Indexed: 03/02/2023]
Abstract
It is commonly believed that solvation effects on the vibrational properties of a solute are easily accounted for by simple rules of thumbs, that is, solvating a polar molecule in a polar medium has the only effect of red shifting all its spectroscopical features and, similarly, solvating a polar molecule in a nonpolar medium has the opposite effect. In this work, we use theoretical vibrational spectroscopy at quasi-classical and quantum approximate semiclassical level to gain atomistic insights about solvent-solute interactions for 2'-deoxyguanosine and the G-quadruplex. We employ the quasi-classical trajectory method to include full anharmonicity into our calculated spectra, and then introduce quantum nuclear effects by means of divide-and-conquer semiclassical spectroscopy calculations. Solvation is treated explicitly leading to a good reproducibility of the available experimental data and reliable predictions when an experimental reference is missing.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Davide Moscato
- Dipartimento di Chimica, Università degli Studi di Milano, Milano, Italy
| | - Fabio Gabas
- Dipartimento di Chimica, Università degli Studi di Milano, Milano, Italy
| | - Riccardo Conte
- Dipartimento di Chimica, Università degli Studi di Milano, Milano, Italy
| | - Michele Ceotto
- Dipartimento di Chimica, Università degli Studi di Milano, Milano, Italy
| |
Collapse
|
7
|
Di Liberto G, Giordano L. Role of solvation model on the stability of oxygenates on Pt(111): A comparison between microsolvation, extended bilayer, and extended metal/water interface. ELECTROCHEMICAL SCIENCE ADVANCES 2023. [DOI: 10.1002/elsa.202100204] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/09/2023] Open
Affiliation(s)
| | - Livia Giordano
- Department of Materials Science University of Milano‐Bicocca Milano Italy
| |
Collapse
|
8
|
Larsson HR, Schröder M, Beckmann R, Brieuc F, Schran C, Marx D, Vendrell O. State-resolved infrared spectrum of the protonated water dimer: revisiting the characteristic proton transfer doublet peak. Chem Sci 2022; 13:11119-11125. [PMID: 36320484 PMCID: PMC9517273 DOI: 10.1039/d2sc03189b] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2022] [Accepted: 08/26/2022] [Indexed: 11/29/2023] Open
Abstract
The infrared (IR) spectra of protonated water clusters encode precise information on the dynamics and structure of the hydrated proton. However, the strong anharmonic coupling and quantum effects of these elusive species remain puzzling up to the present day. Here, we report unequivocal evidence that the interplay between the proton transfer and the water wagging motions in the protonated water dimer (Zundel ion) giving rise to the characteristic doublet peak is both more complex and more sensitive to subtle energetic changes than previously thought. In particular, hitherto overlooked low-intensity satellite peaks in the experimental spectrum are now unveiled and mechanistically assigned. Our findings rely on the comparison of IR spectra obtained using two highly accurate potential energy surfaces in conjunction with highly accurate state-resolved quantum simulations. We demonstrate that these high-accuracy simulations are important for providing definite assignments of the complex IR signals of fluxional molecules.
Collapse
Affiliation(s)
- Henrik R Larsson
- Department of Chemistry and Biochemistry, University of California Merced CA 95343 USA
- Division of Chemistry and Chemical Engineering, California Institute of Technology Pasadena CA 91125 USA
| | - Markus Schröder
- Theoretische Chemie, Physikalisch-Chemisches Institut, Universität Heidelberg Im Neuenheimer Feld 229 D - 69120 Heidelberg Germany
| | - Richard Beckmann
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum 44780 Bochum Germany
| | - Fabien Brieuc
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum 44780 Bochum Germany
| | - Christoph Schran
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum 44780 Bochum Germany
| | - Dominik Marx
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum 44780 Bochum Germany
| | - Oriol Vendrell
- Theoretische Chemie, Physikalisch-Chemisches Institut, Universität Heidelberg Im Neuenheimer Feld 229 D - 69120 Heidelberg Germany
| |
Collapse
|
9
|
Malpathak S, Church MS, Ananth N. A Semiclassical Framework for Mixed Quantum Classical Dynamics. J Phys Chem A 2022; 126:6359-6375. [PMID: 36070472 DOI: 10.1021/acs.jpca.2c03467] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Semiclassical (SC) approximations for quantum dynamic simulations in complex chemical systems range from rigorously accurate methods that are computationally expensive to methods that exhibit near-classical scaling with system size but are limited in their ability to describe quantum effects. In practical studies of high-dimensional reactions, neither extreme is the best choice: frequently a high-level quantum mechanical description is only required for a handful of modes, while the majority of environment modes that do not play a key role in the reactive event of interest are well served with a lower level of theory. In this feature, we introduce modified Filinov filtration as a powerful tool to construct mixed quantum-classical SC theories where different subsystems can be quantized to different extents without introducing ad hoc intersubsystem interaction terms. We demonstrate that these Filinov-based SC methods can systematically tune between quantum and classical limit SC behavior, offering a practical way forward to accurate and computationally efficient simulations of high-dimensional quantum processes.
Collapse
Affiliation(s)
- Shreyas Malpathak
- Department of Chemistry and Chemical Biology, Baker Laboratory, Cornell University Ithaca, New York 14853, United States
| | - Matthew S Church
- Department of Chemistry, Brown University, Providence, Rhode Island 02906, United States
| | - Nandini Ananth
- Department of Chemistry and Chemical Biology, Baker Laboratory, Cornell University Ithaca, New York 14853, United States
| |
Collapse
|
10
|
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
|
11
|
Martínez AG, Gómez PC, de la Moya S, Siehl HU. Structural proton transfer rates in pure water according to Marcus theory and TD-DFT computations. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.119048] [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]
|
12
|
Schwabe B, Niemeyer JC. Deep Zoom-In Simulation of a Fuzzy Dark Matter Galactic Halo. PHYSICAL REVIEW LETTERS 2022; 128:181301. [PMID: 35594112 DOI: 10.1103/physrevlett.128.181301] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/19/2021] [Revised: 03/04/2022] [Accepted: 04/07/2022] [Indexed: 06/15/2023]
Abstract
Fuzzy dark matter (FDM) made of ultralight bosonic particles is a viable alternative to cold dark matter with clearly distinguishable small-scale features in collapsed structures. On large scales, it behaves gravitationally like cold dark matter deviating only by a cutoff in the initial power spectrum and can be studied using N-body methods. In contrast, wave interference effects near the de Broglie scale result in new phenomena unique to FDM. Interfering modes in filaments and halos yield a stochastically oscillating granular structure which condenses into solitonic cores during halo formation. Investigating these highly nonlinear wave phenomena requires the spatially resolved numerical integration of the Schrödinger equation. In previous papers we introduced a hybrid zoom-in scheme that combines N-body methods to model the large-scale gravitational potential around and the mass accretion onto pre-selected halos with simulations of the Schrödinger-Poisson equation to capture wave-like effects inside these halos. In this work, we present a new, substantially improved reconstruction method for the wave function inside of previously collapsed structures. We demonstrate its capabilities with a deep zoom-in simulation of a well-studied sub-L_{*}-sized galactic halo from cosmological initial conditions. With a particle mass of m=2.5×10^{-22} eV and halo mass M_{vir}=1.7×10^{11} M_{⊙} in a (60 h^{-1} comoving Mpc)^{3} cosmological box, it reaches an effective resolution of 20 comoving pc. This pushes the values of m and M accessible to simulations significantly closer to those relevant for studying galaxy evolution in the allowed range of FDM masses.
Collapse
Affiliation(s)
- Bodo Schwabe
- CAPA and Departamento de Física Teórica, Universidad de Zaragoza, 50009 Zaragoza, Spain and Institut für Astrophysik, Georg-August-Universität Göttingen, D-37077 Göttingen, Germany
| | - Jens C Niemeyer
- Institut für Astrophysik, Georg-August-Universität Göttingen, D-37077 Göttingen, Germany
| |
Collapse
|
13
|
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
|
14
|
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
|
15
|
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
|
16
|
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
|
17
|
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
|
18
|
Abstract
Many efforts undertaken to study the solvation process have led to general theories that may describe mean properties, but are unable to provide a detailed understanding at the molecular level. Remarkably, the basic question of how many solvent molecules are necessary to solvate one solute molecule is still open. By exploring several water aggregates of increasing complexity, in this contribution we employ semiclassical spectroscopy to determine on quantum dynamical grounds the minimal network of surrounding water molecules to make the central one display the same vibrational features of liquid water. We find out that double-acceptor double-donor tetrahedral coordination constituting the standard picture is necessary but not sufficient, and that particular care must be reserved for the quantum description of the combination band due to the coupling of the central monomer bending mode with network librations. It is actually our ability to investigate the combination band with a quantum-derived approach that allows us to answer the titular question. The minimal structure eventually responsible for proper solvation is made of a total of 21 water molecules and includes two complete solvation shells, of which the whole first one is tetrahedrally coordinated to the central molecule. How quantum spectroscopic simulations can explain water solvation by comparison with experimental spectra.![]()
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
|
19
|
Aieta C, Bertaina G, Micciarelli M, Ceotto M. Representing molecular ground and excited vibrational eigenstates with nuclear densities obtained from semiclassical initial value representation molecular dynamics. J Chem Phys 2020; 153:214117. [PMID: 33291909 DOI: 10.1063/5.0031391] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present in detail and validate an effective Monte Carlo approach for the calculation of the nuclear vibrational densities via integration of molecular eigenfunctions that we have preliminary employed to calculate the densities of the ground and the excited OH stretch vibrational states in the protonated glycine molecule [Aieta et al., Nat Commun 11, 4348 (2020)]. Here, we first validate and discuss in detail the features of the method on a benchmark water molecule. Then, we apply it to calculate on-the-fly the ab initio anharmonic nuclear densities in the correspondence of the fundamental transitions of NH and CH stretches in protonated glycine. We show how we can gain both qualitative and quantitative physical insight by inspection of different one-nucleus densities and assign a character to spectroscopic absorption peaks using the expansion of vibrational states in terms of harmonic basis functions. The visualization of the nuclear vibrations in a purely quantum picture allows us to observe and quantify the effects of anharmonicity on the molecular structure, also to exploit the effect of IR excitations on specific bonds or functional groups, beyond the harmonic approximation. We also calculate the quantum probability distribution of bond lengths, angles, and dihedrals of the molecule. Notably, we observe how in the case of one type of fundamental NH stretching, the typical harmonic nodal pattern is absent in the anharmonic distribution.
Collapse
Affiliation(s)
- Chiara Aieta
- Dipartimento di Chimica, Università degli Studi di Milano, Via C. Golgi 19, 20133 Milano, Italy
| | - Gianluca Bertaina
- Dipartimento di Chimica, Università degli Studi di Milano, Via C. Golgi 19, 20133 Milano, Italy
| | - Marco Micciarelli
- Dipartimento di Chimica, Università degli Studi di Milano, Via C. Golgi 19, 20133 Milano, Italy
| | - Michele Ceotto
- Dipartimento di Chimica, Università degli Studi di Milano, Via C. Golgi 19, 20133 Milano, Italy
| |
Collapse
|
20
|
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
|
21
|
Aieta C, Micciarelli M, Bertaina G, Ceotto M. Anharmonic quantum nuclear densities from full dimensional vibrational eigenfunctions with application to protonated glycine. Nat Commun 2020; 11:4348. [PMID: 32859910 PMCID: PMC7455743 DOI: 10.1038/s41467-020-18211-3] [Citation(s) in RCA: 23] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/23/2020] [Accepted: 07/29/2020] [Indexed: 11/24/2022] Open
Abstract
The interpretation of molecular vibrational spectroscopic signals in terms of atomic motion is essential to understand molecular mechanisms and for chemical characterization. The signals are usually assigned after harmonic normal mode analysis, even if molecular vibrations are known to be anharmonic. Here we obtain the quantum anharmonic vibrational eigenfunctions of the 11-atom protonated glycine molecule and we calculate the density distribution of its nuclei and its geometry parameters, for both the ground and the O-H stretch excited states, using our semiclassical method based on ab initio molecular dynamics trajectories. Our quantum mechanical results describe a molecule elongated and more flexible with respect to what previously thought. More importantly, our method is able to assign each spectral peak in vibrational spectroscopy by showing quantitatively how normal modes involving different functional groups cooperate to originate that spectroscopic signal. The method will possibly allow for a better rationalization of experimental spectroscopy.
Collapse
Affiliation(s)
- Chiara Aieta
- Dipartimento di Chimica, Università degli Studi di Milano, via C. Golgi 19, 20133, Milano, Italy
| | - Marco Micciarelli
- Dipartimento di Chimica, Università degli Studi di Milano, via C. Golgi 19, 20133, Milano, Italy
| | - Gianluca Bertaina
- Dipartimento di Chimica, Università degli Studi di Milano, via C. Golgi 19, 20133, Milano, Italy
- Istituto Nazionale di Ricerca Metrologica, Strada delle Cacce 91, 10135, Torino, Italy
| | - Michele Ceotto
- Dipartimento di Chimica, Università degli Studi di Milano, via C. Golgi 19, 20133, Milano, Italy.
| |
Collapse
|
22
|
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
|
23
|
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
|
24
|
Conte R, Parma L, Aieta C, Rognoni A, Ceotto M. Improved semiclassical dynamics through adiabatic switching trajectory sampling. J Chem Phys 2019; 151:214107. [PMID: 31822104 DOI: 10.1063/1.5133144] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We introduce an improved semiclassical dynamics approach to quantum vibrational spectroscopy. In this method, a harmonic-based phase space sampling is preliminarily driven toward non-harmonic quantization by slowly switching on the actual potential. The new coordinates and momenta serve as initial conditions for the semiclassical dynamics calculation, leading to a substantial decrease in the number of chaotic trajectories to deal with. Applications are presented for model and molecular systems of increasing dimensionality characterized by moderate or high chaoticity. They include a bidimensional Henon-Heiles potential, water, formaldehyde, and methane. The method improves accuracy and precision of semiclassical results and it can be easily interfaced with all pre-existing semiclassical theories.
Collapse
Affiliation(s)
- Riccardo Conte
- Dipartimento di Chimica, Università degli Studi di Milano, Via Golgi 19, 20133 Milano, Italy
| | - Lorenzo Parma
- 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
| | - Alessandro Rognoni
- 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
|