Toward spectroscopically accurate global ab initio potential energy surface for the acetylene-vinylidene isomerization.

A new full-dimensional global potential energy surface (PES) for the acetylene-vinylidene isomerization on the ground (S0) electronic state has been constructed by fitting ∼37,000 high-level ab initio points using the permutation invariant polynomial-neural network method with a root mean square error of 9.54 cm(-1). The geometries and harmonic vibrational frequencies of acetylene, vinylidene, and all other stationary points (two distinct transition states and one secondary minimum in between) have been determined on this PES. Furthermore, acetylene vibrational energy levels have been calculated using the Lanczos algorithm with an exact (J = 0) Hamiltonian. The vibrational energies up to 12,700 cm(-1) above the zero-point energy are in excellent agreement with the experimentally derived effective Hamiltonians, suggesting that the PES is approaching spectroscopic accuracy. In addition, analyses of the wavefunctions confirm the experimentally observed emergence of the local bending and counter-rotational modes in the highly excited bending vibrational states. The reproduction of the experimentally derived effective Hamiltonians for highly excited bending states signals the coming of age for the ab initio based PES, which can now be trusted for studying the isomerization reaction.

[1]  T. Germann,et al.  Quantum mechanical calculation of resonance tunneling in acetylene isomerization via the vinylidene intermediate , 1998 .

[2]  J. Muckerman,et al.  MODE-SELECTIVE INFRARED EXCITATION OF LINEAR ACETYLENE , 1995 .

[3]  B. Orr.  Spectroscopy and energetics of the acetylene molecule: dynamical complexity alongside structural simplicity , 2006 .

[4]  H. Schaefer.  The 1,2 hydrogen shift: a common vehicle for the disappearance of evanescent molecular species , 1979 .

[5]  J. Rose,et al.  Bending dynamics from acetylene spectra: Normal, local, and precessional modes , 1996 .

[6]  Hua Guo,et al.  Full-dimensional quantum calculation of the vibrational energy / levels of hydrogen peroxide HOOH , 2000 .

[7]  Hans-Joachim Werner,et al.  A simple and efficient CCSD(T)-F12 approximation. , 2007, The Journal of chemical physics.

[8]  W. C. Lineberger,et al.  Observation of X1A1 vinylidene by photoelectron spectroscopy of the C2H2− ion , 1983 .

[9]  Jun Li,et al.  Quantum and quasi-classical dynamics of the OH + CO → H + CO2 reaction on a new permutationally invariant neural network potential energy surface. , 2014, The Journal of chemical physics.

[10]  Hans-Joachim Werner,et al.  Simplified CCSD(T)-F12 methods: theory and benchmarks. , 2009, The Journal of chemical physics.

[11]  Jun Li,et al.  Permutation invariant polynomial neural network approach to fitting potential energy surfaces. II. Four-atom systems. , 2013, The Journal of chemical physics.

[12]  Hua Guo,et al.  Probing highly excited vibrational eigenfunctions using a modified single Lanczos propagation method: Application to acetylene (HCCH) , 2003 .

[13]  Gengxin Chen,et al.  Approximate constants of motion and energy transfer pathways in highly excited acetylene , 1991 .

[14]  S. S. Sidhu,et al.  RRKM model of C2H4 dissociation: Heat of formation of vinylidene , 1989 .

[15]  J. V. Auwera,et al.  Vibrational Spectroscopic Database on Acetylene, X̃ 1Σg+ (12C2H2, 12C2D2, and 13C2H2) , 2003 .

[16]  Guohui Li,et al.  Full-dimensional quantum calculations of vibrational energy levels of acetylene (HCCH) up to 13,000 cm−1 , 2002 .

[17]  Bryan M. Wong,et al.  Electronic signatures of large amplitude motions: dipole moments of vibrationally excited local-bend and local-stretch states of S0 acetylene. , 2006, The journal of physical chemistry. B.

[18]  S. Carter,et al.  A potential energy surface for the ground state of acetylene, H2C2([Xtilde] 1Σ g +) , 1980 .

[19]  G. B. Park,et al.  Full dimensional Franck-Condon factors for the acetylene à (1)Au-X̃ (1)Σ(g)(+) transition. II. Vibrational overlap factors for levels involving excitation in ungerade modes. , 2014, The Journal of chemical physics.

[20]  R. Prosmiti,et al.  Periodic orbits, bifurcation diagrams and the spectroscopy of C2H2 system , 1995 .

[21]  J. L. Kinsey,et al.  Fluorescence and stimulated emission S1 → S0 spectra of acetylene: Regular and ergodic regions , 1985 .

[22]  H. Schaefer,et al.  Vinylidene: a very shallow minimum on the C2H2 potential energy surface , 1981 .

[23]  J. Gauss,et al.  An estimation of the isomerization energy of acetylene , 1999 .

[24]  J Grant Hill,et al.  Correlation consistent basis sets for molecular core-valence effects with explicitly correlated wave functions: the atoms B-Ne and Al-Ar. , 2010, The Journal of chemical physics.

[25]  R. Prosmiti,et al.  Periodic orbits and bifurcation diagrams of acetylene/vinylidene revisited , 2003 .

[26]  Bryan M. Wong,et al.  A new approach toward transition state spectroscopy. , 2013, Faraday discussions.

[27]  M. Jacobson,et al.  The infrared-ultraviolet dispersed fluorescence spectrum of acetylene: New classes of bright states , 2001 .

[28]  Mohammad Bagher Menhaj,et al.  Training feedforward networks with the Marquardt algorithm , 1994, IEEE Trans. Neural Networks.

[29]  H. Köppel,et al.  Barrier recrossing in the vinylidene–acetylene isomerization reaction: A five-dimensional ab initio quantum dynamical investigation , 2001 .

[30]  M. Herman,et al.  THE VIBRATIONAL ENERGY PATTERN IN ACETYLENE (IV) : UPDATED GLOBAL VIBRATION CONSTANTS FOR 12C2H2 , 1999 .

[31]  R. Silbey,et al.  Local mode behavior in the acetylene bending system , 1999 .

[32]  N. Handy,et al.  A Refined Quartic Forcefield for Acetylene: Accurate Calculation of the Vibrational Spectrum , 1993 .

[33]  R. Wyatt,et al.  A finite basis‐discrete variable representation calculation of vibrational levels of planar acetylene , 1992 .

[34]  J. Liévin,et al.  Vinylidene-acetylene cation isomerization investigated by large scale ab initio calculations. , 2006, The Journal of chemical physics.

[35]  J. Muckerman,et al.  Vibrational eigenvalues and eigenfunctions for planar acetylene by wave-packet propagation, and its mode-selective infrared excitation , 1997 .

[36]  B. Li,et al.  Full-dimensional quantum dynamics study of vinylidene-acetylene isomerization: a scheme using the normal mode Hamiltonian. , 2011, Physical chemistry chemical physics : PCCP.

[37]  I. Mills,et al.  Anharmonic force field of acetylene , 1976 .

[38]  A. McCoy,et al.  The bending dynamics of acetylene , 1996 .

[39]  Hua Guo,et al.  Six-dimensional quantum calculations of highly excited vibrational energy levels of hydrogen peroxide and its deuterated isotopomers , 2001 .

[40]  J. Stanton,et al.  High-accuracy estimates for the vinylidene-acetylene isomerization energy and the ground state rotational constants of :C═CH2. , 2013, The journal of physical chemistry. A.

[41]  William F. Polik,et al.  DISPERSED FLUORESCENCE SPECTRUM OF ACETYLENE FROM THE A 1AU ORIGIN : RECOGNITION OF POLYADS AND TEST OF MULTIRESONANT EFFECTIVE HAMILTONIAN MODEL FOR THE X STATE , 1996 .

[42]  T. Carrington,et al.  Calculating vibrational energies and wave functions of vinylidene using a contracted basis with a locally reorthogonalized coupled two-term Lanczos eigensolver. , 2006, The Journal of chemical physics.

[43]  M. E. Kellman,et al.  The dance of molecules: new dynamical perspectives on highly excited molecular vibrations. , 2007, Accounts of chemical research.

[44]  Darin C. Burleigh,et al.  Rotation–vibration interactions in highly excited states of SO2 and H2CO , 1991 .

[45]  J. L. Kinsey,et al.  Evidence of quantum ergodicity in stimulated emission pumping spectra of acetylene , 1985 .

[46]  Joel M. Bowman,et al.  Permutationally invariant potential energy surfaces in high dimensionality , 2009 .

[47]  M. Jacobson,et al.  State-by-state assignment of the bending spectrum of acetylene at 15 000 cm−1: A case study of quantum-classical correspondence , 1999 .

[48]  M. Jacobson,et al.  Anomalously slow intramolecular vibrational redistribution in the acetylene X̃ 1Σg+ state above 10 000 cm−1 of internal energy , 1998 .

[49]  P. Popelier,et al.  Potential energy surfaces fitted by artificial neural networks. , 2010, The journal of physical chemistry. A.

[50]  J. Bowman,et al.  Reduced dimensionality quantum calculations of acetylene↔vinylidene isomerization , 2002 .

[51]  Hua Guo,et al.  A scaled ab initio potential energy surface for acetylene and vinylidene , 2003 .

[52]  Hua Guo,et al.  Isotope effect in normal-to-local transition of acetylene bending modes. , 2012, The Journal of chemical physics.

[53]  R. Silbey,et al.  Intramolecular vibrational redistribution of energy in the stimulated emission pumping spectrum of acetylene , 1993 .

[54]  M. E. Kellman,et al.  Bifurcation phase diagram for C(2)H(2) bending dynamics has a tetracritical point with spectral patterns. , 2010, The journal of physical chemistry. A.

[55]  S. Carter,et al.  Potential models and local mode vibrational eigenvalue calculations for acetylene , 1982 .

[56]  R. Silbey,et al.  Pure bending dynamics in the acetylene X̃ 1Σg+ state up to 15 000 cm−1 of internal energy , 1998 .

[57]  J. Tennyson,et al.  Calculating energy levels of isomerizing tetra-atomic molecules. II. The vibrational states of acetylene and vinylidene. , 2005, The Journal of chemical physics.

[58]  H. Köppel,et al.  Ab initio quantum dynamical study of the vinylidene-acetylene isomerization , 1998 .

[59]  Joel M. Bowman,et al.  A new ab initio potential energy surface describing acetylene/vinylidene isomerization , 2003 .

[60]  J. L. Kinsey,et al.  High resolution spectroscopic detection of acetylene–vinylidene isomerization by spectral cross correlation , 1989 .

[61]  Hua Guo,et al.  Permutation invariant polynomial neural network approach to fitting potential energy surfaces. , 2013, The Journal of chemical physics.

[62]  W. C. Lineberger,et al.  A study of the singlet and triplet states of vinylidene by photoelectron spectroscopy of H2C=C−, D2C=C−, and HDC=C−. Vinylidene–acetylene isomerization , 1989 .

[63]  T. Carrington,et al.  Discrete‐Variable Representations and their Utilization , 2007 .

[64]  Hua Guo,et al.  A single Lanczos propagation method for calculating transition amplitudes. II. Modified QL and symmetry adaptation , 2001 .

[65]  Hua Guo.  Recursive Solutions to Large Eigenproblems in Molecular Spectroscopy and Reaction Dynamics , 2007 .

[66]  E. Sibert,et al.  Highly excited vibrational states of acetylene: A variational calculation , 1993 .

[67]  M. Jacobson,et al.  Acetylene at the Threshold of Isomerization , 2000 .

[68]  E. Sibert,et al.  Quantum, semiclassical and classical dynamics of the bending modes of acetylene , 1996 .

[69]  M. E. Kellman,et al.  Catastrophe map and the role of individual resonances in C2H2 bending dynamics. , 2009, The Journal of chemical physics.

[70]  W. Bian,et al.  Efficient quantum calculations of vibrational states of vinylidene in full dimensionality: a scheme with combination of methods. , 2008, The Journal of chemical physics.

[71]  T. Carrington,et al.  Vinylidene: Potential energy surface and unimolecular reaction dynamics , 1984 .

[72]  B. Li,et al.  Accurate quantum dynamics study on the resonance decay of vinylidene. , 2011, Chemphyschem : a European journal of chemical physics and physical chemistry.

[73]  J. K. Lundberg,et al.  VIBRATIONALLY HIGHLY EXCITED ACETYLENE AS STUDIED BY DISPERSED FLUORESCENCE AND STIMULATED EMISSION PUMPING SPECTROSCOPY , 1991 .

[74]  Zhen Xie,et al.  Permutationally Invariant Polynomial Basis for Molecular Energy Surface Fitting via Monomial Symmetrization. , 2010, Journal of chemical theory and computation.

[75]  A. Varandas,et al.  Accurate MRCI and CC study of the most relevant stationary points and other topographical attributes for the ground-state C(2)H(2) potential energy surface. , 2010, Journal of Physical Chemistry A.

[76]  Z. Vager,et al.  STUDY OF UNIMOLECULAR REACTIONS BY COULOMB EXPLOSION IMAGING : THE NONDECAYING VINYLIDENE , 1998 .

[77]  David Feller,et al.  An expanded calibration study of the explicitly correlated CCSD(T)-F12b method using large basis set standard CCSD(T) atomization energies. , 2013, The Journal of chemical physics.

[78]  Hua Guo,et al.  A nine-dimensional ab initio global potential energy surface for the H₂O⁺ + H₂ → H₃O⁺ + H reaction. , 2014, The Journal of chemical physics.

[79]  John E. Adams,et al.  Reaction path Hamiltonian for polyatomic molecules , 1980 .

[80]  J. Gauss,et al.  Vibrational structure in the vinylidene anion photoelectron spectrum: Closing the gap between theory and experiment , 1999 .

[81]  Robert W. Field,et al.  Numerical pattern recognition analysis of acetylene dispersed fluorescence spectra , 1998 .

[82]  Jun Li,et al.  A nine-dimensional global potential energy surface for NH4(X(2)A(1)) and kinetics studies on the H + NH3↔ H2 + NH2 reaction. , 2014, Physical chemistry chemical physics : PCCP.

[83]  N. Govind,et al.  Long Live Vinylidene! A New View of the H2CC: → HC⋮CH Rearrangement from ab Initio Molecular Dynamics , 2001 .

[84]  J. Behler.  Neural network potential-energy surfaces in chemistry: a tool for large-scale simulations. , 2011, Physical chemistry chemical physics : PCCP.

[85]  Hua Guo,et al.  A single Lanczos propagation method for calculating transition amplitudes , 1999 .