1
|
Dong HC, Hsu PJ, Kuo JL. Searching low-energy conformers of neutral and protonated di-, tri-, and tetra-glycine using first-principles accuracy assisted by the use of neural network potentials. Phys Chem Chem Phys 2024; 26:11126-11139. [PMID: 38530660 DOI: 10.1039/d3cp05659g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/28/2024]
Abstract
In the last ten years, combinations of state-of-the-art gas-phase spectroscopies and quantum chemistry calculations have suggested several intuitive trends in the structure of small polypeptides that may not hold true. For example, the preference for the cis form of the peptide bond and multiple protonated sites was proposed by comparing experimental spectra with low-energy minima obtained from limited structural sampling using various density functional theory methods. For understanding the structures of polypeptides, extensive sampling of their configurational space with high-accuracy computational methods is required. In this work, we demonstrated the use of deep-learning neural network potential (DL-NNP) to assist in exploring the structure and energy landscape of di-, tri-, and tetra-glycine with the accuracy of high-level quantum chemistry methods, and low-energy conformers of small polypeptides can be efficiently located. We hope that the structures of these polypeptides we found and our preliminary analysis will stimulate further experimental investigations.
Collapse
Affiliation(s)
- Hieu Cao Dong
- Institute of Atomic and Molecular Sciences, Academia Sinica, Taipei, 10617, Taiwan.
- Molecular Science and Technology Program, Taiwan International Graduate Program, Academia Sinica, Taipei, 11529, Taiwan
- International Graduate Program of Molecular Science and Technology (NTU-MST), National Taiwan University, Taipei 10617, Taiwan
| | - Po-Jen Hsu
- Institute of Atomic and Molecular Sciences, Academia Sinica, Taipei, 10617, Taiwan.
| | - Jer-Lai Kuo
- Institute of Atomic and Molecular Sciences, Academia Sinica, Taipei, 10617, Taiwan.
- Molecular Science and Technology Program, Taiwan International Graduate Program, Academia Sinica, Taipei, 11529, Taiwan
- International Graduate Program of Molecular Science and Technology (NTU-MST), National Taiwan University, Taipei 10617, Taiwan
- Department of Chemistry, National Tsing Hua University, Hsinchu 30013, Taiwan
| |
Collapse
|
2
|
Abstract
AbstractFormic acid dimer as the prototypical doubly hydrogen-bonded gas-phase species is discussed from the perspective of the three translational and the three rotational degrees of freedom which are lost when two formic acid molecules form a stable complex. The experimental characterisation of these strongly hindered translations and rotations is reviewed, as are attempts to describe the associated fundamental vibrations, their combinations, and their thermal shifts by different electronic structure calculations and vibrational models. A remarkable match is confirmed for the combination of a CCSD(T)-level harmonic treatment and an MP2-level anharmonic VPT2 correction. Qualitatively correct thermal shifts of the vibrational spectra can be obtained from classical molecular dynamics in CCSD(T)-quality force fields. A detailed analysis suggests that this agreement between experiment and composite theoretical treatment is not strongly affected by fortuitous error cancellation but fully converged variational treatments of the six pair or intermolecular modes and their overtones and combinations in this model system would be welcome.
Collapse
|
3
|
Sapunar M, Domcke W, Došlić N. UV absorption spectra of DNA bases in the 350-190 nm range: assignment and state specific analysis of solvation effects. Phys Chem Chem Phys 2019; 21:22782-22793. [PMID: 31595896 DOI: 10.1039/c9cp04662c] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
Abstract
The theoretical assignment of electronic spectra of polyatomic molecules is a challenging problem that requires the specification of the character of a large number of electronic states. We propose a procedure for automatically determining the character of electronic transitions and apply it to the study of UV spectra of DNA bases in the gas phase and in the aqueous environment. The procedure is based on the computation of electronic wave function overlaps and accounts for an extensive sampling of nuclear geometries. Novelties of this work are the theoretical assignment of the electronic spectra of DNA bases up to 190 nm and a state specific analysis of solvation effects. By accounting for different effects contributing to the total solvent shift we obtained a good agreement between the computed and experimental spectra. Effects of vibrational averaging, temperature and solvent-induced structural changes shift excitation energies to lower values. Solvent-solute electrostatic interactions are state specific and strongly destabilize nRyd states, and to lesser extent nπ* and πRyd states. Altogether, this results in the stabilization of ππ* states and destabilization of nπ*, πRyd and nRyd states in solution.
Collapse
Affiliation(s)
- Marin Sapunar
- Department of Physical Chemistry, Ruder Bošković Institute, 10000 Zagreb, Croatia.
| | | | | |
Collapse
|
4
|
Houston PL, Van Hoozen BL, Qu C, Yu Q, Bowman JM. Teaching vibrational spectra to assign themselves. Faraday Discuss 2018; 212:65-82. [PMID: 30259026 DOI: 10.1039/c8fd00075a] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
A new paradigm for assigning vibrational spectra is described. Instead of proceeding from potential to Hamiltonian to eigenvalues/eigenvectors/intensities to spectrum, the new method goes "backwards" directly from spectrum to eigenvectors. The eigenvectors then "assign" the spectrum, in that they identify the basis states that contribute to each eigenvalue. To start, we demonstrate an algorithm that can obtain useful estimates of the eigenvectors connecting a real, symmetric Hamiltonian to its eigenvalues even if the only available information about the Hamiltonian is its diagonal elements. When this algorithm is augmented with information about transition intensities, it can be used to assign a complex vibrational spectrum using only information about (1) eigenvalues (the peak centers of the spectrum) and (2) a harmonic basis set (taken to be the diagonal elements of the Hamiltonian). Examples will be discussed, including application to the vibrationally complex spectral region of the formic acid dimer.
Collapse
Affiliation(s)
- Paul L Houston
- Department of Chemistry and Chemical Biology, Cornell Universtiy, Ithaca, NY 14853, USA. and School of Chemistry and Biochemistry, Georgia Inst. of Technology, Atlanta, GA 30332-0400, USA
| | - Brian L Van Hoozen
- Department of Chemistry, Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| | - Chen Qu
- Department of Chemistry, Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| | - Qi Yu
- Department of Chemistry, Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| | - Joel M Bowman
- Department of Chemistry, Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| |
Collapse
|
5
|
Qu C, Bowman JM. High-dimensional fitting of sparse datasets of CCSD(T) electronic energies and MP2 dipole moments, illustrated for the formic acid dimer and its complex IR spectrum. J Chem Phys 2018; 148:241713. [DOI: 10.1063/1.5017495] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Chen Qu
- Cherry L. Emerson Center for Scientific Computation and Department of Chemistry, Emory University, Atlanta, Georgia 30322, USA
| | - Joel M. Bowman
- Cherry L. Emerson Center for Scientific Computation and Department of Chemistry, Emory University, Atlanta, Georgia 30322, USA
| |
Collapse
|
6
|
Qu C, Bowman JM. IR Spectra of (HCOOH) 2 and (DCOOH) 2: Experiment, VSCF/VCI, and Ab Initio Molecular Dynamics Calculations Using Full-Dimensional Potential and Dipole Moment Surfaces. J Phys Chem Lett 2018; 9:2604-2610. [PMID: 29709189 DOI: 10.1021/acs.jpclett.8b00447] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
We report quantum VSCF/VCI and ab initio molecular dynamics (AIMD) calculations of the IR spectra of (HCOOH)2 and (DCOOH)2, using full-dimensional, ab initio potential energy and dipole moment surfaces (PES and DMS). These surfaces are fits, using permutationally invariant polynomials, to 13 475 ab initio CCSD(T)-F12a electronic energies and MP2 dipole moments. Here "AIMD" means using these ab initio potential and dipole moment surfaces in the MD calculations. The VSCF/VCI calculations use all (24) normal modes for coupling, with a four-mode representation of the potential. The quantum spectra align well with jet-cooled and room-temperature experimental spectra over the spectral range 600-3600 cm-1. Analyses of the complex O-H and C-H stretch bands are made based on the mixing of the VSCF/VCI basis functions. The comparisons of the AIMD IR spectra with both experimental and VSCF/VCI ones provide tests of the accuracy of the AIMD approach. These indicate good accuracy for simple bands but not for the complex O-H stretch band, which is upshifted from experimental and VSCF/VCI bands by roughly 300 cm-1. In addition to testing the AIMD approach, the PES, DMS, and VSCF/VCI calculations for formic acid dimer provide opportunities for testing other methods to represent high-dimensional data and other methods that perform postharmonic vibrational calculations.
Collapse
Affiliation(s)
- Chen Qu
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation , Emory University , Atlanta , Georgia 30322 , United States
| | - Joel M Bowman
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation , Emory University , Atlanta , Georgia 30322 , United States
| |
Collapse
|
7
|
Qu C, Bowman J. Quantum and classical IR spectra of (HCOOH)2, (DCOOH)2 and (DCOOD)2 using ab initio potential energy and dipole moment surfaces. Faraday Discuss 2018; 212:33-49. [DOI: 10.1039/c8fd00077h] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Abstract
Full-dimensional (24 modes) quantum calculation of the IR spectrum of (DCOOD)2, and comparison with classical MD one.
Collapse
Affiliation(s)
- Chen Qu
- Department of Chemistry
- Cherry L. Emerson Center for Scientific Computations
- Emory University
- Atlanta
- USA
| | - Joel M. Bowman
- Department of Chemistry
- Cherry L. Emerson Center for Scientific Computations
- Emory University
- Atlanta
- USA
| |
Collapse
|
8
|
Silber D, Kowalski PM, Traeger F, Buchholz M, Bebensee F, Meyer B, Wöll C. Adsorbate-induced lifting of substrate relaxation is a general mechanism governing titania surface chemistry. Nat Commun 2016; 7:12888. [PMID: 27686286 PMCID: PMC5056433 DOI: 10.1038/ncomms12888] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2016] [Accepted: 08/09/2016] [Indexed: 11/17/2022] Open
Abstract
Under ambient conditions, almost all metals are coated by an oxide. These coatings, the result of a chemical reaction, are not passive. Many of them bind, activate and modify adsorbed molecules, processes that are exploited, for example, in heterogeneous catalysis and photochemistry. Here we report an effect of general importance that governs the bonding, structure formation and dissociation of molecules on oxidic substrates. For a specific example, methanol adsorbed on the rutile TiO2(110) single crystal surface, we demonstrate by using a combination of experimental and theoretical techniques that strongly bonding adsorbates can lift surface relaxations beyond their adsorption site, which leads to a significant substrate-mediated interaction between adsorbates. The result is a complex superstructure consisting of pairs of methanol molecules and unoccupied adsorption sites. Infrared spectroscopy reveals that the paired methanol molecules remain intact and do not deprotonate on the defect-free terraces of the rutile TiO2(110) surface. Molecules on a metal surface may be modified by the presence of oxide layers, but further mechanistic understanding is still required. Here the authors show for methanol on rutile TiO2(110) that strongly bonded adsorbates lift surface relaxations, leading to substrate-mediated interaction between adsorbates.
Collapse
Affiliation(s)
- David Silber
- Institute of Functional Interfaces (IFG), Karlsruhe Institute of Technology (KIT), Hermann-von-Helmholtz-Platz 1, 76344 Eggenstein-Leopoldshafen, Germany
| | - Piotr M Kowalski
- Institute of Energy and Climate Research (IEK-6), Forschungszentrum Jülich, Wilhelm-Johnen-Strasse, 52425 Jülich, Germany
| | - Franziska Traeger
- Westfälische Hochschule, FB 8, August-Schmidt-Ring 10, 45665 Recklinghausen, Germany
| | - Maria Buchholz
- Institute of Functional Interfaces (IFG), Karlsruhe Institute of Technology (KIT), Hermann-von-Helmholtz-Platz 1, 76344 Eggenstein-Leopoldshafen, Germany
| | - Fabian Bebensee
- Institute of Functional Interfaces (IFG), Karlsruhe Institute of Technology (KIT), Hermann-von-Helmholtz-Platz 1, 76344 Eggenstein-Leopoldshafen, Germany
| | - Bernd Meyer
- Interdisciplinary Center for Molecular Materials (ICMM) and Computer-Chemistry-Center (CCC), Friedrich-Alexander-Universität Erlangen-Nürnberg, Nägelsbachstrasse 25, 91052 Erlangen, Germany
| | - Christof Wöll
- Institute of Functional Interfaces (IFG), Karlsruhe Institute of Technology (KIT), Hermann-von-Helmholtz-Platz 1, 76344 Eggenstein-Leopoldshafen, Germany
| |
Collapse
|
9
|
Zhou W, Cao M, Su S, Li N, Zhao X, Wang J, Li X, Hu C. A newly-designed polyoxometalate-based plasmonic visible-light catalyst for the photodegradation of methyl blue. ACTA ACUST UNITED AC 2013. [DOI: 10.1016/j.molcata.2013.01.028] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/01/2022]
|
10
|
Novak J, Mališ M, Prlj A, Ljubić I, Kühn O, Došlić N. Photoinduced Dynamics of Formic Acid Monomers and Dimers: The Role of the Double Hydrogen Bond. J Phys Chem A 2012; 116:11467-75. [DOI: 10.1021/jp3074865] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Jurica Novak
- Rudjer Bošković Institute, 10000 Zagreb, Croatia
| | - Momir Mališ
- Rudjer Bošković Institute, 10000 Zagreb, Croatia
| | - Antonio Prlj
- Faculty of Science, University of Zagreb, 10000 Zagreb, Croatia
| | - Ivan Ljubić
- Rudjer Bošković Institute, 10000 Zagreb, Croatia
| | - Oliver Kühn
- Institut für Physik, Universität Rostock, D-18051 Rostock, Germany
| | - Nađa Došlić
- Rudjer Bošković Institute, 10000 Zagreb, Croatia
| |
Collapse
|
11
|
Affiliation(s)
- Chao Lv
- Key Laboratory of Colloid and Interface Chemistry, Shandong University
| | - Guanchen Xu
- Key Laboratory of Colloid and Interface Chemistry, Shandong University
| | - Xiao Chen
- Key Laboratory of Colloid and Interface Chemistry, Shandong University
| |
Collapse
|
12
|
Saritha B, Durga Prasad M. An Electronic Structure Perspective of the Promoter Modes in Proton Transfer Reactions. J Phys Chem A 2011; 115:2802-10. [DOI: 10.1021/jp108579t] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- B. Saritha
- School of Chemistry, University of Hyderabad, Hyderabad 500046, India
| | - M. Durga Prasad
- School of Chemistry, University of Hyderabad, Hyderabad 500046, India
| |
Collapse
|
13
|
Ghalla H, Rekik N, Michta A, Oujia B, Flakus HT. Theoretical modeling of infrared spectra of the hydrogen and deuterium bond in aspirin crystal. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2010; 75:37-47. [PMID: 19884041 DOI: 10.1016/j.saa.2009.09.029] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/15/2009] [Revised: 09/06/2009] [Accepted: 09/13/2009] [Indexed: 05/28/2023]
Abstract
An extended quantum theoretical approach of the nu(X-H) IR lineshape of cyclic dimers of weakly H-bonded species is proposed. We have extended a previous approach [M.E.-A. Benmalti, P. Blaise, H.T. Flakus, O. Henri-Rousseau, Chem. Phys. 320 (2006) 267] by accounting for the anharmonicity of the slow mode which is described by a "Morse" potential in order to reproduce the polarized infrared spectra of the hydrogen and deuterium bond in acetylsalicylic acid (aspirin) crystals. From comparison of polarized IR spectra of isotopically neat and isotopically diluted aspirin crystals it resulted that centrosymmetric aspirin dimer was the bearer of the crystal main spectral properties. In this approach, the adiabatic approximation is performed for each separate H-bond bridge of the dimer and a strong non-adiabatic correction is introduced into the model via the resonant exchange between the fast mode excited states of the two moieties. Within the strong anharmonic coupling theory, according to which the X-H...Y high-frequency mode is anharmonically coupled to the H-bond bridge, this model incorporated the Davydov coupling between the excited states of the two moieties, the quantum direct and indirect dampings and the anharmonicity for the H-bond bridge. The spectral density is obtained within the linear response theory by Fourier transform of the damped autocorrelation functions. The evaluated spectra are in fairly good agreement with the experimental ones by using a minimum number of independent parameters. The effect of deuteration has been well reproduced by reducing simply the angular frequency of the fast mode and the anharmonic coupling parameter.
Collapse
Affiliation(s)
- Houcine Ghalla
- Laboratoire de Physique Quantique, Faculté des Sciences de Monastir, 5000 route de Kairouan, Tunisia
| | | | | | | | | |
Collapse
|
14
|
Petković M, Novak J, Došlić N. Shaping the infrared spectrum of the acetic acid dimer in the OH-stretching range: Multiple conformers and anharmonic coupling. Chem Phys Lett 2009. [DOI: 10.1016/j.cplett.2009.04.014] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
15
|
Mališ M, Matanović I, Došlić N. A Computational Study of Electronic and Spectroscopic Properties of Formic Acid Dimer Isotopologues. J Phys Chem A 2009; 113:6034-40. [DOI: 10.1021/jp901067u] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Affiliation(s)
- M. Mališ
- Department of Physical Chemistry, R. Bošković Institute, Bijenička 54, 10000, Zagreb, Croatia
| | - I. Matanović
- Department of Physical Chemistry, R. Bošković Institute, Bijenička 54, 10000, Zagreb, Croatia
| | - N. Došlić
- Department of Physical Chemistry, R. Bošković Institute, Bijenička 54, 10000, Zagreb, Croatia
| |
Collapse
|
16
|
Ito F. Infrared spectra of (HCOOH)2 and (DCOOH)2 in rare gas matrices: A comparative study with gas phase spectra. J Chem Phys 2008; 128:114310. [DOI: 10.1063/1.2841078] [Citation(s) in RCA: 36] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022] Open
|
17
|
Matanović I, Došlić N, Johnson BR. Generalized approximation to the reaction path: The formic acid dimer case. J Chem Phys 2008; 128:084103. [DOI: 10.1063/1.2833978] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
|