Cis-->trans, trans-->cis isomerizations and N-O bond dissociation of nitrous acid (HONO) on an ab initio potential surface obtained by novelty sampling and feed-forward neural network fitting.

The isomerization and dissociation dynamics of HONO are investigated on an ab initio potential surface obtained by fitting the results of electronic structure calculations at 21 584 configurations by using previously described novelty sampling and feed-forward neural network (NN) methods. The electronic structure calculations are executed by using GAUSSIAN 98 with a 6-311G(d) basis set at the MP4(SDQ) level of accuracy. The average absolute error of the NN fits varies from 0.012 eV (1.22 kJ mol(-1)) to 0.017 eV (1.64 kJ mol(-1)). The average computation time for a HONO trajectory using a single NN surface is approximately 4.8 s. These computation times compare very favorably with those required by other methods primarily because the NN fitting needs to be executed only one time rather than at every integration point. If the average result obtained from a committee of NNs is employed at each point rather than a single NN, increased fitting accuracy can be achieved at the expense of increased computational requirements. In the present investigation, we find that a committee comprising five NN potentials reduces the average absolute interpolation error to 0.0111 eV (1.07 kJ mol(-1)). Cis-trans isomerization rates with total energy of 1.70 eV (including zero point energy) have been computed for a variety of different initial distributions of the internal energy. In contrast to results previously reported by using an empirical potential, where cis-->trans to trans-->cis rate coefficient ratios at 1.70 eV total energy were found to lie in the range of 2.0-12.9 depending on the vibration mode excited, these ratios on the ab initio NN potential lie in the range of 0.63-1.94. It is suggested that this result is a reflection of much larger intramode coupling terms present in the ab initio potential surface. A direct consequence of this increased coupling is a significant decrease in the mode specific rate enhancement when compared to results obtained by using empirical surfaces. All isomerizations are found to be first order in accordance with the results reported by using empirical potentials. The dissociation rate to NO+OH has been investigated at internal HONO energies of 3.10 and 3.30 eV for different distributions of this energy among the six vibrational modes of HONO. These dissociations are also found to be first order. The computed dissociation rate coefficients exhibit only modest mode specific rate enhancement that is significantly smaller than that obtained on an empirical surface because of the much larger mode couplings present on the ab initio surface.

[1]  Donald L Thompson,et al.  Interpolating moving least-squares methods for fitting potential energy surfaces: Improving efficiency via local approximants. , 2007, The Journal of chemical physics.

[2]  M. Malshe,et al.  Theoretical investigation of the dissociation dynamics of vibrationally excited vinyl bromide on an ab initio potential-energy surface obtained using modified novelty sampling and feedforward neural networks. II. Numerical application of the method. , 2007, The Journal of chemical physics.

[3]  Richard Dawes,et al.  Interpolating moving least-squares methods for fitting potential energy surfaces: computing high-density potential energy surface data from low-density ab initio data points. , 2007, The Journal of chemical physics.

[4]  Donald L Thompson,et al.  Interpolating moving least-squares methods for fitting potential energy surfaces: an application to the H2CN unimolecular reaction. , 2007, The Journal of chemical physics.

[5]  Ranga Komanduri,et al.  Molecular dynamics investigations of the dissociation of SiO2 on an ab initio potential energy surface obtained using neural network methods. , 2006, The Journal of chemical physics.

[6]  R Komanduri,et al.  Theoretical investigation of the dissociation dynamics of vibrationally excited vinyl bromide on an ab initio potential-energy surface obtained using modified novelty sampling and feed-forward neural networks. , 2006, The Journal of chemical physics.

[7]  R Komanduri,et al.  Ab initio potential-energy surfaces for complex, multichannel systems using modified novelty sampling and feedforward neural networks. , 2005, The Journal of chemical physics.

[8]  David S. Watkins,et al.  Fundamentals of Matrix Computations: Watkins/Fundamentals of Matrix Computations , 2005 .

[9]  Akio Kawano,et al.  Interpolating moving least-squares methods for fitting potential energy surfaces: applications to classical dynamics calculations. , 2004, The Journal of chemical physics.

[10]  Rifat Sonmez,et al.  Conceptual cost estimation of building projects with regression analysis and neural networks , 2004 .

[11]  Akio Kawano,et al.  Improving the accuracy of interpolated potential energy surfaces by using an analytical zeroth-order potential function. , 2004, The Journal of chemical physics.

[12]  H. Meyer,et al.  A study of the mode-selective trans--cis isomerization in HONO using ab initio methodology. , 2004, The Journal of chemical physics.

[13]  Donald L. Thompson,et al.  Interpolating moving least-squares methods for fitting potential energy surfaces: Detailed analysis of one-dimensional applications , 2003 .

[14]  D. Thompson,et al.  A theoretical study of cis–trans isomerization in HONO using an empirical valence bond potential , 2003 .

[15]  Michael A. Collins,et al.  Molecular potential-energy surfaces for chemical reaction dynamics , 2002 .

[16]  M. Räsänen,et al.  HONO in solid Kr: Site-selective trans↔cis isomerization with narrow-band infrared radiation , 2000 .

[17]  B. Kuhn,et al.  A new six-dimensional analytical potential up to chemically significant energies for the electronic ground state of hydrogen peroxide , 1999 .

[18]  H. Rabitz,et al.  A global H2O potential energy surface for the reaction O(1D)+H2→OH+H , 1996 .

[19]  Martin T. Hagan,et al.  Neural network design , 1995 .

[20]  Michael A. Collins,et al.  Molecular Potential Energy Surfaces by Interpolation , 1994, International Conference on Computational Science.

[21]  H. W. Schranz,et al.  An efficient microcanonical sampling procedure for molecular systems , 1991 .

[22]  D. Thompson,et al.  Mode specificity and the influence of rotation in cis-trans isomerization and dissociation in HONO , 1989 .

[23]  L. Raff Projection methods for obtaining intramolecular energy transfer rates from classical trajectory results: Application to 1,2‐difluoroethane , 1988 .

[24]  G. A. Petersson,et al.  A complete basis set model chemistry. I. The total energies of closed‐shell atoms and hydrides of the first‐row elements , 1988 .

[25]  Gillian C. Lynch,et al.  Intramolecular energy transfer and cis–trans isomerization in HONO , 1987 .

[26]  M. Räsänen,et al.  Ab initio calculations on hono: Energies, geometries and force fields on different levels of theory , 1985 .

[27]  J. Shirk,et al.  Isomerization of hono in solid nitrogen by selective vibrational excitation , 1983 .

[28]  J. Shirk,et al.  The infrared laser photoisomerization of HONO in solid N2 and Ar , 1982 .

[29]  Mark S. Gordon,et al.  The isomers of silacyclopropane , 1980 .

[30]  A. D. McLean,et al.  Contracted Gaussian basis sets for molecular calculations. I. Second row atoms, Z=11–18 , 1980 .

[31]  S. Nordholm,et al.  An efficient microcanonical sampling method , 1978 .

[32]  J. Pople,et al.  Self—Consistent Molecular Orbital Methods. XII. Further Extensions of Gaussian—Type Basis Sets for Use in Molecular Orbital Studies of Organic Molecules , 1972 .

[33]  J. Pople,et al.  Self‐Consistent Molecular‐Orbital Methods. IX. An Extended Gaussian‐Type Basis for Molecular‐Orbital Studies of Organic Molecules , 1971 .

[34]  M. Plesset,et al.  Note on an Approximation Treatment for Many-Electron Systems , 1934 .

[35]  J. Mason,et al.  Algorithms for approximation , 1987 .

[36]  Michael Baer,et al.  Theory of chemical reaction dynamics , 1985 .

[37]  J. Pople,et al.  Self‐consistent molecular orbital methods. XX. A basis set for correlated wave functions , 1980 .

[38]  P. C. Hariharan,et al.  Accuracy of AH n equilibrium geometries by single determinant molecular orbital theory , 1974 .

[39]  E. Bright Wilson,et al.  Book Reviews: Molecular Vibrations. The Theory of Infrared and Raman Vibrational Spectra , 1955 .

[40]  G. Herzberg,et al.  Spectra of diatomic molecules , 1950 .