Photo-induced desorption of NO from NiO(100): calculation of the four-dimensional potential energy surfaces and systematic wave packet studies

The velocity distributions of the laser-induced desorption of NO molecules from an epitaxially grown film of NiO(100) on Ni(100) have been studied [Mull et al., J. Chem. Phys., 1992, 96, 7108]. A pronounced bimodality of velocity distributions has been found, where the NO molecules desorbing with higher velocities exhibit a coupling to the rotational quantum states J. In this article we present simulations of state resolved velocity distributions on a full ab initio level. As a basis for this quantum mechanical treatment a 4D potential energy surface (PES) was constructed for the electronic ground and a representative excited state, using a NiO5Mg(18+)13 cluster. The PESs of the electronic ground and an excited state were calculated at the CASPT2 and the configuration interaction (CI) level of theory, respectively. Multi-dimensional quantum wave packet simulations on these two surfaces were performed for different sets of degrees of freedom. Our key finding is that at least a 3D wave packet simulation, in which the desorption coordinate Z, polar angle theta and lateral coordinate X are included, is necessary to allow the simulation of experimental velocity distributions. Analysis of the wave packet dynamics demonstrates that essentially the lateral coordinate, which was neglected in previous studies [Klüner et al., Phys. Rev. Lett. 1998, 80, 5208], is responsible for the experimentally observed bimodality. An extensive analysis shows that the bimodality is due to a bifurcation of the wave packet on the excited state PES, where the motion of the molecule parallel to the surface plays a decisive role.     

An eight-degree-of-freedom, time-dependent quantum dynamics study for the H2+C2H reaction on a new modified potential energy surface

An eight-dimensional time-dependent quantum dynamics wave packet approach is performed for the study of the H2+C2H-->H+C2H2 reaction system on a new modified potential energy surface (PES) [L.-P. Ju et al., Chem. Phys. Lett. 409, 249 (2005)]. This new potential energy surface is obtained by modifying Wang and Bowman's old PES [J. Chem. Phys. 101, 8646 (1994)] based on the new ab initio calculation. This new modified PES has a much lower transition state barrier height at 2.29 kcal/mol than Wang and Bowman's old PES at 4.3 kcal/mol. This study shows that the reactivity for this diatom-triatom reaction system is enhanced by vibrational excitations of H2, whereas the vibrational excitations of C2H only have a small effect on the reactivity. Furthermore, the bending excitations of C2H, compared to the ground state reaction probability, hinder the reactivity. The comparison of the rate constant between this calculation and experimental results agrees with each other very well. This comparison indicates that the new modified PES corrects the large barrier height problem in Wang and Bowman's old PES.     

On the role of dynamical barriers in barrierless reactions at low energies: S(1D) + H2

Reaction probabilities as a function of total angular momentum (opacity functions) and the resulting reaction cross sections for the collision of open shell S((1)D) atoms with para-hydrogen have been calculated in the kinetic energy range 0.09-10 meV (1-120 K). The quantum mechanical hyperspherical reactive scattering method and quasi-classical trajectory and statistical quasi-classical trajectory approaches were used. Two different ab initio potential energy surfaces (PESs) have been considered. The widely used reproducing kernel Hilbert space (RKHS) PES by Ho et al. [T.-S. Ho, T. Hollebeek, H. Rabitz, S. D. Chao, R. T. Skodje, A. S. Zyubin, and A. M. Mebel, J. Chem. Phys 116, 4124 (2002)] and the recently published accurate double many-body expansion (DMBE)/complete basis set (CBS) PES by Song and Varandas [Y. Z. Song and A. J. C. Varandas, J. Chem. Phys. 130, 134317 (2009)]. The calculations at low collision energies reveal very different dynamical behaviors on the two PESs. The reactivity on the RKHS PES is found to be considerably larger than that on the DMBE/CBS PES as a result of larger reaction probabilities at low total (here also orbital) angular momentum values and to opacity functions which extend to significantly larger total angular momentum values. The observed differences have their origin in two major distinct topographic features. Although both PESs are essentially barrierless for equilibrium H-H distances, when the H-H bond is compressed the DMBE/CBS PES gives rise to a dynamical barrier which limits the reactivity of the system. This barrier is completely absent in the RHKS PES. In addition, the latter PES exhibits a van der Walls well in the entrance channel which reduces the height of the centrifugal barrier and is able to support resonances. As a result, a significant larger cross section is found on this PES, with marked oscillations attributable to shape resonances and/or to the opening of partial wave contributions. The comparison of the results on both PESs is illustrative of the wealth of the dynamics at low collision energy. It is also illuminating about the difficulties encountered in modeling an all-purpose global potential energy surface.     

Quasiclassical trajectory study of the Cl+CH4 reaction dynamics on a quadratic configuration interaction with single and double excitation interpolated potential energy surface

An ab initio interpolated potential energy surface (PES) for the Cl+CH(4) reactive system has been constructed using the interpolation method of Collins and co-workers [J. Chem. Phys. 102, 5647 (1995); 108, 8302 (1998); 111, 816 (1999); Theor. Chem. Acc. 108, 313 (2002)]. The ab initio calculations have been performed using quadratic configuration interaction with single and double excitation theory to build the PES. A simple scaling all correlation technique has been used to obtain a PES which yields a barrier height and reaction energy in good agreement with high level ab initio calculations and experimental measurements. Using these interpolated PESs, a detailed quasiclassical trajectory study of integral and differential cross sections, product rovibrational populations, and internal energy distributions has been carried out for the Cl+CH(4) and Cl+CD(4) reactions, and the theoretical results have been compared with the available experimental data. It has been shown that the calculated total reaction cross sections versus collision energy for the Cl+CH(4) and Cl+CD(4) reactions is very sensitive to the barrier height. Besides, due to the zero-point energy (ZPE) leakage of the CH(4) molecule to the reaction coordinate in the quasiclassical trajectory (QCT) calculations, the reaction threshold falls below the barrier height of the PES. The ZPE leakage leads to CH(3) and HCl coproducts with internal energy below its corresponding ZPEs. We have shown that a Gaussian binning (GB) analysis of the trajectories yields excitation functions in somehow better agreement with the experimental determinations. The HCl(v'=0) and DCl(v'=0) rotational distributions are as well very sensitive to the ZPE problem. The GB correction narrows and shifts the rotational distributions to lower values of the rotational quantum numbers. However, the present QCT rotational distributions are still hotter than the experimental distributions. In both reactions the angular distributions shift from backward peaked to sideways peaked as collision energy increases, as seen in the experiments and other theoretical calculations.     

A global 12-dimensional ab initio potential energy surface and dynamical studies for the SiH4+H-->SiH3+H2 reaction

A global 12-dimensional ab initio interpolated potential energy surface (PES) for the SiH(4)+H-->SiH(3)+H(2) reaction is presented. The ab initio calculations are based on the unrestricted quadratic configuration interaction treatment with all single and double excitations together with the cc-pVTZ basis set, and the modified Shepard interpolation method of Collins and co-workers [K. C. Thompson et al., J. Chem. Phys. 108, 8302 (1998); M. A. Collins, Theor. Chem. Acc. 108, 313 (2002); R. P. A. Bettens and M. A. Collins, J. Chem. Phys. 111, 816 (1999)] is applied. Using this PES, classical trajectory and variational transition state theory calculations have been carried out, and the computed rate constants are in good agreement with the available experimental data.     

Quantum dynamics of H2, D2, and HD in the small dodecahedral cage of clathrate hydrate: evaluating H2-water nanocage interaction potentials by comparison of theory with inelastic neutron scattering experiments

We have performed rigorous quantum five-dimensional (5D) calculations and analysis of the translation-rotation (T-R) energy levels of one H(2), D(2), and HD molecule inside the small dodecahedral (H(2)O)(20) cage of the structure II clathrate hydrate, which was treated as rigid. The H(2)- cage intermolecular potential energy surface (PES) used previously in the molecular dynamics simulations of the hydrogen hydrates [Alavi et al., J. Chem. Phys. 123, 024507 (2005)] was employed. This PES, denoted here as SPC/E, combines an effective, empirical water-water pair potential [Berendsen et al., J. Phys. Chem. 91, 6269 (1987)] and electrostatic interactions between the partial charges placed on H(2)O and H(2). The 5D T-R eigenstates of HD were calculated also on another 5D H(2)-cage PES denoted PA-D, used by us earlier to investigate the quantum T-R dynamics of H(2) and D(2) in the small cage [Xu et al., J. Phys. Chem. B 110, 24806 (2006)]. In the PA-D PES, the hydrogen-water pair potential is described by the ab initio 5D PES of the isolated H(2)-H(2)O dimer. The quality of the SPC/E and the PA-D H(2)-cage PESs was tested by direct comparison of the T-R excitation energies calculated on them to the results of two recent inelastic neutron scattering (INS) studies of H(2) and HD inside the small clathrate cage. The translational fundamental and overtone excitations, as well as the triplet splittings of the j=0-->j=1 rotational transitions, of H(2) and HD in the small cage calculated on the SPC/E PES agree very well with the INS results and represent a significant improvement over the results computed on the PA-D PES. Our calculations on the SPC/E PES also make predictions about several spectroscopic observables for the encapsulated H(2), D(2), and HD, which have not been measured yet.     

Quasi-classical trajectories study of Ne79Br2(B) vibrational predissociation

A full-dimensional quasi-classical trajectories study on the vibrational predissociation (VP) of the Ne79Br2(B) complex is presented. Following the most recent experiments, the Br2(B) vibrational levels v'=16-29 were explored. The total angular momentum, J, was taken to be zero, and a semiclassical Franck-Condon model to compute initial conditions from quantum distributions was employed. Predissociation lifetimes were extracted from Ne79Br2 population decay by using two different exponential laws. Predicted lifetimes are in excellent agreement with the last experimental results [J. A. Cabrera, C. R. Bieler, B. C. Olbricht, W. E. van der Veer and K. C. Janda, J. Chem. Phys., 2005, 123, 054311]. The Br2 fragment ro-vibrational distributions resulting from the VP of the molecule were obtained from the statistics of classical magnitudes using the standard binning procedure. Computed rotational distributions (for the Deltav'=-1, -2 channels) are also in very good agreement with the experimental results [M. Nejad-Sattari and T. A. Stephenson, J. Chem. Phys., 1997, 106 5454]. The influence of two quantum effects-the closing of the Deltav'=-1 dissociation channel and the intramolecular vibrational relaxation (IVR) mechanism-on the agreement with experimental rotational distributions, is discussed. Due to the classical character of our calculations and the binning procedure we used, the agreement of computed vibrational distributions with experimental and quantum theoretical is qualitative. For instance, for v'=28-for which the Deltav'=-1 channel is experimentally found to be closed-the Deltav'=-2 channel becomes statistically more significant. A discussion on the viability of similar quasi-classical methods to model the VP dynamics of analogous clusters is presented.     

A quantum-classical treatment of non-adiabatic transitions

We have carried out molecular dynamics simulations of non-adiabatic processes with the help of a newly formulated potentially exact quantum-classical approach derived from a method proposed earlier [J. Chem. Phys. 118 (2003) 5302]. In this method, time-dependent Schroedinger equation is solved by representing Ψ on a moving Gauss–Hermite DVR grid, the motion of grid-centre being handled classically, but self consistently with the quantum evolution of the wavefunction. Electronic transitions are allowed anywhere in the configuration space among any number of coupled states. We have tested the method on three model problems proposed by J.C. Tully [J. Chem. Phys. 93 (1990) 1061]. These models are relevant to a wide range of gas-phase and condensed-phase phenomena occurring even at low energies. Excellent agreement of computed transition probabilities with corresponding quantum mechanical (DVR/FFT) results even in the deep quantum regime and its cost-efficiency (computational) are encouraging.     

Adsorption and scattering of H2 and D2 by NiAl(110)

We present quasiclassical dynamics calculations of H2 and D2 scattering by the NiAl(110) surface using a recently proposed six-dimensional potential-energy surface (PES) obtained from density-functional theory calculations. The results for dissociative adsorption confirm several experimental predictions using (rotationally hot) D2 beams, namely, the existence of a dissociation barrier, the small isotopic effect, the importance of vibrational enhancement, and the existence of normal energy scaling. The latter conclusion shows that normal energy scaling is not necessarily associated with weak corrugated surfaces. The results for rotationally elastic and inelastic diffractions are also in reasonable agreement with experiment, but they show that many more diffractive transitions are responsible for the observed structures than previously assumed. This points to the validity of the PES recently proposed [P. Riviere, H. F. Busnengo, and F. Martin, J. Chem. Phys. 121, 751 (2004)] to describe dissociative adsorption as well as rotationally elastic and inelastic diffractions in the H2NiAl(110) system.     

Polarizable and flexible model for ethanol

A polarizable, flexible model for ethanol is obtained based on an extensive series of B3LYP/6-311++G(d,p) calculations and molecular dynamics simulations. The ethanol model includes electric-field dependence in both the atomic charges and the intramolecular degrees of freedom. Field-dependent intramolecular potentials have been attempted only once previously, for OH and HH stretches in water [P. Cicu et al., J. Chem. Phys. 112, 8267 (2000)]. The torsional potential involving the hydrogen-bonding hydrogen in ethanol is found to be particularly field sensitive. The methodology for developing field-dependent potentials can be readily generalized to other molecules and is discussed in detail. Molecular dynamics simulations of bulk ethanol are performed and the results are assessed based on comparisons with the self-diffusion coefficient [N. Karger et al., J. Chem. Phys. 93, 3437 (1990)], dielectric constant [J. T. Kindt and C. A. Schmuttenmaer, J. Phys. Chem. 100, 10373 (1996)], enthalpy of vaporization [R. C. Wilhoit and B. J. Zwolinski, J. Phys. Chem. Ref. Data, Suppl. 2, 2 (1973)], and experimental interatomic distributions [C. J. Benmore and Y. L. Loh, J. Chem. Phys. 112, 5877 (2000)]. The simultaneous variation of the atomic charges and the intramolecular potentials requires modified equations of motion and a multiple time step algorithm has been implemented to solve these equations. The article concludes with a discussion of the bulk structure and properties with an emphasis on the hydrogen bonding network.     

Photodissociation of CH3I: a full-dimensional (9D) quantum dynamics study

The photodissociation of methyl iodide in the A band is studied by full-dimensional (9D) wave packet dynamics calculations using the multiconfigurational time-dependent Hartree approach. The potential energy surfaces employed are based on the diabatic potentials of Xie et al. [J. Phys. Chem. A 2000, 104, 1009] and the vertical excitation energy is taken from recent ab initio calculations [Alekseyev et al. J. Chem. Phys.2007, 126, 234102]. The absorption spectrum calculated for exclusively parallel excitation agrees well with the experimental spectrum of the A band. The electronic population dynamics is found to be strongly dependent on the motion in the torsional coordinate related to the H(3)-C-I bend, which presumably is an artifact of the diabatic model employed. The calculated fully product state-selected partial spectra can be interpreted based on the reflection principle and suggests strong coupling between the C-I stretching and the H(3)-C-I bending motions during the dissociation process. The computed rotational and vibrational product distributions typically reproduce the trends seen in the experiment. In agreement with experiment, a small but significant excitation of the total symmetric stretching and the asymmetric bending modes of the methyl fragment can be seen. In contrast, the umbrella mode of the methyl is found to be too highly excited in the calculated distributions.     

Effect of the overall rotation on the cis-trans isomerization of HONO induced by an external field

Rovibrational eigenenergies of HONO are computed and compared to experimental energies available in the literature. For their computation, we use a previously developed potential energy surface (PES) and a newly derived exact kinetic energy operator (KEO) including the overall rotation for a tetra-atomic molecule in non-orthogonal coordinates. In addition, we use the Heidelberg Multi-Configuration Time-Dependent Hartree (MCTDH) package. We compare the experimental rovibrational eigenvalues of HONO available in the literature with those obtained with MCTDH and a previously developed potential energy surface (PES) [F. Richter et al., J. Chem. Phys., 2004, 120, 1306.] for the cis geometry. The effect of the overall rotation on the process studied in our previous work on HONO [F. Richter et al., J. Chem. Phys., 2007, 127, 164315.] leading to the cis→trans isomerization of HONO is investigated. This effect on this process is found to be weak.     

Investigation of entrance and exit effects on liquid transport through a cylindrical nanopore

The entrance and exit effects on liquid transport through a nano-sized cylindrical pore under different solid wall-liquid interactions were studied by comparing molecular dynamics (MD) results of a finite length nanopore in a membrane with those of an infinite length one. The liquid transport through a finite length nanopore in a membrane was carried out by using a pressure-driven non-equilibrium molecular dynamics (NEMD) method proposed by Huang et al. [C. Huang, K. Nandakumar, P. Choi and L. W. Kostiuk, J. Chem. Phys., 2006, 124, 234701]. The fluid motion through an infinite length nanopore, which had the same cross-stream dimension as the finite length channel in the membrane, but with periodic boundary conditions in the stream-wise direction, was carried out by using the external-field driven NEMD approach [J. Koplik, J. R. Bavanar and J. F. Willemsen, Phys. Rev. Lett., 1988, 60, 1282]. The NEMD results show that the pressure and density distributions averaged over the channel in the radial direction in both finite and infinite length channels are similar, but the radial distributions of the stream-wise velocity were significantly different when the solid wall was repulsive. The entrance and exit effects lead to a decrease in flow rate at about 39% for the repulsive wall and 6% for the neutral-like wall.     

Non-reactive scattering of N2 from the W(110) surface studied with different exchange-correlation functionals

The non-reactive scattering of N(2) from the W(110) surface is studied with six dimensional (6D) classical dynamics and two distinct potential energy surfaces (PES). Here, we use the PESs calculated with density functional theory and two different exchange-correlation functionals, the PW91 [J. E. Perdew et al., Phys. Rev. B, 1992, 46, 6671] and the RPBE [B. Hammer et al., Phys. Rev. B, 1999, 59, 7413]. By analyzing the final rotational state and angular distributions, we extract information on the characteristics of the two PESs in the 6D configurational space. Comparison of the theoretical results with the available experimental data provides detailed information on the validity of each functional. In general, the PW91 PES is more corrugated than the RPBE one in all the configurational space, meaning that there is a stronger dependence of the potential energy on the molecular orientation and position over the surface unit cell. Furthermore, we find that the larger corrugation and the less repulsive character exhibited by the PW91 PES seems to be realistic at distances above the chemisorption well. In contrast, the less corrugated RPBE PES performs better in the region below the chemisorption well.     

Ab initio large-amplitude quantum-tunneling dynamics in vinyl radical: a vibrationally adiabatic approach

Large-amplitude tunneling in vinyl radical over a C2v planar transition state involves CCH bending excitation coupled to all other internal coordinates, resulting in a significant dependence of barrier height and shape on vibrational degrees of freedom at the zero-point level. An ab initio potential surface for vinyl radical has been calculated at the CCSD(T) level (AVnZ; n=2, 3, 4, 5) for vibrationally adiabatic 1D motion along the planar CCH bending tunneling coordinate, extrapolated to the complete basis set (CBS) limit and corrected for anharmonic zero-point effects. The polyatomic reduced moment of inertia is calculated explicitly as a function of tunneling coordinate, with eigenvalues and tunneling splittings obtained from numerical solution of the resulting 1D Schr?dinger equation. Linear scaling of the CBS potential to match predicted and observed tunneling splittings empirically yields an adiabatic barrier height of DeltaEadiab=1696(20) cm(-1) which, when corrected for zero-point energy contributions, translates into an effective barrier of DeltaEeff=1602(20) cm(-1) consistent with estimates (DeltaE=1580(100) cm(-1)) by Tanaka and coworkers [J. Chem. Phys., 2004, 120, 3604-3618]. These zero-point-corrected potential surfaces are used to predict tunneling dynamics in vibrationally excited states of vinyl radical, providing strong support for previous jet-cooled high-resolution infrared studies [Dong et al., J. Phys. Chem. A, 2006, 110, 3059-3070] in the symmetric CH2 stretch mode.     

Hydrogen adsorption strength and sites in the metal organic framework MOF5: Comparing experiment and model calculations

Hydrogen adsorption in porous, high surface area, and stable metal organic frameworks (MOF’s) appears a novel route towards hydrogen storage materials [N.L. Rosi, J. Eckert, M. Eddaoudi, D.T. Vodak, J. Kim, M. O’Keeffe, O.M. Yaghi, Science 300 (2003) 1127; J.L.C. Rowsell, A.R. Millward, K. Sung Park, O.M. Yaghi, J. Am. Chem. Soc. 126 (2004) 5666; G. Ferey, M. Latroche, C. Serre, F. Millange, T. Loiseau, A. Percheron-Guegan, Chem. Commun. (2003) 2976; T. Loiseau, C. Serre, C. Huguenard, G. Fink, F. Taulelle, M. Henry, T. Bataille, G. Férey, Chem. Eur. J. 10 (2004) 1373]. A prerequisite for such materials is sufficient adsorption interaction strength for hydrogen adsorbed on the adsorption sites of the material because this facilitates successful operation under moderate temperature and pressure conditions. Here we report detailed information on the geometry of the hydrogen adsorption sites, based on the analysis of inelastic neutron spectroscopy (INS). The adsorption energies for the metal organic framework MOF5 equal about 800 K for part of the different sites, which is significantly higher than for nanoporous carbon materials (550 K) [H.G. Schimmel, G.J. Kearley, M.G. Nijkamp, C.T. Visser, K.P. de Jong, F.M. Mulder, Chem. Eur. J. 9 (2003) 4764], and is in agreement with what is found in first principles calculations [T. Sagara, J. Klassen, E. Ganz, J. Chem. Phys. 121 (2004) 12543; F.M. Mulder, T.J. Dingemans, M. Wagemaker, G.J. Kearley, Chem. Phys. 317 (2005) 113]. Assignments of the INS spectra is realized using comparison with independently published model calculations [F.M. Mulder, T.J. Dingemans, M. Wagemaker, G.J. Kearley, Chem. Phys. 317 (2005) 113] and structural data [T. Yildirim, M.R. Hartman, Phys. Rev. Lett. 95 (2005) 215504].     

Comparison of model potentials for molecular-dynamics simulations of silica

Structural, thermomechanical, and dynamic properties of pure silica SiO2 are calculated with three different model potentials, namely, the potential suggested by van Beest, Kramer, and van Santen (BKS) [Phys. Rev. Lett. 64, 1955 (1990)], the fluctuating-charge potential with a Morse stretch term for the short-range interactions proposed by Demiralp, Cagin, and Goddard (DCG)[Phys. Rev. Lett. 82, 1708 (1999)], and a polarizable force field proposed by Tangney and Scandolo (TS) [J. Chem. Phys. 117, 8898 (2002)]. The DCG potential had to be modified due to flaws in the original treatment. While BKS reproduces many thermomechanical properties of different polymorphs rather accurately, it also shows qualitatively wrong trends concerning the phononic density of states, an absence of the experimentally observed anomaly in the c/a ratio at the quartz alpha-beta transition, pathological instabilities in the beta-cristobalite phase, and a vastly overestimated transition pressure for the stishovite I --> II transition. These shortcomings are only partially remedied by the modified DCG potential but greatly improved by the TS potential. DCG and TS both reproduce a pressure-induced transition from alpha-quartz to quartz II, predicted theoretically based on the BKS potential.     

Theoretical and experimental study of weakly bound CO2-(pH2)2 trimers

The infrared spectrum of CO(2)-(pH(2))(2) trimers is predicted by performing exact basis-set calculations on a global potential energy surface defined as the sum of accurately known two-body pH(2)-CO(2) (J. Chem. Phys. 2010, 132, 214309) and pH(2)-pH(2) potentials (J. Chem. Phys. 2008, 129, 094304). These results are compared with new spectroscopic measurements for this species, for which 13 transitions are now assigned. A reduced-dimension treatment of the pH(2) rotation has been employed by applying the hindered-rotor averaging technique of Li, Roy, and Le Roy (J. Chem. Phys. 2010, 133, 104305). Three-body effects and the quality of the potential are discussed. A new technique for displaying the three-dimensional pH(2) density in the body-fixed frame is used, and shows that in the ground state the two pH(2) molecules are localized much more closely together than is the case for the two He atoms in the analogous CO(2)-(He)(2) species. A clear tunneling splitting is evident for the torsional motion of the two pH(2) molecules on a ring about the CO(2) molecular axis, in contrast to the case of CO(2)-(He)(2) where a more regular progression of vibrational levels reflects the much lower torsional barrier.