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
|
Domenichini G, Dellago C. Molecular Hessian matrices from a machine learning random forest regression algorithm. J Chem Phys 2023; 159:194111. [PMID: 37982481 DOI: 10.1063/5.0169384] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/25/2023] [Accepted: 10/27/2023] [Indexed: 11/21/2023] Open
Abstract
In this article, we present a machine learning model to obtain fast and accurate estimates of the molecular Hessian matrix. In this model, based on a random forest, the second derivatives of the energy with respect to redundant internal coordinates are learned individually. The internal coordinates together with their specific representation guarantee rotational and translational invariance. The model is trained on a subset of the QM7 dataset but is shown to be applicable to larger molecules picked from the QM9 dataset. From the predicted Hessian, it is also possible to obtain reasonable estimates of the vibrational frequencies, normal modes, and zero point energies of the molecules.
Collapse
Affiliation(s)
- Giorgio Domenichini
- Faculty of Physics, University of Vienna, Kolingasse 14-16, 1090 Vienna, Austria
| | - Christoph Dellago
- Faculty of Physics, University of Vienna, Kolingasse 14-16, 1090 Vienna, Austria
| |
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
|
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
|
6
|
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
|
7
|
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
|
8
|
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
|
9
|
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
|
10
|
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
|
11
|
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
|
12
|
Galimberti DR, Sauer J. Chemically Accurate Vibrational Free Energies of Adsorption from Density Functional Theory Molecular Dynamics: Alkanes in Zeolites. J Chem Theory Comput 2021; 17:5849-5862. [PMID: 34459582 PMCID: PMC8444336 DOI: 10.1021/acs.jctc.1c00519] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
![]()
We present a methodology
to compute, at reduced computational cost,
Gibbs free energies, enthalpies, and entropies of adsorption from
molecular dynamics. We calculate vibrational partition functions from
vibrational energies, which we obtain from the vibrational density
of states by projection on the normal modes. The use of a set of well-chosen
reference structures along the trajectories accounts for the anharmonicities
of the modes. For the adsorption of methane, ethane, and propane in
the H-CHA zeolite, we limit our treatment to a set of vibrational
modes localized at the adsorption site (zeolitic OH group) and the
alkane molecule interacting with it. Only two short trajectories (1–20
ps) are required to reach convergence (<1 kJ/mol) for the thermodynamic
functions. The mean absolute deviations from the experimentally measured
values are 2.6, 2.8, and 4.7 kJ/mol for the Gibbs free energy, the
enthalpy, and the entropy term (−TΔS),
respectively. In particular, the entropy terms show a major improvement
compared to the harmonic approximation and almost reach the accuracy
of the previous use of anharmonic frequencies obtained with curvilinear
distortions of individual modes. The thermodynamic functions so obtained
follow the trend of the experimental values for methane, ethane, and
propane, and the Gibbs free energy of adsorption at experimental conditions
is correctly predicted to change from positive for methane (5.9 kJ/mol)
to negative for ethane (−4.8 kJ/mol) and propane (−7.1
kJ/mol).
Collapse
Affiliation(s)
- Daria Ruth Galimberti
- Institut für Chemie, Humboldt-Universität, Unter den Linden 6, 10117 Berlin, Germany.,Institute for Molecules and Materials, Radboud University, Heyendaalseweg 135, 6525 AJ Nijmegen, The Netherlands
| | - Joachim Sauer
- Institut für Chemie, Humboldt-Universität, Unter den Linden 6, 10117 Berlin, Germany
| |
Collapse
|
13
|
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
|
14
|
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
|
15
|
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
|
16
|
Mino L, Negri C, Santalucia R, Cerrato G, Spoto G, Martra G. Morphology, Surface Structure and Water Adsorption Properties of TiO 2 Nanoparticles: A Comparison of Different Commercial Samples. Molecules 2020; 25:molecules25204605. [PMID: 33050364 PMCID: PMC7587218 DOI: 10.3390/molecules25204605] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/10/2020] [Revised: 10/01/2020] [Accepted: 10/03/2020] [Indexed: 02/06/2023] Open
Abstract
Water is a molecule always present in the reaction environment in photocatalytic and biomedical applications of TiO2 and a better understanding of its interaction with the surface of TiO2 nanoparticles is crucial to develop materials with improved performance. In this contribution, we first studied the nature and the surface structure of the exposed facets of three commercial TiO2 samples (i.e., TiO2 P25, SX001, and PC105) by electron microscopy and IR spectroscopy of adsorbed CO. The morphological information was then correlated with the water adsorption properties, investigated at the molecular level, moving from multilayers of adsorbed H2O to the monolayer, combining medium- and near-IR spectroscopies. Finally, we assessed in a quantitative way the surface hydration state at different water equilibrium pressures by microgravimetric measurements.
Collapse
Affiliation(s)
- Lorenzo Mino
- Department of Chemistry and NIS Centre, University of Torino, via Giuria 7, 10125 Torino, Italy; (C.N.); (R.S.); (G.S.); (G.M.)
- Correspondence: (L.M.); (G.C.)
| | - Chiara Negri
- Department of Chemistry and NIS Centre, University of Torino, via Giuria 7, 10125 Torino, Italy; (C.N.); (R.S.); (G.S.); (G.M.)
- Department of Chemistry, Center for materials science and nanotechnology, University of Oslo, Sem Sælands vei 26, 0371 Oslo, Norway
| | - Rosangela Santalucia
- Department of Chemistry and NIS Centre, University of Torino, via Giuria 7, 10125 Torino, Italy; (C.N.); (R.S.); (G.S.); (G.M.)
| | - Giuseppina Cerrato
- Department of Chemistry and NIS Centre, University of Torino, via Giuria 7, 10125 Torino, Italy; (C.N.); (R.S.); (G.S.); (G.M.)
- Correspondence: (L.M.); (G.C.)
| | - Giuseppe Spoto
- Department of Chemistry and NIS Centre, University of Torino, via Giuria 7, 10125 Torino, Italy; (C.N.); (R.S.); (G.S.); (G.M.)
| | - Gianmario Martra
- Department of Chemistry and NIS Centre, University of Torino, via Giuria 7, 10125 Torino, Italy; (C.N.); (R.S.); (G.S.); (G.M.)
| |
Collapse
|
17
|
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
|