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
|
Shivankar BR, Bhandare VV, Joshi K, Patil VS, Dhotare PS, Sonawane KD, Krishnamurty S. Investigation of cathinone analogs targeting human dopamine transporter using molecular modeling. J Biomol Struct Dyn 2024:1-16. [PMID: 38698732 DOI: 10.1080/07391102.2024.2335303] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/16/2023] [Accepted: 03/20/2024] [Indexed: 05/05/2024]
Abstract
In a step towards understanding the structure-property relationship among Synthetic Cathinones (SCs), a combined methodology based on Density Functional Theory (DFT), Administration, Distribution, Metabolism, Excretion, and Toxicity (ADMET) predictions, docking and molecular dynamics simulations have been applied to correlate physicochemical descriptors of various SCs to their biological activity. The results from DFT and molecular docking studies correlate well with each other explaining the biological activity trends of the studied SCs. Quantum mechanical descriptors viz. polarizability, electron affinity, ionization potential, chemical hardness, electronegativity, molecular electrostatic potential, and ion interaction studies unravel the distinguishingly reactive nature of Group D (pyrrolidine substituted) and Group E (methylenedioxy and pyrrolidine substituted) compounds. According to ADMET analysis, Group D and Group E molecules have a higher probability of permeating through the blood-brain barrier. Molecular docking results indicate that Phe76, Ala77, Asp79, Val152, Tyr156, Phe320, and Phe326 constitute the binding pocket residues of hDAT in which the most active ligands MDPV, MDPBP, and MDPPP are bound. Finally, to validate the derived quantum chemical descriptors and docking results, Molecular Dynamics (MD) simulations are performed with homology-modelled hDAT (human dopamine transporter). The MD simulation results revealed that the majority of SCs remain stable within the hDAT protein's active sites via non-bonded interactions after 100 ns long simulations. The findings from DFT, ADMET analysis, molecular docking, and molecular dynamics simulation studies complement each other suggesting that pyrrolidine-substituted SCs (Group D and E), specifically, MPBP and PVN are proven potent SCs along with MDPV, validating various experimental observations.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Bhavana R Shivankar
- Physical Chemistry Division, CSIR-National Chemical Laboratory, Pune, India
- Academy of Scientific and Innovative Research (AcSIR), Ghaziabad, India
| | | | - Krati Joshi
- Physical Chemistry Division, CSIR-National Chemical Laboratory, Pune, India
| | - Vishal S Patil
- Department of Pharmacology and Toxicology, KLE College of Pharmacy Belagavi, KLE Academy of Higher Education and Research (KAHER), Belagavi, India
- ICMR-National Institute of Traditional Medicine, Belagavi, Karnataka, India
| | | | | | - Sailaja Krishnamurty
- Physical Chemistry Division, CSIR-National Chemical Laboratory, Pune, India
- Academy of Scientific and Innovative Research (AcSIR), Ghaziabad, India
| |
Collapse
|
3
|
Schröder B, Rauhut G. From the Automated Calculation of Potential Energy Surfaces to Accurate Infrared Spectra. J Phys Chem Lett 2024; 15:3159-3169. [PMID: 38478898 PMCID: PMC10961845 DOI: 10.1021/acs.jpclett.4c00186] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/19/2024] [Revised: 02/20/2024] [Accepted: 02/28/2024] [Indexed: 03/22/2024]
Abstract
Advances in the development of quantum chemical methods and progress in multicore architectures in computer science made the simulation of infrared spectra of isolated molecules competitive with respect to established experimental methods. Although it is mainly the multidimensional potential energy surface that controls the accuracy of these calculations, the subsequent vibrational structure calculations need to be carefully converged in order to yield accurate results. As both aspects need to be considered in a balanced way, we focus on approaches for molecules of up to 12-15 atoms with respect to both parts, which have been automated to some extent so that they can be employed in routine applications. Alternatives to machine learning will be discussed, which appear to be attractive, as long as local regions of the potential energy surface are sufficient. The automatization of these methods is still in its infancy, and the generalization to molecules with large amplitude motions or molecular clusters is far from trivial, but many systems relevant for astrophysical studies are already in reach.
Collapse
Affiliation(s)
- Benjamin Schröder
- Institute
of Physical Chemistry, University of Goettingen, Tammannstrasse 6, Göttingen 37077, Germany
| | - Guntram Rauhut
- Institute
for Theoretical Chemistry, University of
Stuttgart, Pfaffenwaldring 55, Stuttgart 70569, Germany
| |
Collapse
|
4
|
Yang Z, Chen G, Shen J, Ma C, Gu J, Zhu C, Li L, Gao H. A hydrogen bonding based SERS method for direct label-free L-hydroxyproline detection. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2023; 299:122834. [PMID: 37178585 DOI: 10.1016/j.saa.2023.122834] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/27/2023] [Revised: 05/03/2023] [Accepted: 05/05/2023] [Indexed: 05/15/2023]
Abstract
The detection of non-protein nitrogen adulterants is a major challenge in dairy testing. As a marker molecule of animal hydrolyzed protein, the presence of non-edible L-hydroxyproline (L-Hyp) molecules can be used to identify low-quality milk containing components of animal hydrolyzed protein. However, it is still difficult to detect L-Hyp directly in milk. The Ag@COF-COOH substrate in this paper can be used to realize label-free L-Hyp sensitive detection based on the hydrogen bond transition mechanism. To explore the mechanism, the binding sites of hydrogen bond interaction have been verified experimentally and computationally, and the charge transfer process was also explained in terms of HOMO/LOMO energy level. In conclusion, the quantitative models for L-Hyp in an aqueous environment and in milk were developed. The limit of detection (LOD) of L-Hyp in an aqueous environment could reach 8.18 ng/mL, with R2 of 0.982. The linear range of quantitative detection in milk was 0.5-1000 μg/mL and the LOD was as low as 0.13 μg/mL. In this work, a hydrogen bond interaction based Surface-enhanced Raman spectroscopy (SERS) method for the label-free detection of L-Hyp was proposed, which complemented the application of SERS technology in the detection of dairy products.
Collapse
Affiliation(s)
- Zichen Yang
- School of Science, Jiangnan University, Wuxi, China; School of Internet of Things Engineering, Jiangnan University, Wuxi, China; Jiangsu Provincial Research Center of Light Industrial Optoelectronic Engineering and Technology, Wuxi, China
| | - Guoqing Chen
- School of Science, Jiangnan University, Wuxi, China; Jiangsu Provincial Research Center of Light Industrial Optoelectronic Engineering and Technology, Wuxi, China.
| | - Jialu Shen
- School of Science, Jiangnan University, Wuxi, China; Jiangsu Provincial Research Center of Light Industrial Optoelectronic Engineering and Technology, Wuxi, China
| | - Chaoqun Ma
- School of Science, Jiangnan University, Wuxi, China; Jiangsu Provincial Research Center of Light Industrial Optoelectronic Engineering and Technology, Wuxi, China
| | - Jiao Gu
- School of Science, Jiangnan University, Wuxi, China; Jiangsu Provincial Research Center of Light Industrial Optoelectronic Engineering and Technology, Wuxi, China
| | - Chun Zhu
- School of Science, Jiangnan University, Wuxi, China; Jiangsu Provincial Research Center of Light Industrial Optoelectronic Engineering and Technology, Wuxi, China
| | - Lei Li
- School of Science, Jiangnan University, Wuxi, China; Jiangsu Provincial Research Center of Light Industrial Optoelectronic Engineering and Technology, Wuxi, China
| | - Hui Gao
- School of Science, Jiangnan University, Wuxi, China; Jiangsu Provincial Research Center of Light Industrial Optoelectronic Engineering and Technology, Wuxi, China
| |
Collapse
|
5
|
Fayaz A, Banik S, Kanchan Roy T. The importance of electron correlations on vibrational anharmonicities and potential energy surfaces. COMPUT THEOR CHEM 2023. [DOI: 10.1016/j.comptc.2023.114059] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
|
6
|
Wang T, Shi H, Yu Z, Liu T. Theoretical insight into the mechanism of palladium-catalyzed oxidative cascade reaction of phenylacetylene and allenes with different assisting groups. J Organomet Chem 2022. [DOI: 10.1016/j.jorganchem.2022.122480] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/31/2022]
|
7
|
Dall’Osto G, Corni S. Time Resolved Raman Scattering of Molecules: A Quantum Mechanics Approach with Stochastic Schroedinger Equation. J Phys Chem A 2022; 126:8088-8100. [PMID: 36278928 PMCID: PMC9639147 DOI: 10.1021/acs.jpca.2c05245] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
Raman scattering is a very powerful tool employed to
characterize
molecular systems. Here we propose a novel theoretical strategy to
calculate the Raman cross-section in time domain, by computing the
cumulative Raman signal emitted during the molecular evolution in
time. Our model is based on a numerical propagation of the vibronic
wave function under the effect of a light pulse of arbitrary shape.
This approach can therefore tackle a variety of experimental setups.
Both resonance and nonresonance Raman scattering can be retrieved,
and also the time-dependent fluorescence emission is computed. The
model has been applied to porphyrin considering both resonance and
nonresonance conditions and varying the incident field duration. Moreover
the effect of the vibrational relaxation, which should be taken into
account when its time scale is similar to that of the Raman emission,
has been included through the stochastic Schroedinger equation approach.
Collapse
Affiliation(s)
- Giulia Dall’Osto
- Department of Chemical Sciences, University of Padova, via Marzolo 1, Padova, 35131, Italy
| | - Stefano Corni
- Department of Chemical Sciences, University of Padova, via Marzolo 1, Padova, 35131, Italy
- CNR Institute of Nanoscience, via Campi 213/A, Modena, 41125, Italy
| |
Collapse
|
8
|
Effects of non-local exchange functionals in the density functional theories for the description of molecular vibrations. J CHEM SCI 2022. [DOI: 10.1007/s12039-022-02061-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
9
|
Trung NQ, Mechler A, Hoa NT, Vo QV. Calculating bond dissociation energies of X-H (X=C, N, O, S) bonds of aromatic systems via density functional theory: a detailed comparison of methods. ROYAL SOCIETY OPEN SCIENCE 2022; 9:220177. [PMID: 35706655 PMCID: PMC9174704 DOI: 10.1098/rsos.220177] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/18/2022] [Accepted: 05/06/2022] [Indexed: 05/03/2023]
Abstract
In this study, the performance of 17 different density functional theory functionals was compared for the calculation of the bond dissociation energy (BDE) values of X-H (X=C, N, O, S) bonds of aromatic compounds. The effect of the size of the basis set (expansions of 6-31(G)) was also assessed for the initial geometry and zero-point energy calculations, followed by the single-point BDE calculations with different model chemistries with the 6-311 + (3df,2p) basis set. It was found that the size of the basis set for geometry optimization has a much smaller effect on the accuracy of BDE than the choice of functional for the following single-point calculations. The M06-2X, M05-2X and M08-HX functionals yielded highly accurate BDE values compared to experimental data (with the average mean unsigned error MUE = 1.2-1.5 kcal mol-1), performing better than any of the other functionals. The results suggest that geometry optimization may be performed with B3LYP functional and a small basis set, whereas the M06-2X, M05-2X and M08-HX density functionals with a suitably large basis set offer the best method for calculating BDEs of ArX-H (X=C, N, O, S) bonds.
Collapse
Affiliation(s)
- Nguyen Quang Trung
- Department of Chemistry, The University of Danang - University of Science and Education, Da Nang 550000, Vietnam
- Quality Assurance and Testing Center 2, Danang 550000, Vietnam
| | - Adam Mechler
- Department of Chemistry and Physics, La Trobe University, Victoria 3086, Australia
| | - Nguyen Thi Hoa
- Academic Affairs, The University of Danang - University of Technology and Education, Danang 550000, Vietnam
| | - Quan V. Vo
- Faculty of Chemical Technology – Environment, The University of Danang - University of Technology and Education, Danang 550000, Vietnam
| |
Collapse
|
10
|
Guo S, Zhu C, Chen G, Gu J, Ma C, Gao H, Li L, Zhang Y, Li X, Wang Z, Wei Y, Wang G, Shen J. A theoretical study on intermolecular hydrogen bonds of isopropanol-water clusters. Theor Chem Acc 2022. [DOI: 10.1007/s00214-022-02865-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
|
11
|
Jakubowska K, Pecul M, Ruud K. Relativistic Four-Component DFT Calculations of Vibrational Frequencies. J Phys Chem A 2021; 125:10315-10320. [PMID: 34843253 PMCID: PMC8667032 DOI: 10.1021/acs.jpca.1c07398] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
We investigate the
effect of relativity on harmonic vibrational
frequencies. Density functional theory (DFT) calculations using the
four-component Dirac–Coulomb Hamiltonian have been performed
for 15 hydrides (H2X, X = O, S, Se, Te, Po; XH3, X = N, P, As, Sb, Bi; and XH4, X = C, Si, Ge, Sn, Pb)
as well as for HC≡CPbH3. The vibrational frequencies
have been calculated using finite differences of the molecular energy
with respect to geometrical distortions of the nuclei. The influences
of the choice of basis set, exchange–correlation functional,
and step length for the numerical differentiation on the calculated
harmonic vibrational frequencies have been tested, and the method
has been found to be numerically robust. Relativistic effects are
noticeable for the heavier congeners H2Te and H2Po, SbH3 and BiH3, and SnH4 and
PbH4 and are much more pronounced for the vibrational modes
with higher frequencies. Spin–orbit effects constitute a very
small fraction of the total relativistic effects, except for H2Te and H2Po. For HC≡CPbH3 we
find that only the frequencies of the modes with large contributions
from Pb displacements are significantly affected by relativity.
Collapse
Affiliation(s)
| | - Magdalena Pecul
- Faculty of Chemistry, University of Warsaw, 02-093 Warsaw, Poland
| | - Kenneth Ruud
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, UiT - The Arctic University of Norway, N-9019 Tromsø, Norway
| |
Collapse
|
12
|
Asif M, Sajid H, Ayub K, Khan AA, Ahmad R, Ans M, Mahmood T. Nano-porous C 4N as a toxic pesticide's scavenger: A quantum chemical approach. J Mol Graph Model 2021; 111:108078. [PMID: 34826716 DOI: 10.1016/j.jmgm.2021.108078] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/20/2021] [Revised: 11/11/2021] [Accepted: 11/12/2021] [Indexed: 11/27/2022]
Abstract
The sensing affinity of C4N is the most fascinating topic of research due to its excellent chemical and electronic properties. Moreover, owing to the highly active porous cavity, C4N can easily accommodate foreign molecules. Herein, we studied the adsorption properties of carbamate insecticides (CMs) namely, Dimetalin (DMT), Carbanolate (CBT), Isolan (ISO) and Propoxur (PRO) using density functional theory calculations. All the results are calculated at widely accepted ωB97XD functional along with 6-31G(d, p) basis set. The calculated counterpoise corrected interaction energy of the reported complexes ranges between -20.05 and -27.04 kcal/mol, however, the interaction distances are found to be higher than 2.00 Å. The values of interacting parameters depict that the carbamate molecules are physisorbed via noncovalent interactions that can easily be reversible. Moreover, the binding of selected insecticides notably changes the electronic structure of C4N. The electronic changes are characterized by the energies of HOMO & LUMO, their energy gaps and CHELPG charge transfer. The charge density difference between C4N surface and carbamate pesticides are characterized by EDD and CDA analysis. Moreover, the ab initio molecular dynamic study reveals that the complexes are stable even at 500 K. The photochemical sensing properties of C4N are estimated by time dependent UV-Vis calculations. The high sensitivity of C4N towards considered analytes enable it to act as a promising sensor for toxic pesticides.
Collapse
Affiliation(s)
- Misbah Asif
- Department of Chemistry, COMSATS University, Abbottabad Campus, Abbottabad, 22060, Pakistan
| | - Hasnain Sajid
- School of Science and Technology, Nottingham Trent University, Clifton Lane, Nottingham, NG11 8NS, UK
| | - Khurshid Ayub
- Department of Chemistry, COMSATS University, Abbottabad Campus, Abbottabad, 22060, Pakistan
| | - Adnan Ali Khan
- Center for Computational Materials Science, University of Malakand, Pakistan; Department of Chemistry, University of Malakand, Pakistan
| | - Rashid Ahmad
- Center for Computational Materials Science, University of Malakand, Pakistan; Department of Chemistry, University of Malakand, Pakistan
| | - Muhammad Ans
- Department of Chemistry, University of Agriculture, Faisalabad, Punjab, Pakistan
| | - Tariq Mahmood
- Department of Chemistry, COMSATS University, Abbottabad Campus, Abbottabad, 22060, Pakistan; Department of Chemistry, College of Science, University of Bahrain, P.O. Box 32038, Bahrain.
| |
Collapse
|
13
|
Sheng M, Silvestrini F, Biczysko M, Puzzarini C. Structural and Vibrational Properties of Amino Acids from Composite Schemes and Double-Hybrid DFT: Hydrogen Bonding in Serine as a Test Case. J Phys Chem A 2021; 125:9099-9114. [PMID: 34623165 DOI: 10.1021/acs.jpca.1c06993] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
The structures, relative stabilities, and vibrational wavenumbers of the two most stable conformers of serine, stabilized by the O-H···N, O-H···O═C and N-H···O-H intramolecular hydrogen bonds, have been evaluated by means of state-of-the-art composite schemes based on coupled-cluster (CC) theory. The so-called "cheap" composite approach (CCSD(T)/(CBS+CV)MP2) allowed determination of accurate equilibrium structures and harmonic vibrational wavenumbers, also pointing out significant corrections beyond the CCSD(T)/cc-pVTZ level. These accurate results stand as a reference for benchmarking selected hybrid and double-hybrid, dispersion-corrected DFT functionals. B2PLYP-D3 and DSDPBEP86 in conjunction with a triple-ζ basis set have been confirmed as effective methodologies for structural and spectroscopic studies of medium-sized flexible biomolecules, also showing intramolecular hydrogen bonding. These best performing double-hybrid functionals have been employed to simulate IR spectra by means of vibrational perturbation theory, also considering hybrid CC/DFT schemes. The best overall agreement with experiment, with mean absolute error of 8 cm-1, has been obtained by combining CCSD(T)/(CBS+CV)MP2 harmonic wavenumbers with B2PLYP-D3/maug-cc-pVTZ anharmonic corrections. Finally, a composite scheme entirely based on CCSD(T) calculations (CCSD(T)/CBS+CV) has been employed for energetics, further confirming that serine II is the most stable conformer, also when zero-point vibrational energy corrections are included.
Collapse
Affiliation(s)
- Mingzhu Sheng
- International Centre for Quantum and Molecular Structures, Physics Department, Shanghai University, 99 Shangda Road, Shanghai, 200444, China
| | - Filippo Silvestrini
- Department of Chemistry "Giacomo Ciamician", University of Bologna, Via F. Selmi 2, 40126 Bologna, Italy
| | - Malgorzata Biczysko
- International Centre for Quantum and Molecular Structures, Physics Department, Shanghai University, 99 Shangda Road, Shanghai, 200444, China
| | - Cristina Puzzarini
- Department of Chemistry "Giacomo Ciamician", University of Bologna, Via F. Selmi 2, 40126 Bologna, Italy
| |
Collapse
|
14
|
Mendolicchio M, Bloino J, Barone V. General Perturb-Then-Diagonalize Model for the Vibrational Frequencies and Intensities of Molecules Belonging to Abelian and Non-Abelian Symmetry Groups. J Chem Theory Comput 2021; 17:4332-4358. [PMID: 34085530 PMCID: PMC8280743 DOI: 10.1021/acs.jctc.1c00240] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/08/2021] [Indexed: 11/29/2022]
Abstract
In this paper, we show that the standard second-order vibrational perturbation theory (VPT2) for Abelian groups can be used also for non-Abelian groups without employing specific equations for two- or threefold degenerate vibrations but rather handling in the proper way all the degeneracy issues and deriving the peculiar spectroscopic signatures of non-Abelian groups (e.g., l -doubling) by a posteriori transformations of the eigenfunctions. Comparison with the results of previous conventional implementations shows a perfect agreement for the vibrational energies of linear and symmetric tops, thus paving the route to the transparent extension of the equations already available for asymmetric tops to the energies of spherical tops and the infrared and Raman intensities of molecules belonging to non-Abelian symmetry groups. The whole procedure has been implemented in our general engine for vibro-rotational computations beyond the rigid rotor/harmonic oscillator model and has been validated on a number of test cases.
Collapse
Affiliation(s)
| | - Julien Bloino
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| | - Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| |
Collapse
|
15
|
Izadnia S, LaForge AC, Stienkemeier F, Cheeseman JR, Bloino J, Cheramy J, Jäger W, Xu Y. Unusual binary aggregates of perylene bisimide revealed by their electronic transitions in helium nanodroplets and DFT calculations. Phys Chem Chem Phys 2021; 23:13862-13872. [PMID: 34159991 DOI: 10.1039/d1cp01923f] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The S1 ← S0 electronic transition of perylene bisimide (PBI) and its binary aggregates were investigated using a combination of helium nanodroplet isolation spectroscopy and computational methods. First, well-resolved vibronic bands of the PBI monomer obtained under the superfluid helium nanodroplet environment were compared to simulated vibronic spectra with anharmonic corrections of the band positions. Second, about ten sets of weaker vibronic bands were observed, which show similar vibronic patterns as that of the PBI monomer and have their band origins red-shifted by about 8 to 218 cm-1. Experimental Poisson curve analyses, performed at the origins of these new sets of bands and the PBI monomer, indicate that the carriers of these weaker red-shifted vibronic bands are binary adducts of PBI. Three types of PBI dimer structures where the electronic transition dipole moments of the two subunits are perpendicular to each other were proposed as possible carriers of these red-shifted vibronic patterns. Extensive vibronic simulations were carried out in a multi-step procedure with TD-DFT, vertical Hessian, and finally adiabatic Hessian approaches. Small red-shifted band origins and very similar vibronic patterns to that of the monomer were predicted for unusual, T-shaped, type I dimer structures and are in close agreement with the experimental data. The combined experimental and theoretical results indicate that the helium nanodroplet environment enables the formation of these unusual T-shaped dimers and stabilizes them.
Collapse
Affiliation(s)
- S Izadnia
- Institute of Physics, University of Freiburg, Hermann-Herder-Str. 3, D-79104 Freiburg, Germany.
| | - A C LaForge
- Institute of Physics, University of Freiburg, Hermann-Herder-Str. 3, D-79104 Freiburg, Germany.
| | - F Stienkemeier
- Institute of Physics, University of Freiburg, Hermann-Herder-Str. 3, D-79104 Freiburg, Germany.
| | - James R Cheeseman
- Gaussian Inc., 340 Quinnipiac St., Bldg. 40, Wallingford, Connecticut 06492-4050, USA
| | - Julien Bloino
- Scuola Normale Superiore, Piazza dei Cavalieri, 7, 56126 Pisa, Italy
| | - J Cheramy
- Department of Chemistry, University of Alberta, Edmonton, Alberta T6G 2G2, Canada.
| | - Wolfgang Jäger
- Department of Chemistry, University of Alberta, Edmonton, Alberta T6G 2G2, Canada.
| | - Yunjie Xu
- Department of Chemistry, University of Alberta, Edmonton, Alberta T6G 2G2, Canada.
| |
Collapse
|
16
|
Rey M, Chizhmakova IS, Nikitin AV, Tyuterev VG. Towards a complete elucidation of the ro-vibrational band structure in the SF 6 infrared spectrum from full quantum-mechanical calculations. Phys Chem Chem Phys 2021; 23:12115-12126. [PMID: 34032236 DOI: 10.1039/d0cp05727d] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
The first accurate and complete theoretical room-temperature rotationally resolved spectra in the range 300-3000 cm-1 are reported for the three most abundant isotopologues (32SF6, 33SF6 and 34SF6) of the sulfur hexafluoride molecule. The literature reports that SF6 is widely used as a prototype molecule for studying the multi-photon excitation processes with powerful lasers in the infrared range. On the other hand, SF6 is an important greenhouse molecule with a very long lifetime in the atmosphere. Because of relatively low vibrational frequencies, the hot bands of this molecule contribute significantly to the absorption infrared spectra even at room temperature. This makes the calculation of complete rovibrational line lists required for fully converged opacity modeling extremely demanding. In order to reduce the computational costs, symmetry was exploited at all stages of the first global variational nuclear motion calculations by means of irreducible tensor operators. More than 2600 new vibrational band centers were predicted using our empirically refined ab initio potential energy surface. Highly excited rotational states were calculated up to J = 121, resulting in 6 billion transitions computed from an ab initio dipole moment surface and distributed over more than 500 cold and hot bands. The final line lists are made available through the TheoReTS information system (http://theorets.univ-reims.fr, http://theorets.tsu.ru). For the first time, the major (ro)vibrational band structures in the wavenumber range corresponding to the strongest absorption in the infra-red are completely elucidated for a seven-atom molecule, providing excellent agreement with the observed spectral patterns. It is shown that the obtained results are more complete than all available line lists, permitting reliable modelling of the temperature dependence of the molecular opacity.
Collapse
Affiliation(s)
- Michaël Rey
- Groupe de Spectrométrie Moléculaire et Atmosphérique, UMR CNRS 7331, BP 1039, F-51687, Reims Cedex 2, France.
| | - Iana S Chizhmakova
- Laboratory of Quantum Mechanics of Molecules and Radiative Processes, Tomsk State University, 36 Lenin Avenue, 634050 Tomsk, Russia
| | - Andrei V Nikitin
- V.E. Zuev Institute of Atmospheric Optics, Russian Academy of Sciences, 1, Akademichesky Avenue, 634055 Tomsk, Russian Federation
| | - Vladimir G Tyuterev
- Groupe de Spectrométrie Moléculaire et Atmosphérique, UMR CNRS 7331, BP 1039, F-51687, Reims Cedex 2, France. and Laboratory of Quantum Mechanics of Molecules and Radiative Processes, Tomsk State University, 36 Lenin Avenue, 634050 Tomsk, Russia
| |
Collapse
|
17
|
Käser S, Boittier ED, Upadhyay M, Meuwly M. Transfer Learning to CCSD(T): Accurate Anharmonic Frequencies from Machine Learning Models. J Chem Theory Comput 2021; 17:3687-3699. [PMID: 33960787 DOI: 10.1021/acs.jctc.1c00249] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
The calculation of the anharmonic modes of small- to medium-sized molecules for assigning experimentally measured frequencies to the corresponding type of molecular motions is computationally challenging at sufficiently high levels of quantum chemical theory. Here, a practical and affordable way to calculate coupled-cluster quality anharmonic frequencies using second-order vibrational perturbation theory (VPT2) from machine-learned models is presented. The approach, referenced as "NN + VPT2", uses a high-dimensional neural network (PhysNet) to learn potential energy surfaces (PESs) at different levels of theory from which harmonic and VPT2 frequencies can be efficiently determined. The NN + VPT2 approach is applied to eight small- to medium-sized molecules (H2CO, trans-HONO, HCOOH, CH3OH, CH3CHO, CH3NO2, CH3COOH, and CH3CONH2) and frequencies are reported from NN-learned models at the MP2/aug-cc-pVTZ, CCSD(T)/aug-cc-pVTZ, and CCSD(T)-F12/aug-cc-pVTZ-F12 levels of theory. For the largest molecules and at the highest levels of theory, transfer learning (TL) is used to determine the necessary full-dimensional, near-equilibrium PESs. Overall, NN + VPT2 yields anharmonic frequencies to within 20 cm-1 of experimentally determined frequencies for close to 90% of the modes for the highest quality PES available and to within 10 cm-1 for more than 60% of the modes. For the MP2 PESs only ∼60% of the NN + VPT2 frequencies were within 20 cm-1 of the experiment, with outliers up to ∼150 cm-1, compared to the experiment. It is also demonstrated that the approach allows to provide correct assignments for strongly interacting modes such as the OH bending and the OH torsional modes in formic acid monomer and the CO-stretch and OH-bend mode in acetic acid.
Collapse
Affiliation(s)
- Silvan Käser
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| | - Eric D Boittier
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| | - Meenu Upadhyay
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| | - Markus Meuwly
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| |
Collapse
|
18
|
Aarjane M, Slassi S, Ghaleb A, Amine A. Synthesis, spectroscopic characterization (FT-IR, NMR) and DFT computational studies of new isoxazoline derived from acridone. J Mol Struct 2021. [DOI: 10.1016/j.molstruc.2021.129921] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
|
19
|
IR spectroscopy of condensed phase systems: Can the environment induce vibrational mode coupling? Chem Phys Lett 2021. [DOI: 10.1016/j.cplett.2020.138168] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
|
20
|
Tourchi AE, Benabdelkrim A, Hammoutène D, Ben Yaghlane S, Abdallah HH, Ben Said R, Linguerri R, Hochlaf M. Theoretical Characterization of the Structure and Spectroscopy of HCNO 2 Isomers and Applications. J Phys Chem A 2020; 124:11061-11071. [PMID: 33347316 DOI: 10.1021/acs.jpca.0c09027] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We carried out a theoretical, fully ab initio, investigation of the stable forms of the [H,C,N,O,O] pentatomic molecular system, whose isomers are involved in fundamental combustion and atmospheric processes and are of potential interest for astrophysics. By adopting the MP2 and CCSD(T) electronic structure methods, combined with extrapolations to the complete basis set (CBS) limit, we characterized 20 low-energy isomers, excluding weak van der Waals complexes. For these molecules, we determined a set of geometrical parameters, relative energies, anharmonic vibrational frequencies, IR intensities, and fragmentation/formation energies from various atomic and/or molecular fragments. We discuss the relevance of the present findings for the search of new molecular species in astrophysical and atmospheric media and give suggestions for their possible detection in laboratory experiments. The set of data provided by the present work should facilitate the identification of these species from their gas-phase and low-temperature solid matrix spectra, whenever measured.
Collapse
Affiliation(s)
- Ala Eddine Tourchi
- Faculty of Chemistry, Laboratory of Thermodynamics and Molecular Modeling, USTHB, BP 32, El Alia, 16111, Bab Ezzouar, Algiers, Algeria.,COSYS/LISIS, Université Gustave Eiffel, 5 Bd Descartes 77454, Champs sur Marne, France
| | - Aicha Benabdelkrim
- Faculty of Chemistry, Laboratory of Thermodynamics and Molecular Modeling, USTHB, BP 32, El Alia, 16111, Bab Ezzouar, Algiers, Algeria.,COSYS/LISIS, Université Gustave Eiffel, 5 Bd Descartes 77454, Champs sur Marne, France
| | - Dalila Hammoutène
- Faculty of Chemistry, Laboratory of Thermodynamics and Molecular Modeling, USTHB, BP 32, El Alia, 16111, Bab Ezzouar, Algiers, Algeria
| | - Saida Ben Yaghlane
- Laboratoire de Spectroscopie Atomique, Moléculaire et Applications LSAMA, Université de Tunis El Manar, Tunis, Tunisia
| | - Hassan H Abdallah
- Department of Chemistry, Salahaddin University-Erbil, Erbil 44001, Iraq
| | - Ridha Ben Said
- Department of Chemistry, College of Science and Arts, Qassim University, Ar Rass, Saudi Arabia
| | - Roberto Linguerri
- COSYS/LISIS, Université Gustave Eiffel, 5 Bd Descartes 77454, Champs sur Marne, France
| | - Majdi Hochlaf
- COSYS/LISIS, Université Gustave Eiffel, 5 Bd Descartes 77454, Champs sur Marne, France
| |
Collapse
|
21
|
Mitra H, Roy TK. Comprehensive Benchmark Results for the Accuracy of Basis Sets for Anharmonic Molecular Vibrations. J Phys Chem A 2020; 124:9203-9221. [DOI: 10.1021/acs.jpca.0c06634] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Hrishit Mitra
- Department of Chemistry and Chemical Sciences, Central University of Jammu, Rahya-Suchani (Bagla), Jammu, Jammu and Kashmir 181143, India
| | - Tapta Kanchan Roy
- Department of Chemistry and Chemical Sciences, Central University of Jammu, Rahya-Suchani (Bagla), Jammu, Jammu and Kashmir 181143, India
| |
Collapse
|
22
|
Benabdelkrim A, Tourchi AE, Hammoutène D, Ben Yaghlane S, Abdallah HH, Linguerri R, Hochlaf M. Characterization of the simplest sulfenyl thiocyanate: isomers, spectroscopy and implications of astrophysical and biological relevance. Phys Chem Chem Phys 2020; 22:17052-17061. [PMID: 32658239 DOI: 10.1039/d0cp02382e] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
Sulfenyl thiocyanate compounds, RSSCN, are involved in the human immune system biochemical processes. They are also the routes for the synthesis of complex S-containing species such as polypeptides, or symmetrical (RSSR) and unsymmetrical disulfides (RSSR'). At present, we have characterized the stable forms of the simplest sulfenyl thiocyanate compound, HSSCN, at the coupled cluster level. We found twenty-three isomers, for which we determined a set of structural parameters, anharmonic frequencies and reaction energies for the formation of the corresponding diatomic + triatomic and atomic + tetratomic fragments. We also discussed the implications of the present findings for biological entities containing a disulfide bridge, where we identified three isomers that may serve as prototypes. Similarities and differences with other S/N hybrid bioactive molecules are also discussed. From an astrophysical point of view, we expect HSSCN isomers to be present in astrophysical media, since several of their molecular fragments have already been detected. In sum, the present set of data can be used for the identification of HSSCN compounds and understanding the physical chemistry of sulfur containing molecules in vivo, in the laboratory and in astrophysical media.
Collapse
Affiliation(s)
- Aicha Benabdelkrim
- USTHB, Faculty of Chemistry, Laboratory of Thermodynamics and Molecular Modeling, BP 32, Al Alia, 16111, Bab Ezzouar, Algiers, Algeria
| | | | | | | | | | | | | |
Collapse
|
23
|
Palivec V, Michal P, Kapitán J, Martinez‐Seara H, Bouř P. Raman Optical Activity of Glucose and Sorbose in Extended Wavenumber Range. Chemphyschem 2020; 21:1272-1279. [DOI: 10.1002/cphc.202000261] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2020] [Revised: 04/22/2020] [Indexed: 12/31/2022]
Affiliation(s)
- Vladimír Palivec
- Institute of Organic Chemistry and BiochemistryAcademy of Sciences Flemingovo náměstí 2 16610 Prague Czech Republic
| | - Pavel Michal
- Department of OpticsPalacký University Olomouc 17. listopadu 12, 77146 Olomouc Czech Republic
| | - Josef Kapitán
- Department of OpticsPalacký University Olomouc 17. listopadu 12, 77146 Olomouc Czech Republic
| | - Hector Martinez‐Seara
- Institute of Organic Chemistry and BiochemistryAcademy of Sciences Flemingovo náměstí 2 16610 Prague Czech Republic
| | - Petr Bouř
- Institute of Organic Chemistry and BiochemistryAcademy of Sciences Flemingovo náměstí 2 16610 Prague Czech Republic
| |
Collapse
|
24
|
Carbonnière P, Erba A, Richter F, Dovesi R, Rerat M. Calculation of Anharmonic IR and Raman Intensities for Periodic Systems from DFT Calculations: Implementation and Validation. J Chem Theory Comput 2020; 16:3343-3351. [PMID: 32275427 DOI: 10.1021/acs.jctc.9b01061] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
An extension of the CRYSTAL program is presented allowing for calculations of anharmonic infrared (IR) intensities and Raman activities for periodic systems. This work is a follow-up of two papers devoted to the computation of anharmonic vibrational states of solids from DFT (density functional theory) calculations (Erba et al. J. Chem. Theory Comput. 2019, 15, 3755-3765 and Erba et al. J. Chem. Theory Comput. 2019, 15, 3766-3777). The approach presented here relies on the evaluation of integrals of the dipole moment and polarizability operators over anharmonic wave functions obtained from either VSCF (vibrational self-consistent field) or VCI (vibrational configuration interaction) calculations. With this extension, the program now allows for a more complete characterization of the vibrational spectroscopic features of solids within the density functional theory. In particular, it is able (i) to provide reliable positions and intensities for the most intense spectral features and (ii) to check whether a first overtone or a combination band has a nonvanishing IR intensity or Raman activity. Therefore, it becomes possible to assign the transition(s) corresponding to satellite peak(s) around a fundamental transition or the overtones or combination bands that may be as intense as their corresponding fundamental transitions through the strongest mode-mode couplings, as in so-called Fermi resonances. The present method is assessed on two molecular systems, H2O and H2CO, as well as on two solid state cases, boron hydrides BH4 and their deuterated species BD4 in a crystalline environment of alkali metals (M = Na, K). The solid state cases are particularly insightful as, in the B-H (or B-D) stretching region here considered, they exhibit many spectral features entirely due to anharmonic effects: two out of three in the IR spectrum and four out of six in the Raman spectrum. All IR and Raman active overtones and combination bands experimentally observed are correctly predicted with our approach. The effect of the adopted quantum-chemical model (DFT exchange-correlation functional/basis set) for the electronic structure calculations on the computed spectra is discussed and found to be significant, which suggests some special care is needed for the analysis of subtle spectral features.
Collapse
Affiliation(s)
- P Carbonnière
- IPREM, Université de Pau et des Pays de l'Adour, IPREM-CAPT UMR CNRS 5254, Hélioparc Pau Pyrénées, 2 avenue du Président Angot, 64053 Pau Cedex 09, Pau, France
| | - A Erba
- Dipartimento di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - F Richter
- IPREM, Université de Pau et des Pays de l'Adour, IPREM-CAPT UMR CNRS 5254, Hélioparc Pau Pyrénées, 2 avenue du Président Angot, 64053 Pau Cedex 09, Pau, France
| | - R Dovesi
- Dipartimento di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - M Rerat
- IPREM, Université de Pau et des Pays de l'Adour, IPREM-CAPT UMR CNRS 5254, Hélioparc Pau Pyrénées, 2 avenue du Président Angot, 64053 Pau Cedex 09, Pau, France
| |
Collapse
|
25
|
Green JA, Improta R. Vibrations of the guanine-cytosine pair in chloroform: an anharmonic computational study. Phys Chem Chem Phys 2020; 22:5509-5522. [PMID: 32104818 DOI: 10.1039/c9cp06373k] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
We compute at the anharmonic level the vibrational spectra of the Watson-Crick dimer formed by guanosine (G) and cytidine (C) in chloroform, together with those of G, C and the most populated GG dimer. The spectra for deuterated and partially deuterated GC are also computed. We use DFT calculations, with B3LYP and CAM-B3LYP as reference functionals. Solvent effects from chloroform are included via the Polarizable Continuum Model (PCM), and by performing tests on models including up two chloroform molecules. Both B3LYP and CAM-B3LYP calculations reproduce the shape of the experimental spectra well in the fingerprint region (1500-1700 cm-1) and in the N-H stretching region (2800-3600 cm-1), with B3LYP providing better quantitative agreement with experiments. According to our calculations, the N-H amido streching mode of G falls at ∼2900 cm-1, while the N-H amino of G and C falls at ∼3100 cm-1 when hydrogen-bonded, or ∼3500 cm-1 when free. Overtone and combination bands strongly contribute to the absorption band at ∼3300 cm-1. Inclusion of bulk solvent effects significantly increases the accuracy of the computed spectra, while solute-solvent interactions have a smaller, though still noticeable, effect. Some key aspects of the anharmonic treatment of strongly vibrationally coupled supermolecular systems and the related methodological issues are also discussed.
Collapse
Affiliation(s)
- James A Green
- Istituto di Biostrutture e Bioimmagini-CNR, Via Mezzocannone 16, I-80134 Napoli, Italy.
| | - Roberto Improta
- Istituto di Biostrutture e Bioimmagini-CNR, Via Mezzocannone 16, I-80134 Napoli, Italy.
| |
Collapse
|
26
|
Pouchan C. Quantum molecular dynamic approach to modeling the vibrational spectra of propynenitrile and its 1:1 complexes with water. COMPUT THEOR CHEM 2019. [DOI: 10.1016/j.comptc.2019.06.004] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
27
|
Puzzarini C, Bloino J, Tasinato N, Barone V. Accuracy and Interpretability: The Devil and the Holy Grail. New Routes across Old Boundaries in Computational Spectroscopy. Chem Rev 2019; 119:8131-8191. [DOI: 10.1021/acs.chemrev.9b00007] [Citation(s) in RCA: 114] [Impact Index Per Article: 22.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Affiliation(s)
- Cristina Puzzarini
- Dipartimento di Chimica “Giacomo Ciamician”, Università di Bologna, Via F. Selmi 2, I-40126 Bologna, Italy
| | - Julien Bloino
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| | - Nicola Tasinato
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| | - Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| |
Collapse
|
28
|
Erba A, Maul J, Ferrabone M, Dovesi R, Rérat M, Carbonnière P. Anharmonic Vibrational States of Solids from DFT Calculations. Part II: Implementation of the VSCF and VCI Methods. J Chem Theory Comput 2019; 15:3766-3777. [DOI: 10.1021/acs.jctc.9b00294] [Citation(s) in RCA: 26] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Alessandro Erba
- Dipartimento di Chimica, Universitá di Torino, via Giuria 5, 10125 Torino, Italy
| | - Jefferson Maul
- Dipartimento di Chimica, Universitá di Torino, via Giuria 5, 10125 Torino, Italy
| | - Matteo Ferrabone
- Dipartimento di Chimica, Universitá di Torino, via Giuria 5, 10125 Torino, Italy
| | - Roberto Dovesi
- Dipartimento di Chimica, Universitá di Torino, via Giuria 5, 10125 Torino, Italy
| | - Michel Rérat
- IPREM, Université de Pau et des Pays de l’Adour, IPREM-CAPT UMR CNRS 5254, Hélioparc Pau Pyrénées, 2 avenue du Président Angot, 64053 PAU CEDEX 9, Pau, France
| | - Philippe Carbonnière
- IPREM, Université de Pau et des Pays de l’Adour, IPREM-CAPT UMR CNRS 5254, Hélioparc Pau Pyrénées, 2 avenue du Président Angot, 64053 PAU CEDEX 9, Pau, France
| |
Collapse
|
29
|
DFT modelling of the infrared spectra for the isolated and the micro-hydrated forms of uracil. Theor Chem Acc 2019. [DOI: 10.1007/s00214-019-2431-4] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
|
30
|
Ionic reaction products of iodine with pyridine, 4-methylpyridine, and 4-tert-butylpyridine in a polyethylene matrix. A FTIR polarization spectroscopic investigation. Chem Phys Lett 2019. [DOI: 10.1016/j.cplett.2018.12.010] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|
31
|
Krasnoshchekov SV, Schutski RS, Craig NC, Sibaev M, Crittenden DL. Comparing the accuracy of perturbative and variational calculations for predicting fundamental vibrational frequencies of dihalomethanes. J Chem Phys 2018; 148:084102. [DOI: 10.1063/1.5020295] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Sergey V. Krasnoshchekov
- Department of Chemistry, Lomonosov Moscow State University, Leninskiye Gory 1-3, Moscow 119991, Russian Federation
| | | | - Norman C. Craig
- Department of Chemistry and Biochemistry, Oberlin College, Oberlin, Ohio 44074, USA
| | - Marat Sibaev
- Research School of Chemistry, Australian National University, Canberra, ACT 2601, Australia
| | | |
Collapse
|
32
|
Ravichandran L, Banik S. Performance of different density functionals for the calculation of vibrational frequencies with vibrational coupled cluster method in bosonic representation. Theor Chem Acc 2017. [DOI: 10.1007/s00214-017-2177-9] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
33
|
Jungwirth J, Šebestík J, Šafařík M, Kapitán J, Bouř P. Quantitative Determination of Ala-Ala Conformer Ratios in Solution by Decomposition of Raman Optical Activity Spectra. J Phys Chem B 2017; 121:8956-8964. [DOI: 10.1021/acs.jpcb.7b07154] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Jakub Jungwirth
- Institute
of Organic Chemistry and Biochemistry, Academy of Sciences, Flemingovo
náměstí 2, 16610 Prague, Czech Republic
- Faculty
of Mathematics and Physics, Charles University, Ke Karlovu 3, 12116 Prague, Czech Republic
| | - Jaroslav Šebestík
- Institute
of Organic Chemistry and Biochemistry, Academy of Sciences, Flemingovo
náměstí 2, 16610 Prague, Czech Republic
| | - Martin Šafařík
- Institute
of Organic Chemistry and Biochemistry, Academy of Sciences, Flemingovo
náměstí 2, 16610 Prague, Czech Republic
| | - Josef Kapitán
- Department
of Optics, Palacký University, 17. listopadu 12, 77146 Olomouc, Czech Republic
| | - Petr Bouř
- Institute
of Organic Chemistry and Biochemistry, Academy of Sciences, Flemingovo
náměstí 2, 16610 Prague, Czech Republic
| |
Collapse
|
34
|
Spectroscopic Characterization of Key Aromatic and Heterocyclic Molecules: A Route toward the Origin of Life. ACTA ACUST UNITED AC 2017; 154. [DOI: 10.3847/1538-3881/aa7d54] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
|
35
|
Hansen PE, Spanget-Larsen J. NMR and IR Investigations of Strong Intramolecular Hydrogen Bonds. Molecules 2017; 22:E552. [PMID: 28353675 PMCID: PMC6154318 DOI: 10.3390/molecules22040552] [Citation(s) in RCA: 70] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2017] [Revised: 03/20/2017] [Accepted: 03/24/2017] [Indexed: 11/24/2022] Open
Abstract
For the purpose of this review, strong hydrogen bonds have been defined on the basis of experimental data, such as OH stretching wavenumbers, νOH, and OH chemical shifts, δOH (in the latter case, after correction for ring current effects). Limits for O-H···Y systems are taken as 2800 > νOH > 1800 cm-1, and 19 ppm > δOH > 15 ppm. Recent results as well as an account of theoretical advances are presented for a series of important classes of compounds such as β-diketone enols, β-thioxoketone enols, Mannich bases, proton sponges, quinoline N-oxides and diacid anions. The O···O distance has long been used as a parameter for hydrogen bond strength in O-H···O systems. On a broad scale, a correlation between OH stretching wavenumbers and O···O distances is observed, as demonstrated experimentally as well as theoretically, but for substituted β-diketone enols this correlation is relatively weak.
Collapse
Affiliation(s)
- Poul Erik Hansen
- Department of Science and Environment, Roskilde University, Universitetsvej 1, P.O. Box 260, DK-4000 Roskilde, Denmark.
| | - Jens Spanget-Larsen
- Department of Science and Environment, Roskilde University, Universitetsvej 1, P.O. Box 260, DK-4000 Roskilde, Denmark.
| |
Collapse
|
36
|
Thicoipe S, Carbonniere P, Pouchan C. DFT modelling of the infrared spectra for isolated and aqueous forms of adenine. Theor Chem Acc 2017. [DOI: 10.1007/s00214-017-2076-0] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
37
|
Intramolecular hydrogen bonding in myricetin and myricitrin. Quantum chemical calculations and vibrational spectroscopy. J Mol Struct 2017. [DOI: 10.1016/j.molstruc.2016.11.069] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
38
|
Infrared Investigations of the Neutral-Ionic Phase Transition in TTF-CA and Its Dynamics. CRYSTALS 2017. [DOI: 10.3390/cryst7010017] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
39
|
Cannelli O, Giovannini T, Baiardi A, Carlotti B, Elisei F, Cappelli C. Understanding the interplay between the solvent and nuclear rearrangements in the negative solvatochromism of a push–pull flexible quinolinium cation. Phys Chem Chem Phys 2017; 19:32544-32555. [DOI: 10.1039/c7cp07104c] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The main effects (solvation, vibronic progression) affecting the band position and shape of a push–pull flexible quinolinium cation OPA are highlighted.
Collapse
Affiliation(s)
- Oliviero Cannelli
- Scuola Normale Superiore
- Pisa
- Italy
- Laboratoire de Spectroscopie Ultrarapide
- ISIC-FSB
| | | | | | - Benedetta Carlotti
- Department of Chemistry
- Biology and Biotechnology and Centro di Eccellenza sui Materiali Innovativi Nanostrutturati (CEMIN)
- University of Perugia
- via Elce di Sotto 8
- 06123 Perugia
| | - Fausto Elisei
- Department of Chemistry
- Biology and Biotechnology and Centro di Eccellenza sui Materiali Innovativi Nanostrutturati (CEMIN)
- University of Perugia
- via Elce di Sotto 8
- 06123 Perugia
| | | |
Collapse
|
40
|
Ruggiero MT, Zeitler JA. Resolving the Origins of Crystalline Anharmonicity Using Terahertz Time-Domain Spectroscopy and ab Initio Simulations. J Phys Chem B 2016; 120:11733-11739. [PMID: 27766874 DOI: 10.1021/acs.jpcb.6b10248] [Citation(s) in RCA: 35] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Anharmonicity has been shown to be an important piece of the fundamental framework that dictates numerous observable phenomena. In particular, anharmonicity is the driving force of vibrational relaxation processes, mechanisms that are integral to the proper function of numerous chemical processes. However, elucidating its origins has proven difficult due to experimental and theoretical challenges, specifically related to separating the anharmonic contributions from other unrelated effects. While no one technique is particularly suited for providing a complete picture of anharmonicity, by combining multiple complementary methods such a characterization can be made. In this study the role of individual atomic interactions on the anharmonic properties of crystalline purine, the building block of many DNA and RNA nucleobases, is studied by experimental terahertz time-domain spectroscopy and first-principles density functional theory (DFT) and ab initio molecular dynamics simulations (AIMD). In particular, the detailed vibrational information provided by the DFT calculations is used to interpret the atomic origins of anharmonic-related effects as determined by the AIMD calculations, which are in good agreement with the experimental data. The results highlight that anharmonicity is especially pronounced in the intermolecular interactions, particularly along the amine hydrogen bond coordinate, and yields valuable insight into what is similarly observed complex biosystems and crystalline solids.
Collapse
Affiliation(s)
- Michael T Ruggiero
- Department of Chemical Engineering and Biotechnology, University of Cambridge , Philippa Fawcett Drive, Cambridge CB3 0AS, United Kingdom
| | - J Axel Zeitler
- Department of Chemical Engineering and Biotechnology, University of Cambridge , Philippa Fawcett Drive, Cambridge CB3 0AS, United Kingdom
| |
Collapse
|
41
|
Barnes L, Schindler B, Compagnon I, Allouche AR. Fast and accurate hybrid QM//MM approach for computing anharmonic corrections to vibrational frequencies. J Mol Model 2016; 22:285. [DOI: 10.1007/s00894-016-3135-5] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/03/2016] [Accepted: 10/09/2016] [Indexed: 10/20/2022]
|
42
|
Shiekh BA, Kaur D, Seth B, Mahajan S. The theoretical-cum-statistical approach for the investigation of reaction NO2+ O(3P) → NO + O2 using SCTST and a full anharmonic VPT2 model. Chem Phys Lett 2016. [DOI: 10.1016/j.cplett.2016.08.058] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
|
43
|
On the choice electronic structure method to calculate the quartic potential energy surface for the vibrational calculation of polyatomic molecules. Theor Chem Acc 2016. [DOI: 10.1007/s00214-016-1962-1] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
44
|
Fornaro T, Biczysko M, Bloino J, Barone V. Reliable vibrational wavenumbers for C=O and N-H stretchings of isolated and hydrogen-bonded nucleic acid bases. Phys Chem Chem Phys 2016; 18:8479-90. [PMID: 26940362 PMCID: PMC5612391 DOI: 10.1039/c5cp07386c] [Citation(s) in RCA: 37] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
The accurate prediction of vibrational wavenumbers for functional groups involved in hydrogen-bonded bridges remains an important challenge for computational spectroscopy. For the specific case of the C=O and N-H stretching modes of nucleobases and their oligomers, the paucity of experimental reference values needs to be compensated by reliable computational data, which require the use of approaches going beyond the standard harmonic oscillator model. Test computations performed for model systems (formamide, acetamide and their cyclic homodimers) in the framework of the second order vibrational perturbation theory (VPT2) confirmed that anharmonic corrections can be safely computed by global hybrid (GHF) or double hybrid (DHF) functionals, whereas the harmonic part is particularly challenging. As a matter of fact, GHFs perform quite poorly and even DHFs, while fully satisfactory for C=O stretchings, face unexpected difficulties when dealing with N-H stretchings. On these grounds, a linear regression for N-H stretchings has been obtained and validated for the heterodimers formed by 4-aminopyrimidine with 6-methyl-4-pyrimidinone (4APM-M4PMN) and by uracil with water. In view of the good performance of this computational model, we have built a training set of B2PLYP-D3/maug-cc-pVTZ harmonic wavenumbers (including linear regression scaling for N-H) for six-different uracil dimers and a validation set including 4APM-M4PMN, one of the most stable hydrogen-bonded adenine homodimers, as well as the adenine-uracil, adenine-thymine, guanine-cytosine and adenine-4-thiouracil heterodimers. Because of the unfavourable scaling of DHF harmonic wavenumbers with the dimensions of the investigated systems, we have optimized a linear regression of B3LYP-D3/N07D harmonic wavenumbers for the training set, which has been next checked against the validation set. This relatively cheap model, which shows very good agreement with experimental data (average errors of about 10 cm(-1)), paves the route toward a reliable analysis of spectroscopic signatures for larger polynucleotides.
Collapse
Affiliation(s)
- Teresa Fornaro
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy.
| | - Malgorzata Biczysko
- International Centre for Quantum and Molecular Structures, College of Sciences, Shanghai University, 99 Shangda Road, Shanghai, 200444, China.
| | - Julien Bloino
- Consiglio Nazionale delle Ricerche, Istituto di Chimica dei Composti OrganoMetallici (ICCOM-CNR), UOS di Pisa, Area della Ricerca CNR, Via G. Moruzzi 1, I-56124 Pisa, Italy
| | - Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy.
| |
Collapse
|
45
|
Cornaton Y, Ringholm M, Louant O, Ruud K. Analytic calculations of anharmonic infrared and Raman vibrational spectra. Phys Chem Chem Phys 2016; 18:4201-15. [PMID: 26784673 PMCID: PMC5063043 DOI: 10.1039/c5cp06657c] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2015] [Accepted: 01/05/2016] [Indexed: 11/21/2022]
Abstract
Using a recently developed recursive scheme for the calculation of high-order geometric derivatives of frequency-dependent molecular properties [Ringholm et al., J. Comp. Chem., 2014, 35, 622], we present the first analytic calculations of anharmonic infrared (IR) and Raman spectra including anharmonicity both in the vibrational frequencies and in the IR and Raman intensities. In the case of anharmonic corrections to the Raman intensities, this involves the calculation of fifth-order energy derivatives-that is, the third-order geometric derivatives of the frequency-dependent polarizability. The approach is applicable to both Hartree-Fock and Kohn-Sham density functional theory. Using generalized vibrational perturbation theory to second order, we have calculated the anharmonic infrared and Raman spectra of the non- and partially deuterated isotopomers of nitromethane, where the inclusion of anharmonic effects introduces combination and overtone bands that are observed in the experimental spectra. For the major features of the spectra, the inclusion of anharmonicities in the calculation of the vibrational frequencies is more important than anharmonic effects in the calculated infrared and Raman intensities. Using methanimine as a trial system, we demonstrate that the analytic approach avoids errors in the calculated spectra that may arise if numerical differentiation schemes are used.
Collapse
Affiliation(s)
- Yann Cornaton
- Centre for Theoretical and Computational Chemistry, Department of Chemistry, University of Tromsø-The Arctic University of Norway, N-9037 Tromsø, Norway.
| | - Magnus Ringholm
- Centre for Theoretical and Computational Chemistry, Department of Chemistry, University of Tromsø-The Arctic University of Norway, N-9037 Tromsø, Norway.
| | - Orian Louant
- Laboratory of Theoretical Chemistry, UCPTS, University of Namur, Rue de Bruxelles 61, B-5000 Namur, Belgium
| | - Kenneth Ruud
- Centre for Theoretical and Computational Chemistry, Department of Chemistry, University of Tromsø-The Arctic University of Norway, N-9037 Tromsø, Norway.
| |
Collapse
|
46
|
The role of torsional motion on the properties of propiolic acid and its H/D isotopic analogs: A density functional study using SCTST and a full anharmonic VPT2 model. Chem Phys Lett 2016. [DOI: 10.1016/j.cplett.2016.01.015] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
|
47
|
Al-Shalalfeh MM, Saleh TA, Al-Saadi AA. Silver colloid and film substrates in surface-enhanced Raman scattering for 2-thiouracil detection. RSC Adv 2016. [DOI: 10.1039/c6ra14832h] [Citation(s) in RCA: 72] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
Abstract
Surface-enhanced Raman scattering-based silver substrates were designed and fabricated for the detection of 2-thiouracil (2-TU).
Collapse
Affiliation(s)
- Mutasem M. Al-Shalalfeh
- Department of Chemistry
- King Fahd University of Petroleum & Minerals
- Dhahran 31261
- Saudi Arabia
| | - Tawfik A. Saleh
- Department of Chemistry
- King Fahd University of Petroleum & Minerals
- Dhahran 31261
- Saudi Arabia
| | - Abdulaziz A. Al-Saadi
- Department of Chemistry
- King Fahd University of Petroleum & Minerals
- Dhahran 31261
- Saudi Arabia
| |
Collapse
|
48
|
Jedidi A, Li R, Fornasiero P, Cavallo L, Carbonniere P. Vibrational Fingerprints of Low-Lying Pt(n)P(2n) (n = 1-5) Cluster Structures from Global Optimization Based on Density Functional Theory Potential Energy Surfaces. J Phys Chem A 2015; 119:11711-8. [PMID: 26566005 DOI: 10.1021/acs.jpca.5b08495] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Abstract
Vibrational fingerprints of small Pt(n)P(2n) (n = 1-5) clusters were computed from their low-lying structures located from a global exploration of their DFT potential energy surfaces with the GSAM code. Five DFT methods were assessed from the CCSD(T) wavenumbers of PtP2 species and CCSD relative energies of Pt2P4 structures. The eight first Pt(n)P(2n) isomers found are reported. The vibrational computations reveal (i) the absence of clear signatures made by overtone or combination bands due to very weak mechanical and electrical anharmonicities and (ii) some significant and recurrent vibrational fingerprints in correlation with the different PP bonding situations in the Pt(n)P(2n) structures.
Collapse
Affiliation(s)
- Abdesslem Jedidi
- Division of Physical Sciences and Engineering, KAUST Catalysis Center, King Abdullah University of Science and Technology (KAUST) , 4700 KAUST, Thuwal 23955-6900, Saudi Arabia
| | - Rui Li
- Groupe de Chimie Théorique et Réactivité, IPREM/ECP UMR CNRS 5254, Université de Pau et des Pays de l'Adour , F-64000 Pau, France
| | - Paolo Fornasiero
- Department of Chemical and Pharmaceutical Sciences, INSTM, ICCOM-CNR, University of Trieste , Via L. Giorgieri 1, 34127 Trieste, Italy
| | - Luigi Cavallo
- Division of Physical Sciences and Engineering, KAUST Catalysis Center, King Abdullah University of Science and Technology (KAUST) , 4700 KAUST, Thuwal 23955-6900, Saudi Arabia
| | - Philippe Carbonniere
- Groupe de Chimie Théorique et Réactivité, IPREM/ECP UMR CNRS 5254, Université de Pau et des Pays de l'Adour , F-64000 Pau, France
| |
Collapse
|
49
|
Puzzarini C, Barone V. Theoretical Study of the X2NO Systems (X = F, Cl, Br, I): Effects of Halogen Substitution on Structural and Spectroscopic Properties. J Chem Theory Comput 2015; 5:2378-87. [PMID: 26616619 DOI: 10.1021/ct9001762] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Structural and spectroscopic properties of the X2 NO series of radicals, with X = F, Cl, Br, I, have been computed by the coupled cluster ansatz in conjunction with hierarchical series of basis sets, accounting for, in most cases, core correlation effects and extrapolation to the complete basis set limit. Namely, equilibrium structures, vibrational frequencies, and hyperfine coupling constants have been considered. Methods rooted into the density functional theory have been used to estimate anharmonic and, in conjunction with the polarizable continuum model, environmental effects. The remarkable agreement with the available experimental data, limited to the lighter member of the series, confirms the reliability of our computational approach and suggests that the data for heavier species represent reliable benchmarks for future experimental data and/or cheaper computational methods devised for larger systems.
Collapse
Affiliation(s)
- Cristina Puzzarini
- Dipartimento di Chimica "G. Ciamician", Università di Bologna, Via F. Selmi 2, 40126 Bologna, Italy, and Scuola Normale Superiore di Pisa, Piazza dei Cavalieri 7, 56126 Pisa, Italy, IPCF-CNR, via Moruzzi 1, 56124 Pisa, Italy
| | - Vincenzo Barone
- Dipartimento di Chimica "G. Ciamician", Università di Bologna, Via F. Selmi 2, 40126 Bologna, Italy, and Scuola Normale Superiore di Pisa, Piazza dei Cavalieri 7, 56126 Pisa, Italy, IPCF-CNR, via Moruzzi 1, 56124 Pisa, Italy
| |
Collapse
|
50
|
Cappelli C, Monti S, Scalmani G, Barone V. On the Calculation of Vibrational Frequencies for Molecules in Solution Beyond the Harmonic Approximation. J Chem Theory Comput 2015; 6:1660-9. [PMID: 26615698 DOI: 10.1021/ct100048g] [Citation(s) in RCA: 49] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We report some results on the calculation of vibrational spectra of molecules in condensed phase with accounting simultaneously for anharmonicity and solute-solvent interactions, the latter being described by means of the polarizable continuum model (PCM). Density functional theory force fields are employed as well as a new implementation of the PCM cavity and its derivatives. The results obtained for formaldehyde and simple peptide prototypes show that our approach is able to yield a quantitative agreement with experiments for vacuo-to-solvent harmonic and anharmonic frequency shifts.
Collapse
Affiliation(s)
- Chiara Cappelli
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via Risorgimento 35, I-56126 Pisa, Italy, Istituto per i Processi Chimico-Fisici del Consiglio Nazionale delle Ricerche, via Moruzzi,1 I-56124 Pisa, Italy, Gaussian, Inc., 340 Quinnipiac Street, Building 40, Wallingford, Connecticut 06492, and Scuola Normale Superiore, Piazza dei Cavalieri, 7 I-56126 Pisa, Italy
| | - Susanna Monti
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via Risorgimento 35, I-56126 Pisa, Italy, Istituto per i Processi Chimico-Fisici del Consiglio Nazionale delle Ricerche, via Moruzzi,1 I-56124 Pisa, Italy, Gaussian, Inc., 340 Quinnipiac Street, Building 40, Wallingford, Connecticut 06492, and Scuola Normale Superiore, Piazza dei Cavalieri, 7 I-56126 Pisa, Italy
| | - Giovanni Scalmani
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via Risorgimento 35, I-56126 Pisa, Italy, Istituto per i Processi Chimico-Fisici del Consiglio Nazionale delle Ricerche, via Moruzzi,1 I-56124 Pisa, Italy, Gaussian, Inc., 340 Quinnipiac Street, Building 40, Wallingford, Connecticut 06492, and Scuola Normale Superiore, Piazza dei Cavalieri, 7 I-56126 Pisa, Italy
| | - Vincenzo Barone
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via Risorgimento 35, I-56126 Pisa, Italy, Istituto per i Processi Chimico-Fisici del Consiglio Nazionale delle Ricerche, via Moruzzi,1 I-56124 Pisa, Italy, Gaussian, Inc., 340 Quinnipiac Street, Building 40, Wallingford, Connecticut 06492, and Scuola Normale Superiore, Piazza dei Cavalieri, 7 I-56126 Pisa, Italy
| |
Collapse
|