1
|
Panek PT, Hoeske AA, Jacob CR. On the choice of coordinates in anharmonic theoretical vibrational spectroscopy: Harmonic vs. anharmonic coupling in vibrational configuration interaction. J Chem Phys 2019; 150:054107. [PMID: 30736699 DOI: 10.1063/1.5083186] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023] Open
Abstract
By a suitable choice of coordinates, the computational effort required for calculations of anharmonic vibrational spectra can be reduced significantly. By using suitable localized-mode coordinates obtained from an orthogonal transformation of the conventionally used normal-mode coordinates, anharmonic couplings between modes can be significantly reduced. However, such a transformation introduces harmonic couplings between the localized modes. To elucidate the role of these harmonic couplings, we consider the vibrational self-consistent field (VSCF)/vibrational configuration interaction (VCI) calculations for both few-mode model systems and for ethene as a molecular test case. We show that large harmonic couplings can result in significant errors in localized-mode L-VSCF/L-VCI calculations and study the convergence with respect to the size of the VCI excitation space. To further elucidate the errors introduced by harmonic couplings, we discuss the connection between L-VSCF/L-VCI and vibrational exciton models. With the help of our results, we propose an algorithm for the localization of normal modes in suitable subsets that are chosen to strictly limit the errors introduced by the harmonic couplings while still leading to maximally localized modes.
Collapse
Affiliation(s)
- Paweł T Panek
- Technische Universität Braunschweig, Institute of Physical and Theoretical Chemistry, Gaußstraße 17, 38106 Braunschweig, Germany
| | - Adrian A Hoeske
- Technische Universität Braunschweig, Institute of Physical and Theoretical Chemistry, Gaußstraße 17, 38106 Braunschweig, Germany
| | - Christoph R Jacob
- Technische Universität Braunschweig, Institute of Physical and Theoretical Chemistry, Gaußstraße 17, 38106 Braunschweig, Germany
| |
Collapse
|
2
|
Bulik IW, Frisch MJ, Vaccaro PH. Fixed-Node, Importance-Sampling Diffusion Monte Carlo for Vibrational Structure with Accurate and Compact Trial States. J Chem Theory Comput 2018; 14:1554-1563. [DOI: 10.1021/acs.jctc.8b00016] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Ireneusz W. Bulik
- Department of Chemistry, Yale University, New Haven, Connecticut 06511, United States
| | - Michael J. Frisch
- Gaussian, Inc., 340 Quinnipiac St. Bldg. 40, Wallingford, Connecticut 06492, United States
| | - Patrick H. Vaccaro
- Department of Chemistry, Yale University, New Haven, Connecticut 06511, United States
| |
Collapse
|
3
|
Bulik IW, Frisch MJ, Vaccaro PH. Vibrational self-consistent field theory using optimized curvilinear coordinates. J Chem Phys 2017; 147:044110. [DOI: 10.1063/1.4995440] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Ireneusz W. Bulik
- Department of Chemistry, Yale University, New Haven, Connecticut 06511, USA
| | - Michael J. Frisch
- Gaussian, Inc., 340 Quinnipiac St. Bldg. 40, Wallingford, Connecticut 06492, USA
| | - Patrick H. Vaccaro
- Department of Chemistry, Yale University, New Haven, Connecticut 06511, USA
| |
Collapse
|
4
|
Sibaev M, Crittenden DL. Balancing accuracy and efficiency in selecting vibrational configuration interaction basis states using vibrational perturbation theory. J Chem Phys 2016. [DOI: 10.1063/1.4960600] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
5
|
Sibaev M, Crittenden DL. An efficient and numerically stable procedure for generating sextic force fields in normal mode coordinates. J Chem Phys 2016; 144:214107. [DOI: 10.1063/1.4953080] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- M. Sibaev
- Department of Chemistry, University of Canterbury, Christchurch, New Zealand
| | - D. L. Crittenden
- Department of Chemistry, University of Canterbury, Christchurch, New Zealand
| |
Collapse
|
6
|
Mancera LA, Benoit DM. Vibrational anharmonicity of small gold and silver clusters using the VSCF method. Phys Chem Chem Phys 2016; 18:529-49. [DOI: 10.1039/c5cp05283a] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We study the vibrational spectra of small neutral gold (Au2–Au10) and silver (Ag2–Au5) clusters using the vibrational self-consistent field method (VSCF) in order to account for anharmonicity.
Collapse
Affiliation(s)
- Luis A. Mancera
- Institute of Theoretical Chemistry
- University of Ulm
- D-89069 Ulm
- Germany
| | | |
Collapse
|
7
|
An alternative methodology to assess the quality of empirical potentials for small gold clusters. COMPUT THEOR CHEM 2015. [DOI: 10.1016/j.comptc.2015.05.026] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
8
|
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
|
9
|
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
|
10
|
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
|
11
|
Panek PT, Jacob CR. Efficient calculation of anharmonic vibrational spectra of large molecules with localized modes. Chemphyschem 2014; 15:3365-77. [PMID: 25080397 DOI: 10.1002/cphc.201402251] [Citation(s) in RCA: 54] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/17/2014] [Indexed: 11/11/2022]
Abstract
The analysis and interpretation of the vibrational spectra of complex (bio)molecular systems, such as polypeptides and proteins, requires support from quantum-chemical calculations. Such calculations are currently restricted to the harmonic approximation. Here, we show how one of the main bottlenecks in such calculations, the evaluation of the potential energy surface, can be overcome by using localized modes instead of the commonly employed normal modes. We apply such local vibrational self-consistent field (L-VSCF) and vibrational configuration interaction (L-VCI) calculations to a cyclic water tetramer and a helical hexa-alanine peptide. The results show that the use of localized modes is equivalent to the commonly used normal modes, but offers several advantages. First, a faster convergence with respect to the excitation level is observed in L-VCI calculations. Second, the localized modes provide a reduced representation of the couplings between modes that show a regular coupling pattern. This can be used to disregard a significant number of small two-mode potentials a priori. Several such reduced coupling approximations are explored, and we show that the number of single-point calculations required to evaluate the potential energy surface can be significantly reduced without introducing noticeable errors in the resulting vibrational spectra.
Collapse
Affiliation(s)
- Paweł T Panek
- Karlsruhe Institute of Technology (KIT), Center for Functional Nanostructures and Institute of Physical Chemistry, Wolfgang-Gaede-Str. 1a, 76131 Karlsruhe (Germany) http://www.christophjacob.eu
| | | |
Collapse
|
12
|
Wang XG, Carrington T. Rovibrational levels and wavefunctions of Cl−H2O. J Chem Phys 2014; 140:204306. [DOI: 10.1063/1.4875798] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
13
|
Mancera LA, Benoit DM. Towards an understanding of the vibrational spectrum of the neutral Au7cluster. Phys Chem Chem Phys 2013; 15:1929-43. [DOI: 10.1039/c2cp43211k] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
|
14
|
Roy TK, Gerber RB. Vibrational self-consistent field calculations for spectroscopy of biological molecules: new algorithmic developments and applications. Phys Chem Chem Phys 2013; 15:9468-92. [DOI: 10.1039/c3cp50739d] [Citation(s) in RCA: 141] [Impact Index Per Article: 12.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
15
|
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
|
16
|
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
|
17
|
Petković M. O–H Stretch in Phenol and Its Hydrogen-Bonded Complexes: Band Position and Relaxation Pathways. J Phys Chem A 2011; 116:364-71. [DOI: 10.1021/jp209897y] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
Affiliation(s)
- Milena Petković
- Faculty of Physical Chemistry, University of Belgrade, Studentski trg 12, 11 158 Belgrade, Serbia
| |
Collapse
|
18
|
Jin L, Simons JP, Gerber RB. Structures of the xylose–water complex: Energetics, transitions between conformers and spectroscopy. Chem Phys Lett 2011. [DOI: 10.1016/j.cplett.2011.11.008] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/15/2022]
|
19
|
Meng K, Wang J. Anharmonic overtone and combination states of glycine and two model peptides examined by vibrational self-consistent field theory. Phys Chem Chem Phys 2011; 13:2001-13. [DOI: 10.1039/c0cp01177k] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
20
|
Suwan I, Gerber R. VSCF in internal coordinates and the calculation of anharmonic torsional mode transitions. Chem Phys 2010. [DOI: 10.1016/j.chemphys.2010.05.025] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
21
|
Gerber RB, Sebek J. Dynamics simulations of atmospherically relevant molecular reactions. INT REV PHYS CHEM 2009. [DOI: 10.1080/01442350903016684] [Citation(s) in RCA: 29] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
22
|
Mielke SL, Truhlar DG. Improved Methods for Feynman Path Integral Calculations of Vibrational−Rotational Free Energies and Application to Isotopic Fractionation of Hydrated Chloride Ions. J Phys Chem A 2009; 113:4817-27. [DOI: 10.1021/jp900834u] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Steven L. Mielke
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431
| | - Donald G. Truhlar
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431
| |
Collapse
|
23
|
Manzhos S, Carrington T. Using neural networks, optimized coordinates, and high-dimensional model representations to obtain a vinyl bromide potential surface. J Chem Phys 2009; 129:224104. [PMID: 19071904 DOI: 10.1063/1.3021471] [Citation(s) in RCA: 81] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We demonstrate that it is possible to obtain good potentials using high-dimensional model representations (HDMRs) fitted with neural networks (NNs) from data in 12 dimensions and 15 dimensions. The HDMR represents the potential as a sum of lower-dimensional functions and our NN-based approach makes it possible to obtain all of these functions from one set of fitting points. To reduce the number of terms in the HDMR, we use optimized redundant coordinates. By using exponential neurons, one obtains a potential in sum-of-products form, which greatly facilitates quantum dynamics calculations. A 12-dimensional (reference) potential surface for vinyl bromide is first refitted to show that it can be represented as a sum of two-dimensional functions. To fit 3d functions of the original coordinates, to improve the potential, a huge amount of data would be required. Redundant coordinates avoid this problem. They enable us to bypass the combinatorial explosion of the number of terms which plagues all HDMR and multimode-type methods. We also fit to a set of approximately 70,000 ab initio points for vinyl bromide in 15 dimensions [M. Malshe et al., J. Chem. Phys. 127, 134105 (2007)] and show that it is possible to obtain a surface in sum-of-products form of quality similar to the quality of the full-dimensional fit. Although we obtain a full-dimensional surface, we limit the cost of the fitting by building it from fits of six-dimensional functions, each of which requires only a small NN.
Collapse
Affiliation(s)
- Sergei Manzhos
- Département de chimie, Université de Montréal, Case postale 6128, succursale Centre-ville Montréal, (Québec) H3C 3J7 Canada.
| | | |
Collapse
|
24
|
Benoit DM. Fast vibrational calculation of anharmonic OH-stretch frequencies for two low-energy noradrenaline conformers. J Chem Phys 2008; 129:234304. [DOI: 10.1063/1.3040427] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
25
|
|
26
|
Avila G, Kroes GJ, van Hemert MC. The photodissociation of the water dimer in the à band: A twelve-dimensional quasiclassical study. J Chem Phys 2008; 128:144313. [DOI: 10.1063/1.2868775] [Citation(s) in RCA: 15] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
27
|
Scribano Y, Benoit DM. Calculation of vibrational frequencies through a variational reduced-coupling approach. J Chem Phys 2007; 127:164118. [DOI: 10.1063/1.2798104] [Citation(s) in RCA: 36] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
28
|
Miller Y, Chaban GM, Zhou J, Asmis KR, Neumark DM, Gerber RB. Vibrational spectroscopy of (SO42−)∙(H2O)n clusters, n=1–5: Harmonic and anharmonic calculations and experiment. J Chem Phys 2007; 127:094305. [PMID: 17824737 DOI: 10.1063/1.2764074] [Citation(s) in RCA: 50] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The vibrational spectroscopy of (SO4(2-)).(H2O)n is studied by theoretical calculations for n=1-5, and the results are compared with experiments for n=3-5. The calculations use both ab initio MP2 and DFT/B3LYP potential energy surfaces. Both harmonic and anharmonic calculations are reported, the latter with the CC-VSCF method. The main findings are the following: (1) With one exception (H2O bending mode), the anharmonicity of the observed transitions, all in the experimental window of 540-1850 cm(-1), is negligible. The computed anharmonic coupling suggests that intramolecular vibrational redistribution does not play any role for the observed linewidths. (2) Comparison with experiment at the harmonic level of computed fundamental frequencies indicates that MP2 is significantly more accurate than DFT/B3LYP for these systems. (3) Strong anharmonic effects are, however, calculated for numerous transitions of these systems, which are outside the present observation window. These include fundamentals as well as combination modes. (4) Combination modes for the n=1 and n=2 clusters are computed. Several relatively strong combination transitions are predicted. These show strong anharmonic effects. (5) An interesting effect of the zero point energy (ZPE) on structure is found for (SO4(2-)).(H2O)(5): The global minimum of the potential energy corresponds to a C(s) structure, but with incorporation of ZPE the lowest energy structure is C2v, in accordance with experiment. (6) No stable structures were found for (OH-).(HSO4-).(H2O)n, for n<or=5.
Collapse
Affiliation(s)
- Yifat Miller
- Department of Physical Chemistry, The Hebrew University, Jerusalem 91904, Israel
| | | | | | | | | | | |
Collapse
|
29
|
Rheinecker J, Bowman JM. The calculated infrared spectrum of Cl-H2O using a new full dimensional ab initio potential surface and dipole moment surface. J Chem Phys 2007; 125:133206. [PMID: 17029453 DOI: 10.1063/1.2209675] [Citation(s) in RCA: 43] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We report a full dimensional, ab initio-based global potential energy surface (PES) and dipole moment surface for Cl-H2O. Both surfaces are symmetric with respect to interchange of the H atoms. The PES is a fit to thousands of electronic energies calculated using the coupled-cluster method [CCSD(T)] with a moderately large basis (aug-cc-pVTZ). Vibrational energies and wave functions are accurately obtained using MULTIMODE. The wave function and dipole moment surface are used to calculate and analyze the pure infrared spectrum at 0 K which is compared with experiment. Vibrational energies and the infrared spectra for DOD and HOD/DOH are also presented.
Collapse
Affiliation(s)
- Jaime Rheinecker
- Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| | | |
Collapse
|
30
|
Manzhos S, Carrington T. Using redundant coordinates to represent potential energy surfaces with lower-dimensional functions. J Chem Phys 2007; 127:014103. [PMID: 17627333 DOI: 10.1063/1.2746846] [Citation(s) in RCA: 74] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We propose a method for fitting potential energy surfaces with a sum of component functions of lower dimensionality. This form facilitates quantum dynamics calculations. We show that it is possible to reduce the dimensionality of the component functions by introducing new and redundant coordinates obtained with linear transformations. The transformations are obtained from a neural network. Different coordinates are used for different component functions and the new coordinates are determined as the potential is fitted. The quality of the fits and the generality of the method are illustrated by fitting reference potential surfaces of hydrogen peroxide and of the reaction OH+H(2)-->H(2)O+H.
Collapse
Affiliation(s)
- Sergei Manzhos
- Département de Chimie, Université de Montréal, CP 6128, succursale Centre-ville, Montréal (Québec) H3C 3J7, Canada.
| | | |
Collapse
|
31
|
Anharmonic vibrational spectroscopy calculations with electronic structure potentials: comparison of MP2 and DFT for organic molecules. Theor Chem Acc 2007. [DOI: 10.1007/s00214-007-0299-1] [Citation(s) in RCA: 44] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
32
|
|
33
|
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
|
34
|
Benoit DM. Efficient correlation-corrected vibrational self-consistent field computation of OH-stretch frequencies using a low-scaling algorithm. J Chem Phys 2006; 125:244110. [PMID: 17199343 DOI: 10.1063/1.2423006] [Citation(s) in RCA: 49] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
The authors present a new computational scheme to perform accurate and fast direct correlation-corrected vibrational self-consistent field (CC-VSCF) computations for a selected number of vibrational modes, which is aimed at predicting a few vibrations in large molecular systems. The method is based on a systematic selection of vibrational mode-mode coupling terms, leading to the direct ab initio construction of a sparse potential energy surface. The computational scaling of the CC-VSCF computation on the generated surface is then further reduced by using a screening procedure for the correlation-correction contributions. The proposed method is applied to the computation of the OH-stretch frequency of five aliphatic alcohols. The authors investigate the influence of different pseudopotential and all-electron basis sets on the quality of the correlated potential energy surfaces computed and on the OH-stretch frequencies calculated for each surface. With the help of these test systems, the authors show that their method offers a computational scaling that is two orders of magnitude lower than a standard CC-VSCF method and that it is of equal accuracy.
Collapse
Affiliation(s)
- David M Benoit
- Nachwuchsgruppe Theorie-SFB 569, Albert-Einstein-Allee 11, University of Ulm, D-89081 Ulm, Germany.
| |
Collapse
|
35
|
Malaspina T, Fileti EE, Riveros JM, Canuto S. Ab Initio Study of the Isomeric Equilibrium of the HCN···H2O and H2O···HCN Hydrogen-Bonded Clusters. J Phys Chem A 2006; 110:10303-8. [PMID: 16928122 DOI: 10.1021/jp062780p] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
An ab initio study of the stability, spectroscopic properties, and isomeric equilibrium of the hydrogen-bonded HCN...H2O and H2O...HCN isomers is presented. Density functional theory and perturbative second-order MP2 and coupled-cluster CCSD(T) calculations were carried out and binding energies obtained with correlation-consistent basis sets including extrapolation to the infinity basis set level. At the best theoretical level, CCSD(T), the H2O...HCN complex is more stable than the HCN...H2O complex by ca. 6.3 kJ mol(-1). Rotational and vibrational spectra, including anharmonic corrections, are calculated. These calculated spectroscopic data are used to obtain thermochemical contributions to the thermodynamic functions and hence the Gibbs free energy. The relative free energies are used to estimate the equilibrium constant for isomerism. We find that under typical conditions of supersonic expansion experiments (T < 150 K) H2O...HCN is essentially the only isomer present. Furthermore, our calculations indicate that the hydrogen-bonded cluster becomes favorable over the separated moieties at temperatures below 200 K.
Collapse
Affiliation(s)
- T Malaspina
- Instituto de Física, Universidade de São Paulo, CP 66318, 05315-970, São Paulo, SP, Brazil
| | | | | | | |
Collapse
|
36
|
Chakraborty A, Truhlar DG. Converged vibrational energy levels and quantum mechanical vibrational partition function of ethane. J Chem Phys 2006; 124:184310. [PMID: 16709109 DOI: 10.1063/1.2193521] [Citation(s) in RCA: 15] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
Abstract
The vibrational partition function of ethane is calculated in the temperature range of 200-600 K using well-converged energy levels that were calculated by vibrational configuration interaction, and the results are compared to the harmonic oscillator partition function. This provides the first test of the harmonic oscillator approximation for a molecule with more than five atoms. The absolute free energies computed by the harmonic oscillator approximation are in error by 0.59-0.62 kcal/mol over the 200-600 K temperature range.
Collapse
Affiliation(s)
- Arindam Chakraborty
- Department of Chemistry, University of Minnesota, Minneapolis, Minnesota 55455-0431, USA
| | | |
Collapse
|
37
|
|
38
|
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
|
39
|
Fridgen TD, McMahon TB, Maître P, Lemaire J. Experimental infrared spectra of Cl−(ROH) (R = H, CH3, CH3CH2) complexes in the gas-phase. Phys Chem Chem Phys 2006; 8:2483-90. [PMID: 16721432 DOI: 10.1039/b603102a] [Citation(s) in RCA: 15] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Infrared multiple photon dissociation spectra for the chloride ion solvated by either water, methanol or ethanol have been recorded using an FTICR spectrometer coupled to a free-electron laser, and are presented here along with assignments to the observed bands. The assignments made to the Cl(-)/H(2)O, Cl(-)(CH(3)OH), and Cl(-)(CH(3)CH(2)OH) spectra are based on comparison with the neutral H(2)O, CH(3)OH, and CH(3)CH(2)OH spectra, respectively. This work confirms that a band observed around 1400 cm(-1) in the Cl(-)(H(2)O) spectrum is not due to the Ar tag in Ar predissociation spectra. The carrier of this band is, most likely, the first overtone of the OHCl bend. Based on the position of the overtone in the IRMPD spectrum, 1375 cm(-1), the fundamental must occur very close to 700 cm(-1) and observation of this band should aid theoretical treatments of the spectrum of this complex. B3LYP/6-311++G(2df,2pd) calculations are shown to reproduce the IRMPD spectra of all three solvated chloride species. They also predict that attaching one or two Ar atoms to the Cl(-)(H(2)O) complex results in a shift of no more than a few wavenumbers in the fundamental bands for the bare complex, in agreement with previous experiment. For both alcohol-Cl(-) complexes, the S(N)2 "backside attack" isomers are not observed and Cl(-) is predicted theoretically, and confirmed experimentally, to be bound to the hydroxyl hydrogen. For Cl(-)(CH(3)CH(2)OH), the trans and gauche conformers are similar in energy, with the gauche conformer predicted to be thermodynamically favoured. The experimental infrared spectrum agrees well with that predicted for the gauche conformer but a mixture of gauche and anti conformers cannot be ruled out based on the experimental spectra nor on the computed thermochemistry.
Collapse
Affiliation(s)
- Travis D Fridgen
- Department of Chemistry, Memorial University of Newfoundland, St. John's, NL, Canada.
| | | | | | | |
Collapse
|
40
|
Day MB, Kirschner KN, Shields GC. Global Search for Minimum Energy (H2O)n Clusters, n = 3−5. J Phys Chem A 2005; 109:6773-8. [PMID: 16834031 DOI: 10.1021/jp0513317] [Citation(s) in RCA: 73] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
The Gaussian-3 (G3) model chemistry method has been used to calculate the relative deltaG(o) values for all possible conformers of neutral clusters of water, (H2O)n, where n = 3-5. A complete 12-fold conformational search around each hydrogen bond produced 144, 1728, and 20,736 initial starting structures of the water trimer, tetramer, and pentamer. These structures were optimized with PM3, followed by HF/6-31G* optimization, and then with the G3 model chemistry. Only two trimers are present on the G3 potential energy hypersurface. We identified 5 tetramers and 10 pentamers on the potential energy and free-energy hypersurfaces at 298 K. None of these 17 structures were linear; all linear starting models folded into cyclic or three-dimensional structures. The cyclic pentamer is the most stable isomer at 298 K. On the basis of this and previous studies, we expect the cyclic tetramers and pentamers to be the most significant cyclic water clusters in the atmosphere.
Collapse
Affiliation(s)
- Mary Beth Day
- Department of Chemistry, Hamilton College, 198 College Hill Road, Clinton, New York 13323, USA
| | | | | |
Collapse
|
41
|
Miller Y, Chaban GM, Gerber RB. Ab Initio Vibrational Calculations for H2SO4 and H2SO4·H2O: Spectroscopy and the Nature of the Anharmonic Couplings. J Phys Chem A 2005; 109:6565-74. [PMID: 16834003 DOI: 10.1021/jp058110l] [Citation(s) in RCA: 85] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Vibrational frequencies for fundamental, overtone, and combination excitations of sulfuric acid (H2SO4) and of sulfuric acid monohydrate cluster (H2SO4 x H2O) are computed directly from ab initio MP2/TZP potential surface points using the correlation-corrected vibrational self-consistent field (CC-VSCF) method, which includes anharmonic effects. The results are compared with experiment. The computed transitions show in nearly all cases good agreement with experimental data and consistent improvement over the harmonic approximation. The CC-VSCF improvements over the harmonic approximation are largest for the overtone and combination excitations and for the OH stretching fundamental. The agreement between the calculations and experiment also supports the validity of the MP2/TZP potential surfaces. Anharmonic coupling between different vibrational modes is found to significantly affect the vibrational frequencies. Analysis of the mean magnitude of the anharmonic coupling interactions between different pairs of normal modes is carried out. The results suggest possible mechanisms for the internal flow of vibrational energy in H2SO4 and H2SO4 x H2O.
Collapse
Affiliation(s)
- Y Miller
- Department of Physical Chemistry and Fritz Haber Research Center, The Hebrew University, Jerusalem 91904, Israel
| | | | | |
Collapse
|
42
|
Theoretical study of anharmonic vibrational spectra of HNO3, HNO3–H2O, HNO4: Fundamental, overtone and combination excitations. Chem Phys 2005. [DOI: 10.1016/j.chemphys.2005.01.012] [Citation(s) in RCA: 52] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|
43
|
Begue D, Carbonniere P, Pouchan C. Calculations of Vibrational Energy Levels by Using a Hybrid ab Initio and DFT Quartic Force Field: Application to Acetonitrile. J Phys Chem A 2005; 109:4611-6. [PMID: 16833799 DOI: 10.1021/jp0406114] [Citation(s) in RCA: 98] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
A hybrid quartic force field with quadratic force constants calculated at the CCSD(T)/cc-pVTZ level and cubic and quartic terms determined by a B3LYP/cc-pVTZ treatment is proposed to compute the vibrational energy levels of acetonitrile from a variational method. Fundamentals and overtones calculated in the range of 300-3200 cm(-1) are in fair agreement with the 31 observed data, with an absolute mean deviation of less than 0.8%. These results allow us to explain several Fermi resonances.
Collapse
Affiliation(s)
- D Begue
- Laboratoire de Chimie Théorique et Physico-Chimie Moléculaire - UMR 5624, Fédération de Recherche IPREM 2606, Université de Pau et des Pays de l'Adour, IFR - Rue Jules Ferry, 64000 PAU, France.
| | | | | |
Collapse
|
44
|
Diri K, Myshakin EM, Jordan KD. On the Contribution of Vibrational Anharmonicity to the Binding Energies of Water Clusters. J Phys Chem A 2005; 109:4005-9. [PMID: 16833722 DOI: 10.1021/jp050004w] [Citation(s) in RCA: 39] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
The second-order vibrational perturbation theory method has been used together with the B3LYP and MP2 electronic structure methods to investigate the effects of anharmonicity on the vibrational zero-point energy (ZPE) contributions to the binding energies of (H2O)n, n = 2-6, clusters. For the low-lying isomers of (H2O)6, the anharmonicity correction to the binding energy is calculated to range from -248 to -355 cm(-1). It is also demonstrated that although high-order electron correlation effects are important for the individual vibrational frequencies, they are relatively unimportant for the net ZPE contributions to the binding energies of water clusters.
Collapse
Affiliation(s)
- Kadir Diri
- Department of Chemistry and Center for Molecular and Materials Simulations, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, USA
| | | | | |
Collapse
|
45
|
Daniel Boese A, Klopper W, Martin * JML. Anharmonic force fields and thermodynamic functions using density functional theory. Mol Phys 2005. [DOI: 10.1080/00268970512331339369] [Citation(s) in RCA: 45] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
46
|
Czarnik-Matusewicz B, Rospenk M, Koll A, Mavri J. Influence of Substituents on the Anharmonicity of νs(OH) Vibration in Phenol Derivatives Explored by Experimental and Theoretical Approach. J Phys Chem A 2005; 109:2317-24. [PMID: 16839002 DOI: 10.1021/jp0455273] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Very good reproducibility of the first five vibrational transitions of phenol in the gas phase by the MP2/ 6-31G potential for O-H bond stretching was found. The vibrational levels were calculated by a program for variational solving of the time-independent Schrödinger equation in one dimension. Relative intensities of particular transitions were determined on the basis of the function of the dipole moment. The substituent effects on the nu(s)(OH) transitions and on the intensity of these transitions, as well as on the structure of eleven phenols, was analyzed as a function of the pK(a) values.
Collapse
|
47
|
Lynch VA, Mielke SL, Truhlar DG. Accurate vibrational-rotational partition functions and standard-state free energy values for H2O2 from Monte Carlo path-integral calculations. J Chem Phys 2004; 121:5148-62. [PMID: 15352807 DOI: 10.1063/1.1782511] [Citation(s) in RCA: 40] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Accurate quantum mechanical partition functions and absolute free energies of H(2)O(2) are determined using a realistic potential energy surface [J. Koput, S. Carter, and N. C. Handy, J. Phys. Chem. A 102, 6325 (1998)] for temperatures ranging from 300 to 2,400 K by using Monte Carlo path integral calculations with new, efficient polyatomic importance sampling methods. The path centroids are sampled in Jacobi coordinates via a set of independent ziggurat schemes. The calculations employed enhanced-same-path extrapolation of trapezoidal Trotter Fourier path integrals, and the paths were constructed using fast Fourier sine transforms. Importance sampling was also used in Fourier coefficient space, and adaptively optimized stratified sampling was used in configuration space. The free energy values obtained from the path-integral calculations are compared to separable-mode approximations, to the Pitzer-Gwinn approximation, and to values in thermodynamic tables. Our calculations support the recently proposed revisions to the JANAF tables.
Collapse
Affiliation(s)
- Vanessa Audette Lynch
- Department of Chemistry and Supercomputer Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, USA
| | | | | |
Collapse
|
48
|
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
|
49
|
Barone V. Vibrational zero-point energies and thermodynamic functions beyond the harmonic approximation. J Chem Phys 2004; 120:3059-65. [PMID: 15268458 DOI: 10.1063/1.1637580] [Citation(s) in RCA: 457] [Impact Index Per Article: 22.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
This paper compares harmonic and anharmonic zero-point energies and thermodynamic functions for a number of molecules of small and medium size. Anharmonic corrections cannot be neglected for quantitative studies, but can be obtained quite effectively by a perturbative treatment including cubic force constants to the second order and semidiagonal quartic constants to the first order. Simple finite difference equations provide all the necessary terms by at most 6N-11 Hessian evaluations, where N is the number of atoms in the system. Accurate values are obtained by this method using the Becke three parameter Lee-Yang-Parr functional, medium size basis sets, and, when needed, proper treatment of internal rotations. The whole model has been completely automated in the Gaussian package.
Collapse
Affiliation(s)
- Vincenzo Barone
- Laboratorio di Struttura e Dinamica Molecolare, Dipartimento di Chimica, Universita "Federico II," Complesso Universitario Monte Sant' Angelo, Via Cintia, I-80126 Napoli, Italy.
| |
Collapse
|
50
|
Benoit DM. Fast vibrational self-consistent field calculations through a reduced mode–mode coupling scheme. J Chem Phys 2004; 120:562-73. [PMID: 15267890 DOI: 10.1063/1.1631817] [Citation(s) in RCA: 96] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present a new methodology to perform fast correlation-corrected vibrational self-consistent field (CC-VSCF) calculations using ab initio potential energy points calculated on the fly. Our method is based on the replacement of all-electron basis sets with a pseudo-potential basis for heavy atoms, and on an efficient reduction of the number of pair-coupling elements used in the CC-VSCF procedure. The method is applied to several test systems: H2O, NH3, and CH4, where it proves to be efficient, providing a speedup factor of 2 compared to a standard CC-VSCF calculation. We also apply our technique to the simulation of the vibrational spectrum of ethane and show that very accurate results can be obtained with a substantial speedup for this system.
Collapse
Affiliation(s)
- David M Benoit
- Department of Chemistry, Queen Mary, University of London, Mile End Road, London, E1 4NS, United Kingdom.
| |
Collapse
|