1
|
Mendolicchio M, Bloino J, Barone V. Perturb-Then-Diagonalize Vibrational Engine Exploiting Curvilinear Internal Coordinates. J Chem Theory Comput 2022; 18:7603-7619. [PMID: 36322968 DOI: 10.1021/acs.jctc.2c00773] [Citation(s) in RCA: 14] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
The present paper is devoted to the implementation and validation of a second-order perturbative approach to anharmonic vibrations, followed by variational treatment of strong couplings (GVPT2) based on curvilinear internal coordinates. The main difference with respect to the customary Cartesian-based formulation is that the kinetic energy operator is no longer diagonal, and has to be expanded as well, leading to additional terms which have to be taken into proper account. It is, however, possible to recast all the equations as well-defined generalizations of the corresponding Cartesian-based counterparts, thus achieving a remarkable simplification of the new implementation. Particular attention is paid to the treatment of Fermi resonances with significant number of test cases analyzed fully, validating the new implementation. The results obtained in this work confirm that curvilinear coordinates strongly reduce the strength of inter-mode couplings compared to their Cartesian counterparts. This increases the reliability of low-order perturbative treatments for semi-rigid molecules and paves the way toward the reliable representation of more flexible molecules where small- and large-amplitude motions can be safely decoupled and treated at different levels of theory.
Collapse
Affiliation(s)
- Marco Mendolicchio
- Scuola Superiore Meridionale, Largo S. Marcellino 10, Napoli I-80138, Italy
| | - Julien Bloino
- Scuola Normale Superiore, Piazza dei Cavalieri 7, Pisa I-56126, Italy
| | - Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, Pisa I-56126, Italy
| |
Collapse
|
2
|
Hino K, Kurashige Y. Matrix Product State Formulation of the MCTDH Theory in Local Mode Representations for Anharmonic Potentials. J Chem Theory Comput 2022; 18:3347-3356. [PMID: 35606892 DOI: 10.1021/acs.jctc.2c00243] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The matrix product state formulation of the multiconfiguration time-dependent Hartree theory, MPS-MCTDH, reported previously [Kurashige, J. Chem. Phys. 2018, 19, 194114] is extended to realistic anharmonic potentials with n-mode representations beyond the linear vibronic coupling model. For realistic vibrational potentials, the local mode representation should give a more compact representation of the potentials, i.e., lowering the dimensionality of the entanglements, than the normal coordinates, and the MPS-MCTDH formulation should work more efficiently and maintain the accuracy with a small bond dimension of the MPS ansatz. In fact, it was confirmed that the use of the local coordinates made the interaction matrices diagonal dominant and the number of terms in the n-body expansion of the potentials was significantly reduced. The method was applied to the IR spectrum of the CH2O molecule, the zero-point energies, and the vibrational energy redistribution dynamics of polyenes C2nH2n+2. The results showed that the efficiency of the MPS-MCTDH method is significantly accelerated by the use of local coordinates even if the long-range interactions are included in the potential.
Collapse
Affiliation(s)
- Kentaro Hino
- Department of Chemistry, Graduate School of Science, Kyoto University, Kitashirakawa Oiwake-cho, Sakyo-ku, Kyoto 606-8502, Japan
| | - Yuki Kurashige
- Department of Chemistry, Graduate School of Science, Kyoto University, Kitashirakawa Oiwake-cho, Sakyo-ku, Kyoto 606-8502, Japan
| |
Collapse
|
3
|
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
|
4
|
Donati G, Petrone A, Rega N. Multiresolution continuous wavelet transform for studying coupled solute-solvent vibrations via ab initio molecular dynamics. Phys Chem Chem Phys 2020; 22:22645-22661. [PMID: 33015693 DOI: 10.1039/d0cp02495c] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
Vibrational analysis in solution and the theoretical determination of infrared and Raman spectra are of key importance in many fields of chemical interest. Vibrational band dynamics of molecules and their sensitivity to the environment can also be captured by these spectroscopies in their time dependent version. However, it is often difficult to provide an interpretation of the experimental data at the molecular scale, such as molecular mechanisms or the processes hidden behind them. In this work, we present a theoretical-computational protocol based on ab initio molecular dynamics simulations and a combination of normal-like (generalized) mode analysis of solute-solvent clusters with a wavelet transform, for the first time. The case study is the vibrational dynamics of N-methyl-acetamide (NMA) in water solution, a well-known model of hydration of peptides and proteins. Amide modes are typical bands of peptide and protein backbone, and their couplings with the environment are very challenging in terms of the accurate prediction of solvent induced intensity and frequency shifts. The contribution of water molecules surrounding NMA to the composition of generalized and time resolved modes is introduced in our vibrational analysis, showing unequivocally its influence on the amide mode spectra. It is also shown that such mode compositions need the inclusion of the first shell solvent molecules to be accurately described. The wavelet analysis is proven to be strongly recommended to follow the time evolution of the spectra, and to capture vibrational band couplings and frequency shifts over time, preserving at the same time a well-balanced time-frequency resolution. This peculiar feature also allows one to perform a combined structural-vibrational analysis, where the different strengths of hydrogen bond interactions can quantitatively affect the amide bands over time at finite temperature. The proposed method allows for the direct connection between vibrational modes and local structural changes, providing a link from the spectroscopic observable to the structure, in this case the peptide backbone, and its hydration layouts.
Collapse
Affiliation(s)
- Greta Donati
- Department of Chemical Sciences, University of Napoli Federico II, Complesso Universitario di M. S. Angelo, Via Cintia, I-80126 Napoli, Italy.
| | | | | |
Collapse
|
5
|
Matthews DA, Cheng L, Harding ME, Lipparini F, Stopkowicz S, Jagau TC, Szalay PG, Gauss J, Stanton JF. Coupled-cluster techniques for computational chemistry: The CFOUR program package. J Chem Phys 2020; 152:214108. [DOI: 10.1063/5.0004837] [Citation(s) in RCA: 214] [Impact Index Per Article: 53.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022] Open
Affiliation(s)
- Devin A. Matthews
- Department of Chemistry, Southern Methodist University, Dallas, Texas 75275, USA
| | - Lan Cheng
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, USA
| | - Michael E. Harding
- Institut für Physikalische Chemie, Karlsruher Institut für Technologie (KIT), Kaiserstr. 12, D-76131 Karlsruhe, Germany
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, I-56124 Pisa, Italy
| | - Stella Stopkowicz
- Department Chemie, Johannes Gutenberg-Universität Mainz, Duesbergweg 10-14, D-55128 Mainz, Germany
| | - Thomas-C. Jagau
- Department of Chemistry, University of Munich (LMU), Butenandtstr. 5-13, D-81377 Munich, Germany
| | - Péter G. Szalay
- ELTE Eötvös Loránd University, Institute of Chemistry, Pázmány Péter sétány 1/A, H-1117 Budapest, Hungary
| | - Jürgen Gauss
- Department Chemie, Johannes Gutenberg-Universität Mainz, Duesbergweg 10-14, D-55128 Mainz, Germany
| | - John F. Stanton
- Quantum Theory Project, Departments of Chemistry and Physics, University of Florida, Gainesville, Florida 32611, USA
| |
Collapse
|
6
|
Patel AHG, Mohammed AAK, Limacher PA, Ayers PW. Finite Field Method for Nonlinear Optical Property Prediction Using Rational Function Approximants. J Phys Chem A 2017. [PMID: 28632384 DOI: 10.1021/acs.jpca.7b04049] [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/28/2022]
Abstract
The finite field (FF) method is a quick, easy-to-implement tool for the prediction of nonlinear optical properties. Here, we present and explore a novel variant of the FF method, which uses a rational function to fit a molecule's energy with respect to an electric field. Similarly to previous FF methods, factors crucial for the method's accuracy were tuned. These factors include the number of terms in the function, the distribution of fields used to construct the approximation, and the initial field in the approximation. It was found that the approximant form that best fits the energy has four numerator terms and three denominator terms. To determine a reasonable field distribution, the common ratio of a geometric progression was optimized to √2. Finally, an algorithm for determining a good initial field guess was devised. The optimized FF method was used to compute the polarizability and second hyperpolarizability for a set of 121 molecules and the first hyperpolarizability for a set of 91 molecules. The results from this were compared to a previous polynomial-based FF method. It was found that using a rational function gives higher errors compared to the polynomial model. However, unlike the polynomial model, no subsequent refinement steps were needed to obtain usable results. An overall comparison of the behavior of the two methods also shows that the rational function is less sensitive to the chosen initial field, making it a good choice for new quantum chemistry codes.
Collapse
Affiliation(s)
- Anand H G Patel
- Department of Chemistry & Chemical Biology, McMaster University 1280 Main Street West, Hamilton, Ontario, Canada L8S 4L8
| | - Ahmed A K Mohammed
- Department of Chemistry & Chemical Biology, McMaster University 1280 Main Street West, Hamilton, Ontario, Canada L8S 4L8
| | - Peter A Limacher
- Institute of Physical Chemistry, Karlsruhe Institute of Technology , 76131 Karlsruhe, Germany
| | - Paul W Ayers
- Department of Chemistry & Chemical Biology, McMaster University 1280 Main Street West, Hamilton, Ontario, Canada L8S 4L8
| |
Collapse
|
7
|
Faber R, Kaminsky J, Sauer SPA. Rovibrational and Temperature Effects in Theoretical Studies of NMR Parameters. GAS PHASE NMR 2016. [DOI: 10.1039/9781782623816-00218] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Abstract
The demand for high precision calculations of NMR shieldings (or their related values, chemical shifts δ) and spin–spin coupling constants facilitating and supporting detailed interpretations of NMR spectra increases hand in hand with the development of computational techniques and hardware resources. Highly sophisticated calculations including even relativistic effects are nowadays possible for these properties. However, NMR parameters depend not only on molecular structure and environment but also on molecular flexibility and temperature and the apparent success of theoretical predictions for molecular equilibrium geometries creates a demand for zero-point vibrational and temperature corrections. In this chapter we describe briefly the theory behind rovibrational corrections and review then some important contributions to this field.
Collapse
Affiliation(s)
- Rasmus Faber
- Department of Chemistry, University of Copenhagen Universitetsparken 5 DK-2100 Copenhagen Ø Denmark
| | - Jakub Kaminsky
- Department of Molecular Spectroscopy, Institute of Organic Chemistry and Biochemistry 166 10 Prague Czech Republic
| | - Stephan P. A. Sauer
- Department of Chemistry, University of Copenhagen Universitetsparken 5 DK-2100 Copenhagen Ø Denmark
| |
Collapse
|
8
|
Barone V, Biczysko M, Bloino J, Cimino P, Penocchio E, Puzzarini C. CC/DFT Route toward Accurate Structures and Spectroscopic Features for Observed and Elusive Conformers of Flexible Molecules: Pyruvic Acid as a Case Study. J Chem Theory Comput 2015; 11:4342-63. [PMID: 26575928 PMCID: PMC5905675 DOI: 10.1021/acs.jctc.5b00580] [Citation(s) in RCA: 64] [Impact Index Per Article: 7.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The structures and relative stabilities as well as the rotational and vibrational spectra of the three low-energy conformers of pyruvic acid (PA) have been characterized using a state-of-the-art quantum-mechanical approach designed for flexible molecules. By making use of the available experimental rotational constants for several isotopologues of the most stable PA conformer, Tc-PA, the semiexperimental equilibrium structure has been derived. The latter provides a reference for the pure theoretical determination of the equilibrium geometries for all conformers, thus confirming for these structures an accuracy of 0.001 Å and 0.1 deg for bond lengths and angles, respectively. Highly accurate relative energies of all conformers (Tc-, Tt-, and Ct-PA) and of the transition states connecting them are provided along with the thermodynamic properties at low and high temperatures, thus leading to conformational enthalpies accurate to 1 kJ mol(-1). Concerning microwave spectroscopy, rotational constants accurate to about 20 MHz are provided for the Tt- and Ct-PA conformers, together with the computed centrifugal-distortion constants and dipole moments required to simulate their rotational spectra. For Ct-PA, vibrational frequencies in the mid-infrared region accurate to 10 cm(-1) are reported along with theoretical estimates for the transitions in the near-infrared range, and the corresponding infrared spectrum including fundamental transitions, overtones, and combination bands has been simulated. In addition to the new data described above, theoretical results for the Tc- and Tt-PA conformers are compared with all available experimental data to further confirm the accuracy of the hybrid coupled-cluster/density functional theory (CC/DFT) protocol applied in the present study. Finally, we discuss in detail the accuracy of computational models fully based on double-hybrid DFT functionals (mainly at the B2PLYP/aug-cc-pVTZ level) that avoid the use of very expensive CC calculations.
Collapse
Affiliation(s)
- Vincenzo Barone
- Scuola Normale Superiore,
Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| | - Malgorzata Biczysko
- 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
| | - 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
| | - Paola Cimino
- Dipartimento di Scienze
Farmaceutiche, Università degli Studi di
Salerno, via Ponte don Melillo, I-84084 Fisciano
(SA), Italy
| | | | - Cristina Puzzarini
- Dipartimento di Chimica
“Giacomo Ciamician”,
Università di Bologna, Via Selmi 2, I-40126
Bologna, Italy
| |
Collapse
|
9
|
Piccardo M, Bloino J, Barone V. Generalized Vibrational Perturbation Theory for Rotovibrational Energies of Linear, Symmetric and Asymmetric Tops: Theory, Approximations, and Automated Approaches to Deal with Medium-to-Large Molecular Systems. INTERNATIONAL JOURNAL OF QUANTUM CHEMISTRY 2015; 115:948-982. [PMID: 26345131 PMCID: PMC4553754 DOI: 10.1002/qua.24931] [Citation(s) in RCA: 72] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 03/02/2015] [Revised: 04/03/2015] [Accepted: 04/08/2015] [Indexed: 06/05/2023]
Abstract
Models going beyond the rigid-rotor and the harmonic oscillator levels are mandatory for providing accurate theoretical predictions for several spectroscopic properties. Different strategies have been devised for this purpose. Among them, the treatment by perturbation theory of the molecular Hamiltonian after its expansion in power series of products of vibrational and rotational operators, also referred to as vibrational perturbation theory (VPT), is particularly appealing for its computational efficiency to treat medium-to-large systems. Moreover, generalized (GVPT) strategies combining the use of perturbative and variational formalisms can be adopted to further improve the accuracy of the results, with the first approach used for weakly coupled terms, and the second one to handle tightly coupled ones. In this context, the GVPT formulation for asymmetric, symmetric, and linear tops is revisited and fully generalized to both minima and first-order saddle points of the molecular potential energy surface. The computational strategies and approximations that can be adopted in dealing with GVPT computations are pointed out, with a particular attention devoted to the treatment of symmetry and degeneracies. A number of tests and applications are discussed, to show the possibilities of the developments, as regards both the variety of treatable systems and eligible methods.
Collapse
Affiliation(s)
- Matteo Piccardo
- Scuola Normale Superiore Piazza dei Cavalieri 7, I-56126, Pisa, Italy E-mail:
| | - Julien Bloino
- Scuola Normale Superiore Piazza dei Cavalieri 7, I-56126, Pisa, Italy E-mail: ; Consiglio Nazionale delle Ricerche, Istituto di Chimica dei Composti Organometallici (CNR-ICCOM) UOS di Pisa Via G. Moruzzi, 1, I-56124, Italy
| | - Vincenzo Barone
- Scuola Normale Superiore Piazza dei Cavalieri 7, I-56126, Pisa, Italy E-mail:
| |
Collapse
|
10
|
Fornaro T, Burini D, Biczysko M, Barone V. Hydrogen-Bonding Effects on Infrared Spectra from Anharmonic Computations: Uracil–Water Complexes and Uracil Dimers. J Phys Chem A 2015; 119:4224-36. [DOI: 10.1021/acs.jpca.5b01561] [Citation(s) in RCA: 99] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Affiliation(s)
- Teresa Fornaro
- Scuola Normale Superiore, Piazza
dei Cavalieri 7, I-56126 Pisa, Italy
| | - Diletta Burini
- Dipartimento
di Matematica e Informatica, Università di Perugia, INFN Sezione
Perugia Via Vanvitelli, I-106123 Perugia, Italy
| | - Malgorzata Biczysko
- Physics
Department, and International Centre for Quantum and Molecular Structures, Shanghai University, 99 Shangda Road, Shanghai, 200444 China
- 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
|
11
|
Bloino J. A VPT2 Route to Near-Infrared Spectroscopy: The Role of Mechanical and Electrical Anharmonicity. J Phys Chem A 2015; 119:5269-87. [DOI: 10.1021/jp509985u] [Citation(s) in RCA: 76] [Impact Index Per Article: 8.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Affiliation(s)
- 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, 56124 Pisa, Italy
- Scuola Normale Superiore, piazza dei Cavalieri 7, I-56126 Pisa, Italy
| |
Collapse
|
12
|
Fornaro T, Carnimeo I, Biczysko M. Toward Feasible and Comprehensive Computational Protocol for Simulation of the Spectroscopic Properties of Large Molecular Systems: The Anharmonic Infrared Spectrum of Uracil in the Solid State by the Reduced Dimensionality/Hybrid VPT2 Approach. J Phys Chem A 2014; 119:5313-26. [DOI: 10.1021/jp510101y] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Affiliation(s)
- Teresa Fornaro
- Scuola Normale Superiore, Piazza
dei Cavalieri 7, I-56126 Pisa, Italy
| | - Ivan Carnimeo
- Scuola Normale Superiore, Piazza
dei Cavalieri 7, I-56126 Pisa, Italy
- Compunet, Istituto Italiano di Tecnologia (IIT), Via Morego 30, I-16163 Genova, Italy
| | - Malgorzata Biczysko
- Scuola Normale Superiore, Piazza
dei Cavalieri 7, I-56126 Pisa, Italy
- Istituto di Chimica
dei Composti OrganoMetallici (ICCOM-CNR), UOS di Pisa, Area della
Ricerca CNR, Consiglio Nazionale delle Ricerche, Via G. Moruzzi 1, I-56124 Pisa, Italy
| |
Collapse
|
13
|
Maeda S, Taketsugu T, Morokuma K, Ohno K. Anharmonic Downward Distortion Following for Automated Exploration of Quantum Chemical Potential Energy Surfaces. BULLETIN OF THE CHEMICAL SOCIETY OF JAPAN 2014. [DOI: 10.1246/bcsj.20140189] [Citation(s) in RCA: 32] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Affiliation(s)
- Satoshi Maeda
- Department of Chemistry, Faculty of Science, Hokkaido University
| | | | - Keiji Morokuma
- Fukui Institute for Fundamental Chemistry, Kyoto University
- Cherry L. Emerson Center for Scientific Computation and Department of Chemistry, Emory University
| | - Koichi Ohno
- Graduate School of Science, Tohoku University
- Institute for Quantum Chemical Exploration
| |
Collapse
|
14
|
Barone V, Baiardi A, Bloino J. New developments of a multifrequency virtual spectrometer: stereo-electronic, dynamical, and environmental effects on chiroptical spectra. Chirality 2014; 26:588-600. [PMID: 24839096 PMCID: PMC4604663 DOI: 10.1002/chir.22325] [Citation(s) in RCA: 32] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/10/2013] [Revised: 01/29/2014] [Accepted: 03/03/2014] [Indexed: 11/09/2022]
Abstract
Computational spectroscopy has recently evolved from a field reserved for specialists toward a general tool allowing interpretations and analyses of experimental results. However, the current practice of providing tables of transitions for rigid geometries, possibly tuned by phenomenological broadening, is by far too naive. In order to improve this situation in the last few years we have been developing a general, robust, and user-friendly virtual spectrometer (VS) able to complement experimental studies for complex systems in condensed phases. The VS is based on flexible graphical pre- and postprocessing tools interfaced with general number-crunching software. This last tool is rooted in several electronic structure methodologies (DFT, TD-DFT, post-Hartree-Fock), powerful discrete/continuum models for describing environmental effects, and general vibrational and vibronic models. These last topics are the main focus of this work, which sketches our latest developments related to effective inclusion of anharmonic contributions, together with time-independent and/or time-dependent descriptions of vibronic transitions including Franck-Condon, Herzberg-Teller, and Duschinsky effects. Some test cases are described in some detail with the aim of showing the role of different effects in ruling vibrational (VCD) and electronic (ECD, CPL) chiral spectroscopies.
Collapse
Affiliation(s)
- Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Alberto Baiardi
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Julien Bloino
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
- Consiglio Nazionale delle Ricerche, Istituto di Chimica dei Composti Organometallici, UOS di Pisa, Via G. Moruzzi 1, 56124 Pisa, Italy
| |
Collapse
|
15
|
Fornaro T, Biczysko M, Monti S, Barone V. Dispersion corrected DFT approaches for anharmonic vibrational frequency calculations: nucleobases and their dimers. Phys Chem Chem Phys 2014; 16:10112-28. [PMID: 24531740 PMCID: PMC4612423 DOI: 10.1039/c3cp54724h] [Citation(s) in RCA: 68] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Abstract
Computational spectroscopy techniques have become in the last few years an effective means to analyze and assign infrared (IR) spectra of molecular systems of increasing dimensions and in different environments. However, transition from compilation of harmonic data to fully anharmonic simulations of spectra is still underway. The most promising results for large systems have been obtained, in our opinion, by perturbative vibrational approaches based on potential energy surfaces computed by hybrid (especially B3LYP) density functionals and medium size (e.g. SNSD) basis sets. In this framework, we are actively developing a comprehensive and robust computational protocol aimed at quantitative reproduction of the spectra of nucleic acid base complexes and their adsorption on solid supports (organic/inorganic). In this contribution we report the essential results of the first step devoted to isolated monomers and dimers. It is well known that in order to model the vibrational spectra of weakly bound molecular complexes dispersion interactions should be taken into proper account. In this work we have chosen two popular and inexpensive approaches to model dispersion interactions, namely the semi-empirical dispersion correction (D3) and pseudopotential based (DCP) methodologies both in conjunction with the B3LYP functional. These have been used for simulating fully anharmonic IR spectra of nucleobases and their dimers through generalized second order vibrational perturbation theory (GVPT2). We have studied, in particular, isolated adenine, hypoxanthine, uracil, thymine and cytosine, the hydrogen-bonded and stacked adenine and uracil dimers, and the stacked adenine-naphthalene heterodimer. Anharmonic frequencies are compared with standard B3LYP results and experimental findings, while the computed interaction energies and structures of complexes are compared to the best available theoretical estimates.
Collapse
Affiliation(s)
- Teresa Fornaro
- Scuola Normale Superiore, piazza dei Cavalieri 7, I-56126 Pisa, Italy.
| | | | | | | |
Collapse
|
16
|
Barone V, Biczysko M, Bloino J. Fully anharmonic IR and Raman spectra of medium-size molecular systems: accuracy and interpretation. Phys Chem Chem Phys 2014; 16:1759-87. [PMID: 24346191 PMCID: PMC4604664 DOI: 10.1039/c3cp53413h] [Citation(s) in RCA: 295] [Impact Index Per Article: 29.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
Computation of full infrared (IR) and Raman spectra (including absolute intensities and transition energies) for medium- and large-sized molecular systems beyond the harmonic approximation is one of the most interesting challenges of contemporary computational chemistry. Contrary to common beliefs, low-order perturbation theory is able to deliver results of high accuracy (actually often better than those issuing from current direct dynamics approaches) provided that anharmonic resonances are properly managed. This perspective sketches the recent developments in our research group toward the development of a robust and user-friendly virtual spectrometer rooted in second-order vibrational perturbation theory (VPT2) and usable also by non-specialists essentially as a black-box procedure. Several examples are explicitly worked out in order to illustrate the features of our computational tool together with the most important ongoing developments.
Collapse
Affiliation(s)
- Vincenzo Barone
- Scuola Normale Superiore, piazza dei Cavalieri 7, I-56126 Pisa, Italy.
| | | | | |
Collapse
|
17
|
Aidas K, Angeli C, Bak KL, Bakken V, Bast R, Boman L, Christiansen O, Cimiraglia R, Coriani S, Dahle P, Dalskov EK, Ekström U, Enevoldsen T, Eriksen JJ, Ettenhuber P, Fernández B, Ferrighi L, Fliegl H, Frediani L, Hald K, Halkier A, Hättig C, Heiberg H, Helgaker T, Hennum AC, Hettema H, Hjertenæs E, Høst S, Høyvik IM, Iozzi MF, Jansík B, Jensen HJA, Jonsson D, Jørgensen P, Kauczor J, Kirpekar S, Kjærgaard T, Klopper W, Knecht S, Kobayashi R, Koch H, Kongsted J, Krapp A, Kristensen K, Ligabue A, Lutnæs OB, Melo JI, Mikkelsen KV, Myhre RH, Neiss C, Nielsen CB, Norman P, Olsen J, Olsen JMH, Osted A, Packer MJ, Pawlowski F, Pedersen TB, Provasi PF, Reine S, Rinkevicius Z, Ruden TA, Ruud K, Rybkin VV, Sałek P, Samson CCM, de Merás AS, Saue T, Sauer SPA, Schimmelpfennig B, Sneskov K, Steindal AH, Sylvester-Hvid KO, Taylor PR, Teale AM, Tellgren EI, Tew DP, Thorvaldsen AJ, Thøgersen L, Vahtras O, Watson MA, Wilson DJD, Ziolkowski M, Agren H. The Dalton quantum chemistry program system. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2013; 4:269-284. [PMID: 25309629 PMCID: PMC4171759 DOI: 10.1002/wcms.1172] [Citation(s) in RCA: 870] [Impact Index Per Article: 79.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
Dalton is a powerful general-purpose program system for the study of molecular electronic structure at the Hartree-Fock, Kohn-Sham, multiconfigurational self-consistent-field, Møller-Plesset, configuration-interaction, and coupled-cluster levels of theory. Apart from the total energy, a wide variety of molecular properties may be calculated using these electronic-structure models. Molecular gradients and Hessians are available for geometry optimizations, molecular dynamics, and vibrational studies, whereas magnetic resonance and optical activity can be studied in a gauge-origin-invariant manner. Frequency-dependent molecular properties can be calculated using linear, quadratic, and cubic response theory. A large number of singlet and triplet perturbation operators are available for the study of one-, two-, and three-photon processes. Environmental effects may be included using various dielectric-medium and quantum-mechanics/molecular-mechanics models. Large molecules may be studied using linear-scaling and massively parallel algorithms. Dalton is distributed at no cost from http://www.daltonprogram.org for a number of UNIX platforms.
Collapse
Affiliation(s)
- Kestutis Aidas
- Department of General Physics and Spectroscopy, Faculty of Physics, Vilnius University Vilnius, Lithuania
| | | | - Keld L Bak
- Aarhus University School of Engineering Aarhus, Denmark
| | - Vebjørn Bakken
- Faculty of Mathematics and Natural Sciences, University of Oslo Oslo, Norway
| | - Radovan Bast
- Department of Theoretical Chemistry and Biology, School of Biotechnology, KTH Royal Institute of Technology Stockholm, Sweden
| | | | | | | | - Sonia Coriani
- Department of Chemical and Pharmaceutical Sciences, University of Trieste Trieste, Italy
| | - Pål Dahle
- Norwegian Computing Center Oslo, Norway
| | | | - Ulf Ekström
- CTCC, Department of Chemistry, University of Oslo Oslo, Norway
| | - Thomas Enevoldsen
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark Odense, Denmark
| | | | | | - Berta Fernández
- Department of Physical Chemistry and Center for Research in Biological Chemistry and Molecular Materials (CIQUS), University of Santiago de Compostela Santiago de Compostela, Spain
| | - Lara Ferrighi
- CTCC, Department of Chemistry, UiT The Arctic University of Norway, Tromsø Norway
| | - Heike Fliegl
- CTCC, Department of Chemistry, University of Oslo Oslo, Norway
| | - Luca Frediani
- CTCC, Department of Chemistry, UiT The Arctic University of Norway, Tromsø Norway
| | | | | | - Christof Hättig
- Department of Theoretical Chemistry, Ruhr-University Bochum Bochum, Germany
| | | | - Trygve Helgaker
- CTCC, Department of Chemistry, University of Oslo Oslo, Norway
| | | | - Hinne Hettema
- Department of Philosophy, The University of Auckland Auckland, New Zealand
| | - Eirik Hjertenæs
- Department of Chemistry, Norwegian University of Science and Technology Trondheim, Norway
| | - Stinne Høst
- Department of Geoscience, Aarhus University Aarhus, Denmark
| | | | | | | | - Hans Jørgen Aa Jensen
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark Odense, Denmark
| | - Dan Jonsson
- High-Performance Computing Group, UiT The Arctic University of Norway, Tromsø Norway
| | - Poul Jørgensen
- Department of Chemistry, Aarhus University Aarhus, Denmark
| | - Joanna Kauczor
- Department of Physics, Chemistry and Biology, Linköping University Linköping, Sweden
| | | | | | - Wim Klopper
- Institute of Physical Chemistry, Karlsruhe Institute of Technology Karlsruhe, Germany
| | - Stefan Knecht
- Laboratory of Physical Chemistry, ETH Zürich Zürich, Switzerland
| | - Rika Kobayashi
- Australian National University Supercomputer Facility Canberra, Australia
| | - Henrik Koch
- Department of Chemistry, Norwegian University of Science and Technology Trondheim, Norway
| | - Jacob Kongsted
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark Odense, Denmark
| | | | | | - Andrea Ligabue
- Computer Services: Networks and Systems, University of Modena and Reggio Emilia Modena, Italy
| | | | - Juan I Melo
- Physics Department, FCEyN-UBA and IFIBA-CONICET, Universidad de Buenos Aires Buenos Aires, Argentina
| | - Kurt V Mikkelsen
- Department of Chemistry, University of Copenhagen, Copenhagen Denmark
| | - Rolf H Myhre
- Department of Chemistry, Norwegian University of Science and Technology Trondheim, Norway
| | - Christian Neiss
- Department of Chemistry and Pharmacy, Friedrich-Alexander University Erlangen-Nürnberg Erlangen, Germany
| | | | - Patrick Norman
- Department of Physics, Chemistry and Biology, Linköping University Linköping, Sweden
| | - Jeppe Olsen
- Department of Chemistry, Aarhus University Aarhus, Denmark
| | - Jógvan Magnus H Olsen
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark Odense, Denmark
| | | | - Martin J Packer
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark Odense, Denmark
| | - Filip Pawlowski
- Institute of Physics, Kazimierz Wielki University Bydgoszcz, Poland
| | | | - Patricio F Provasi
- Department of Physics, University of Northeastern and IMIT-CONICET Corrientes, Argentina
| | - Simen Reine
- CTCC, Department of Chemistry, University of Oslo Oslo, Norway
| | - Zilvinas Rinkevicius
- Department of Theoretical Chemistry and Biology, School of Biotechnology and Swedish e-Science Research Center (SeRC), KTH Royal Institute of Technology Stockholm, Sweden
| | | | - Kenneth Ruud
- CTCC, Department of Chemistry, UiT The Arctic University of Norway, Tromsø Norway
| | - Vladimir V Rybkin
- Institute of Physical Chemistry, Karlsruhe Institute of Technology Karlsruhe, Germany
| | | | - Claire C M Samson
- Institute of Physical Chemistry, Karlsruhe Institute of Technology Karlsruhe, Germany
| | | | - Trond Saue
- Paul Sabatier University Toulouse, France
| | - Stephan P A Sauer
- Department of Chemistry, University of Copenhagen, Copenhagen Denmark
| | - Bernd Schimmelpfennig
- Institute for Nuclear Waste Disposal, Karlsruhe Institute of Technology Karlsruhe, Germany
| | | | - Arnfinn H Steindal
- CTCC, Department of Chemistry, UiT The Arctic University of Norway, Tromsø Norway
| | | | - Peter R Taylor
- VLSCI and School of Chemistry, University of Melbourne Parkville, Australia
| | - Andrew M Teale
- School of Chemistry, University of Nottingham Nottingham, UK
| | - Erik I Tellgren
- CTCC, Department of Chemistry, University of Oslo Oslo, Norway
| | - David P Tew
- School of Chemistry, University of Bristol Bristol, UK
| | | | | | - Olav Vahtras
- Department of Theoretical Chemistry and Biology, School of Biotechnology, KTH Royal Institute of Technology Stockholm, Sweden
| | - Mark A Watson
- Department of Chemistry, Princeton University Princeton, New Jersey
| | - David J D Wilson
- Department of Chemistry and La Trobe Institute for Molecular Sciences, La Trobe University Melbourne, Australia
| | - Marcin Ziolkowski
- CoE for Next Generation Computing, Clemson University Clemson, South Carolina
| | - Hans Agren
- Department of Theoretical Chemistry and Biology, School of Biotechnology, KTH Royal Institute of Technology Stockholm, Sweden
| |
Collapse
|
18
|
Charmet AP, Quartarone G, Ronchin L, Tortato C, Vavasori A. Quantum Chemical Investigation on Indole: Vibrational Force Field and Theoretical Determination of Its Aqueous pKa Value. J Phys Chem A 2013; 117:6846-58. [DOI: 10.1021/jp4049692] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Andrea Pietropolli Charmet
- Dipartimento
di Scienze Molecolari e Nanosistemi, Università Ca’ Foscari Venezia, Calle
Larga Santa Marta 2137, I-30123, Venezia, Italy
| | - Giuseppe Quartarone
- Dipartimento
di Scienze Molecolari e Nanosistemi, Università Ca’ Foscari Venezia, Calle
Larga Santa Marta 2137, I-30123, Venezia, Italy
| | - Lucio Ronchin
- Dipartimento
di Scienze Molecolari e Nanosistemi, Università Ca’ Foscari Venezia, Calle
Larga Santa Marta 2137, I-30123, Venezia, Italy
| | - Claudio Tortato
- Dipartimento
di Scienze Molecolari e Nanosistemi, Università Ca’ Foscari Venezia, Calle
Larga Santa Marta 2137, I-30123, Venezia, Italy
| | - Andrea Vavasori
- Dipartimento
di Scienze Molecolari e Nanosistemi, Università Ca’ Foscari Venezia, Calle
Larga Santa Marta 2137, I-30123, Venezia, Italy
| |
Collapse
|
19
|
Zen A, Zhelyazov D, Guidoni L. Optimized Structure and Vibrational Properties by Error Affected Potential Energy Surfaces. J Chem Theory Comput 2012; 8:4204-4215. [PMID: 24093004 DOI: 10.1021/ct300576n] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The precise theoretical determination of the geometrical parameters of molecules at the minima of their potential energy surface and of the corresponding vibrational properties are of fundamental importance for the interpretation of vibrational spectroscopy experiments. Quantum Monte Carlo techniques are correlated electronic structure methods promising for large molecules, which are intrinsically affected by stochastic errors on both energy and force calculations, making the mentioned calculations more challenging with respect to other more traditional quantum chemistry tools. To circumvent this drawback in the present work, we formulate the general problem of evaluating the molecular equilibrium structures, the harmonic frequencies, and the anharmonic coefficients of an error affected potential energy surface. The proposed approach, based on a multidimensional fitting procedure, is illustrated together with a critical evaluation of systematic and statistical errors. We observe that the use of forces instead of energies in the fitting procedure reduces the statistical uncertainty of the vibrational parameters by 1 order of magnitude. Preliminary results based on variational Monte Carlo calculations on the water molecule demonstrate the possibility to evaluate geometrical parameters and harmonic and anharmonic coefficients at this level of theory with an affordable computational cost and a small stochastic uncertainty (<0.07% for geometries and <0.7% for vibrational properties).
Collapse
Affiliation(s)
- Andrea Zen
- Dipartimento di Fisica, La Sapienza-Università di Roma, P.le Aldo Moro 2, 00185 Roma, Italy
| | | | | |
Collapse
|
20
|
Bloino J, Barone V. A second-order perturbation theory route to vibrational averages and transition properties of molecules: general formulation and application to infrared and vibrational circular dichroism spectroscopies. J Chem Phys 2012; 136:124108. [PMID: 22462836 DOI: 10.1063/1.3695210] [Citation(s) in RCA: 259] [Impact Index Per Article: 21.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023] Open
Abstract
A general formulation to compute anharmonic vibrational averages and transition properties at the second-order of perturbation theory is derived from the Rayleigh-Schrödinger development. This approach is intended to be applicable to any property expanded as a Taylor series up to the third order with respect to normal coordinates or their associated momenta. The equations are straightforward to implement and can be easily adapted to various properties, as illustrated for the case of electric and magnetic dipole moments. From those, infrared and vibrational circular dichroism spectra can be readily obtained. This fully automatic procedure has been applied to several chiral molecules of small-to-medium sizes and compared to the standard double harmonic approximation and to experimental data.
Collapse
Affiliation(s)
- 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, 56124 Pisa, Italy.
| | | |
Collapse
|
21
|
Bloino J, Biczysko M, Barone V. General Perturbative Approach for Spectroscopy, Thermodynamics, and Kinetics: Methodological Background and Benchmark Studies. J Chem Theory Comput 2012; 8:1015-36. [DOI: 10.1021/ct200814m] [Citation(s) in RCA: 213] [Impact Index Per Article: 17.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Affiliation(s)
- Julien Bloino
- Scuola Normale Superiore, piazza
dei Cavalieri 7, I-56126 Pisa, Italy
- Consiglio Nazionale delle Ricerche, Istituto di Chimica dei Composti Organometallici (CNR-ICCOM), UOS di Pisa, Via G. Moruzzi, 1 I-56124 Italy
| | - Malgorzata Biczysko
- Center for Nanotechnology Innovation
@NEST, Istituto Italiano di Tecnologia, Piazza San Silvestro 12, I-56127
Pisa, Italy
| | - Vincenzo Barone
- Scuola Normale Superiore, piazza
dei Cavalieri 7, I-56126 Pisa, Italy
| |
Collapse
|
22
|
Biczysko M, Bloino J, Carnimeo I, Panek P, Barone V. Fully ab initio IR spectra for complex molecular systems from perturbative vibrational approaches: Glycine as a test case. J Mol Struct 2012. [DOI: 10.1016/j.molstruc.2011.10.012] [Citation(s) in RCA: 43] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
23
|
Barone V, Baiardi A, Biczysko M, Bloino J, Cappelli C, Lipparini F. Implementation and validation of a multi-purpose virtual spectrometer for large systems in complex environments. Phys Chem Chem Phys 2012; 14:12404-22. [DOI: 10.1039/c2cp41006k] [Citation(s) in RCA: 122] [Impact Index Per Article: 10.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
|
24
|
Puzzarini C, Biczysko M, Barone V. Accurate Anharmonic Vibrational Frequencies for Uracil: The Performance of Composite Schemes and Hybrid CC/DFT Model. J Chem Theory Comput 2011; 7:3702-10. [DOI: 10.1021/ct200552m] [Citation(s) in RCA: 98] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Cristina Puzzarini
- Dipartimento di Chimica “G. Ciamician”, Università di Bologna, Via F. Selmi 2, 40126 Bologna, Italy
| | - Malgorzata Biczysko
- Center for Nanotechnology Innovation @NEST, Istituto Italiano di Tecnologia, Piazza San Silvestro, 12 - 56127 Pisa, Italy
| | - Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| |
Collapse
|
25
|
Vázquez J, Harding ME, Stanton JF, Gauss J. Vibrational Energy Levels via Finite-Basis Calculations Using a Quasi-Analytic Form of the Kinetic Energy. J Chem Theory Comput 2011; 7:1428-42. [DOI: 10.1021/ct100711u] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Juana Vázquez
- Center for Theoretical Chemistry, Department of Chemistry and Biochemistry, University of Texas at Austin, Austin, Texas 78712, United States
- Institut für Physikalische Chemie, Universität Mainz, D-55099 Mainz, Germany
| | - Michael E. Harding
- Center for Theoretical Chemistry, Department of Chemistry and Biochemistry, University of Texas at Austin, Austin, Texas 78712, United States
- Institut für Physikalische Chemie, Universität Mainz, D-55099 Mainz, Germany
| | - John F. Stanton
- Center for Theoretical Chemistry, Department of Chemistry and Biochemistry, University of Texas at Austin, Austin, Texas 78712, United States
| | - Jürgen Gauss
- Institut für Physikalische Chemie, Universität Mainz, D-55099 Mainz, Germany
| |
Collapse
|
26
|
Pedone A, Biczysko M, Barone V. Environmental effects in computational spectroscopy: accuracy and interpretation. Chemphyschem 2010; 11:1812-32. [PMID: 20358575 DOI: 10.1002/cphc.200900976] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Abstract
Spectroscopic techniques are valuable tools for understanding the structure and dynamics of complex systems, such as biomolecules or nanomaterials. Most of the current research is devoted to the development of new experimental techniques for improving the intrinsic resolution of different spectra. However, the subtle interplay of several different effects acting at different length and time scales still makes the interpretation and analysis of such spectra a very difficult task. In this respect, computational spectroscopy is becoming a needful and versatile tool for the assignment and interpretation of experimental spectra. It is in fact possible nowadays to model with relatively high accuracy the physical-chemical properties of complex molecules in different environments, and to link spectroscopic evidence directly to the structural and dynamical properties of optically or magnetically active solvated probes. In this Review, significant steps toward the simulation of entire spectra in condensed phases are presented together with some basic aspects of computational spectroscopy, which highlight how intramolecular and intermolecular degrees of freedom influence several spectroscopic parameters.
Collapse
Affiliation(s)
- Alfonso Pedone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, Pisa, Italy
| | | | | |
Collapse
|
27
|
Biczysko M, Panek P, Scalmani G, Bloino J, Barone V. Harmonic and Anharmonic Vibrational Frequency Calculations with the Double-Hybrid B2PLYP Method: Analytic Second Derivatives and Benchmark Studies. J Chem Theory Comput 2010; 6:2115-25. [PMID: 26615939 DOI: 10.1021/ct100212p] [Citation(s) in RCA: 221] [Impact Index Per Article: 15.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
This work aims to provide reliable benchmark data on the accuracy of harmonic and anharmonic vibrational frequencies computed with the B2PLYP double-hybrid density functional method. The exchange-correlation contributions required for the B2PLYP analytical second derivatives are presented here, which allow for the effective calculation of harmonic frequency as well as cubic and semidiagonal quartic force fields. The latter, in turn, are necessary to compute the anharmonic vibrational frequencies with the perturbative approach (VPT2). The quality of harmonic vibrational frequencies computed in conjunction with basis sets of double- to quadruple-ζ quality has been checked against reference data from the F38 benchmark set. Then, for an additional set of small closed- and open-shell systems, both harmonic frequencies and anharmonic contributions computed at the B2PLYP/N07D and the B2PLYP/aug-cc-pVTZ levels have been compared to their CCSD(T) counterparts. Moreover, for selected medium-size molecules (furan, pyrrole, thiophene, uracil, anisole, phenol, and pyridine), anharmonic frequencies have been compared to well established experimental results. Such benchmark studies have shown that the B2PLYP/N07D model provides good quality harmonic frequencies and describes correctly anharmonic contributions, the latter being of similar accuracy to their B3LYP/N07D counterparts, but obtained at significantly larger computational cost. Additionally, increased accuracy can be obtained by adopting hybrid models where the B2PLYP/N07D anharmonic contributions are combined with harmonic frequencies computed with more accurate quantum mechanical (QM) approaches or by B2PLYP with larger basis sets. This work confirmed also that most of the recently developed density functionals are significantly less suited for vibrational computations, while the B2PLYP method can be recommended for spectroscopic studies where a good accuracy of vibrational properties is required.
Collapse
Affiliation(s)
- Malgorzata Biczysko
- Scuola Normale Superiore, piazza dei Cavalieri 7, 56126 Pisa, Italy, Dipartimento di Chimica "Paolo Corradini" and CR-INSTM Village Università di Napoli Federico II, Complesso Univ. Monte S. Angelo, via Cintia, 80126 Napoli, Italy, Faculty of Chemistry, University of Wroclaw, ul. Joliot-Curie 14, 50-383 Wroclaw, Poland, Gaussian, Inc., 340 Quinnipiac St., Bldg 40, Wallingford, Connecticut 06492, and Istituto Nazionale di Fisica Nucleare (INFN), Sezione di Pisa, Pisa, Italy
| | - Pawel Panek
- Scuola Normale Superiore, piazza dei Cavalieri 7, 56126 Pisa, Italy, Dipartimento di Chimica "Paolo Corradini" and CR-INSTM Village Università di Napoli Federico II, Complesso Univ. Monte S. Angelo, via Cintia, 80126 Napoli, Italy, Faculty of Chemistry, University of Wroclaw, ul. Joliot-Curie 14, 50-383 Wroclaw, Poland, Gaussian, Inc., 340 Quinnipiac St., Bldg 40, Wallingford, Connecticut 06492, and Istituto Nazionale di Fisica Nucleare (INFN), Sezione di Pisa, Pisa, Italy
| | - Giovanni Scalmani
- Scuola Normale Superiore, piazza dei Cavalieri 7, 56126 Pisa, Italy, Dipartimento di Chimica "Paolo Corradini" and CR-INSTM Village Università di Napoli Federico II, Complesso Univ. Monte S. Angelo, via Cintia, 80126 Napoli, Italy, Faculty of Chemistry, University of Wroclaw, ul. Joliot-Curie 14, 50-383 Wroclaw, Poland, Gaussian, Inc., 340 Quinnipiac St., Bldg 40, Wallingford, Connecticut 06492, and Istituto Nazionale di Fisica Nucleare (INFN), Sezione di Pisa, Pisa, Italy
| | - Julien Bloino
- Scuola Normale Superiore, piazza dei Cavalieri 7, 56126 Pisa, Italy, Dipartimento di Chimica "Paolo Corradini" and CR-INSTM Village Università di Napoli Federico II, Complesso Univ. Monte S. Angelo, via Cintia, 80126 Napoli, Italy, Faculty of Chemistry, University of Wroclaw, ul. Joliot-Curie 14, 50-383 Wroclaw, Poland, Gaussian, Inc., 340 Quinnipiac St., Bldg 40, Wallingford, Connecticut 06492, and Istituto Nazionale di Fisica Nucleare (INFN), Sezione di Pisa, Pisa, Italy
| | - Vincenzo Barone
- Scuola Normale Superiore, piazza dei Cavalieri 7, 56126 Pisa, Italy, Dipartimento di Chimica "Paolo Corradini" and CR-INSTM Village Università di Napoli Federico II, Complesso Univ. Monte S. Angelo, via Cintia, 80126 Napoli, Italy, Faculty of Chemistry, University of Wroclaw, ul. Joliot-Curie 14, 50-383 Wroclaw, Poland, Gaussian, Inc., 340 Quinnipiac St., Bldg 40, Wallingford, Connecticut 06492, and Istituto Nazionale di Fisica Nucleare (INFN), Sezione di Pisa, Pisa, Italy
| |
Collapse
|
28
|
Esteban AL, Díez E, Galache MP, San Fabián J, Casanueva J, Contreras RH. Vibrational contributions to vicinal proton–proton coupling constants3JHH. Mol Phys 2010. [DOI: 10.1080/00268971003630687] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
29
|
Puzzarini C, Biczysko M, Barone V. Accurate Harmonic/Anharmonic Vibrational Frequencies for Open-Shell Systems: Performances of the B3LYP/N07D Model for Semirigid Free Radicals Benchmarked by CCSD(T) Computations. J Chem Theory Comput 2010; 6:828-38. [DOI: 10.1021/ct900594h] [Citation(s) in RCA: 109] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
Affiliation(s)
- Cristina Puzzarini
- Dipartimento di Chimica “G. Ciamician”, Università di Bologna, Via F. Selmi 2, 40126 Bologna, Italy, Dipartimento di Chimica “Paolo Corradini” and CR-INSTM Village, Università di Napoli Federico II, Complesso Univ. Monte S. Angelo, via Cintia, 80126 Napoli, Italy, and Scuola Normale Superiore, piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Malgorzata Biczysko
- Dipartimento di Chimica “G. Ciamician”, Università di Bologna, Via F. Selmi 2, 40126 Bologna, Italy, Dipartimento di Chimica “Paolo Corradini” and CR-INSTM Village, Università di Napoli Federico II, Complesso Univ. Monte S. Angelo, via Cintia, 80126 Napoli, Italy, and Scuola Normale Superiore, piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Vincenzo Barone
- Dipartimento di Chimica “G. Ciamician”, Università di Bologna, Via F. Selmi 2, 40126 Bologna, Italy, Dipartimento di Chimica “Paolo Corradini” and CR-INSTM Village, Università di Napoli Federico II, Complesso Univ. Monte S. Angelo, via Cintia, 80126 Napoli, Italy, and Scuola Normale Superiore, piazza dei Cavalieri 7, 56126 Pisa, Italy
| |
Collapse
|
30
|
Barone V, Biczysko M, Brancato G. Extending the Range of Computational Spectroscopy by QM/MM Approaches: Time-Dependent and Time-Independent Routes. ADVANCES IN QUANTUM CHEMISTRY 2010. [DOI: 10.1016/s0065-3276(10)59002-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/20/2023]
|
31
|
Matito E, Toffoli D, Christiansen O. A hierarchy of potential energy surfaces constructed from energies and energy derivatives calculated on grids. J Chem Phys 2009; 130:134104. [DOI: 10.1063/1.3092921] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
32
|
|
33
|
Maeda S, Watanabe Y, Ohno K. Finding important anharmonic terms in the sixth-order potential energy function by the scaled hypersphere search method: An application to vibrational analyses of molecules and clusters. J Chem Phys 2008; 128:144111. [DOI: 10.1063/1.2884348] [Citation(s) in RCA: 31] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
34
|
Beste A, Buchanan AC, Britt PF, Hathorn BC, Harrison RJ. Kinetic Analysis of the Pyrolysis of Phenethyl Phenyl Ether: Computational Prediction of α/β-Selectivities. J Phys Chem A 2007; 111:12118-26. [DOI: 10.1021/jp075861+] [Citation(s) in RCA: 56] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
35
|
Harding ME, Gauss J, Pflüger K, Werner HJ. High-Accuracy Extrapolated Ab Initio Thermochemistry of Vinyl Chloride. J Phys Chem A 2007; 111:13623-8. [DOI: 10.1021/jp0756348] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Affiliation(s)
- Michael E. Harding
- Institut für Physikalische Chemie, Universität Mainz, Jakob-Welder-Weg 11, D-55099 Mainz, Germany, and Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| | - Jürgen Gauss
- Institut für Physikalische Chemie, Universität Mainz, Jakob-Welder-Weg 11, D-55099 Mainz, Germany, and Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| | - Klaus Pflüger
- Institut für Physikalische Chemie, Universität Mainz, Jakob-Welder-Weg 11, D-55099 Mainz, Germany, and Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| | - Hans-Joachim Werner
- Institut für Physikalische Chemie, Universität Mainz, Jakob-Welder-Weg 11, D-55099 Mainz, Germany, and Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| |
Collapse
|
36
|
Yagi K, Hirata S, Hirao K. Efficient configuration selection scheme for vibrational second-order perturbation theory. J Chem Phys 2007; 127:034111. [PMID: 17655435 DOI: 10.1063/1.2748774] [Citation(s) in RCA: 53] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A fast algorithm of vibrational second-order Moller-Plesset perturbation theory is proposed, enabling a substantial reduction in the number of vibrational self-consistent-field (VSCF) configurations that need to be summed in the calculations. Important configurations are identified a priori by assuming that a reference VSCF wave function is approximated well by harmonic oscillator wave functions and that fifth- and higher-order anharmonicities are negligible. The proposed scheme has reduced the number of VSCF configurations by more than 100 times for formaldehyde, ethylene, and furazan with an error in computed frequencies being not more than a few cm(-1).
Collapse
Affiliation(s)
- Kiyoshi Yagi
- Department of Applied Chemistry, School of Engineering, The University of Tokyo, Tokyo 113-8656, Japan
| | | | | |
Collapse
|
37
|
Barone V. Anharmonic vibrational properties by a fully automated second-order perturbative approach. J Chem Phys 2007; 122:14108. [PMID: 15638643 DOI: 10.1063/1.1824881] [Citation(s) in RCA: 1142] [Impact Index Per Article: 67.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
This paper describes the implementation of a fully automated code for the building of anharmonic force constants and their use in a second-order perturbative evaluation of vibrorotational parameters. Next, a number of test applications are discussed, which show the strengths and limits of various computational levels.
Collapse
Affiliation(s)
- Vincenzo Barone
- Laboratorio di Struttura e Dinamica Molecolare, Dipartimento di Chimica, Università Federico II, Complesso Universitario Monte Sant'Angelo, Via Cintia, I-80126 Naples, Italy.
| |
Collapse
|
38
|
Kongsted J, Christiansen O. Automatic generation of force fields and property surfaces for use in variational vibrational calculations of anharmonic vibrational energies and zero-point vibrational averaged properties. J Chem Phys 2006; 125:124108. [PMID: 17014167 DOI: 10.1063/1.2352734] [Citation(s) in RCA: 65] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
An automatic and general procedure for the calculation of geometrical derivatives of the energy and general property surfaces for molecular systems is developed and implemented. General expressions for an n-mode representation are derived, where the n-mode representation includes only the couplings between n or less degrees of freedom. The general expressions are specialized to derivative force fields and property surfaces, and a scheme for calculation of the numerical derivatives is implemented. The implementation is interfaced to electronic structure programs and may be used for both ground and excited electronic states. The implementation is done in the context of a vibrational structure program and can be used in combination with vibrational self-consistent field (VSCF), vibrational configuration interaction (VCI), vibrational Moller-Plesset, and vibrational coupled cluster calculations of anharmonic wave functions and calculation of vibrational averaged properties at the VSCF and VCI levels. Sample calculations are presented for fundamental vibrational energies and vibrationally averaged dipole moments and frequency dependent polarizabilities and hyperpolarizabilities of water and formaldehyde.
Collapse
Affiliation(s)
- Jacob Kongsted
- Department of Chemistry, University of Aarhus, DK-8000 Aarhus C, Denmark.
| | | |
Collapse
|
39
|
|
40
|
Rodriguez-Garcia V, Yagi K, Hirao K, Iwata S, Hirata S. Franck-Condon factors based on anharmonic vibrational wave functions of polyatomic molecules. J Chem Phys 2006; 125:014109. [PMID: 16863289 DOI: 10.1063/1.2209676] [Citation(s) in RCA: 36] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
Abstract
Franck-Condon (FC) integrals of polyatomic molecules are computed on the basis of vibrational self-consistent-field (VSCF) or configuration-interaction (VCI) calculations capable of including vibrational anharmonicity to any desired extent (within certain molecular size limits). The anharmonic vibrational wave functions of the initial and final states are expanded unambiguously by harmonic oscillator basis functions of normal coordinates of the respective electronic states. The anharmonic FC integrals are then obtained as linear combinations of harmonic counterparts, which can, in turn, be evaluated by established techniques taking account of the Duschinsky rotations, geometry displacements, and frequency changes. Alternatively, anharmonic wave functions of both states are expanded by basis functions of just one electronic state, permitting the FC integral to be evaluated directly by the Gauss-Hermite quadrature used in the VSCF and VCI steps [Bowman et al., Mol. Phys. 104, 33 (2006)]. These methods in conjunction with the VCI and coupled-cluster with singles, doubles, and perturbative triples [CCSD(T)] method have predicted the peak positions and intensities of the vibrational manifold in the X 2B1 photoelectron band of H2O with quantitative accuracy. It has revealed that two weakly visible peaks are the result of intensity borrowing from nearby states through anharmonic couplings, an effect explained qualitatively by VSCF and quantitatively by VCI, but not by the harmonic approximation. The X 2B2 photoelectron band of H2CO is less accurately reproduced by this method, likely because of the inability of CCSD(T)/cc-pVTZ to describe the potential energy surface of open-shell H2CO+ with the same high accuracy as in H2O+.
Collapse
Affiliation(s)
- Valerie Rodriguez-Garcia
- Quantum Theory Project, Department of Chemistry, University of Florida, Gainesville, Florida 32611-8435, USA
| | | | | | | | | |
Collapse
|
41
|
|
42
|
Oyanagi C, Yagi K, Taketsugu T, Hirao K. Highly accurate potential-energy and dipole moment surfaces for vibrational state calculations of methane. J Chem Phys 2006; 124:64311. [PMID: 16483211 DOI: 10.1063/1.2162891] [Citation(s) in RCA: 37] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Abstract
Full-dimensional ab initio potential-energy surface (PES) and dipole moment surface are constructed for a methane molecule at the CCSD(T)/cc-pVTZ and MP2/cc-pVTZ levels of theory, respectively, by the modified Shepard interpolation method based on the fourth-order Taylor expansion [MSI(4th)]. The reference points for the interpolation have been set in the coupling region of CH symmetric and antisymmetric stretching modes so as to reproduce the vibrational energy levels related to CH stretching vibrations. The vibrational configuration-interaction calculations have been performed to obtain the energy levels and the absorption intensities up to 9000 cm(-1) with the use of MSI(4th)-PES. The calculated fundamental frequencies and low-lying vibrational energy levels show that MSI(4th) is superior to the widely employed quartic force field, giving a better agreement with the experimental values. The absorption bands of overtones as well as combination bands, which are caused by purely anharmonic effects, have been obtained up to 9000 cm(-1). Strongly coupled states with visible intensity have been found in the 6500-9000 cm(-1) region where the experimental data are still lacking.
Collapse
Affiliation(s)
- Chikako Oyanagi
- Department of Chemistry, Ochanomizu University, Tokyo, Japan
| | | | | | | |
Collapse
|
43
|
Dunn ME, Evans TM, Kirschner KN, Shields GC. Prediction of accurate anharmonic experimental vibrational frequencies for water clusters, (H2O)n, n=2-5. J Phys Chem A 2006; 110:303-9. [PMID: 16392869 PMCID: PMC2548414 DOI: 10.1021/jp054958y] [Citation(s) in RCA: 84] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Accurate anharmonic experimental vibrational frequencies for water clusters consisting of 2-5 water molecules have been predicted on the basis of comparing different methods with MP2/aug-cc-pVTZ calculated and experimental anharmonic frequencies. The combination of using HF/6-31G* scaled frequencies for intramolecular modes and anharmonic frequencies for intermolecular modes gives excellent agreement with experiment for the water dimer and trimer and are as good as the expensive anharmonic MP2 calculations. The water trimer, the cyclic Ci and S4 tetramers, and the cyclic pentamer all have unique peaks in the infrared spectrum between 500 and 800 cm-1 and between 3400 and 3700 cm-1. Under the right experimental conditions these different clusters can be uniquely identified using high-resolution IR spectroscopy.
Collapse
Affiliation(s)
- Meghan E Dunn
- Department of Chemistry, Hamilton College, 198 College Hill Road, Clinton, New York 13323, USA
| | | | | | | |
Collapse
|
44
|
A scaled hypersphere interpolation technique for efficient construction of multidimensional potential energy surfaces. Chem Phys Lett 2005. [DOI: 10.1016/j.cplett.2005.08.063] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
45
|
Barone V, Carbonniere P, Pouchan C. Accurate vibrational spectra and magnetic properties of organic free radicals: The case of H2CN. J Chem Phys 2005; 122:224308. [PMID: 15974669 DOI: 10.1063/1.1924592] [Citation(s) in RCA: 44] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present the structural, magnetic, and vibrational properties of H2CN computed using a second-order perturbative approach in which equilibrium values and harmonic frequencies evaluated at the coupled-cluster level are combined with anharmonic and vibrational averaging contributions obtained by hybrid Hartree-Fock/Kohn-Sham methods. Our computations lead to remarkably accurate results and suggest some revision of the experimental vibrational assignments.
Collapse
Affiliation(s)
- Vincenzo Barone
- Dipartimento di Chimica, Università Federico II, Complesso Monte S. Angelo, via Cintia, 1-80126 Napoli, Italy.
| | | | | |
Collapse
|
46
|
Abstract
Using higher levels of wave-function-based electronic structure theory than previously applied, as well as density functional theory (B-LYP and B3-LYP functionals), all theoretical models conclude that three ONOOH conformers are stationary point minima, in disagreement with some of the previous studies that we survey. In order of increasing energy, these are the cis-cis, cis-perp, and trans-perp conformers. Basis sets including diffuse functions seem to be needed to obtain a qualitatively correct representation of the internal rotation potential energy surface at higher levels of theory. Internal rotation about the peroxide bond involving the cis-cis, cis-gauche transition structure (TS), cis-perp, and cis-trans TS conformers is studied in detail. To help ascertain the relative stability of the cis-perp conformer, multireference configuration interaction energy calculations are carried out, and rule of thumb estimates of multireference character in the ground-state wave functions of the ONOOH conformers are considered. CCSD(T)/aug-cc-pVTZ physical properties (geometries, rotational constants, electric dipole moments, harmonic vibrational frequencies, and infrared intensities) are compared with the analogous experimental data wherever possible, and also with density functional theory. Where such experimental data are nonexistent, the CCSD(T) and B3-LYP results are useful representations. For example, the electric dipole moment |mu(e)| of the cis-cis conformer is predicted to be 0.97+/-0.03 D. CCSD(T) energies, extrapolated to the aug-cc-pVNZ limit, are employed in isodesmic reaction schemes to derive zero Kelvin heats of formation and bond dissociation energies of the ONOOH stationary point minima. In agreement with recent gas-phase experiments, the peroxide bond dissociation energies of the cis-cis and trans-perp conformers are calculated as 19.3+/-0.4 and 16.0+/-0.4 kcalmol, respectively. The lowest energy cis-cis conformer is less stable than nitric acid by 28.1+/-0.4 kcalmol at 0 K.
Collapse
Affiliation(s)
- Mark P McGrath
- Department of Chemistry, University of California, Irvine, CA 92697, USA.
| | | |
Collapse
|
47
|
Ruden TA, Helgaker T, Jørgensen P, Olsen J. Coupled-cluster connected quadruples and quintuples corrections to the harmonic vibrational frequencies and equilibrium bond distances of HF, N2, F2, and CO. J Chem Phys 2004; 121:5874-84. [PMID: 15367015 DOI: 10.1063/1.1780155] [Citation(s) in RCA: 112] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Using the coupled-cluster method, we have examined the contributions from the connected quadruple and quintuple virtual excitations to the harmonic vibrational frequencies and equilibrium bond distances of HF, N(2), F(2), and CO. Whereas the largest quadruples contributions are -18.8 cm(-1) to the harmonic frequency of N(2) and 0.43 pm to the bond distance of F(2), the largest quintuples contributions are -3.9 cm(-1) to the harmonic frequency of N(2) and 0.03 pm to the bond distances of N(2) and F(2). As we improve the description by going from the coupled-cluster singles-and-doubles model with a perturbative triples correction [CCSD(T)] to the coupled-cluster singles-doubles-triples-and-quadruples model, the mean and maximum absolute errors in the calculated frequencies relative to experiment are reduced from 11.3 and 15.8 cm(-1), respectively, to 3.2 and 4.7 cm(-1); for the bond distances, the mean and maximum absolute errors are reduced from 0.16 and 0.47 pm, respectively, to 0.04 and 0.13 pm. The calculations presented here confirm previous observations that, to some extent, the relatively small errors in the CCSD(T) equilibrium bond distances and harmonic frequencies arise from a cancellation of errors in the approximate (perturbative) treatment of the connected triples and the neglect of higher-order connected excitations. Further inclusion of quintuples contributions, relativistic corrections, and adiabatic corrections reduces the mean and maximum absolute errors to 1.1 and 2.3 cm(-1), respectively, for the harmonic frequencies and to 0.02 and 0.05 pm, respectively, for the bond distances.
Collapse
Affiliation(s)
- Torgeir A Ruden
- Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, N-0315 Oslo, Norway
| | | | | | | |
Collapse
|
48
|
Yagi K, Hirao K, Taketsugu T, Schmidt MW, Gordon MS. Ab initiovibrational state calculations with a quartic force field: Applications to H2CO, C2H4, CH3OH, CH3CCH, and C6H6. J Chem Phys 2004; 121:1383-9. [PMID: 15260682 DOI: 10.1063/1.1764501] [Citation(s) in RCA: 158] [Impact Index Per Article: 7.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Abstract
For polyatomic molecules, n-mode coupling representations of the quartic force field (nMR-QFF) are presented, which include terms up to n normal coordinate couplings in a fourth-order polynomial potential energy function. The computational scheme to evaluate third-and fourth-order derivatives by finite differentiations of the energy is fully described. The code to generate the nMR-QFF has been implemented into GAMESS program package and interfaced with the vibrational self-consistent field (VSCF) and correlation corrected VSCF (cc-VSCF) methods. As a demonstration, fundamental frequencies have been calculated by the cc-VSCF method based on 2MR-QFF for formaldehyde, ethylene, methanol, propyne, and benzene. The applications show that 2MR-QFF is a highly accurate potential energy function, with errors of 1.0-1.9% relative to the experimental value in fundamental frequencies. This approach will help quantitative evaluations of vibrational energies of a general molecule with a reasonable computational cost.
Collapse
Affiliation(s)
- Kiyoshi Yagi
- Department of Applied Chemistry, Graduate School of Engineering, The University of Tokyo, Tokyo 113-8656, Japan
| | | | | | | | | |
Collapse
|
49
|
Barone V. Accurate Vibrational Spectra of Large Molecules by Density Functional Computations beyond the Harmonic Approximation: The Case of Azabenzenes. J Phys Chem A 2004. [DOI: 10.1021/jp0497662] [Citation(s) in RCA: 130] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Vincenzo Barone
- Laboratorio di Struttura e Dinamica Molecolare, Dipartimento di Chimica, Università “Federico II”, Complesso Universitario Monte Sant'Angelo, Via Cintia, I-80126 Napoli, Italy
| |
Collapse
|
50
|
Barone V, Festa G, Grandi A, Rega N, Sanna N. Accurate vibrational spectra of large molecules by density functional computations beyond the harmonic approximation: the case of uracil and 2-thiouracil. Chem Phys Lett 2004. [DOI: 10.1016/j.cplett.2004.03.024] [Citation(s) in RCA: 96] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|