1
|
Xu R, Jiang Z, Yang Q, Bloino J, Biczysko M. Harmonic and anharmonic vibrational computations for biomolecular building blocks: Benchmarking DFT and basis sets by theoretical and experimental IR spectrum of glycine conformers. J Comput Chem 2024; 45:1846-1869. [PMID: 38682874 DOI: 10.1002/jcc.27377] [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: 02/12/2024] [Revised: 04/01/2024] [Accepted: 04/02/2024] [Indexed: 05/01/2024]
Abstract
Advanced vibrational spectroscopic experiments have reached a level of sophistication that can only be matched by numerical simulations in order to provide an unequivocal analysis, a crucial step to understand the structure-function relationship of biomolecules. While density functional theory (DFT) has become the standard method when targeting medium-size or larger systems, the problem of its reliability and accuracy are well-known and have been abundantly documented. To establish a reliable computational protocol, especially when accuracy is critical, a tailored benchmark is usually required. This is generally done over a short list of known candidates, with the basis set often fixed a priori. In this work, we present a systematic study of the performance of DFT-based hybrid and double-hybrid functionals in the prediction of vibrational energies and infrared intensities at the harmonic level and beyond, considering anharmonic effects through vibrational perturbation theory at the second order. The study is performed for the six-lowest energy glycine conformers, utilizing available "state-of-the-art" accurate theoretical and experimental data as reference. Focusing on the most intense fundamental vibrations in the mid-infrared range of glycine conformers, the role of the basis sets is also investigated considering the balance between computational cost and accuracy. Targeting larger systems, a broad range of hybrid schemes with different computational costs is also tested.
Collapse
Affiliation(s)
- Ruiqin Xu
- Department of Physics, College of Sciences, Shanghai University, Shanghai, China
| | | | - Qin Yang
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Science, Prague, Czechia
| | - Julien Bloino
- Classe di Scienze, Scuola Normale Superiore, Pisa, Italy
| | - Malgorzata Biczysko
- Department of Physics, College of Sciences, Shanghai University, Shanghai, China
| |
Collapse
|
2
|
Capone M, Parisse G, Narzi D, Guidoni L. Unravelling Mn 4Ca cluster vibrations in the S 1, S 2 and S 3 states of the Kok-Joliot cycle of photosystem II. Phys Chem Chem Phys 2024. [PMID: 39037338 DOI: 10.1039/d4cp01307g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/23/2024]
Abstract
Vibrational spectroscopy serves as a powerful tool for characterizing intermediate states within the Kok-Joliot cycle. In this study, we employ a QM/MM molecular dynamics framework to calculate the room temperature infrared absorption spectra of the S1, S2, and S3 states via the Fourier transform of the dipole time auto-correlation function. To better analyze the computational data and assign spectral peaks, we introduce an approach based on dipole-dipole correlation function of cluster moieties of the reaction center. Our analysis reveals variation in the infrared signature of the Mn4Ca cluster along the Kok-Joliot cycle, attributed to its increasing symmetry and rigidity resulting from the rising oxidation state of the Mn ions. Furthermore, we successfully assign the debated contributions in the frequency range around 600 cm-1. This computational methodology provides valuable insights for deciphering experimental infrared spectra and understanding the water oxidation process in both biological and artificial systems.
Collapse
Affiliation(s)
- Matteo Capone
- Università degli studi dell'Aquila, Dipartimento di Scienze Fisiche e Chimiche, L'Aquila, Italy.
| | - Gianluca Parisse
- Università degli studi dell'Aquila, Dipartimento di Scienze Fisiche e Chimiche, L'Aquila, Italy.
| | - Daniele Narzi
- Università degli studi dell'Aquila, Dipartimento di Scienze Fisiche e Chimiche, L'Aquila, Italy.
| | - Leonardo Guidoni
- Università degli studi dell'Aquila, Dipartimento di Scienze Fisiche e Chimiche, L'Aquila, Italy.
| |
Collapse
|
3
|
Chen Y, Pios SV, Gelin MF, Chen L. Accelerating Molecular Vibrational Spectra Simulations with a Physically Informed Deep Learning Model. J Chem Theory Comput 2024; 20:4703-4710. [PMID: 38825857 DOI: 10.1021/acs.jctc.4c00173] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/04/2024]
Abstract
In recent years, machine learning (ML) surrogate models have emerged as an indispensable tool to accelerate simulations of physical and chemical processes. However, there is still a lack of ML models that can accurately predict molecular vibrational spectra. Here, we present a highly efficient multitask ML surrogate model termed Vibrational Spectra Neural Network (VSpecNN), to accurately calculate infrared (IR) and Raman spectra based on dipole moments and polarizabilities obtained on-the-fly via ML-enhanced molecular dynamics simulations. The methodology is applied to pyrazine, a prototypical polyatomic chromophore. The VSpecNN-predicted energies are well within the chemical accuracy (1 kcal/mol), and the errors for VSpecNN-predicted forces are only half of those obtained from a popular high-performance ML model. Compared to the ab initio reference, the VSpecNN-predicted frequencies of IR and Raman spectra differ only by less than 5.87 cm-1, and the intensities of IR spectra and the depolarization ratios of Raman spectra are well reproduced. The VSpecNN model developed in this work highlights the importance of constructing highly accurate neural network potentials for predicting molecular vibrational spectra.
Collapse
Affiliation(s)
| | | | - Maxim F Gelin
- School of Science, Hangzhou Dianzi University, Hangzhou 310018, China
| | | |
Collapse
|
4
|
Zhang Y, Wang Y, Xu X, Chen Z, Yang Y. Vibrational Spectra of Highly Anharmonic Water Clusters: Molecular Dynamics and Harmonic Analysis Revisited with Constrained Nuclear-Electronic Orbital Methods. J Chem Theory Comput 2023; 19:9358-9368. [PMID: 38096546 DOI: 10.1021/acs.jctc.3c01037] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2023]
Abstract
Vibrational spectroscopy is widely used to gain insights into structural and dynamic properties of chemical, biological, and materials systems. Thus, an efficient and accurate method to simulate vibrational spectra is desired. In this paper, we justify and employ a microcanonical molecular simulation scheme to calculate the vibrational spectra of three challenging water clusters: the neutral water dimer (H4O2), the protonated water trimer (H7O3+), and the protonated water tetramer (H9O4+). We find that with the accurate description of quantum nuclear delocalization effects through the constrained nuclear-electronic orbital framework, including vibrational mode coupling effects through molecular dynamics simulations can additionally improve the vibrational spectrum calculations. In contrast, without the quantum nuclear delocalization picture, conventional ab initio molecular dynamics may even lead to less accurate results than harmonic analysis.
Collapse
Affiliation(s)
- Yuzhe Zhang
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin-Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Yiwen Wang
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin-Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Xi Xu
- Center for Advanced Materials Research, Beijing Normal University, Zhuhai 519087, China
| | - Zehua Chen
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin-Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Yang Yang
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin-Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| |
Collapse
|
5
|
Paschoal VH, Ribeiro MCC. DFT and ab initio molecular dynamics simulation study of the infrared spectrum of the protic ionic liquid 2-hydroxyethylammonium formate. Phys Chem Chem Phys 2023; 25:26475-26485. [PMID: 37753589 DOI: 10.1039/d3cp02914j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/28/2023]
Abstract
Protic ionic liquids (PILs) typically show a complex band shape in their infrared (IR) spectra in the high-frequency range due to the hydrogen stretching vibrations of functional groups forming rather strong hydrogen bonds (H-bonds). In the low-frequency range, the intermolecular stretching mode of the H-bond leaves a mark in the far-IR spectrum of PILs. In this study, the IR spectrum of the PIL 2-hydroxyethylammonium formate, [HOCH2CH2NH3][HCOO], is investigated in order to identify the different modes that contribute to the high-frequency band shape, i.e. the cation ν(NH), ν(OH), and ν(CH) modes, and the anion ν(CH) mode, as well as the intermolecular mode of the strongest H-bond in the far-IR spectrum. The assignment is validated by quantum chemistry calculations of clusters at the density functional theory (DFT) level for four ionic pairs and by ab initio molecular dynamics (AIMD) simulations of ten ionic pairs. There is good agreement between the vibrational frequencies obtained from DFT and AIMD simulations for both the high- and low-frequency ranges. Based on the calculations, the strong H-bond interaction between the cation -NH3 group and [HCOO]- gives a broad band envelope associated with the ν(NH) mode in the high-frequency range of the IR spectrum on which there are narrower peaks corresponding to the ν(OH) and ν(CH) modes. In the far-IR (FIR) spectrum, the anions' rattling motion gives a broad feature with a maximum at 160 cm-1, while the H-bond's intermolecular NH⋯O stretching mode appears as a peak at 255 cm-1.
Collapse
Affiliation(s)
- Vitor Hugo Paschoal
- Laboratório de Espectroscopia Molecular, Departamento de Química Fundamental, Instituto de Química, Universidade de São Paulo, 05513-970, São Paulo, SP, Brazil.
| | - Mauro C C Ribeiro
- Laboratório de Espectroscopia Molecular, Departamento de Química Fundamental, Instituto de Química, Universidade de São Paulo, 05513-970, São Paulo, SP, Brazil.
| |
Collapse
|
6
|
Ge F, Zhang L, Hou YF, Chen Y, Ullah A, Dral PO. Four-Dimensional-Spacetime Atomistic Artificial Intelligence Models. J Phys Chem Lett 2023; 14:7732-7743. [PMID: 37606602 DOI: 10.1021/acs.jpclett.3c01592] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/23/2023]
Abstract
We demonstrate that AI can learn atomistic systems in the four-dimensional (4D) spacetime. For this, we introduce the 4D-spacetime GICnet model, which for the given initial conditions (nuclear positions and velocities at time zero) can predict nuclear positions and velocities as a continuous function of time up to the distant future. Such models of molecules can be unrolled in the time dimension to yield long-time high-resolution molecular dynamics trajectories with high efficiency and accuracy. 4D-spacetime models can make predictions for different times in any order and do not need a stepwise evaluation of forces and integration of the equations of motions at discretized time steps, which is a major advance over traditional, cost-inefficient molecular dynamics. These models can be used to speed up dynamics, simulate vibrational spectra, and obtain deeper insight into nuclear motions, as we demonstrate for a series of organic molecules.
Collapse
Affiliation(s)
- Fuchun Ge
- State Key Laboratory of Physical Chemistry of Solid Surfaces, College of Chemistry and Chemical Engineering, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and Innovation Laboratory for Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen, Fujian 361005, China
| | - Lina Zhang
- State Key Laboratory of Physical Chemistry of Solid Surfaces, College of Chemistry and Chemical Engineering, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and Innovation Laboratory for Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen, Fujian 361005, China
| | - Yi-Fan Hou
- State Key Laboratory of Physical Chemistry of Solid Surfaces, College of Chemistry and Chemical Engineering, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and Innovation Laboratory for Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen, Fujian 361005, China
| | - Yuxinxin Chen
- State Key Laboratory of Physical Chemistry of Solid Surfaces, College of Chemistry and Chemical Engineering, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and Innovation Laboratory for Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen, Fujian 361005, China
| | - Arif Ullah
- School of Physics and Optoelectronic Engineering, Anhui University, Hefei 230601, China
| | - Pavlo O Dral
- State Key Laboratory of Physical Chemistry of Solid Surfaces, College of Chemistry and Chemical Engineering, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and Innovation Laboratory for Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen, Fujian 361005, China
| |
Collapse
|
7
|
Bowles J, Jähnigen S, Vuilleumier R, Calvo F, Clavaguéra C, Agostini F. Influence of the environment on the infrared spectrum of alanine: An effective mode analysis. J Chem Phys 2023; 158:094305. [PMID: 36889973 DOI: 10.1063/5.0135608] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/10/2023] Open
Abstract
The vibrational spectrum of the alanine amino acid was computationally determined in the infrared range 1000-2000 cm-1, under various environments encompassing the gas, hydrated, and crystalline phases, by means of classical molecular dynamics trajectories, carried out with the Atomic Multipole Optimized Energetics for Biomolecular Simulation polarizable force field. An effective mode analysis was performed, in which the spectra are optimally decomposed into different absorption bands arising from well-defined internal modes. In the gas phase, this analysis allows us to unravel the significant differences between the spectra obtained for the neutral and zwitterionic forms of alanine. In condensed phases, the method provides invaluable insight into the molecular origins of the vibrational bands and further shows that peaks with similar positions can be traced to rather different molecular motions.
Collapse
Affiliation(s)
- Jessica Bowles
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, 91405 Orsay, France
| | - Sascha Jähnigen
- Institut des Sciences Moléculaires d'Orsay (ISMO), CNRS, Université Paris-Saclay, 91405 Orsay, France
| | - Rodolphe Vuilleumier
- PASTEUR Laboratory, Département de Chimie, Ecole Normale Supérieure, PSL University, Sorbonne Université, CNRS, 75005 Paris, France
| | - Florent Calvo
- Université Grenoble Alpes, CNRS, LIPhy, F-38000 Grenoble, France
| | - Carine Clavaguéra
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, 91405 Orsay, France
| | - Federica Agostini
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, 91405 Orsay, France
| |
Collapse
|
8
|
Andersson Å, Yatsyna V, Linares M, Rijs A, Zhaunerchyk V. Indication of 3 10-Helix Structure in Gas-Phase Neutral Pentaalanine. J Phys Chem A 2023; 127:938-945. [PMID: 36669091 PMCID: PMC9900583 DOI: 10.1021/acs.jpca.2c07863] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Abstract
We investigate the gas-phase structure of the neutral pentaalanine peptide. The IR spectrum in the 340-1820 cm-1 frequency range is obtained by employing supersonic jet cooling, infrared multiphoton dissociation, and vacuum-ultraviolet action spectroscopy. Comparison with quantum chemical spectral calculations suggests that the molecule assumes multiple stable conformations, mainly of two structure types. In the most stable conformation theoretically found, the N-terminus forms a C5 ring and the backbone resembles that of an 310-helix with two β-turns. Additionally, the conformational preferences of pentaalanine have been evaluated using Born-Oppenheimer molecular dynamics, showing that a nonzero simulation time step causes a systematic frequency shift.
Collapse
Affiliation(s)
- Åke Andersson
- Department
of Physics, University of Gothenburg, 41296 Gothenburg, Sweden
| | - Vasyl Yatsyna
- Department
of Physics, University of Gothenburg, 41296 Gothenburg, Sweden,FELIX
Laboratory, Institute for Molecules and Materials, Radboud University, Toernooiveld 7, 6525 ED Nijmegen, The Netherlands,Laboratoire
de Chimie Physique Moléculaire, École
Polytechnique Fédérale de Lausanne, EPFL SB ISIC LCPM, Station 6, CH-1015 Lausanne, Switzerland
| | - Mathieu Linares
- Laboratory
of Organic Electronics and Group of Scientific Visualization Department
of Science and Technology (ITN), Linköping
University, 601 74 Norrköping, Sweden
| | - Anouk Rijs
- FELIX
Laboratory, Institute for Molecules and Materials, Radboud University, Toernooiveld 7, 6525 ED Nijmegen, The Netherlands,Division
of BioAnalytical Chemistry, AIMMS Amsterdam Institute of Molecular
and Life Sciences, Vrije Universiteit Amsterdam, De Boelelaan 1108, 1081 HV, Amsterdam, The Netherlands,
| | - Vitali Zhaunerchyk
- Department
of Physics, University of Gothenburg, 41296 Gothenburg, Sweden,
| |
Collapse
|
9
|
Lykhin AO, Truhlar DG, Gagliardi L. Dipole Moment Calculations Using Multiconfiguration Pair-Density Functional Theory and Hybrid Multiconfiguration Pair-Density Functional Theory. J Chem Theory Comput 2021; 17:7586-7601. [PMID: 34793166 DOI: 10.1021/acs.jctc.1c00915] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The dipole moment is the molecular property that most directly indicates molecular polarity. The accuracy of computed dipole moments depends strongly on the quality of the calculated electron density, and the breakdown of single-reference methods for strongly correlated systems can lead to poor predictions of the dipole moments in those cases. Here, we derive the analytical expression for obtaining the electric dipole moment by multiconfiguration pair-density functional theory (MC-PDFT), and we assess the accuracy of MC-PDFT for predicting dipole moments at equilibrium and nonequilibrium geometries. We show that MC-PDFT dipole moment curves have reasonable behavior even for stretched geometries, and they significantly improve upon the CASSCF results by capturing more electron correlation. The analysis of a dataset consisting of 18 first-row transition-metal diatomics and 6 main-group polyatomic molecules with a multireference character suggests that MC-PDFT and its hybrid extension (HMC-PDFT) perform comparably to CASPT2 and MRCISD+Q methods and have a mean unsigned deviation of 0.2-0.3 D with respect to the best available dipole moment reference values. We explored the dependence of the predicted dipole moments upon the choice of the on-top density functional and active space, and we recommend the tPBE and hybrid tPBE0 on-top choices for the functionals combined with the moderate correlated-participating-orbitals scheme for selecting the active space. With these choices, the mean unsigned deviations (in debyes) of the calculated equilibrium dipole moments from the best estimates are 0.77 for CASSCF, 0.29 for MC-PDFT, 0.24 for HMC-PDFT, 0.28 for CASPT2, and 0.25 for MRCISD+Q. These results are encouraging because the computational cost of MC-PDFT or HMC-PDFT is largely reduced compared to the CASPT2 and MRCISD+Q methods.
Collapse
Affiliation(s)
- Aleksandr O Lykhin
- Department of Chemistry, Pritzker School of Molecular Engineering, The James Franck Institute and Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, United States
| | - Laura Gagliardi
- Department of Chemistry, Pritzker School of Molecular Engineering, The James Franck Institute and Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, Illinois 60637, United States.,Argonne National Laboratory, Lemont, Illinois 60439, United States
| |
Collapse
|
10
|
Gaigeot MP. Some opinions on MD-based vibrational spectroscopy of gas phase molecules and their assembly: An overview of what has been achieved and where to go. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2021; 260:119864. [PMID: 34052762 DOI: 10.1016/j.saa.2021.119864] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/14/2021] [Revised: 04/13/2021] [Accepted: 04/18/2021] [Indexed: 06/12/2023]
Abstract
We hereby review molecular dynamics simulations for anharmonic gas phase spectroscopy and provide some of our opinions of where the field is heading. With these new directions, the theoretical IR/Raman spectroscopy of large (bio)-molecular systems will be more easily achievable over longer time-scale MD trajectories for an increase in accuracy of the MD-IR and MD-Raman calculated spectra. With the new directions presented here, the high throughput 'decoding' of experimental IR/Raman spectra into 3D-structures should thus be possible, hence advancing e.g. the field of MS-IR for structural characterization by spectroscopy. We also review the assignment of vibrational spectra in terms of anharmonic molecular modes from the MD trajectories, and especially introduce our recent developments based on Graph Theory algorithms. Graph Theory algorithmic is also introduced in this review for the identification of the molecular 3D-structures sampled over MD trajectories.
Collapse
Affiliation(s)
- Marie-Pierre Gaigeot
- Université Paris-Saclay, Univ Evry, CNRS, LAMBE UMR8587, 91025 Evry-Courcouronnes, France.
| |
Collapse
|
11
|
Barone V, Alessandrini S, Biczysko M, Cheeseman JR, Clary DC, McCoy AB, DiRisio RJ, Neese F, Melosso M, Puzzarini C. Computational molecular spectroscopy. ACTA ACUST UNITED AC 2021. [DOI: 10.1038/s43586-021-00034-1] [Citation(s) in RCA: 23] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
|
12
|
Gerber RB. My Trajectory in Molecular Reaction Dynamics and Spectroscopy. Annu Rev Phys Chem 2021; 72:1-34. [PMID: 33276702 DOI: 10.1146/annurev-physchem-090519-124238] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
This is the story of a career in theoretical chemistry during a time of dramatic changes in the field due to phenomenal growth in the availability of computational power. It is likewise the story of the highly gifted graduate students and postdoctoral fellows that I was fortunate to mentor throughout my career. It includes reminiscences of the great mentors that I had and of the exciting collaborations with both experimentalists and theorists on which I built much of my research. This is an account of the developments of exciting scientific disciplines in which I was involved: vibrational spectroscopy, molecular reaction mechanisms and dynamics, e.g., in atmospheric chemistry, and the prediction of new, exotic molecules, in particular noble gas molecules. From my very first project to my current work, my career in science has brought me the excitement and fascination of research. What a wonderful pursuit!
Collapse
Affiliation(s)
- Robert Benny Gerber
- The Fritz Haber Research Center and Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem 91904, Israel; .,Department of Chemistry, University of California, Irvine, California 92697, USA
| |
Collapse
|
13
|
Hassan I, Ferraro F, Imhof P. Effect of the Hydration Shell on the Carbonyl Vibration in the Ala-Leu-Ala-Leu Peptide. Molecules 2021; 26:2148. [PMID: 33917998 PMCID: PMC8068333 DOI: 10.3390/molecules26082148] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/01/2021] [Revised: 03/30/2021] [Accepted: 03/30/2021] [Indexed: 11/16/2022] Open
Abstract
The vibrational spectrum of the Ala-Leu-Ala-Leu peptide in solution, computed from first-principles simulations, shows a prominent band in the amide I region that is assigned to stretching of carbonyl groups. Close inspection reveals combined but slightly different contributions by the three carbonyl groups of the peptide. The shift in their exact vibrational signature is in agreement with the different probabilities of these groups to form hydrogen bonds with the solvent. The central carbonyl group has a hydrogen bond probability intermediate to the other two groups due to interchanges between different hydrogen-bonded states. Analysis of the interaction energies of individual water molecules with that group shows that shifts in its frequency are directly related to the interactions with the water molecules in the first hydration shell. The interaction strength is well correlated with the hydrogen bond distance and hydrogen bond angle, though there is no perfect match, allowing geometrical criteria for hydrogen bonds to be used as long as the sampling is sufficient to consider averages. The hydrogen bond state of a carbonyl group can therefore serve as an indicator of the solvent's effect on the vibrational frequency.
Collapse
Affiliation(s)
- Irtaza Hassan
- Institute for Theoretical Physics, Freie Universtiät Berlin, Arnimallee 14, 14195 Berlin, Germany;
| | - Federica Ferraro
- Computer Chemistry Center, Friedrich-Alexander University (FAU) Erlangen-Nürnberg, Nägelsbachstrasse 25, 91052 Erlangen, Germany;
| | - Petra Imhof
- Institute for Theoretical Physics, Freie Universtiät Berlin, Arnimallee 14, 14195 Berlin, Germany;
- Computer Chemistry Center, Friedrich-Alexander University (FAU) Erlangen-Nürnberg, Nägelsbachstrasse 25, 91052 Erlangen, Germany;
| |
Collapse
|
14
|
Roy TK, Gerber RB. Dual Basis Approach for Ab Initio Anharmonic Calculations of Vibrational Spectroscopy: Application to Microsolvated Biomolecules. J Chem Theory Comput 2020; 16:7005-7016. [PMID: 32991804 DOI: 10.1021/acs.jctc.0c00725] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
A dual electronic basis set approach is introduced for more efficient but accurate calculations of the anharmonic vibrational spectra in the framework of the vibrational self-consistent field (VSCF) theory. In this approach, an accurate basis set is used to compute the vibrational spectra at the harmonic level. The results are used to scale the potential surface from a more modest but much more efficient basis set. The scaling is such that at the harmonic level the new, scaled potential agrees with one of the accurate basis sets. The approach is tested in the application of the microsolvated, protected amino acid Ac-Phe-OMe, using the scaled anharmonic hybrid potential in the VSCF and VSCF-PT2 algorithms. The hybrid potential method yields results that are in good accord with the experiment and very close to those obtained in calculations with the high-level, very costly potential from the large basis set. At the same time, the hybrid potential calculations are considerably less expensive. The results of the hybrid calculations are much more accurate than those computed from the potential surface corresponding to the modest basis set. The results are very encouraging for using the hybrid potential method for inexpensive yet sufficiently accurate anharmonic calculations for the spectra of large biomolecules.
Collapse
Affiliation(s)
- Tapta Kanchan Roy
- Department of Chemistry and Chemical Sciences, Central University of Jammu, Rahya-Suchani (Bagla), Dist: Samba Jammu-181143, India
| | - R Benny Gerber
- Institute of Chemistry, The Hebrew University, Jerusalem 91904, Israel.,Department of Chemistry, University of California, Irvine, California 92697, United States
| |
Collapse
|
15
|
Abstract
Recently, AIMD (ab initio molecular dynamics) has been extensively employed to explore the dynamical information of electronic systems. However, it remains extremely challenging to reliably predict the properties of nanosystems with a radical nature using conventional electronic structure methods (e.g., Kohn-Sham density functional theory) due to the presence of static correlation. To address this challenge, we combine the recently formulated TAO-DFT (thermally-assisted-occupation density functional theory) with AIMD. The resulting TAO-AIMD method is employed to investigate the instantaneous/average radical nature and infrared spectra of n-acenes containing n linearly fused benzene rings (n = 2-8) at 300 K. According to the TAO-AIMD simulations, on average, the smaller n-acenes (up to n = 5) possess a nonradical nature, and the larger n-acenes (n = 6-8) possess an increasing radical nature, showing remarkable similarities to the ground-state counterparts at 0 K. Besides, the infrared spectra of n-acenes obtained with the TAO-AIMD simulations are in qualitative agreement with the existing experimental data.
Collapse
Affiliation(s)
- Shaozhi Li
- Department of Physics, National Taiwan University, Taipei, Taiwan
| | - Jeng-Da Chai
- Department of Physics, National Taiwan University, Taipei, Taiwan
- Center for Theoretical Physics and Center for Quantum Science and Engineering, National Taiwan University, Taipei, Taiwan
| |
Collapse
|
16
|
Prentice JCA, Aarons J, Womack JC, Allen AEA, Andrinopoulos L, Anton L, Bell RA, Bhandari A, Bramley GA, Charlton RJ, Clements RJ, Cole DJ, Constantinescu G, Corsetti F, Dubois SMM, Duff KKB, Escartín JM, Greco A, Hill Q, Lee LP, Linscott E, O'Regan DD, Phipps MJS, Ratcliff LE, Serrano ÁR, Tait EW, Teobaldi G, Vitale V, Yeung N, Zuehlsdorff TJ, Dziedzic J, Haynes PD, Hine NDM, Mostofi AA, Payne MC, Skylaris CK. The ONETEP linear-scaling density functional theory program. J Chem Phys 2020; 152:174111. [PMID: 32384832 DOI: 10.1063/5.0004445] [Citation(s) in RCA: 51] [Impact Index Per Article: 12.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Abstract
We present an overview of the onetep program for linear-scaling density functional theory (DFT) calculations with large basis set (plane-wave) accuracy on parallel computers. The DFT energy is computed from the density matrix, which is constructed from spatially localized orbitals we call Non-orthogonal Generalized Wannier Functions (NGWFs), expressed in terms of periodic sinc (psinc) functions. During the calculation, both the density matrix and the NGWFs are optimized with localization constraints. By taking advantage of localization, onetep is able to perform calculations including thousands of atoms with computational effort, which scales linearly with the number or atoms. The code has a large and diverse range of capabilities, explored in this paper, including different boundary conditions, various exchange-correlation functionals (with and without exact exchange), finite electronic temperature methods for metallic systems, methods for strongly correlated systems, molecular dynamics, vibrational calculations, time-dependent DFT, electronic transport, core loss spectroscopy, implicit solvation, quantum mechanical (QM)/molecular mechanical and QM-in-QM embedding, density of states calculations, distributed multipole analysis, and methods for partitioning charges and interactions between fragments. Calculations with onetep provide unique insights into large and complex systems that require an accurate atomic-level description, ranging from biomolecular to chemical, to materials, and to physical problems, as we show with a small selection of illustrative examples. onetep has always aimed to be at the cutting edge of method and software developments, and it serves as a platform for developing new methods of electronic structure simulation. We therefore conclude by describing some of the challenges and directions for its future developments and applications.
Collapse
Affiliation(s)
- Joseph C A Prentice
- Department of Materials, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - Jolyon Aarons
- Department of Physics, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, United Kingdom
| | - James C Womack
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Alice E A Allen
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - Lampros Andrinopoulos
- Department of Physics, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - Lucian Anton
- UKAEA, Culham Science Centre, Abingdon OX14 3DB, United Kingdom
| | - Robert A Bell
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - Arihant Bhandari
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Gabriel A Bramley
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Robert J Charlton
- Department of Materials, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - Rebecca J Clements
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Daniel J Cole
- School of Natural and Environmental Sciences, Newcastle University, Newcastle upon Tyne NE1 7RU, United Kingdom
| | - Gabriel Constantinescu
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - Fabiano Corsetti
- Department of Materials, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - Simon M-M Dubois
- Institute of Condensed Matter and Nanosciences, Université Catholique de Louvain, Louvain-la-Neuve, Belgium
| | - Kevin K B Duff
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - José María Escartín
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - Andrea Greco
- Department of Physics, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - Quintin Hill
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Louis P Lee
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - Edward Linscott
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - David D O'Regan
- School of Physics, AMBER, and CRANN Institute, Trinity College Dublin, The University of Dublin, Dublin 2, Ireland
| | - Maximillian J S Phipps
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Laura E Ratcliff
- Department of Materials, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - Álvaro Ruiz Serrano
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Edward W Tait
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - Gilberto Teobaldi
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Valerio Vitale
- Department of Materials, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - Nelson Yeung
- Department of Physics, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, United Kingdom
| | - Tim J Zuehlsdorff
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA
| | - Jacek Dziedzic
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Peter D Haynes
- Department of Materials, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - Nicholas D M Hine
- Department of Physics, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, United Kingdom
| | - Arash A Mostofi
- Department of Materials, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - Mike C Payne
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - Chris-Kriton Skylaris
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| |
Collapse
|
17
|
Arroyave JM, Ambrusi RE, Pronsato ME, Juan A, Pistonesi MF, Centurión ME. Experimental and DFT Studies of Hybrid Silver/Cdots Nanoparticles. J Phys Chem B 2020; 124:2425-2435. [PMID: 32134662 DOI: 10.1021/acs.jpcb.9b10430] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
In this work, the combination of experimental and theoretical results was employed to confirm an interaction between Cdots and AgNPs in the silver/Cdots hybrid nanoparticles. The experimental data obtained by UV-vis, IR, ζ potential, and TGA techniques were correlated and interpreted by calculations obtained by DFT. In particular, an interaction between the -COO- functional group of the Cdots with AgNPs was revealed. As consequence of this interaction, a frequency shift and a higher absorption intensity in the IR of the -OH group in the Cdots was theoretically predicted and also observed in the experimental IR spectra. Moreover, a bonding and charge distribution analysis was also carried out. These results constitute new physical insight for the Ag@Cdots system. Additionally, based in this type of interaction, energy calculations explained the negative charge surrounding the AgNPs, which was detected by ζ potential measurements. This systematic methodology not only is useful for this nanoparticles system but also could be used to analyze the interaction between the components that constitute other types of hybrid nanoparticles.
Collapse
Affiliation(s)
- Jeison Manuel Arroyave
- INQUISUR, Department of Chemistry, Universidad Nacional del Sur, Avenida Alem 1253, B8000CPB, Bahía Blanca, Argentina
| | - Rubén E Ambrusi
- IFISUR, Department of Physics, Universidad Nacional del Sur, Avenida Alem 1253, B8000CPB, Bahía Blanca, Argentina
| | - María Estela Pronsato
- IFISUR, Department of Physics, Universidad Nacional del Sur, Avenida Alem 1253, B8000CPB, Bahía Blanca, Argentina
| | - Alfredo Juan
- IFISUR, Department of Physics, Universidad Nacional del Sur, Avenida Alem 1253, B8000CPB, Bahía Blanca, Argentina
| | - Marcelo Fabian Pistonesi
- INQUISUR, Department of Chemistry, Universidad Nacional del Sur, Avenida Alem 1253, B8000CPB, Bahía Blanca, Argentina
| | - María Eugenia Centurión
- INQUISUR, Department of Chemistry, Universidad Nacional del Sur, Avenida Alem 1253, B8000CPB, Bahía Blanca, Argentina
| |
Collapse
|
18
|
Dvores MP, Çarçabal P, Maître P, Simons JP, Gerber RB. Gas phase dynamics, conformational transitions and spectroscopy of charged saccharides: the oxocarbenium ion, protonated anhydrogalactose and protonated methyl galactopyranoside. Phys Chem Chem Phys 2020; 22:4144-4157. [PMID: 32039431 DOI: 10.1039/c9cp06572e] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
Protonated intermediates are postulated to be involved in the rate determining step of many sugar reactions. This paper presents a study of protonated sugar species, isolated in the gas phase, using a combination of infrared multiple photon dissociation (IRMPD) spectroscopy, classical ab initio molecular dynamics (AIMD) and quantum mechanical vibrational self-consistent field (VSCF) calculations. It provides a likely identification of the reactive intermediate oxocarbenium ion structure in a d-galactosyl system as well as the saccharide pyrolysis product anhydrogalactose (that suggests oxocarbenium ion stabilization), along with the spectrum of the protonated parent species: methyl d-galactopyranoside-H+. Its vibrational fingerprint indicates intramolecular proton sharing. Classical AIMD simulations for galactosyl oxocarbenium ions, conducted in the temperature range ∼300-350 K (using B3LYP potentials on-the-fly) reveal efficient transitions on the picosecond timescale. Multiple conformers are likely to exist under the experimental conditions and along with static VSCF calculations, they have facilitated the identification of the individual structural motifs of the galactosyl oxocarbenium ion and protonated anhydrogalactose ion conformers that contribute to the observed experimental spectra. These results demonstrate the power of experimental IRMPD spectroscopy combined with dynamics simulations and with computational spectroscopy at the anharmonic level to unravel conformer structures of protonated saccharides, and to provide information on their lifetimes.
Collapse
Affiliation(s)
- M P Dvores
- Institute of Chemistry and the Fritz Haber Research Center, The Hebrew University, Jerusalem 91904, Israel.
| | - P Çarçabal
- Université Paris-Saclay, CNRS, Institut des Sciences Moléculaires d'Orsay, 91405, Orsay, France
| | - P Maître
- Université Paris-Saclay, CNRS, Institut de Chimie Physique, 91405, Orsay, France
| | - J P Simons
- Department of Chemistry, Physical and Theoretical Chemistry Laboratory, Oxford University, South Parks Road, Oxford OX1 3QZ, UK
| | - R B Gerber
- Institute of Chemistry and the Fritz Haber Research Center, The Hebrew University, Jerusalem 91904, Israel. and Department of Chemistry, University of California Irvine, CA 92697, USA
| |
Collapse
|
19
|
Bakels S, Gaigeot MP, Rijs AM. Gas-Phase Infrared Spectroscopy of Neutral Peptides: Insights from the Far-IR and THz Domain. Chem Rev 2020; 120:3233-3260. [PMID: 32073261 PMCID: PMC7146864 DOI: 10.1021/acs.chemrev.9b00547] [Citation(s) in RCA: 34] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
![]()
Gas-phase, double
resonance IR spectroscopy has proven to be an
excellent approach to obtain structural information on peptides ranging
from single amino acids to large peptides and peptide clusters. In
this review, we discuss the state-of-the-art of infrared action spectroscopy
of peptides in the far-IR and THz regime. An introduction to the field
of far-IR spectroscopy is given, thereby highlighting the opportunities
that are provided for gas-phase research on neutral peptides. Current
experimental methods, including spectroscopic schemes, have been reviewed.
Structural information from the experimental far-IR spectra can be
obtained with the help of suitable theoretical approaches such as
dynamical DFT techniques and the recently developed Graph Theory.
The aim of this review is to underline how the synergy between far-IR
spectroscopy and theory can provide an unprecedented picture of the
structure of neutral biomolecules in the gas phase. The far-IR signatures
of the discussed studies are summarized in a far-IR map, in order
to gain insight into the origin of the far-IR localized and delocalized
motions present in peptides and where they can be found in the electromagnetic
spectrum.
Collapse
Affiliation(s)
- Sjors Bakels
- Radboud University, Institute for Molecules and Materials, FELIX Laboratory, Toernooiveld 7-c, 6525 ED Nijmegen, The Netherlands
| | - Marie-Pierre Gaigeot
- LAMBE CNRS UMR8587, Université d'Evry val d'Essonne, Blvd F. Mitterrand, Bât Maupertuis, 91025 Evry, France
| | - Anouk M Rijs
- Radboud University, Institute for Molecules and Materials, FELIX Laboratory, Toernooiveld 7-c, 6525 ED Nijmegen, The Netherlands
| |
Collapse
|
20
|
Djouama R, Hamza Reguig F, Bounaceur B, Dauchez M, Krallafa AM. Structure, Dynamics, and Photophysical Properties of a Series of [Pt(NH 2R) 4]–[PtX 4] Complexes. J Phys Chem A 2020; 124:911-923. [DOI: 10.1021/acs.jpca.9b10300] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Rabie Djouama
- LCPM Laboratory, Faculty of Sciences, Chemistry Department, University of Oran 1 Ahmed BenBella, 31000 Oran, Algeria
| | - Farouk Hamza Reguig
- LCPM Laboratory, Faculty of Sciences, Chemistry Department, University of Oran 1 Ahmed BenBella, 31000 Oran, Algeria
| | - Boumediene Bounaceur
- LCPM Laboratory, Faculty of Sciences, Chemistry Department, University of Oran 1 Ahmed BenBella, 31000 Oran, Algeria
| | - Manuel Dauchez
- SirMa CNRS UMR 7369, MEDyC, University of Reims Champagne Ardenne, Reims 51097, France
| | - Abdelghani Mohamed Krallafa
- LCPM Laboratory, Faculty of Sciences, Chemistry Department, University of Oran 1 Ahmed BenBella, 31000 Oran, Algeria
- SirMa CNRS UMR 7369, MEDyC, University of Reims Champagne Ardenne, Reims 51097, France
| |
Collapse
|
21
|
Samala NR, Agmon N. Temperature Dependence of Intramolecular Vibrational Bands in Small Water Clusters. J Phys Chem B 2019; 123:9428-9442. [PMID: 31553613 DOI: 10.1021/acs.jpcb.9b07777] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
Cyclic water clusters are pivotal for understanding atmospheric reactions as well as liquid water, yet the temperature (T) dependence of their dynamics and spectroscopy is poorly studied. The development of highly accurate water potentials, such as MB-pol, partly rectifies this. It remains to account for the quantum nuclear effects (NQE), because quantum nuclear dynamics become increasingly inaccurate at low temperatures. From a practical point of view, we find that NQE can be accounted for simply by subtracting a constant from the frequencies obtained from the velocity autocorrelation functions (VACF) of classical nuclear dynamics, resulting in unprecedented agreement with experiment, mostly within 5 cm-1. We have performed classical simulations of (H2O)n clusters (n = 2-5) from 20 K and up to their melting temperature, calculating both all-atom and partial VACF, thus generating the temperature dependence of the vibrational frequencies (IR and Raman bands). Focusing on the hydrogen-bonded (HBed) OH stretch and HOH bend, we find opposing T dependencies. The HBed OH modes blue shift linearly with T, attributed to ring expansion rather than any specific conformational change. The lowest-frequency Raman concerted mode is predicted to show the largest such shift. In contrast, the HOH bend undergoes a red-shift, with the highest frequency concerted band undergoing the largest red-shift. These results can be explained by a coupled-oscillator model for n hydrogen atoms on a ring, constrained to move either tangentially (stretch) or perpendicularly (bend) to the ring. With increasing temperature and weakening of HBs, the intrinsic force constant increases (stretch) or remains constant (bend), while the nearest-neighbor coupling constant decreases, and this results in the interesting behavior revealed herein. T-dependent Raman studies are required for testing some of these predictions.
Collapse
Affiliation(s)
- Nagaprasad Reddy Samala
- The Fritz Haber Research Center, Institute of Chemistry , The Hebrew University of Jerusalem , Jerusalem 91904 , Israel
| | - Noam Agmon
- The Fritz Haber Research Center, Institute of Chemistry , The Hebrew University of Jerusalem , Jerusalem 91904 , Israel
| |
Collapse
|
22
|
Tychengulova A, Capone M, Pitari F, Guidoni L. Molecular Vibrations of an Oxygen-Evolving Complex and Its Synthetic Mimic. Chemistry 2019; 25:13385-13395. [PMID: 31340068 DOI: 10.1002/chem.201902621] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2019] [Revised: 07/18/2019] [Indexed: 11/07/2022]
Abstract
Bio-inspired catalysis for artificial photosynthesis has been widely studied for decades, in particular, with the purpose of using bio-disposable and non-toxic metals as building blocks. The characterisation of such catalysts has been achieved by using different kinds of spectroscopic methods, from X-ray crystallography to NMR spectroscopy. An artificial Mn4 CaO4 cubane cluster with dangling Mn4 was synthesised in 2015 [Zhang et al. Science 2015, 348, 690-693]; this cluster showed many structural similarities to that of the natural oxygen-evolving complex. An accurate structural and spectroscopic comparison between the natural and artificial systems is highly relevant to understand the catalytic mechanism. Among data from different techniques, the differential FTIR spectra (Sn+1 -Sn ) of photosystem II are still lacking a complete interpretation. The availability of IR data of the artificial cluster offers a unique opportunity to assign absolute absorption spectra on a well-defined and easier to interpret analogous moiety. The present work aims to investigate the novel inorganic compound as a model system for an oxygen-evolving complex through measurement of its spectroscopic properties. The experimental results are compared with calculations by using a variety of theoretical methods (normal mode analysis, effective normal mode analysis) in the S1 state. We underline the similarities and the differences in the computational spectra based on atomistic models of Mn4 CaO5 and Mn4 CaO4 complexes.
Collapse
Affiliation(s)
- Aliya Tychengulova
- Department of Basic and Applied Sciences for Engineering, Sapienza University of Rome, Via Scarpa 16, 00161, Rome, Italy
| | - Mateo Capone
- Department of Engineering, Computer Science and Mathematics, University of L'Aquila, Via Vetoio Coppito, 67100, L'Aquila, Italy
| | - Fabio Pitari
- Department of Engineering, Computer Science and Mathematics, University of L'Aquila, Via Vetoio Coppito, 67100, L'Aquila, Italy
- Current address: CINECA High Performance Computing Department, Via Magnanelli, 40033, Casalecchio di Reno, Italy
| | - Leonardo Guidoni
- Department of Physical and Chemical Science, University of L'Aquila, Via Vetoio Coppito, 67100, L'Aquila, Italy
| |
Collapse
|
23
|
Conte R, Gabas F, Botti G, Zhuang Y, Ceotto M. Semiclassical vibrational spectroscopy with Hessian databases. J Chem Phys 2019; 150:244118. [PMID: 31255076 DOI: 10.1063/1.5109086] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Abstract
We report on a new approach to ease the computational overhead of ab initio "on-the-fly" semiclassical dynamics simulations for vibrational spectroscopy. The well known bottleneck of such computations lies in the necessity to estimate the Hessian matrix for propagating the semiclassical pre-exponential factor at each step along the dynamics. The procedure proposed here is based on the creation of a dynamical database of Hessians and associated molecular geometries able to speed up calculations while preserving the accuracy of results at a satisfactory level. This new approach can be interfaced to both analytical potential energy surfaces and on-the-fly dynamics, allowing one to study even large systems previously not achievable. We present results obtained for semiclassical vibrational power spectra of methane, glycine, and N-acetyl-L-phenylalaninyl-L-methionine-amide, a molecule of biological interest made of 46 atoms.
Collapse
Affiliation(s)
- Riccardo Conte
- 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
| | - Giacomo Botti
- Dipartimento di Chimica, Università degli Studi di Milano, via Golgi 19, 20133 Milano, Italy
| | - Yu Zhuang
- Department of Computer Science, Texas Tech University, Lubbock, Texas 79409-3104, USA
| | - Michele Ceotto
- Dipartimento di Chimica, Università degli Studi di Milano, via Golgi 19, 20133 Milano, Italy
| |
Collapse
|
24
|
Micciarelli M, Gabas F, Conte R, Ceotto M. An effective semiclassical approach to IR spectroscopy. J Chem Phys 2019; 150:184113. [DOI: 10.1063/1.5096968] [Citation(s) in RCA: 28] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Marco Micciarelli
- Dipartimento di Chimica, Università degli Studi di Milano, Via C. Golgi 19, 20133 Milano, Italy
| | - Fabio Gabas
- Dipartimento di Chimica, Università degli Studi di Milano, Via C. Golgi 19, 20133 Milano, Italy
| | - Riccardo Conte
- 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
|
25
|
Calvo F. Conformational diversity in deprotonated water clusters and anharmonic infrared spectra. MOLECULAR SIMULATION 2019. [DOI: 10.1080/08927022.2018.1513653] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
Affiliation(s)
- F. Calvo
- Université Grenoble Alpes, CNRS, LIPhy, Grenoble, France
| |
Collapse
|
26
|
Panuszko A, Nowak M, Bruździak P, Stasiulewicz M, Stangret J. Amides as models to study the hydration of proteins and peptides — spectroscopic and theoretical approach on hydration in various temperatures. J Mol Liq 2019. [DOI: 10.1016/j.molliq.2019.01.086] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
|
27
|
Imoto S, Marx D. Pressure response of the THz spectrum of bulk liquid water revealed by intermolecular instantaneous normal mode analysis. J Chem Phys 2019; 150:084502. [PMID: 30823759 DOI: 10.1063/1.5080381] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The radial distribution functions of liquid water are known to change significantly their shape upon hydrostatic compression from ambient conditions deep into the kbar pressure regime. It has been shown that despite their eye-catching changes, the fundamental locally tetrahedral fourfold H-bonding pattern that characterizes ambient water is preserved up to about 10 kbar (1 GPa), which is the stability limit of liquid water at 300 K. The observed increase in coordination number comes from pushing water molecules into the first coordination sphere without establishing an H-bond, resulting in roughly two such additional interstitial molecules at 10 kbar. THz spectroscopy has been firmly established as a powerful experimental technique to analyze H-bonding in aqueous solutions given that it directly probes the far-infrared lineshape and thus the prominent H-bond network mode around 180 cm-1. We, therefore, set out to assess pressure effects on the THz response of liquid water at 10 kbar in comparison to the 1 bar (0.1 MPa) reference, both at 300 K, with the aim to trace back the related lineshape changes to the structural level. To this end, we employ the instantaneous normal mode approximation to rigorously separate the H-bonding peak from the large background arising from the pronounced librational tail. By exactly decomposing the total molecular dynamics into hindered translations, hindered rotations, and intramolecular vibrations, we find that the H-bonding peak arises from translation-translation and translation-rotation correlations, which are successively decomposed down to the level of distinct local H-bond environments. Our utmost detailed analysis based on molecular pair classifications unveils that H-bonded double-donor water pairs contribute most to the THz response around 180 cm-1, whereas interstitial waters are negligible. Moreover, short double-donor H-bonds have their peak maximum significantly shifted toward higher frequencies with respect to such long H-bonds. In conjunction with an increasing relative population of these short H-bonds versus the long ones (while the population of other water pair classes is essentially pressure insensitive), this explains not only the blue-shift of the H-bonding peak by about 20-30 cm-1 in total from 1 bar to 10 kbar but also the filling of the shallow local minimum of the THz lineshape located in between the network peak and the red-wing of the librational band at 1 bar. Based on the changing populations as a function of pressure, we are also able to roughly estimate the pressure-dependence of the H-bond network mode and find that its pressure response and thus the blue-shifting are most pronounced at low kbar pressures.
Collapse
Affiliation(s)
- Sho Imoto
- 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
| |
Collapse
|
28
|
Gupta PK, Esser A, Forbert H, Marx D. Toward theoretical terahertz spectroscopy of glassy aqueous solutions: partially frozen solute-solvent couplings of glycine in water. Phys Chem Chem Phys 2019; 21:4975-4987. [PMID: 30758388 DOI: 10.1039/c8cp07489e] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The molecular-level understanding of THz spectra of aqueous solutions under ambient conditions has been greatly advanced in recent years. Here, we go beyond previous analyses by performing ab initio molecular dynamics simulations of glycine in water with artificially frozen solute or solvent molecules, respectively, while computing the total THz response as well as its decomposition into mode-specific resonances based on the "supermolecular solvation complex" technique. Clamping the water molecules and keeping glycine moving breaks the coupling of glycine to the structural dynamics of the solvent, however, the polarization and dielectric solvation effects in the static solvation cage are still at work since the full electronic structure of the quenched solvent is taken into account. The complementary approach of fixing glycine reveals both the dynamical and electronic response of the solvation cage at the level of its THz response. Moreover, to quantitatively account for the electronic contribution solely due to solvent embedding, the solute species is "vertically desolvated", thus preserving the fully coupled solute-solvent motion in terms of the solute's structural dynamics in solution, while its electronic structure is no longer subject to solute-solvent polarization and charge transfer effects. When referenced to the free simulation of Gly(aq), this three-fold approach allows us to decompose the THz spectral contributions due to the correlated solute-solvent dynamics into entirely structural and purely electronic effects. Beyond providing hitherto unknown insights, the observed systematic changes of THz spectra in terms of peak shifts and lineshape modulations due to conformational freezing and frozen solvation cages might be useful to investigate the solvation of molecules in highly viscous H-bonding solvents such as ionic liquids and even in cryogenic ices as relevant to polar stratospheric and dark interstellar clouds.
Collapse
Affiliation(s)
- Prashant Kumar Gupta
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany
| | | | | | | |
Collapse
|
29
|
Holroyd LF, Bühl M, Gaigeot MP, van Mourik T. Thermodynamics of 5-Bromouracil Tautomerization From First-Principles Molecular Dynamics Simulations. ADVANCES IN QUANTUM CHEMISTRY 2019. [DOI: 10.1016/bs.aiq.2018.05.001] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
|
30
|
Galimberti DR, Bougueroua S, Mahé J, Tommasini M, Rijs AM, Gaigeot MP. Conformational assignment of gas phase peptides and their H-bonded complexes using far-IR/THz: IR-UV ion dip experiment, DFT-MD spectroscopy, and graph theory for mode assignment. Faraday Discuss 2019; 217:67-97. [DOI: 10.1039/c8fd00211h] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Graph theory based vibrational modes as new entities for vibrational THz spectroscopy.
Collapse
Affiliation(s)
| | | | - Jérôme Mahé
- LAMBE UMR8587
- Univ Evry
- Université Paris-Saclay
- CNRS
- 91025 Evry
| | - Matteo Tommasini
- Department of Chemistry, Materials, Chemical Engineering “G. Natta” Politecnico di Milano
- 20133 Milano
- Italy
| | - Anouk M. Rijs
- Radboud University
- Institute for Molecules and Materials
- FELIX Laboratory
- 6525 ED Nijmegen
- The Netherlands
| | | |
Collapse
|
31
|
Roy TK, Kopysov V, Pereverzev A, Šebek J, Gerber RB, Boyarkin OV. Intrinsic structure of pentapeptide Leu-enkephalin: geometry optimization and validation by comparison of VSCF-PT2 calculations with cold ion spectroscopy. Phys Chem Chem Phys 2018; 20:24894-24901. [PMID: 30234204 DOI: 10.1039/c8cp03989e] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
The intrinsic structure of an opioid peptide [Ala2, Leu5]-leucine enkephalin (ALE) has been investigated using first-principles based vibrational self-consistent field (VSCF) theory and cold ion spectroscopy. IR-UV double resonance spectroscopy revealed the presence of only one highly abundant conformer of the singly protonated ALE, isolated and cryogenically cooled in the gas phase. High-level quantum mechanical calculations of electronic structures in conjunction with a systematic conformational search allowed for finding a few low-energy candidate structures. In order to identify the observed structure, we computed vibrational spectra of the candidate structures and employed the theory at the semi-empirically scaled harmonic level and at the first-principles based anharmonic VSCF levels. The best match between the calculated "anharmonic" and the measured spectra appeared, indeed, for the most stable candidate. An average of two spectra calculated with different quantum mechanical potentials is proposed for the best match with experiment. The match thus validates the calculated intrinsic structure of ALE and demonstrates the predictive power of first-principles theory for solving structures of such large molecules.
Collapse
Affiliation(s)
- Tapta Kanchan Roy
- Department of Chemistry and Chemical Sciences, Central University of Jammu, Rahya-Suchani (Bagla), Jammu 181143, India
| | | | | | | | | | | |
Collapse
|
32
|
Anharmonic vibrational spectroscopy calculations using the ab initio CSP method: Applications to H2CO3, (H2CO3)2, H2CO3-H2O and isotopologues. Chem Phys 2018. [DOI: 10.1016/j.chemphys.2017.12.015] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
33
|
Gabas F, Di Liberto G, Conte R, Ceotto M. Protonated glycine supramolecular systems: the need for quantum dynamics. Chem Sci 2018; 9:7894-7901. [PMID: 30542548 PMCID: PMC6237109 DOI: 10.1039/c8sc03041c] [Citation(s) in RCA: 26] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/09/2018] [Accepted: 09/16/2018] [Indexed: 02/02/2023] Open
Abstract
Quantum mechanical simulations unequivocally explain experimental IR spectra of protonated supramolecular systems.
IR spectroscopy is one of the most commonly employed techniques to study molecular vibrations and interactions. However, characterization of experimental IR spectra is not always straightforward. This is the case of protonated glycine supramolecular systems like Gly2H+ and (GlyH + nH2), whose IR spectra raise questions which have still to find definitive answers even after theoretical spectroscopy investigations. Specifically, the assignment of the conformer responsible for the spectrum of the protonated glycine dimer (Gly2H+) has led to much controversy and it is still debated, while structural hypotheses formulated to explain the main experimental spectral features of (GlyH + nH2) systems have not been theoretically confirmed. We demonstrate that simulations must account for quantum dynamical effects in order to resolve these open issues. This is achieved by means of our divide-and-conquer semiclassical initial value representation technique, which approximates the quantum dynamics of high dimensional systems with remarkable accuracy and outperforms not only the commonly employed but unfit scaled-harmonic approaches, but also pure classical dynamics simulations. Besides the specific insights concerning the two particular cases presented here, the general conclusion is that, due to the widespread presence of protonated systems in chemistry, quantum dynamics may play a prominent role and should not be totally overlooked even when dealing with large systems including biological structures.
Collapse
Affiliation(s)
- Fabio Gabas
- Dipartimento di Chimica , Università degli Studi di Milano , via Golgi 19 , 20133 Milano , Italy . ;
| | - Giovanni Di Liberto
- 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
|
34
|
Drużbicki K, Krzystyniak M, Hollas D, Kapil V, Slavíček P, Romanelli G, Fernandez-Alonso F. Hydrogen dynamics in solid formic acid: insights from simulations with quantum colored-noise thermostats. ACTA ACUST UNITED AC 2018. [DOI: 10.1088/1742-6596/1055/1/012003] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
|
35
|
Di Liberto G, Conte R, Ceotto M. "Divide and conquer" semiclassical molecular dynamics: A practical method for spectroscopic calculations of high dimensional molecular systems. J Chem Phys 2018; 148:014307. [PMID: 29306274 DOI: 10.1063/1.5010388] [Citation(s) in RCA: 30] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023] Open
Abstract
We extensively describe our recently established "divide-and-conquer" semiclassical method [M. Ceotto, G. Di Liberto, and R. Conte, Phys. Rev. Lett. 119, 010401 (2017)] and propose a new implementation of it to increase the accuracy of results. The technique permits us to perform spectroscopic calculations of high-dimensional systems by dividing the full-dimensional problem into a set of smaller dimensional ones. The partition procedure, originally based on a dynamical analysis of the Hessian matrix, is here more rigorously achieved through a hierarchical subspace-separation criterion based on Liouville's theorem. Comparisons of calculated vibrational frequencies to exact quantum ones for a set of molecules including benzene show that the new implementation performs better than the original one and that, on average, the loss in accuracy with respect to full-dimensional semiclassical calculations is reduced to only 10 wavenumbers. Furthermore, by investigating the challenging Zundel cation, we also demonstrate that the "divide-and-conquer" approach allows us to deal with complex strongly anharmonic molecular systems. Overall the method very much helps the assignment and physical interpretation of experimental IR spectra by providing accurate vibrational fundamentals and overtones decomposed into reduced dimensionality spectra.
Collapse
Affiliation(s)
- Giovanni Di Liberto
- Dipartimento di Chimica, Università degli Studi di Milano, Via C. Golgi 19, 20133 Milano, Italy
| | - Riccardo Conte
- 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
|
36
|
Frederix PWJM, Patmanidis I, Marrink SJ. Molecular simulations of self-assembling bio-inspired supramolecular systems and their connection to experiments. Chem Soc Rev 2018; 47:3470-3489. [PMID: 29688238 PMCID: PMC5961611 DOI: 10.1039/c8cs00040a] [Citation(s) in RCA: 92] [Impact Index Per Article: 15.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/13/2018] [Indexed: 01/01/2023]
Abstract
In bionanotechnology, the field of creating functional materials consisting of bio-inspired molecules, the function and shape of a nanostructure only appear through the assembly of many small molecules together. The large number of building blocks required to define a nanostructure combined with the many degrees of freedom in packing small molecules has long precluded molecular simulations, but recent advances in computational hardware as well as software have made classical simulations available to this strongly expanding field. Here, we review the state of the art in simulations of self-assembling bio-inspired supramolecular systems. We will first discuss progress in force fields, simulation protocols and enhanced sampling techniques using recent examples. Secondly, we will focus on efforts to enable the comparison of experimentally accessible observables and computational results. Experimental quantities that can be measured by microscopy, spectroscopy and scattering can be linked to simulation output either directly or indirectly, via quantum mechanical or semi-empirical techniques. Overall, we aim to provide an overview of the various computational approaches to understand not only the molecular architecture of nanostructures, but also the mechanism of their formation.
Collapse
Affiliation(s)
- Pim W. J. M. Frederix
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Materials , University of Groningen , Groningen , The Netherlands . ;
| | - Ilias Patmanidis
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Materials , University of Groningen , Groningen , The Netherlands . ;
| | - Siewert J. Marrink
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Materials , University of Groningen , Groningen , The Netherlands . ;
| |
Collapse
|
37
|
Hassan I, Donati L, Stensitzki T, Keller BG, Heyne K, Imhof P. The vibrational spectrum of the hydrated alanine-leucine peptide in the amide region from IR experiments and first principles calculations. Chem Phys Lett 2018. [DOI: 10.1016/j.cplett.2018.03.026] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
|
38
|
Levina E, Penkov NV, Rodionova NN, Tarasov SA, Barykina DV, Vener MV. Hydration of the Carboxylate Group in Anti-Inflammatory Drugs: ATR-IR and Computational Studies of Aqueous Solution of Sodium Diclofenac. ACS OMEGA 2018; 3:302-313. [PMID: 30023777 PMCID: PMC6044930 DOI: 10.1021/acsomega.7b01034] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/20/2017] [Accepted: 12/27/2017] [Indexed: 05/12/2023]
Abstract
Diclofenac (active ingredient of Voltaren) has a significant, multifaceted role in medicine, pharmacy, and biochemistry. Its physical properties and impact on biomolecular structures still attract essential scientific interest. However, its interaction with water has not been described yet at the molecular level. In the present study, we shed light on the interaction between the steric hindrance (the intramolecular N-H···O bond, etc.) carboxylate group (-CO2-) with water. Aqueous solution of sodium declofenac is investigated using attenuated total reflection-infrared (ATR-IR) and computational approaches, i.e., classical molecular dynamics (MD) simulations and density functional theory (DFT). Our coupled classical MD simulations, DFT calculations, and ATR-IR spectroscopy results indicated that the -CO2- group of the diclofenac anion undergoes strong specific interactions with the water molecules. The combined experimental and theoretical techniques provide significant insights into the spectroscopic manifestation of these interactions and the structure of the hydration shell of the -CO2- group. Moreover, the developed methodology for the theoretical analysis of the ATR-IR spectrum could serve as a template for the future IR/Raman studies of the strong interaction between the steric hindrance -CO2- group of bioactive molecules with the water molecules in dilute aqueous solutions.
Collapse
Affiliation(s)
- Elena
O. Levina
- Department
of Molecular and Chemical Physics, Moscow
Institute of Physics and Technology, 7 Institutskiy per., 141700 Dolgoprudny, Russia
| | - Nikita V. Penkov
- Department
of Methods of Optical and Spectral Analysis, Institute of Cell Biophysics, Russian Academy of Sciences, 3 Institutskaya Street, 142292 Pushchino, Russia
| | - Natalia N. Rodionova
- OOO
“NPF” Materia Medica Holding, 47-1 Trifonovskaya Street, 129272 Moscow, Russia
| | - Sergey A. Tarasov
- OOO
“NPF” Materia Medica Holding, 47-1 Trifonovskaya Street, 129272 Moscow, Russia
| | - Daria V. Barykina
- OOO
“NPF” Materia Medica Holding, 47-1 Trifonovskaya Street, 129272 Moscow, Russia
| | - Mikhail V. Vener
- Department
of Quantum Chemistry, Mendeleev University
of Chemical Technology, 9 Miusskaya Square, 125047 Moscow, Russia
- E-mail:
| |
Collapse
|
39
|
Imoto S, Forbert H, Marx D. Aqueous TMAO solutions as seen by theoretical THz spectroscopy: hydrophilic versus hydrophobic water. Phys Chem Chem Phys 2018; 20:6146-6158. [DOI: 10.1039/c7cp07003a] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
All THz resonances of aqueous TMAO solutions are computed and assigned based on ab initio molecular dynamics simulations.
Collapse
Affiliation(s)
- Sho Imoto
- Lehrstuhl für Theoretische Chemie
- Ruhr-Universität Bochum
- 44780 Bochum
- Germany
| | - Harald Forbert
- Center for Solvation Science ZEMOS
- Ruhr-Universität Bochum
- 44780 Bochum
- Germany
| | - Dominik Marx
- Lehrstuhl für Theoretische Chemie
- Ruhr-Universität Bochum
- 44780 Bochum
- Germany
| |
Collapse
|
40
|
Gastegger M, Behler J, Marquetand P. Machine learning molecular dynamics for the simulation of infrared spectra. Chem Sci 2017; 8:6924-6935. [PMID: 29147518 PMCID: PMC5636952 DOI: 10.1039/c7sc02267k] [Citation(s) in RCA: 255] [Impact Index Per Article: 36.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2017] [Accepted: 08/08/2017] [Indexed: 12/28/2022] Open
Abstract
Machine learning has emerged as an invaluable tool in many research areas. In the present work, we harness this power to predict highly accurate molecular infrared spectra with unprecedented computational efficiency. To account for vibrational anharmonic and dynamical effects - typically neglected by conventional quantum chemistry approaches - we base our machine learning strategy on ab initio molecular dynamics simulations. While these simulations are usually extremely time consuming even for small molecules, we overcome these limitations by leveraging the power of a variety of machine learning techniques, not only accelerating simulations by several orders of magnitude, but also greatly extending the size of systems that can be treated. To this end, we develop a molecular dipole moment model based on environment dependent neural network charges and combine it with the neural network potential approach of Behler and Parrinello. Contrary to the prevalent big data philosophy, we are able to obtain very accurate machine learning models for the prediction of infrared spectra based on only a few hundreds of electronic structure reference points. This is made possible through the use of molecular forces during neural network potential training and the introduction of a fully automated sampling scheme. We demonstrate the power of our machine learning approach by applying it to model the infrared spectra of a methanol molecule, n-alkanes containing up to 200 atoms and the protonated alanine tripeptide, which at the same time represents the first application of machine learning techniques to simulate the dynamics of a peptide. In all of these case studies we find an excellent agreement between the infrared spectra predicted via machine learning models and the respective theoretical and experimental spectra.
Collapse
Affiliation(s)
- Michael Gastegger
- University of Vienna , Faculty of Chemistry , Institute of Theoretical Chemistry , Währinger Str. 17 , 1090 Vienna , Austria . ; ; Tel: +43 1 4277 52764
| | - Jörg Behler
- Universität Göttingen , Institut für Physikalische Chemie , Theoretische Chemie , Tammannstr. 6 , 37077 Göttingen , Germany
| | - Philipp Marquetand
- University of Vienna , Faculty of Chemistry , Institute of Theoretical Chemistry , Währinger Str. 17 , 1090 Vienna , Austria . ; ; Tel: +43 1 4277 52764
| |
Collapse
|
41
|
Alata I, Pérez-Mellor A, Ben Nasr F, Scuderi D, Steinmetz V, Gobert F, Jaïdane NE, Zehnacker-Rentien A. Does the Residues Chirality Modify the Conformation of a Cyclo-Dipeptide? Vibrational Spectroscopy of Protonated Cyclo-diphenylalanine in the Gas Phase. J Phys Chem A 2017; 121:7130-7138. [DOI: 10.1021/acs.jpca.7b06159] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Ivan Alata
- Institut
des Sciences Moléculaires d’Orsay, CNRS, Univ. Paris-Sud, Université Paris-Saclay, F-91405 Orsay, France
| | - Ariel Pérez-Mellor
- Institut
des Sciences Moléculaires d’Orsay, CNRS, Univ. Paris-Sud, Université Paris-Saclay, F-91405 Orsay, France
| | - Feriel Ben Nasr
- Institut
des Sciences Moléculaires d’Orsay, CNRS, Univ. Paris-Sud, Université Paris-Saclay, F-91405 Orsay, France
- Laboratoire
de Spectroscopie Atomique Moléculaire et Applications, Université de Tunis El Manar, Tunis 1060, Tunisia
| | - Debora Scuderi
- Laboratoire
de Chimie Physique, CNRS, UMR8000, Univ. Paris-Sud, Orsay F-91405, France
| | - Vincent Steinmetz
- Laboratoire
de Chimie Physique, CNRS, UMR8000, Univ. Paris-Sud, Orsay F-91405, France
| | - Fabrice Gobert
- Laboratoire
de Chimie Physique, CNRS, UMR8000, Univ. Paris-Sud, Orsay F-91405, France
| | - Nejm-Eddine Jaïdane
- Laboratoire
de Spectroscopie Atomique Moléculaire et Applications, Université de Tunis El Manar, Tunis 1060, Tunisia
| | - Anne Zehnacker-Rentien
- Institut
des Sciences Moléculaires d’Orsay, CNRS, Univ. Paris-Sud, Université Paris-Saclay, F-91405 Orsay, France
| |
Collapse
|
42
|
Galimberti DR, Milani A, Gaigeot MP, Radice S, Tonelli C, Picozzi R, Castiglioni C. Static vs dynamic DFT prediction of IR spectra of flexible molecules in the condensed phase: The (ClCF 2CF(CF 3)OCF 2CH 3) liquid as a test case. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2017; 183:195-203. [PMID: 28448957 DOI: 10.1016/j.saa.2017.04.013] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/10/2016] [Revised: 04/07/2017] [Accepted: 04/14/2017] [Indexed: 06/07/2023]
Abstract
First-principles molecular dynamics (FPMD) simulations in the framework of Density Functional Theory (DFT) are carried out for the prediction of the infrared spectrum of the fluorinated molecule ClCF2CF(CF3)OCF2CH3 in liquid and gas phase. This molecule is characterized by a flexible structure, allowing the co-existence of several stable conformers, that differ by values of the torsional angles. FPMD computed spectra are compared to the experimental ones, and to Boltzmann weighted IR spectra based on gas phase calculations.
Collapse
Affiliation(s)
- Daria Ruth Galimberti
- Politecnico di Milano - Dip. Chimica, Materiali, Ing. Chimica "G. Natta", Piazza Leonardo da Vinci 32, 20133 Milan, Italy.
| | - Alberto Milani
- Politecnico di Milano - Dip. Chimica, Materiali, Ing. Chimica "G. Natta", Piazza Leonardo da Vinci 32, 20133 Milan, Italy
| | - Marie-Pierre Gaigeot
- LAMBE CNRS UMR8587, Laboratoire Analyse et Modélisation pour la Biologie et l'Environnement, Université d'Evry val d'Essonne, Boulevard F. Mitterrand, Bat Maupertuis, 91025 Evry, France
| | - Stefano Radice
- Solvay Specialty Polymers RD&T Center, Viale Lombardia 20, 20021 Bollate, MI, Italy
| | - Claudio Tonelli
- Solvay Specialty Polymers RD&T Center, Viale Lombardia 20, 20021 Bollate, MI, Italy
| | - Rosaldo Picozzi
- Solvay Specialty Polymers RD&T Center, Viale Lombardia 20, 20021 Bollate, MI, Italy
| | - Chiara Castiglioni
- Politecnico di Milano - Dip. Chimica, Materiali, Ing. Chimica "G. Natta", Piazza Leonardo da Vinci 32, 20133 Milan, Italy
| |
Collapse
|
43
|
Galimberti DR, Milani A, Tommasini M, Castiglioni C, Gaigeot MP. Combining Static and Dynamical Approaches for Infrared Spectra Calculations of Gas Phase Molecules and Clusters. J Chem Theory Comput 2017; 13:3802-3813. [DOI: 10.1021/acs.jctc.7b00471] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Affiliation(s)
- Daria R. Galimberti
- Dip.
Chimica, Materiali, Ing. Chimica “G. Natta”, Politecnico di Milano, 20133 Milan, Italy
- LAMBE
CNRS UMR8587, Université d’Evry val d’Essonne, 91025 Evry, France
- Université Paris-Saclay, 91190 Saint-Aubin, France
| | - Alberto Milani
- Dip.
Chimica, Materiali, Ing. Chimica “G. Natta”, Politecnico di Milano, 20133 Milan, Italy
| | - Matteo Tommasini
- Dip.
Chimica, Materiali, Ing. Chimica “G. Natta”, Politecnico di Milano, 20133 Milan, Italy
| | - Chiara Castiglioni
- Dip.
Chimica, Materiali, Ing. Chimica “G. Natta”, Politecnico di Milano, 20133 Milan, Italy
| | - Marie-Pierre Gaigeot
- LAMBE
CNRS UMR8587, Université d’Evry val d’Essonne, 91025 Evry, France
- Université Paris-Saclay, 91190 Saint-Aubin, France
| |
Collapse
|
44
|
Wang J. Ultrafast two-dimensional infrared spectroscopy for molecular structures and dynamics with expanding wavelength range and increasing sensitivities: from experimental and computational perspectives. INT REV PHYS CHEM 2017. [DOI: 10.1080/0144235x.2017.1321856] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Affiliation(s)
- Jianping Wang
- Beijing National Laboratory for Molecular Sciences, Molecular Reaction Dynamics Laboratory, CAS Research/Education Center for Excellence in Molecular Sciences, Institute of Chemistry, The Chinese Academy of Sciences, Beijing, P.R. China
- College of Chemistry, University of Chinese Academy of Sciences, Beijing, P.R. China
| |
Collapse
|
45
|
Schienbein P, Schwaab G, Forbert H, Havenith M, Marx D. Correlations in the Solute-Solvent Dynamics Reach Beyond the First Hydration Shell of Ions. J Phys Chem Lett 2017; 8:2373-2380. [PMID: 28488865 DOI: 10.1021/acs.jpclett.7b00713] [Citation(s) in RCA: 47] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
While the real-space structure of solvation shells has been explored for decades, a dynamical perspective that directly relies on changes in the H-bond network became accessible more recently mainly via far-infrared (THz) spectroscopies. A remaining key question is how many hydration shells are affected by ion-induced network perturbations. We disclose that theoretical THz difference spectra of aqueous salt solutions can be deciphered in terms of only a handful of dipolar auto- and cross-correlations, including the second solvation shell. This emphasizes the importance of cross-correlations being often neglected in multicomponent models. Analogously, experimental THz responses of simple ions can be deciphered in a similar way. Dramatic intensity cancellations due to large positive and negative contributions are found to effectively shift intensity maxima. Thus, THz spectroscopy provides an unprecedented view on the details of hydration dynamics, which can be understood by a combination of experiment and theory.
Collapse
Affiliation(s)
- Philipp Schienbein
- Lehrstuhl für Theoretische Chemie, ‡Lehrstuhl für Physikalische Chemie II, and §Center for Solvation Science ZEMOS, Ruhr-Universität Bochum , 44780 Bochum, Germany
| | - Gerhard Schwaab
- Lehrstuhl für Theoretische Chemie, ‡Lehrstuhl für Physikalische Chemie II, and §Center for Solvation Science ZEMOS, Ruhr-Universität Bochum , 44780 Bochum, Germany
| | - Harald Forbert
- Lehrstuhl für Theoretische Chemie, ‡Lehrstuhl für Physikalische Chemie II, and §Center for Solvation Science ZEMOS, Ruhr-Universität Bochum , 44780 Bochum, Germany
| | - Martina Havenith
- Lehrstuhl für Theoretische Chemie, ‡Lehrstuhl für Physikalische Chemie II, and §Center for Solvation Science ZEMOS, Ruhr-Universität Bochum , 44780 Bochum, Germany
| | - Dominik Marx
- Lehrstuhl für Theoretische Chemie, ‡Lehrstuhl für Physikalische Chemie II, and §Center for Solvation Science ZEMOS, Ruhr-Universität Bochum , 44780 Bochum, Germany
| |
Collapse
|
46
|
Krzystyniak M, Drużbicki K, Romanelli G, Gutmann MJ, Rudić S, Imberti S, Fernandez-Alonso F. Nuclear dynamics and phase polymorphism in solid formic acid. Phys Chem Chem Phys 2017; 19:9064-9074. [PMID: 28304035 DOI: 10.1039/c7cp00997f] [Citation(s) in RCA: 29] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
We apply a unique sequence of structural and dynamical neutron-scattering techniques, augmented with density-functional electronic-structure calculations, to establish the degree of polymorphism in an archetypal hydrogen-bonded system - crystalline formic acid. Using this combination of experimental and theoretical techniques, the hypothesis by Zelsmann on the coexistence of the β1 and β2 phases above 220 K is tested. Contrary to the postulated scenario of proton-transfer-driven phase coexistence, the emerging picture is one of a quantitatively different structural change over this temperature range, whereby the loosening of crystal packing promotes temperature-induced shearing of the hydrogen-bonded chains. The presented work, therefore, solves a fifty-year-old puzzle and provides a suitable framework for the use neutron-Compton-scattering techniques in the exploration of phase polymorphism in condensed matter.
Collapse
Affiliation(s)
- Maciej Krzystyniak
- ISIS Facility, Rutherford Appleton Laboratory, Chilton, Didcot, Oxfordshire OX11 0QX, UK. and School of Science and Technology, Nottingham Trent University, Clifton Campus, Nottingham, NG11 8NS, UK
| | - Kacper Drużbicki
- Faculty of Physics, Adam Mickiewicz University, Umultowska St. 85, 61-614 Poznan, Poland
| | - Giovanni Romanelli
- ISIS Facility, Rutherford Appleton Laboratory, Chilton, Didcot, Oxfordshire OX11 0QX, UK.
| | - Matthias J Gutmann
- ISIS Facility, Rutherford Appleton Laboratory, Chilton, Didcot, Oxfordshire OX11 0QX, UK.
| | - Svemir Rudić
- ISIS Facility, Rutherford Appleton Laboratory, Chilton, Didcot, Oxfordshire OX11 0QX, UK.
| | - Silvia Imberti
- ISIS Facility, Rutherford Appleton Laboratory, Chilton, Didcot, Oxfordshire OX11 0QX, UK.
| | - Felix Fernandez-Alonso
- ISIS Facility, Rutherford Appleton Laboratory, Chilton, Didcot, Oxfordshire OX11 0QX, UK. and Department of Physics and Astronomy, University College London, Gower Street, London, WC1E 6BT, UK
| |
Collapse
|
47
|
Yatsyna V, Bakker DJ, Feifel R, Rijs AM, Zhaunerchyk V. Far-infrared amide IV-VI spectroscopy of isolated 2- and 4-Methylacetanilide. J Chem Phys 2017; 145:104309. [PMID: 27634262 DOI: 10.1063/1.4962360] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
Abstract
Delocalized molecular vibrations in the far-infrared and THz ranges are highly sensitive to the molecular structure, as well as to intra- and inter-molecular interactions. Thus, spectroscopic studies of biomolecular structures can greatly benefit from an extension of the conventional mid-infrared to the far-infrared wavelength range. In this work, the conformer-specific gas-phase far-infrared spectra of two aromatic molecules containing the peptide -CO-NH- link, namely, 2- and 4-Methylacetanilide, are investigated. The planar conformations with trans configuration of the peptide link have only been observed in the supersonic-jet expansion. The corresponding far-infrared signatures associated with the vibrations of the peptide -CO-NH- moiety, the so-called amide IV-VI bands, have been assigned and compared with the results of density functional theory frequency calculations based on the anharmonic vibrational second-order perturbation theory approach. The analysis of the experimental and theoretical data shows that the amide IV-VI bands are highly diagnostic for the geometry of the peptide moiety and the molecular backbone. They are also strongly blue-shifted upon formation of the NH⋯O-C hydrogen bonding, which is, for example, responsible for the formation of secondary protein structures. Furthermore, the amide IV-VI bands are also diagnostic for the cis configuration of the peptide link, which can be present in cyclic peptides. The experimental gas-phase data presented in this work can assist the vibrational assignment of similar biologically important systems, either isolated or in natural environments.
Collapse
Affiliation(s)
- Vasyl Yatsyna
- Department of Physics, University of Gothenburg, 412 96 Gothenburg, Sweden
| | - Daniël J Bakker
- Institute for Molecules and Materials, FELIX Laboratory, Radboud University, Toernoovield 7-c, 6525 ED Nijmegen, The Netherlands
| | - Raimund Feifel
- Department of Physics, University of Gothenburg, 412 96 Gothenburg, Sweden
| | - Anouk M Rijs
- Institute for Molecules and Materials, FELIX Laboratory, Radboud University, Toernoovield 7-c, 6525 ED Nijmegen, The Netherlands
| | - Vitali Zhaunerchyk
- Department of Physics, University of Gothenburg, 412 96 Gothenburg, Sweden
| |
Collapse
|
48
|
Mahé J, Bakker DJ, Jaeqx S, Rijs AM, Gaigeot MP. Mapping gas phase dipeptide motions in the far-infrared and terahertz domain. Phys Chem Chem Phys 2017; 19:13778-13787. [DOI: 10.1039/c7cp00369b] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
Vibrational signatures of Ac-Phe-AA-NH2 dipeptides are recorded and analysed in the far IR/THz spectral domain (100–800 cm−1, 3–24 THz), with the ‘AA’ amino acid chosen within the series ‘AA’ = Gly, Ala, Pro, Cys, Ser, Val. Phe stands for phenylalanine.
Collapse
Affiliation(s)
- Jérôme Mahé
- LAMBE CNRS UMR8587
- Université d'Evry val d'Essonne
- 91025 Evry
- France
- Université Paris-Saclay
| | - Daniël J. Bakker
- Radboud University
- Institute for Molecules and Materials
- FELIX Laboratory
- 6525 ED Nijmegen
- The Netherlands
| | - Sander Jaeqx
- Radboud University
- Institute for Molecules and Materials
- FELIX Laboratory
- 6525 ED Nijmegen
- The Netherlands
| | - Anouk M. Rijs
- Radboud University
- Institute for Molecules and Materials
- FELIX Laboratory
- 6525 ED Nijmegen
- The Netherlands
| | - Marie-Pierre Gaigeot
- LAMBE CNRS UMR8587
- Université d'Evry val d'Essonne
- 91025 Evry
- France
- Université Paris-Saclay
| |
Collapse
|
49
|
Lanza G, Chiacchio MA. Effects of Hydration on the Zwitterion Trialanine Conformation by Electronic Structure Theory. J Phys Chem B 2016; 120:11705-11719. [DOI: 10.1021/acs.jpcb.6b08108] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Affiliation(s)
- Giuseppe Lanza
- Dipartimento
di Scienze del
Farmaco, Università di Catania, Viale A. Doria 6, 95125 Catania, Italy
| | - Maria A. Chiacchio
- Dipartimento
di Scienze del
Farmaco, Università di Catania, Viale A. Doria 6, 95125 Catania, Italy
| |
Collapse
|
50
|
Schwing K, Gerhards M. Investigations on isolated peptides by combined IR/UV spectroscopy in a molecular beam – structure, aggregation, solvation and molecular recognition. INT REV PHYS CHEM 2016. [DOI: 10.1080/0144235x.2016.1229331] [Citation(s) in RCA: 39] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
|