Biomedical subjects

Kirk A Peterson

Publications and source records attributed to Kirk A Peterson.

At least 19 recordsLinked to original sources

PubMed2006

On the spectroscopic and thermochemical properties of ClO, BrO, IO, and their anions.

A coupled cluster composite approach has been used to accurately determine the spectroscopic constants, bond dissociation energies, and heats of formation for the X1(2)II(3/2) states of the halogen oxides ClO, BrO, and IO, as well as their negative ions ClO-, BrO-, and IO-. After determining the frozen core, complete basis set (CBS) limit CCSD(T) values, corrections were added for core-valence correlation, relativistic effects (scalar and spin-orbit), the pseudopotential approximation (BrO and IO), iterative connected triple excitations (CCSDT), and iterative quadruples (CCSDTQ). The final ab initio equilibrium bond lengths and harmonic frequencies for ClO and BrO differ from their accurate experimental values by an average of just 0.0005 A and 0.8 cm-1, respectively. The bond length of IO is overestimated by 0.0047 A, presumably due to an underestimation of molecular spin-orbit coupling effects. Spectroscopic constants for the spin-orbit excited X2(2)III(1/2) states are also reported for each species. The predicted bond lengths and harmonic frequencies for the closed-shell anions are expected to be accurate to within about 0.001 A and 2 cm-1, respectively. The dissociation energies of the radicals have been determined by both direct calculation and through use of negative ion thermochemical cycles, which made use of a small amount of accurate experimental data. The resulting values of D0, 63.5, 55.8, and 54.2 kcal/mol for ClO, BrO, and IO, respectively, are the most accurate ab initio values to date, and those for ClO and BrO differ from their experimental values by just 0.1 kcal/mol. These dissociation energies lead to heats of formation, DeltaH(f) (298 K), of 24.2 +/- 0.3, 29.6 +/- 0.4, and 29.9 +/- 0.6 kcal/mol for ClO, BrO, and IO, respectively. Also, the final calculated electron affinities are all within 0.2 kcal/mol of their experimental values. Improved pseudopotential parameters for the iodine atom are also reported, together with revised correlation consistent basis sets for this atom.

Journal Article↗

PubMed2006

The nature of halogen...halogen synthons: crystallographic and theoretical studies.

A study of the halogen...halogen contacts in organic compounds using ab initio calculations and the results of previously reported crystallographic studies show that these interactions are controlled by electrostatics. These contacts can be represented by the geometric parameters of the C--X1...X2--C moieties (where theta1=C--X1...X2 and theta2=X1...X2--C; ri=X1...X2 distance). The distributions of the contacts within the sum of van der Waals radii (rvdW) versus thetai (theta1=theta2) show a maximum at theta approximately 150 degrees for X=Cl, Br, and I. This maximum is not seen in the distribution of F...F contacts. These results are in good agreement with our ab initio calculations. The theoretical results show that the position of the maximum depends on three factors: 1) The type of halogen atom, 2) the hybridization of the ipso carbon atom, and 3) the nature of the other atoms that are bonded to the ipso carbon atom apart from the halogen atom. Calculations show that the strength of these contacts decreases in the following order: I...I>Br...Br>Cl...Cl. Their relative strengths decrease as a function of the hybridization of the ipso carbon atom in the following order: sp2>sp>sp3. Attaching an electronegative atom to the carbon atom strengthens the halogen...halogen contacts. An electrostatic model is proposed based on two assumptions: 1) The presence of a positive electrostatic end cap on the halogen atom (except for fluorine) and 2) the electronic charge is anisotropically distributed around the halogen atom.

Algorithms↗

PubMed2006

Chemically accurate thermochemistry of cadmium: An ab initio study of Cd + XY (X = H, O, Cl, Br; Y = Cl, Br).

Using a composite coupled cluster method employing sequences of correlation consistent basis sets for complete basis set (CBS) extrapolations and with explicit treatment of core-valence correlation and scalar and spin-orbit relativistic effects, the 0 K enthalpies of a wide range of cadmium-halide reactions, namely, Cd + (HCl, HBr, ClO, BrO, Cl2, BrCl, Br2) have been determined to an estimated accuracy of +/-1 kcal/mol. In addition, accurate equilibrium geometries, harmonic frequencies, and dissociation energies have been calculated at the same level of theory for all the diatomic (e.g., CdH, CdO, CdCl, CdBr) and triatomic (CdHCl, CdHBr, CdClO, CdBrO, CdCl2, CdBrCl, CdBr2) species involved in these reactions, some for the very first time. Like their mercury analogues, all of the abstraction reactions are predicted to be endothermic, while the insertion reactions are strongly exothermic with the formation of stable linear, Cd-centric complexes. With the exception of CdH and the reactions involving this species, the present results for the remaining Cd-containing systems are believed to be the most accurate to date.

Journal Article↗

PubMed2006

An ab initio study of the low-lying electronic states of S3.

Accurate calculations of the low-lying singlet and triplet electronic states of thiozone, S(3), have been carried out using large multireference configuration interaction wave functions. Cuts of the full potential energy surfaces along the stretching and bending coordinates have been presented, together with the vertical excitation spectra. The strong experimentally observed absorption around 395 nm is assigned to the 1 (1)B(2) state, which correlates to ground state products. Absorption at wavelengths shorter than 260 nm is predicted to lead to singlet excited state products, S(2) (a (1)Delta(g))+S((1)D). The spectroscopic properties of the X (3)Sigma(g) (-), a (1)Delta(g), and b (1)Sigma(g) (+) electronic states of the S(2) radical have also been accurately characterized in this work. The investigations of the low-lying electronic states were accompanied by accurate ground state coupled cluster calculations of the thermochemistry of both S(2) and S(3) using large correlation consistent basis sets with corrections for core-valence correlation, scalar relativity, and atomic spin-orbit effects. Resulting values for D(0)(S(2)+S) and SigmaD(0) for S(3) are predicted to be 61.3 and 162.7 kcal/mol, respectively, with conservative uncertainties of +/-1 kcal/mol. Analogous calculations predict the C(2v)-D(3h) (open-cyclic) isomerization energy of S(3) to be 4.4+/-0.5 kcal/mol.

Journal Article↗

PubMed2006

Basis set limit electronic excitation energies, ionization potentials, and electron affinities for the 3d transition metal atoms: Coupled cluster and multireference methods.

Recently developed correlation consistent basis sets for the first row transition metal elements Sc-Zn have been utilized to determine complete basis set (CBS) scalar relativistic electron affinities, ionization potentials, and 4s(2)3d(n-2)-4s(1)d(n-1) electronic excitation energies with single reference coupled cluster methods [CCSD(T), CCSDT, and CCSDTQ] and multireference configuration interaction with three reference spaces: 3d4s, 3d4s4p, and 3d4s4p3d'. The theoretical values calculated with the highest order coupled cluster techniques at the CBS limit, including extrapolations to full configuration interaction, are well within 1 kcal/mol of the corresponding experimental data. For the early transition metal elements (Sc-Mn) the internally contracted multireference averaged coupled pair functional method yielded excellent agreement with experiment; however, the atomic properties for the late transition metals (Mn-Zn) proved to be much more difficult to describe with this level of theory, even with the largest reference function of the present work.

Journal Article↗

PubMed2006

Ab initio prediction of the potential energy surface and vibration-rotation energy levels of BeH2.

The equilibrium structure and potential energy surface of beryllium dihydride BeH(2) in its ground electronic state have been determined from highly accurate ab initio calculations. The vibration-rotation energy levels of three isotopomers BeH(2), BeD(2), and BeHD were predicted using the variational method. The calculated spectroscopic constants are in remarkably good agreement with the existing experimental data (sub-cm(-1) accuracy) and should be useful in a further analysis of high-resolution vibration-rotation spectra of all three isotopomers.

Journal Article↗

PubMed2006

Sources of error in electronic structure calculations on small chemical systems.

The sources of error in electronic structure calculations arising from the truncation of the one-particle and n-particle expansions are examined with very large correlation consistent basis sets, in some cases up through valence 10-zeta quality, and coupled cluster methods, up through connected quadruple excitations. A limited number of full configuration interaction corrections are also considered. For cases where full configuration interaction calculations were unavailable or prohibitively expensive, a continued fraction approximation was used. In addition, errors arising from corevalence and relativistic corrections are also probed for a number of small chemical systems. The accuracies of several formulas for estimating total energies and atomization energies in the complete basis set limit are compared in light of the present large basis set findings. In agreement with previous work, the CCSD(T) method is found to provide results that are closer to the CCSDTQ and full configuration-interaction results than the less approximate CCSDT method.

Journal Article↗

PubMed2006

The low-lying electronic states of nickel cyanide and isocyanide: A theoretical investigation.

At different levels of coupled cluster theory optimum structures, energetics, and harmonic vibrational frequencies for several low-lying doublet and quartet electronic states of linear NiCN and NiNC were studied using four contracted Gaussian basis sets, ranging from Ni[6s5p4d2f], CN[4s3p2d] to Ni[8s7p5d3f2g1h], CN[5s4p3d2f1g]. The most reliable predictions were obtained with a relativistic Douglas-Kroll restricted open-shell-based coupled cluster method including singles, doubles, and perturbative triple excitations [DK-R/UCCSD(T)]. This level of theory was used in conjunction with correlation-consistent polarized valence Douglas-Kroll recontracted quadruple-zeta basis sets (cc-pVQZDK). The energetic ordering of the electronic states of NiCN is predicted to be 2delta < 2sigma+ < 2pi < 4delta < 4pi and that of NiNC is 2delta approximately 2sigma+ < 2pi < 4delta < 4pi < 4sigma-. Our theoretical investigation supports the assignment of the ground-state term symbol, the Ni-C stretching frequency, and the bending frequency for the ground electronic state of NiCN by Kingston et al. [J. Mol. Spectrosc. 215, 106 (2002)] and by Sheridan and Ziurys [J. Chem. Phys. 118, 6370 (2003)]. The predicted structure of the 2delta ground state of NiCN, r(e)(Ni-C) = 1.822 angstroms and r(e)(C-N) = 1.167 angstroms, at DK-R/UCCSD(T)/cc-pVQZDK shows excellent agreement with the experimentally determined Ni-C bond length of 1.826 A and less satisfactory agreement for the C-N bond length of 1.153 angstroms [J. Chem. Phys. 118, 6370 (2003)]. It is also concluded that the metal-to-ligand pi back donation is weak or negligible. Additionally, we found that on the 2delta surface the linear cyanide isomer lies lower in energy than the linear isocyanide isomer by 12.2 kcal mol(-1).

Journal Article↗

PubMed2005

Ab initio thermochemistry involving heavy atoms: an investigation of the reactions Hg + IX (X = I, Br, Cl, O).

Accurate 0 K enthalpies have been calculated for reactions of mercury with a series of small iodine-containing molecules (I2, IBr, ICl, and IO). The calculations have been carried out with the coupled cluster singles and doubles method with a perturbative correction for connected triple excitations [CCSD(T)] using sequences of correlation consistent basis sets and accurate relativistic pseudopotentials. Corrections have been included to account for core-valence correlation, spin-orbit coupling, scalar relativity, and the Lamb shift. In a few cases coupled cluster calculations with iterative triple (CCSDT) and quadruple (CCSDTQ) excitations have been carried out to estimate the effects of higher order electron correlation. The pseudopotential calculations have also been compared to all electron calculations using second- and third-order Douglas-Kroll-Hess Hamiltonians. In addition to the reaction enthalpies, heats of formation, bond lengths, and harmonic vibrational frequencies have been calculated for the stable triatomic products HgI2, HgIBr, HgICl, and HgIO. Accurate dissociation energies, equilibrium bond lengths, and harmonic vibrational frequencies have also been calculated for each of the diatomic molecules involved in this study (HgI, HgBr, HgCl, HgO, I2, IBr, ICl, and IO). The reported enthalpies are expected to have accuracies of 1 kcal/mol or better.

Journal Article↗

PubMed2005

Accurate global potential energy surface and reaction dynamics for the ground state of HgBr2.

A global potential energy surface (PES) for the (1)A' ground state of HgBr(2) has been constructed in order to determine the rate constants for atmospherically important reactions involving mercury and bromine. The total energy of HgBr(2) was calculated by the multireference configuration interaction level of theory with series of correlation consistent basis sets up to quadruple-zeta quality with subsequent extrapolation to the complete basis set limit. An additive correction for spin-orbit coupling was also included. The global PES was represented piecewise by interpolating three separate parts of the surface with the reproducing kernel Hilbert space method and connecting them smoothly by switch functions. Quasiclassical trajectory calculations carried out on the surface yielded 298 K thermal rate constants of 3.89 x 10(-11) cm(3)/(mol.s) for the abstraction reaction HgBr + Br --> Hg + Br(2), 2.98 x 10(-11) cm(3)/(mol.s) for the recombination reaction Br + HgBr --> HgBr(2), and 3.97 x 10(-11) cm(3)/(mol.s) for the exchange reaction Br + HgBr --> BrHg + Br. The insertion reaction Hg + Br(2) --> HgBr(2) was found to have a high barrier of 27.2 kcal/mol and a very small rate constant of just 2.74 x 10(-31) cm(3)/(mol.s) determined by the microcanonical variational transition state theory method. The implications of the obtained results to the description of the mechanism of recently observed polar tropospheric mercury depletion events are briefly discussed.

Journal Article↗

PubMed2005

A hierarchical family of three-dimensional potential energy surfaces for He-CO.

A hierarchical family of five three-dimensional potential energy surfaces has been developed for the benchmark He-CO system. Four surfaces were obtained at the coupled cluster singles and doubles level of theory with a perturbational estimate of triple excitations, CCSD(T), and range in quality from the doubly augmented double-zeta basis set to the complete basis set (CBS) limit. The fifth corresponds to an approximate CCSDT/CBS surface (CCSD with iterative triples/CBS, denoted CBS+corr). The CBS limit results were obtained by pointwise basis set extrapolations of the individual counterpoise-corrected interaction energies. For each surface, over 1000 interaction energies were accurately interpolated using a reproducing kernel Hilbert space approach with an R-6+R-7 asymptotic form. In each case, both three-dimensional and effective two-dimensional surfaces were developed. In standard Jacobi coordinates, the final CBS+corr surface has a global minimum at rCO=2.1322a0,R=6.418a0, and gamma=70.84 degrees with a well depth of -22.34 cm-1. The other four surfaces have well depths ranging from -14.83 cm-1 [CCSD(T)/d-aug-cc-pVDZ] to -22.02 cm-1 [CCSD(T)/CBS]. For each of these surfaces the infrared spectrum has been accurately calculated and compared to experiment, as well as to previous theoretical and empirical surfaces. The final CBS+corr surface exhibits root-mean-square and maximum errors compared to experiment (4He) of just 0.03 and 0.04 cm-1, respectively, for all 42 transitions and is the most accurate ab initio surface to date for this system. Other quantities investigated include the interaction second virial coefficient, the integral cross sections, and thermal rate coefficients for rotational relaxation of CO by He, and rate coefficients for CO vibrational relaxation by He. All the observable quantities showed a smooth convergence with respect to the quality of the underlying interaction surface.

Journal Article↗

PubMed2005

Systematically convergent basis sets for transition metals. I. All-electron correlation consistent basis sets for the 3d elements Sc-Zn.

Sequences of basis sets that systematically converge towards the complete basis set (CBS) limit have been developed for the first-row transition metal elements Sc-Zn. Two families of basis sets, nonrelativistic and Douglas-Kroll-Hess (-DK) relativistic, are presented that range in quality from triple-zeta to quintuple-zeta. Separate sets are developed for the description of valence (3d4s) electron correlation (cc-pVnZ and cc-pVnZ-DK; n = T,Q, 5) and valence plus outer-core (3s3p3d4s) correlation (cc-pwCVnZ and cc-pwCVnZ-DK; n = T,Q, 5), as well as these sets augmented by additional diffuse functions for the description of negative ions and weak interactions (aug-cc-pVnZ and aug-cc-pVnZ-DK). Extensive benchmark calculations at the coupled cluster level of theory are presented for atomic excitation energies, ionization potentials, and electron affinities, as well as molecular calculations on selected hydrides (TiH, MnH, CuH) and other diatomics (TiF, Cu2). In addition to observing systematic convergence towards the CBS limits, both 3s3p electron correlation and scalar relativity are calculated to strongly impact many of the atomic and molecular properties investigated for these first-row transition metal species.

Journal Article↗

PubMed2005

Heats of formation of xenon fluorides and the fluxionality of XeF(6) from high level electronic structure calculations.

Atomization energies at 0 K and heats of formation at 0 and 298 K are predicted for XeF(+), XeF(-), XeF(2), XeF(4), XeF(5)(-), and XeF(6) from coupled cluster theory (CCSD(T)) calculations with new correlation-consistent basis sets for Xe. To achieve near chemical accuracy (+/-1 kcal/mol), up to four corrections were added to the complete basis set binding energies based on frozen core coupled cluster theory energies: a correction for core-valence effects, a correction for scalar relativistic effects, a correction for first-order atomic spin-orbit effects, and in some cases, a second-order spin-orbit correction. Vibrational zero-point energies were computed at the coupled cluster level of theory. The structure of XeF(6) is difficult to obtain with the C(3)(v)() and O(h)() structures having essentially the same energy. The O(h)() structure is only 0.19 kcal/mol below the C(3)(v)() one at the CCSD(T)/CBS level using an approximate geometry for the C(3)(v)() structure. With an optimized C(3)(v)() geometry, the C(3)(v)() structure would probably become slightly lower in energy than the O(h)() one. The calculated heats of formation for the neutral XeF(n)() fluorides are less negative than the experimental values from the equilibrium measurements by 2.0, 7.7, and 12.2 kcal/mol for n = 2, 4, and 6, respectively. For the experimental values, derived from the photoionization measurements, this discrepancy becomes even larger, suggesting a need for a redetermination of the experimental values. Evidence is presented for the fluxionality of XeF(6) caused by the presence of a sterically active, free valence electron pair on Xe.

Journal Article↗

PubMed2005

Application of equation-of-motion coupled-cluster methods to low-lying singlet and triplet electronic states of HBO and BOH.

The equilibrium structures and physical properties of the X (1)sigma(+) linear electronic states, linear excited singlet and triplet electronic states of hydroboron monoxide (HBO) (A (1)sigma(-), B (1)delta, a (3)sigma(+), and b (3)delta) and boron hydroxide (BOH) (A (1)sigma(+), B (1)Pi, and b (3)Pi), and their bent counterparts (HBO a (3)A('), b (3)A("), A (1)A("), B (1)A(') and BOH X (1)A('), b (3)A('), c (3)A("), A (1)A('), B (1)A('), C (1)A(")) are investigated using excited electronic state ab initio equation-of-motion coupled-cluster (EOM-CC) methods. A new implementation of open-shell EOM-CC including iterative partial triple excitations (EOM-CC3) was tested. Coupled-cluster wave functions with single and double excitations (CCSD), single, double, and iterative partial triple excitations (CC3), and single, double, and full triple excitations (CCSDT) are employed with the correlation-consistent quadruple and quintuple zeta basis sets. The linear HBO X (1)sigma(+) state is predicted to lie 48.3 kcal mol(-1) (2.09 eV) lower in energy than the BOH X (1)sigma(+) linear stationary point at the CCSDT level of theory. The CCSDT BOH barrier to linearity is predicted to lie 3.7 kcal mol(-1) (0.16 eV). With a harmonic zero-point vibrational energy correction, the HBO X (1)sigma(+)-BOH X (1)A(') energy difference is 45.2 kcal mol(-1) (1.96 eV). The lowest triplet excited electronic state of HBO, a (3)A('), has a predicted excitation energy (T(e)) of 115 kcal mol(-1) (4.97 eV) from the HBO ground state minimum, while the lowest-bound BOH excited electronic state, b (3)A('), has a T(e) of 70.2 kcal mol(-1) (3.04 eV) with respect to BOH X (1)A('). The T(e) values predicted for the lowest singlet excited states are A (1)A(")<--X (1)sigma(+)=139 kcal mol(-1) (6.01 eV) for HBO and A (1)A(')<--X (1)A(')=102 kcal mol(-1) (4.42 eV) for BOH. Also for BOH, the triplet vertical transition energies are b (3)A(')<--X (1)A(')=71.4 kcal mol(-1) (3.10 eV) and c (3)A(")<--X (1)A(')=87.2 kcal mol(-1) (3.78 eV).

Journal Article↗

PubMed2005

Benchmark calculations of the complete configuration-interaction limit of Born-Oppenheimer diagonal corrections to the saddle points of isotopomers of the H + H2 reaction.

We present a detailed ab initio study of the effect that the Born-Oppenheimer diagonal correction (BODC) has on the saddle-point properties of the H3 system and its isotopomers. Benchmark values are presented that are estimated to be within 0.1 cm(-1) of the complete configuration-interaction limit. We consider the basis set and correlation treatment requirements for accurate BODC calculations, and both are observed to be more favorable than for the Born-Oppenheimer energies. The BODC raises the H + H2 barrier height by 0.1532 kcal/mol and slightly narrows the barrier--with the imaginary frequency increasing by approximately 2%.

Journal Article↗

PubMed2005

Methyl cation affinities of rare gases and nitrogen and the heat of formation of diazomethane.

The methyl cation affinities of the rare gases have been calculated at 0 and 298 K by using coupled cluster theory including noniterative, quasiperturbative triple excitations with the new correlation-consistent basis sets for Xe up through aug-cc-pV5Z in some cases. To achieve near chemical accuracy (+/-1 kcal/mol) in the thermodynamic properties, we add to the estimated complete basis set valence binding energies, based on frozen core coupled cluster theory energies, two corrections: (1) a core/valence correction and (2) a scalar relativistic correction. Vibrational zero-point energies were computed at the coupled cluster level of theory at the CCSD(T)/aug-cc-pVDZ level. The calculated rare gas methyl cation affinities (MCA in kcal/mol) at 298 K are the following: MCA(He) = 1.7, MCA(Ne) = 2.5, MCA(Ar) = 16.9, MCA(Kr) = 25.5, and MCA(Xe) = 36.6. Because of the importance of the MCA(N(2)) in the experimental measurements of the MCA scale, we calculated a number of quantities associated with CH(3)N(2)(+) and CH(2)N(2). The calculated values for diazomethane at 298 K are: DeltaH(f)(CH(2)N(2)) = 65.3 kcal/mol, PA(CH(2)N(2)) = 211.9 kcal/mol, and MCA(N(2)) = 43.2 kcal/mol.

Journal Article↗

PubMed2005

An ab initio study of the lowest electronic states of yttrium dicarbide, YC2.

The low-lying electronic states of yttrium dicarbide have been calculated using highly correlated wave functions and systematic sequences of correlation consistent basis sets. For the (2)A(1) ground electronic state, the near-equilibrium potential energy surface (PES) has been calculated using the coupled cluster method in conjunction with basis sets ranging in size from double to quintuple zeta. The relativistic effects have been taken into account by using pseudopotentials for the Y atom. After extrapolation to the complete basis set limit, additional corrections due to core-valence correlation and spin-orbit effects have also been included. The same approach has been followed for the (2)B(1), (2)B(2), and (2)A(2) states but only the C(2V) PESs have been considered in these cases. For the two (2)A(1) electronic excited states and, for comparison purposes, for the ground state, the multireference configuration interaction (MRCI) approach has been used in conjunction with double-zeta and triple-zeta basis sets for the construction of the PES. The molecular and spectroscopic properties predicted for the ground and excited states investigated in this work compare well with the available experimental data, particularly for the ground electronic state. The 0 K dissociation enthalpy of YC(2), DeltaH(Y-C2)(0 K), and its atomization enthalpy, SigmaD(0), are predicted to be 148.4 and 291.5 kcal/mol, respectively.

Journal Article↗

PubMed2005

Effect of the Ar-Ni(s) potential on the cross section for Ar+CH4/Ni{111} collision-induced desorption and the need for a more accurate CH4/Ni{111} potential.

In a previous paper [L. Sun, P. de Sainte Claire, O. Meroueh, and W. L Hase, J. Chem. Phys. 114, 535 (2001)], a classical trajectory simulation was reported of CH(4) desorption from Ni{111} by Ar-atom collisions. At an incident angle theta(i) of 60 degrees (with respect to the surface normal), the calculated collision-induced desorption (CID) cross sections are in excellent agreement with experiment. However, for smaller incident angles the calculated cross sections are larger than the experimental values and for normal collisions, theta(i)=0 degrees , the calculated cross sections are approximately a factor of 2 larger. This trajectory study used an analytic function for the Ar+Ni(s) intermolecular potential which gives an Ar-Ni{111} potential energy minimum which is an order of magnitude too deep. In the work reported here, the previous trajectory study is repeated with an Ar+Ni(s) analytic intermolecular potential which gives an accurate Ar-Ni{111} potential energy minimum and also has a different surface corrugation than the previous potential. Though there are significant differences between the two Ar+Ni(s) analytic potentials, they have no important effects on the CID dynamics and the cross sections reported here are nearly identical to the previous values. Zero-point energy motions of the surface and the CH(4)-Ni(s) intermolecular modes are considered in the simulation and they are found to have a negligible effect on the CID cross sections. Calculations of the intermolecular potential between CH(4) and a Ni atom, at various levels of theory, suggest that there are substantial approximations in the ab initio calculation used to develop the CH(4)+Ni{111} potential. The implication is that the differences between the trajectory and experimental CID cross sections may arise from an inaccurate CH(4)+Ni{111} potential used in the trajectory simulation.

Journal Article↗