1
|
Wang XG, Carrington T. Computing excited OH stretch states of water dimer in 12D using contracted intermolecular and intramolecular basis functions. J Chem Phys 2023; 158:084107. [PMID: 36859104 DOI: 10.1063/5.0139586] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023] Open
Abstract
Due to the ubiquity and importance of water, water dimer has been intensively studied. Computing the (ro-)vibrational spectrum of water dimer is challenging. The potential has eight wells separated by low barriers, which makes harmonic approximations of limited utility. A variational approach is imperative, but difficult because there are 12 coupled vibrational coordinates. In this paper, we use a product contracted basis whose functions are products of intramolecular and intermolecular functions computed using an iterative eigensolver. An intermediate matrix F facilitates calculating matrix elements. Using F, it is possible to do calculations on a general potential without storing the potential on the full quadrature grid. We find that surprisingly many intermolecular functions are required. This is due to the importance of coupling between inter- and intra-molecular coordinates. The full G16 symmetry of water dimer is exploited. We calculate, for the first time, monomer excited stretch states and compare P(1) transition frequencies with their experimental counterparts. We also compare with experimental vibrational shifts and tunneling splittings. Surprisingly, we find that the largest tunneling splitting, which does not involve the interchange of the two monomers, is smaller in the asymmetric stretch excited state than in the ground state. Differences between levels we compute and those obtained with a [6+6]D adiabatic approximation [Leforestier et al. J. Chem. Phys. 137 014305 (2012)] are ∼0.6 cm-1 for states without monomer excitation, ∼4 cm-1 for monomer excited bend states, and as large as ∼10 cm-1 for monomer excited stretch states.
Collapse
Affiliation(s)
- Xiao-Gang Wang
- Chemistry Department, Queen's University, Kingston, Ontario K7L 3N6, Canada
| | - Tucker Carrington
- Chemistry Department, Queen's University, Kingston, Ontario K7L 3N6, Canada
| |
Collapse
|
2
|
Felker PM, Bačić Z. Erratum: "Intermolecular vibrational states of HF trimer from rigorous nine-dimensional quantum calculations: Strong coupling between intermolecular bending and stretching vibrations and the importance of the three-body interactions" [J. Chem. Phys. 157, 194103 (2022)]. J Chem Phys 2022; 157:239902. [PMID: 36550034 DOI: 10.1063/5.0137134] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Affiliation(s)
- Peter M Felker
- Department of Chemistry and Biochemistry, University of California, Los Angeles, California 90095-1569, USA
| | - Zlatko Bačić
- Department of Chemistry, New York University, New York, New York 10003, USA
| |
Collapse
|
3
|
Wang XG, Carrington T. Using nondirect product Wigner D basis functions and the symmetry-adapted Lanczos algorithm to compute the ro-vibrational spectrum of CH 4-H 2O. J Chem Phys 2021; 154:124112. [PMID: 33810654 DOI: 10.1063/5.0044010] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
By doing calculations on the methane-water van der Waals complex, we demonstrate that highly converged energy levels and wavefunctions can be obtained using Wigner D basis functions and the Symmetry-Adapted Lanczos (SAL) method. The Wigner D basis is a nondirect product basis and, therefore, efficient when the kinetic energy operator has accessible singularities. The SAL method makes it possible to exploit symmetry to label energy levels and reduce the cost of the calculation, without explicitly using symmetry-adapted basis functions. Line strengths are computed, and new bands are identified. In particular, we find unusually strong transitions between states associated with the isomers of the global minimum and the secondary minimum.
Collapse
Affiliation(s)
- Xiao-Gang Wang
- Chemistry Department, Queen's University, Kingston, Ontario K7L 3N6, Canada
| | - Tucker Carrington
- Chemistry Department, Queen's University, Kingston, Ontario K7L 3N6, Canada
| |
Collapse
|
4
|
Quintas-Sánchez E, Dawes R, Wang XG, Carrington T. Computational study of the rovibrational spectrum of CO 2-N 2. Phys Chem Chem Phys 2020; 22:22674-22683. [PMID: 33016299 DOI: 10.1039/d0cp04186f] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The CO2-N2 complex is formed from two key components of Earth's atmosphere, and as such, has received some attention from both experimental and theoretical studies. On the theory side, a potential energy surface (PES) based on high level ab initio data was reported [Nasri et al., J. Chem. Phys., 2015, 142, 174301] and then used in more recently reported rovibrational calculations [Lara-Moreno et al., Phys. Chem. Chem. Phys., 2019, 21, 3550]. Accuracy of about 1 percent was achieved for calculated rotational transitions of the ground vibrational state of the complex, compared with previously reported microwave spectra. However, a very recent measurement of the geared bending mode frequency [Barclay et al., J. Chem. Phys., 2020, 153, 014303] recorded a value of 21.4 cm-1, which is wildly different from the corresponding calculated value of 45.9 cm-1. To provide some insight into this discrepancy, we have constructed a new more accurate PES, and used it to perform highly converged variational rovibrational calculations. Our new results yield a value of 21.1 cm-1 for that bending frequency, in close agreement with the experiment. We also obtain significantly improved predicted rotational transitions. Finally, we note that a very shallow well, previously reported as a distinct second isomer, is not found on our new PES, but rather a transition structure is seen in that location.
Collapse
|
5
|
Simmons J, Wang XG, Carrington T. Computational Study of the Rovibrational Spectra of CH 2D + and CHD 2. J Phys Chem A 2019; 123:10281-10289. [PMID: 31657568 DOI: 10.1021/acs.jpca.9b09045] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
In this paper, we present rovibrational energy levels of CH2D+ and CHD2+. They are computed with a large basis and the Lanczos algorithm. CH2D+ and CHD2+ are believed to play an important role in interstellar space, but so far, there are no definitive observations. The predictions of this paper should facilitate detection. For CH2D+, two CH stretch bands have been studied at high resolution. Compared to our calculated energies, the root-mean-square error is 0.08 cm-1. For CHD2+, one CH stretch band has been studied at high resolution. Compared to our calculated energies, the root-mean-square error is 0.5 cm-1. Errors are larger, for both isotopologues, for bend states. We attribute these errors to the potential energy surface. Wave function and probability distribution plots are used to make assignments. The ν1 band of CHD2+ is significantly perturbed, and according to our calculations, the 3ν3 state is closest and might be the most important perturber.
Collapse
Affiliation(s)
- Jesse Simmons
- Department of Chemistry , and Department of Physics , Queen's University , Kingston , Ontario K7L 3N6 , Canada
| | - Xiao-Gang Wang
- Department of Chemistry , and Department of Physics , Queen's University , Kingston , Ontario K7L 3N6 , Canada
| | - Tucker Carrington
- Department of Chemistry , and Department of Physics , Queen's University , Kingston , Ontario K7L 3N6 , Canada
| |
Collapse
|
6
|
Castro-Juárez E, Wang XG, Carrington T, Quintas-Sánchez E, Dawes R. Computational study of the ro-vibrational spectrum of CO-CO 2. J Chem Phys 2019; 151:084307. [PMID: 31470713 DOI: 10.1063/1.5119762] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
An accurate ab initio ground-state intermolecular potential energy surface (PES) was determined for the CO-CO2 van der Waals dimer. The Lanczos algorithm was used to compute rovibrational energies on this PES. For both the C-in and O-in T-shaped isomers, the fundamental transition frequencies agree well with previous experimental results. We confirm that the in-plane states previously observed are geared states. In addition, we have computed and assigned many other vibrational states. The rotational constants we determine from J = 1 energy levels agree well with their experimental counterparts. Planar and out-of-plane cuts of some of the wavefunctions we compute are quite different, indicating strong coupling between the bend and torsional modes. Because the stable isomers are T-shaped, vibration along the out-of-plane coordinates is very floppy. In CO-CO2, when the molecule is out-of-plane, interconversion of the isomers is possible, but the barrier height is higher than the in-plane geared barrier height.
Collapse
Affiliation(s)
| | - Xiao-Gang Wang
- Chemistry Department, Queen's University, Kingston, Ontario K7L 3N6, Canada
| | - Tucker Carrington
- Chemistry Department, Queen's University, Kingston, Ontario K7L 3N6, Canada
| | - Ernesto Quintas-Sánchez
- Department of Chemistry, Missouri University of Science and Technology, Rolla, Missouri 65409, USA
| | - Richard Dawes
- Department of Chemistry, Missouri University of Science and Technology, Rolla, Missouri 65409, USA
| |
Collapse
|
7
|
Toombs E, Carrington T. Redundant coordinates in quantum mechanics. Chem Phys 2018. [DOI: 10.1016/j.chemphys.2018.07.047] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
8
|
Wang XG, Carrington T. Calculated rotation-bending energy levels of CH5+ and a comparison with experiment. J Chem Phys 2016; 144:204304. [DOI: 10.1063/1.4948549] [Citation(s) in RCA: 29] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Xiao-Gang Wang
- Chemistry Department, Queen’s University, Kingston, Ontario K7L 3N6, Canada
| | - Tucker Carrington
- Chemistry Department, Queen’s University, Kingston, Ontario K7L 3N6, Canada
| |
Collapse
|
9
|
Knaanie R, Šebek J, Tsuge M, Myllys N, Khriachtchev L, Räsänen M, Albee B, Potma EO, Gerber RB. Infrared Spectrum of Toluene: Comparison of Anharmonic Isolated-Molecule Calculations and Experiments in Liquid Phase and in a Ne Matrix. J Phys Chem A 2016; 120:3380-9. [DOI: 10.1021/acs.jpca.6b01604] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Roie Knaanie
- Institute
of Chemistry and The Fritz Haber Research Center, The Hebrew University, Jerusalem 91904, Israel
| | - Jiří Šebek
- Institute
of Chemistry and The Fritz Haber Research Center, The Hebrew University, Jerusalem 91904, Israel
- Department
of Physical Chemistry, University of Chemistry and Technology, Prague 16628, Czech Republic
| | - Masashi Tsuge
- Department
of Chemistry, University of Helsinki, P.O. Box 55, FI-00014 Helsinki, Finland
| | - Nanna Myllys
- Department
of Chemistry, University of Helsinki, P.O. Box 55, FI-00014 Helsinki, Finland
| | - Leonid Khriachtchev
- Department
of Chemistry, University of Helsinki, P.O. Box 55, FI-00014 Helsinki, Finland
| | - Markku Räsänen
- Department
of Chemistry, University of Helsinki, P.O. Box 55, FI-00014 Helsinki, Finland
| | - Brian Albee
- Department
of Chemistry, University of California, Irvine, California 92697, United States
| | - Eric O. Potma
- Department
of Chemistry, University of California, Irvine, California 92697, United States
| | - R. Benny Gerber
- Institute
of Chemistry and The Fritz Haber Research Center, The Hebrew University, Jerusalem 91904, Israel
- Department
of Chemistry, University of Helsinki, P.O. Box 55, FI-00014 Helsinki, Finland
- Department
of Chemistry, University of California, Irvine, California 92697, United States
| |
Collapse
|
10
|
Samanta AK, Wang Y, Mancini JS, Bowman JM, Reisler H. Energetics and Predissociation Dynamics of Small Water, HCl, and Mixed HCl–Water Clusters. Chem Rev 2016; 116:4913-36. [DOI: 10.1021/acs.chemrev.5b00506] [Citation(s) in RCA: 41] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Amit K. Samanta
- Department
of Chemistry, University of Southern California, Los Angeles, California 90089-0482, United States
| | - Yimin Wang
- Department
of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - John S. Mancini
- Department
of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - Joel M. Bowman
- Department
of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - Hanna Reisler
- Department
of Chemistry, University of Southern California, Los Angeles, California 90089-0482, United States
| |
Collapse
|
11
|
Wang XG, Carrington T. The vibration-rotation-tunneling levels of N2-H2O and N2-D2O. J Chem Phys 2015; 143:024303. [PMID: 26178101 DOI: 10.1063/1.4923339] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022] Open
Abstract
In this paper, we report vibration-rotation-tunneling levels of the van der Waals clusters N2-H2O and N2-D2O computed from an ab initio potential energy surface. The only dynamical approximation is that the monomers are rigid. We use a symmetry adapted Lanczos algorithm and an uncoupled product basis set. The pattern of the cluster's levels is complicated by splittings caused by H-H exchange tunneling (larger splitting) and N-N exchange tunneling (smaller splitting). An interesting result that emerges from our calculation is that whereas in N2-H2O, the symmetric H-H tunnelling state is below the anti-symmetric H-H tunnelling state for both K = 0 and K = 1, the order is reversed in N2-D2O for K = 1. The only experimental splitting measurements are the D-D exchange tunneling splittings reported by Zhu et al. [J. Chem. Phys. 139, 214309 (2013)] for N2-D2O in the v2 = 1 region of D2O. Due to the inverted order of the split levels, they measure the sum of the K = 0 and K = 1 tunneling splittings, which is in excellent agreement with our calculated result. Other splittings we predict, in particular those of N2-H2O, may guide future experiments.
Collapse
Affiliation(s)
- Xiao-Gang Wang
- Chemistry Department, Queen's University, Kingston, Ontario K7L 3N6, Canada
| | - Tucker Carrington
- Chemistry Department, Queen's University, Kingston, Ontario K7L 3N6, Canada
| |
Collapse
|
12
|
Abstract
I review two new ideas for coping with the size of large product basis sets and large product grids when one computes vibrational energy levels. The first is based on a tensor reduction scheme. It exploits advantages of a sum-of-products potential. The key idea is to use a basis each of whose function is a sum of optimized products and to compress the number of terms in each basis function. When the potential does not have the sum-of-products form, it is usually necessary to use quadrature. The second idea uses a nondirect product grid that has structure and is therefore compatible with efficient matrix–vector products.
Collapse
Affiliation(s)
- Tucker Carrington
- Chemistry Department, Queen’s University, Kingston, ON K7L 3N6, Canada
- Chemistry Department, Queen’s University, Kingston, ON K7L 3N6, Canada
| |
Collapse
|
13
|
Wang XG, Carrington T. Using experimental data and a contracted basis Lanczos method to determine an accurate methane potential energy surface from a least squares optimization. J Chem Phys 2014; 141:154106. [DOI: 10.1063/1.4896569] [Citation(s) in RCA: 49] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Xiao-Gang Wang
- Chemistry Department, Queen's University, Kingston, Ontario K7L 3N6, Canada
| | - Tucker Carrington
- Chemistry Department, Queen's University, Kingston, Ontario K7L 3N6, Canada
| |
Collapse
|
14
|
Roy TK, Carrington T, Gerber RB. Approximate first-principles anharmonic calculations of polyatomic spectra using MP2 and B3LYP potentials: comparisons with experiment. J Phys Chem A 2014; 118:6730-9. [PMID: 24977304 DOI: 10.1021/jp5060155] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Anharmonic vibrational spectroscopy calculations using MP2 and B3LYP computed potential surfaces are carried out for a series of molecules, and frequencies and intensities are compared with those from experiment. The vibrational self-consistent field with second-order perturbation correction (VSCF-PT2) is used in computing the spectra. The test calculations have been performed for the molecules HNO3, C2H4, C2H4O, H2SO4, CH3COOH, glycine, and alanine. Both MP2 and B3LYP give results in good accord with experimental frequencies, though, on the whole, MP2 gives very slightly better agreement. A statistical analysis of deviations in frequencies from experiment is carried out that gives interesting insights. The most probable percentage deviation from experimental frequencies is about -2% (to the red of the experiment) for B3LYP and +2% (to the blue of the experiment) for MP2. There is a higher probability for relatively large percentage deviations when B3LYP is used. The calculated intensities are also found to be in good accord with experiment, but the percentage deviations are much larger than those for frequencies. The results show that both MP2 and B3LYP potentials, used in VSCF-PT2 calculations, account well for anharmonic effects in the spectroscopy of molecules of the types considered.
Collapse
Affiliation(s)
- Tapta Kanchan Roy
- Institute of Chemistry and The Fritz Haber Research Center, The Hebrew University , Jerusalem 91904, Israel
| | | | | |
Collapse
|
15
|
Brown J, Wang XG, Carrington T, Grubbs GS, Dawes R. Computational study of the rovibrational spectrum of CO2–CS2. J Chem Phys 2014; 140:114303. [DOI: 10.1063/1.4867792] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
16
|
Mancini JS, Bowman JM. A New Many-Body Potential Energy Surface for HCl Clusters and Its Application to Anharmonic Spectroscopy and Vibration–Vibration Energy Transfer in the HCl Trimer. J Phys Chem A 2014; 118:7367-74. [DOI: 10.1021/jp412264t] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- John S. Mancini
- Cherry
L. Emerson Center
for Scientific Computation and Department of Chemistry, Emory University, Atlanta, Georgia 30322, United States
| | - Joel M. Bowman
- Cherry
L. Emerson Center
for Scientific Computation and Department of Chemistry, Emory University, Atlanta, Georgia 30322, United States
| |
Collapse
|
17
|
Asselin P, Soulard P, Madebène B, Goubet M, Huet TR, Georges R, Pirali O, Roy P. The cyclic ground state structure of the HF trimer revealed by far infrared jet-cooled Fourier transform spectroscopy. Phys Chem Chem Phys 2014; 16:4797-806. [DOI: 10.1039/c3cp55047h] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
|
18
|
Brown J, Wang XG, Carrington T. Calculating and assigning rovibrational energy levels of (15N2O)2, (15N14NO)2, 14N2O-15N2O and 15N14NO-15N2O. Phys Chem Chem Phys 2013; 15:19159-68. [PMID: 24104969 DOI: 10.1039/c3cp52548a] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
In this paper we report transition frequencies and rotational constants computed for several isotopologues of the nitrous oxide dimer. A previously reported intermolecular potential, the symmetry adapted Lanczos algorithm and an uncoupled product basis set are used to do the calculations. Rotational transition frequencies and rotational constants are in good agreement with experiment. We calculate states localized in both polar and nonpolar wells on the potential surface. Two of the four isotopologues we study have inequivalent monomers. They have wavefunctions localized over a single polar well.
Collapse
Affiliation(s)
- James Brown
- Queen's University, Kingston, Ontario K7L 3N6, Canada
| | | | | |
Collapse
|
19
|
Pesonen J. Constrained molecular vibration-rotation Hamiltonians: contravariant metric tensor. J Chem Phys 2013; 139:144310. [PMID: 24116622 DOI: 10.1063/1.4823835] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Here, I present a practical recipe for obtaining contravariant vibration-rotation metric tensors, and thus the kinetic energy operators, when some degrees of freedom are constrained rigidly. An element of the contravariant metric tensor is obtained as a sum of dot products of contravariant measuring vectors, which are obtained from their unconstrained counterparts by adding a frozen mode correction. The present method applies in principle for any choice of shape coordinates and a body-frame for which the contravariant measuring vectors can be evaluated. In contrast to the existing methods, the present method does not involve evaluation of covariant metric tensors, matrix inversions, chain rules of derivation, or numerical differentiation. It is applied in the sequel paper [L. Partanen, J. Pesonen, E. Sjöholm, and L. Halonen, J. Chem. Phys. 139, 144311 (2013)] to study the effects of several different approximations to the kinetic energy operator, when the two large-amplitude OH-torsional motions in H2SO4 are of interest.
Collapse
Affiliation(s)
- Janne Pesonen
- Department of Chemistry, University of Helsinki, P.O. Box 55 (A. I. Virtasen aukio 1), Helsinki FIN-00014, Finland
| |
Collapse
|
20
|
Avila G, Carrington T. Solving the Schroedinger equation using Smolyak interpolants. J Chem Phys 2013; 139:134114. [DOI: 10.1063/1.4821348] [Citation(s) in RCA: 37] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
21
|
Wang XG, Carrington T. Computing rovibrational levels of polyatomic molecules with polyspherical coordinates and a contracted basis built with a K-independent vibrational primitive basis. Mol Phys 2013. [DOI: 10.1080/00268976.2013.808387] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Affiliation(s)
- Xiao-Gang Wang
- Chemistry Department, Queen’s University , Kingston, Ontario, Canada
| | - Tucker Carrington
- Chemistry Department, Queen’s University , Kingston, Ontario, Canada
| |
Collapse
|
22
|
Wang XG, Carrington T. Computing rovibrational levels of methane with curvilinear internal vibrational coordinates and an Eckart frame. J Chem Phys 2013; 138:104106. [DOI: 10.1063/1.4793474] [Citation(s) in RCA: 66] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
23
|
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]
|
24
|
Avila G, Carrington T. Solving the vibrational Schrödinger equation using bases pruned to include strongly coupled functions and compatible quadratures. J Chem Phys 2012; 137:174108. [DOI: 10.1063/1.4764099] [Citation(s) in RCA: 54] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
25
|
|
26
|
Brown J, Wang XG, Dawes R, Carrington T. Computational study of the rovibrational spectrum of (OCS)2. J Chem Phys 2012; 136:134306. [DOI: 10.1063/1.3697679] [Citation(s) in RCA: 50] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
|
27
|
|
28
|
GUO HUA, CHEN RONGQING, XIE DAIQIAN. CALCULATION OF TRANSITION AMPLITUDES WITH A SINGLE LANCZOS PROPAGATION. JOURNAL OF THEORETICAL & COMPUTATIONAL CHEMISTRY 2012. [DOI: 10.1142/s0219633602000129] [Citation(s) in RCA: 36] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
We review in this article a recently proposed energy-global method that is capable of calculating the entire transition amplitude matrix with a single Lanczos propagation. This method requires neither explicit computation nor storage of the eigenfunctions, rendering it extremely memory efficient. Procedures are proposed to handle situations where "spurious" eigenvalues aggregate around true eigenvalues due to round-off errors. This method is amenable to both real-symmetric and complex-symmetric Hamiltonians. Applications to molecular spectra and reactive scattering are presented. Its relationships with other methods are also discussed.
Collapse
Affiliation(s)
- HUA GUO
- Department of Chemistry, University of New Mexico, Albuquerque, New Mexico 87131, USA
| | - RONGQING CHEN
- Department of Chemistry, University of New Mexico, Albuquerque, New Mexico 87131, USA
| | - DAIQIAN XIE
- Department of Chemistry and Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing, Jiangsu, 210093, P. R. China
| |
Collapse
|
29
|
GATTI FABIEN, IUNG CHRISTOPHE. EXACT AND CONSTRAINED KINETIC ENERGY OPERATORS IN POLYSPHERICAL COORDINATES. JOURNAL OF THEORETICAL & COMPUTATIONAL CHEMISTRY 2011. [DOI: 10.1142/s0219633603000653] [Citation(s) in RCA: 35] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
In this paper the derivation of exact and constrained kinetic energy operators in curvilinear coordinates is presented in a very general context. The polyspherical approach allows us to provide simple and explicit expressions of these operators which are well adapted to the numerical approaches used to solve the Schrödinger equation. These expressions can be exploited to treat a lot of problems such as the calculation of infrared or photoabsorption spectra as well as the study of intramolecular vibrational relaxation of molecules or the study of reactive scattering systems. Special emphasis is placed on concrete applications.
Collapse
Affiliation(s)
- FABIEN GATTI
- Laboratoire Structure et Dynamique des Systèmes Moléculaires et Solides, (UMR 5636), CC 014, Université des Sciences et Techniques du Languedoc, 34095 Montpellier Cedex 05, France
| | - CHRISTOPHE IUNG
- Laboratoire Structure et Dynamique des Systèmes Moléculaires et Solides, (UMR 5636), CC 014, Université des Sciences et Techniques du Languedoc, 34095 Montpellier Cedex 05, France
| |
Collapse
|
30
|
Carrington T, Wang XG. Computing ro-vibrational spectra of van der Waals molecules. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2011. [DOI: 10.1002/wcms.73] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
|
31
|
Avila G, Carrington T. Using a pruned basis, a non-product quadrature grid, and the exact Watson normal-coordinate kinetic energy operator to solve the vibrational Schrödinger equation for C2H4. J Chem Phys 2011; 135:064101. [DOI: 10.1063/1.3617249] [Citation(s) in RCA: 109] [Impact Index Per Article: 8.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
32
|
Avila G, Carrington T. Using nonproduct quadrature grids to solve the vibrational Schrödinger equation in 12D. J Chem Phys 2011; 134:054126. [DOI: 10.1063/1.3549817] [Citation(s) in RCA: 99] [Impact Index Per Article: 7.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
33
|
Wang XG, Carrington T. Theoretical study of the rovibrational spectrum of H2O–H2. J Chem Phys 2011; 134:044313. [DOI: 10.1063/1.3533230] [Citation(s) in RCA: 50] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
34
|
|
35
|
Dawes R, Wang XG, Jasper AW, Carrington T. Nitrous oxide dimer: A new potential energy surface and rovibrational spectrum of the nonpolar isomer. J Chem Phys 2010; 133:134304. [DOI: 10.1063/1.3494542] [Citation(s) in RCA: 85] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
36
|
Abstract
We report calculated microwave and infrared rovibrational transitions of the van der Waals complex He2–OCS. The calculations were done using a product basis, a Lanczos eigensolver, and potentials built from He–OCS, and He–He potential functions taken from the literature. All five of the large amplitude coordinates are treated exactly and calculations are done for J values up to five. All rovibrational levels are converged to 0.001 cm–1 by using basis sets with as many as 87 million funcions. Good agreement is found with previously reported experimental results. Although we assume that the dipole moment is along the OCS axis, we find transitions with appreciable intensity between different torsion states.
Collapse
Affiliation(s)
- Xiao-Gang Wang
- Chemistry Department, Queen’s University, Kingston, ON K7L 3N6, Canada
| | | |
Collapse
|
37
|
Wang XG, Carrington T, McKellar ARW. Theoretical and experimental study of the rovibrational spectrum of He(2)-CO. J Phys Chem A 2009; 113:13331-41. [PMID: 19842654 DOI: 10.1021/jp904778f] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We report calculated microwave and infrared ro-vibrational transitions of the van der Waals complex He(2)-CO. The calculations were done using a product basis and a Lanczos eigensolver, together with He-CO and He-He potential functions taken from the literature. The results are found to be in good agreement with previously reported experimental results, and they enable the experimental assignments to be clarified, augmented, and (in one case) corrected. Unlike some other van der Waals complexes with two He atoms such as He(2)-N(2)O, it is not possible to associate a set of energy levels with the "torsional" motion of the two He atoms on a ring encircling the dopant (in this case CO). Although we assume that the dipole moment is along the CO axis we find nonetheless that many transitions have appreciable intensity due to ro-vibrational coupling.
Collapse
Affiliation(s)
- Xiao-Gang Wang
- Chemistry Department, Queen's University, Kingston, Ontario K7L 3N6, Canada
| | | | | |
Collapse
|
38
|
Avila G, Carrington T. Nonproduct quadrature grids for solving the vibrational Schrödinger equation. J Chem Phys 2009; 131:174103. [DOI: 10.1063/1.3246593] [Citation(s) in RCA: 91] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
39
|
|
40
|
|
41
|
Bowman JM, Carrington T, Meyer HD. Variational quantum approaches for computing vibrational energies of polyatomic molecules. Mol Phys 2008. [DOI: 10.1080/00268970802258609] [Citation(s) in RCA: 303] [Impact Index Per Article: 18.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
|
42
|
Tremblay JC, Beyvers S, Saalfrank P. Selective excitation of coupled CO vibrations on a dissipative Cu(100) surface by shaped infrared laser pulses. J Chem Phys 2008; 128:194709. [DOI: 10.1063/1.2916710] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
43
|
Wang XG, Carrington T. Using a nondirect product discrete variable representation for angular coordinates to compute vibrational levels of polyatomic molecules. J Chem Phys 2008; 128:194109. [DOI: 10.1063/1.2918498] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
44
|
Recursive Solutions to Large Eigenproblems in Molecular Spectroscopy and Reaction Dynamics. REVIEWS IN COMPUTATIONAL CHEMISTRY 2007. [DOI: 10.1002/9780470189078.ch7] [Citation(s) in RCA: 49] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/21/2023]
|
45
|
Manzhos S, Carrington T. Using neural networks to represent potential surfaces as sums of products. J Chem Phys 2007; 125:194105. [PMID: 17129087 DOI: 10.1063/1.2387950] [Citation(s) in RCA: 129] [Impact Index Per Article: 7.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
By using exponential activation functions with a neural network (NN) method we show that it is possible to fit potentials to a sum-of-products form. The sum-of-products form is desirable because it reduces the cost of doing the quadratures required for quantum dynamics calculations. It also greatly facilitates the use of the multiconfiguration time dependent Hartree method. Unlike potfit product representation algorithm, the new NN approach does not require using a grid of points. It also produces sum-of-products potentials with fewer terms. As the number of dimensions is increased, we expect the advantages of the exponential NN idea to become more significant.
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
|
46
|
Wugt Larsen R, Suhm MA. Cooperative organic hydrogen bonds: The librational modes of cyclic methanol clusters. J Chem Phys 2006; 125:154314. [PMID: 17059262 DOI: 10.1063/1.2358349] [Citation(s) in RCA: 38] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022] Open
Abstract
Intermolecular hydrogen bond libration modes of isolated cyclic methanol trimers (approximately 613 cm(-1)) and tetramers (695 and 760 cm(-1)) are observed in pulsed jet Fourier transform infrared spectra and found to exhibit sizeable anharmonicity and mode coupling effects, opening the way for a microscopic interpretation of the broad librational bands of alcohols. The correlation of experimental OH stretching and OH libration band intensities provides important constraints for theoretical band strengths, cluster densities, and size assignments.
Collapse
Affiliation(s)
- R Wugt Larsen
- Institut für Physikalische Chemie, Universität Göttingen, Tammannstrasse 6, 37077 Göttingen, Germany
| | | |
Collapse
|
47
|
Tremblay JC, Carrington T. Calculating vibrational energies and wave functions of vinylidene using a contracted basis with a locally reorthogonalized coupled two-term Lanczos eigensolver. J Chem Phys 2006; 125:094311. [PMID: 16965084 DOI: 10.1063/1.2234480] [Citation(s) in RCA: 83] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We use a contracted basis+Lanczos eigensolver approach to compute vinylidene-like vibrational states of the acetylene-vinylidene system. To overcome problems caused by loss of orthogonality of the Lanczos vectors we reorthogonalize Lanczos vector and use a coupled two-term approach. The calculations are done in CC-HH diatom-diatom Jacobi coordinates which make it easy to compute states one irreducible representation at a time. The most costly parts of the calculation are parallelized and scale well. We estimate that the vinylidene energies we compute are converged to approximately 1 cm(-1).
Collapse
Affiliation(s)
- Jean Christophe Tremblay
- Département de chimie, Université de Montréal, C.P. 6128, succursale Centre-ville, Montréal (Québec) H3C 3J7, Canada.
| | | |
Collapse
|
48
|
Abstract
This article reviews new methods for computing vibrational energy levels of small polyatomic molecules. The principal impediment to the calculation of energy levels is the size of the required basis set. If one uses a product basis the Hamiltonian matrix for a four-atom molecule is too large to store in core memory. We discuss iterative methods that enable one to use a product basis to compute energy levels (and spectra) without storing a Hamiltonian matrix. Despite the advantages of iterative methods it is not possible, using product basis functions, to calculate vibrational spectra of molecules with more than four atoms. A very recent method combining contracted basis functions and the Lanczos algorithm with which vibrational energy levels of methane have been computed is described. New ideas, based on exploiting preconditioning, for reducing the number of matrix-vector products required to converge energy levels of interest are also summarized.Key words: vibrational energy levels, kinetic energy operators, Lanczos algorithm, contracted basis functions, preconditioning.
Collapse
|
49
|
Yu HG. Full-dimensional quantum calculations of vibrational spectra of six-atom molecules. I. Theory and numerical results. J Chem Phys 2004; 120:2270-84. [PMID: 15268366 DOI: 10.1063/1.1636456] [Citation(s) in RCA: 41] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Two quantum mechanical Hamiltonians have been derived in orthogonal polyspherical coordinates, which can be formed by Jacobi and/or Radau vectors etc., for the study of the vibrational spectra of six-atom molecules. The Hamiltonians are expressed in an explicit Hermitian form in the spatial representation. Their matrix representations are described in both full discrete variable representation (DVR) and mixed DVR/nondirect product finite basis representation (FBR) bases. The two-layer Lanczos iteration algorithm [H.-G. Yu, J. Chem. Phys. 117, 8190 (2002)] is employed to solve the eigenvalue problem of the system. A strategy regarding how to carry out the Hamiltonian-vector products for a high-dimensional problem is discussed. By exploiting the inversion symmetry of molecules, a unitary sequential 1D matrix-vector multiplication algorithm is proposed to perform the action of the Hamiltonian on the wavefunction in a symmetrically adapted DVR or FBR basis in the azimuthal angular variables. An application to the vibrational energy levels of the molecular hydrogen trimer (H2)3 in full dimension (12D) is presented. Results show that the rigid-H2 approximation can underestimate the binding energy of the trimer by 27%. Finally, it is demonstrated that the two-layer Lanczos algorithm is also capable of computing the eigenvectors of the system with minor effort.
Collapse
Affiliation(s)
- Hua-Gen Yu
- Department of Chemistry, Brookhaven National Laboratory, Upton, New York 11973-5000, USA.
| |
Collapse
|
50
|
Montgomery J, Poirier B. Eigenspectra calculations using Cartesian coordinates and a rotational symmetry adapted Lanczos method. J Chem Phys 2003. [DOI: 10.1063/1.1599352] [Citation(s) in RCA: 15] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023] Open
|