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.     

Vibrationally inelastic collisions in H+ +CO system: comparing quantum calculations with experiments

State-resolved cross beam experiments [H. Udseth et al., J. Chem. Phys. 60, 3051 (1974); J. Krutein and F. Linder, J. Chem. Phys. 71, 599 (1979); G. Niedner-Schatteburg and J. P. Toennies, Adv. Chem. Phys. LXXXII, 553 (1992)], coupled with proton energy loss spectroscopy for the inelastic scattering of H(+) from CO in the collision range of 10-30 eV show very low vibrational excitation of the target molecule. Stimulated by the experimentally observed low vibrational inelasticity in the system the ground and the first two low-lying excited electronic potential-energy surfaces have been computed using the ab initio multireference configuration interaction method. Quantum dynamics has been performed on the ground potential energy surface in the framework of vibrational close-coupling rotational infinite-order sudden approximation. The various computed dynamical attributes such as differential and integral cross sections, and average vibrational energy transfer are analyzed in detail, and compared successfully with the available experimental results.     

Quantum control of molecular vibrational and rotational excitations in a homonuclear diatomic molecule: a full three-dimensional treatment with polarization forces

The optimal control of the vibrational excitation of the hydrogen molecule [Balint-Kurti et al., J. Chem. Phys. 122, 084110 (2005)] utilizing polarization forces is extended to three dimensions. The polarizability of the molecule, to first and higher orders, is accounted for using explicit ab initio calculations of the molecular electronic energy in the presence of an electric field. Optimal control theory is then used to design infrared laser pulses that selectively excite the molecule to preselected vibrational-rotational states. The amplitude of the electric field of the optimized pulses is restricted so that there is no significant ionization during the process, and a new frequency sifting method is used to simplify the frequency spectrum of the pulse. The frequency spectra of the optimized laser pulses for processes involving rotational excitation are more complex than those relating to processes involving only vibrational excitation.     

A pulsed-field ionization photoelectron secondary ion coincidence study of the H2+ (X,upsilon+=0-15,N+=1)+He proton transfer reaction

The endothermic proton transfer reaction, H2+(upsilon+)+He-->HeH+ + H(DeltaE=0.806 eV), is investigated over a broad range of reactant vibrational levels using high-resolution vacuum ultraviolet to prepare reactant ions either through excitation of autoionization resonances, or using the pulsed-field ionization-photoelectron-secondary ion coincidence (PFI-PESICO) approach. In the former case, the translational energy dependence of the integral reaction cross sections are measured for upsilon+=0-3 with high signal-to-noise using the guided-ion beam technique. PFI-PESICO cross sections are reported for upsilon+=1-15 and upsilon+=0-12 at center-of-mass collision energies of 0.6 and 3.1 eV, respectively. All ion reactant states selected by the PFI-PESICO scheme are in the N+=1 rotational level. The experimental cross sections are complemented with quasiclassical trajectory (QCT) calculations performed on the ab initio potential energy surface provided by Palmieri et al. [Mol. Phys. 98, 1839 (2000)]. The QCT cross sections are significantly lower than the experimental results near threshold, consistent with important contributions due to resonances observed in quantum scattering studies. At total energies above 2 eV, the QCT calculations are in excellent agreement with the present results. PFI-PESICO time-of-flight (TOF) measurements are also reported for upsilon+=3 and 4 at a collision energy of 0.6 eV. The velocity inverted TOF spectra are consistent with the prevalence of a spectator-stripping mechanism.     

Three-dimensional potential energy surface of the Ar-OH(2Pi i) complex

Pure rotational transitions in the ground state for Ar-OH and Ar-OD [Y. Ohshima et al., J. Chem. Phys. 95, 7001 (1991) and Y. Endo et al., Faraday Discuss. 97, 341 (1994)], those in the excited states of the OH vibration, nu(s)=1 and 2, observed by Fourier-transform microwave spectroscopy in the present study, rotation-vibration transitions observed by infrared-ultraviolet double-resonance spectroscopy [K. M. Beck et al., Chem. Phys. Lett. 162, 203 (1989) and R. T. Bonn et al., J. Chem. Phys. 112, 4942 (2000)], and the P-level structure observed by stimulated emission pumping spectroscopy [M. T. Berry et al., Chem. Phys. Lett. 178, 301 (1991)] have been simultaneously analyzed to determine the potential energy surface of Ar-OH in the ground state. A Schrodinger equation, considering all the freedom of motions for an atom-diatom system in the Jacobi coordinate, R, theta, and r, was numerically solved to obtain energies of the rovibrational energy levels using the discrete variable representation method. A three-dimensional potential energy surface is determined by a least-squares fitting. In the analysis the potential parameters, obtained by ab initio calculations at the RCCSD(T) level of theory with a set of basis functions of aug-cc-pVTZ and midbond functions, are used as initial values. The determined intermolecular potential energy surface and its dependence on the OH monomer bond length are compared with those of an isovalent radical complex, Ar-SH.     

Laser excitation spectrum and spectroscopic potential parameters of Cd2 molecule in the 1u(5(3)P2) energy state

Excitation spectrum of Cd2 van der Waals complex was observed in a supersonic free-jet expansion beam. Cadmium dimers were seeded in argon environment (a carrier gas) while a pulsed dye laser served as an excitation light source. Well resolved structures of vibrational bands arising from transition between 1u(3pi(u)) and X0g+ (1sigmag+) states were observed. Using a Birge-Sponer method for analyzing the vibrational excitation spectrum, the first-time determination of the spectroscopical potential parameters of the 1u molecular state was performed. We obtained D'e = 723 +/- 10 cm(-1), omega'e = 28.9 +/- 1.0 cm(-1), omega'e x'e = 0.260 +/- 0.002 cm(-1) and R'e = 3.93 +/- 0.05 A. The latter value was estimated with the help of the computer-simulation of the spectrum. The results are compared with recent results of ab initio calculation of Czuchaj et al. (Chem. Phys. Lett. 225 (1994) 233) and demonstrate a reasonable agreement.     

A dual-level state-specific time-dependent density-functional theory

A highly efficient new algorithm for time-dependent density-functional theory (TDDFT) calculations is presented. In this algorithm, a dual-level approach to speed up DFT calculations (Nakajima and Hirao, J Chem Phys 2006, 124, 184108) is combined with a state-specific (SS) algorithm for TDDFT (Chiba et al., Chem Phys Lett 2006, 420, 391). The dual-level SS-TDDFT algorithm was applied to excitation energy calculations of typical small molecules, the Q bands of the chlorophyll A molecule, the charge-transfer energy of the zincbacteriochlorin-bacteriochlorin model system, and the lowest-lying excitation of the circumcoronene molecule. As a result, it was found that the dual-level SS-TDDFT gave correct excitation energies with errors of 0.2-0.3 eV from the standard TDDFT approach, with much lower CPU times for various types of excitation energies of large-scale molecules.     

Three dimensional quantum dynamics of (H-, H2) and its isotopic variants

We present the results of a time-dependent quantum mechanical investigation using centrifugal sudden approximation in the form of reaction probability as a function of collision energy (E(trans)) in the range 0.3-3.0 eV for a range of total angular momentum (J) values and the excitation function sigma(E(trans)) for the exchange reaction H(-) + H(2) (v = 0, j = 0) --> H(2) + H(-) and its isotopic variants in three dimensions on an accurate ab initio potential energy surface published recently (J. Chem. Phys. 2004, 121, 9343). The excitation function results are shown to be in excellent agreement with those obtained from crossed beam measurements by Zimmer and Linder for H(-) + D(2) collisions for energies below the threshold for electron detachment channel and somewhat larger than the most recent results of Haufler et al. for (H(-), D(2)) and (D(-), H(2)) collisions.     

Pulse shaping effect on two-photon excitation efficiency of alpha-perylene crystals and perylene in chloroform solution

We demonstrated that the two-photon excitation efficiency of perylene in chloroform solution as well as that of crystalline perylene was dramatically increased by optimizing the shape of intense femtosecond laser pulses of a regenerative amplifier output. The efficiency was three times higher than for an unshaped single femtosecond pulse with the pulse width of shorter than 50 fs. The pulse shape optimized for the solution sample was a pulse train with a repetition frequency of about 340 cm(-1), and the pulse shape optimized for crystalline perylene was very similar. These results supported our previous findings on alpha-perylene crystals using weak femtosecond pulses from a mode-locked laser oscillator [T. Okada et al. J. Chem. Phys. 121, 6385 (2004)]. Furthermore, it was confirmed that the shaped pulse optimized for the liquid sample could also increase the two-photon excitation efficiency of alpha-perylene crystals and vice versa. We concluded that the mechanism for the increase in excitation efficiency of perylene in chloroform was almost the same as that for alpha-perylene crystal, and that the efficiency increased mainly through intramolecular dynamical processes. Processes involving intermolecular interactions and/or energy states delocalized over the crystal cannot play the major role.     

Quantum dynamics of H + LiH reaction and its isotopic variants

Time-dependent quantum wave packet dynamics study is carried out to investigate the initial state selected channel specific reactivity of H + LiH collisional system on a new and more accurate ab initio potential energy surface developed by Wernli et al. [J. Phys. Chem. A 113, 1121 (2009)]. The H + LiH reaction proceeds through LiH depletion and H-exchange paths. While the former path is highly exoergic (by ～2.258 eV), the latter path is thermoneutral. State selected and energy resolved integral reaction cross sections and thermal rate constants are reported and compared with the literature data. The reactivity of the LiH depletion channel is found to be greater than the H-exchange channel. Rotational excitation of the reagent LiH molecule causes a decrease of reactivity of both the channels. On the other hand, the vibrational excitation of the reagent LiH decreases the reactivity of the LiH depletion channel and increases the reactivity of the H-exchange channel. The effect of isotopic substitution (H by D) on the reaction dynamics is also examined.     

Accuracy of recent potential energy surfaces for the He-N2 interaction. I. Virial and bulk transport coefficients

A new exchange-Coulomb semiempirical model potential energy surface for the He-N2 interaction has been developed. Together with two recent high-level ab initio potential energy surfaces, it has been tested for the reliability of its predictions of second-virial coefficients and bulk transport phenomena in binary mixtures of He and N2. The agreement with the relevant available measurements is generally within experimental uncertainty for the exchange-Coulomb surface and the ab initio surface of Patel et al. [J. Chem. Phys. 119, 909 (2003)], but with slightly poorer agreement for the earlier ab initio surface of Hu and Thakkar [J. Chem. Phys. 104, 2541 (1996)].     

Constant pressure ab initio molecular dynamics with discrete variable representation basis sets

The use of discrete variable representation (DVR) basis sets within ab initio molecular dynamics calculations allows the latter to be performed with converged energies and, more importantly, converged forces. In this paper, we show how to carry out ab initio molecular dynamics calculations in the isothermal-isobaric ensemble with fully flexible simulation boxes within the DVR basis set framework. In particular, we derive the appropriate DVR based expression for the pressure tensor when the electronic structure is represented using Kohn-Sham density functional theory, and we examine the convergence of this expression as a function of the basis set size. An illustrative example using 64 silicon atoms in a fully flexible box using a combination of the Martyna-Tobias-Klein [Martyna et al., J. Chem. Phys. 101, 4177 (1994)] and Car-Parrinello [Car and Parinello, Phys. Rev. Lett. 55, 2471 (1985)] algorithms is presented to demonstrate the efficacy of the approach.     

Energy-flow dynamics in the molecular channel of propanal photodissociation, C2H5CHO --> C2H6 + CO: direct ab initio molecular dynamics study

Direct ab initio molecular dynamics calculations have been carried out for the molecular channel of the photodissociation of propanal, C2H5CHO --> C2H6 + CO, at the RMP2(full)/cc-pVDZ level of ab initio molecular orbital theory. The initial conditions were generated using the microcanonical sampling to put the excess energy randomly into all vibrational modes of the TS. Starting from the TS, a total of approximately 700 trajectories were numerically integrated for 100 fs. The obtained final energy distributions for the C2H6 and CO fragments and their relative translational motion were found to be quite similar to those obtained for the acetaldehyde reaction, CH3CHO --> CH4 + CO, in our previous study (Chem. Phys. Lett. 2006, 421, 549) despite the fact that the number of degree of freedom for C2H6 is larger than that for CH4. The coupling between the intrinsic reaction coordinate and one of the generalized normal modes orthogonal to it was predicted substantially strong around s = 1.4 amu(1/2) bohr, and it is expected that the energy flow out of C2H6 proceeds through this coupling. However, the obtained energy distributions strongly suggest that the coupling among the modes in C2H6 is quite small and the intramolecular energy redistribution does not occur efficiently in this molecule.     

A one-electron model for the aqueous electron that includes many-body electron-water polarization: Bulk equilibrium structure, vertical electron binding energy, and optical absorption spectrum

Previously, we reported an electron-water pseudopotential designed to be used in conjunction with a polarizable water model, in order to describe the hydrated electron [L. D. Jacobson et al., J. Chem. Phys. 130, 124115 (2009)]. Subsequently, we found this model to be inadequate for the aqueous electron in bulk water, and here we report a reparametrization of the model. Unlike the previous model, the current version is not fit directly to any observables; rather, we use an ab initio exchange-correlation potential, along with a repulsive potential that is fit to reproduce the density maximum of the excess electron's wave function within the static-exchange approximation. The new parametrization performs at least as well as the previous model, as compared to ab initio benchmarks for (H(2)O)(n) (-) clusters, and also predicts reasonable values for the diffusion coefficient, radius of gyration, and absorption maximum of the bulk species. The new model predicts a vertical electron binding energy of 3.7 eV in bulk water, which is 1.4 eV smaller than the value obtained using nonpolarizable models; the difference represents the solvent's electronic reorganization energy following electron detachment. We find that the electron's first solvation shell is quite loose, which may be responsible for the electron's large, positive entropy of hydration. Many-body polarization alters the electronic absorption line shape in a qualitative way, giving rise to a high-energy tail that is observed experimentally but is absent in previous simulations. In our model, this feature arises from spatially diffuse excited states that are bound only by electronic reorganization (i.e., solvent polarization) following electronic excitation.     

High resolution electron attachment to CO₂ clusters

Electron attachment to CO? clusters performed at high energy resolution (0.1 eV) is studied for the first time in the extended electron energy range from threshold (0 eV) to about 10 eV. Dissociative electron attachment (DEA) to single molecules yields O(-) as the only fragment ion arising from the well known (2)Π(u) shape resonance (ion yield centered at 4.4 eV) and a core excited resonance (at 8.2 eV). On proceeding to CO? clusters, non-dissociated complexes of the form (CO?)(n)(-) including the monomer CO?(-) are generated as well as solvated fragment ions of the form (CO?)(n)O(-). The non-decomposed complexes appear already within a resonant feature near threshold (0 eV) and also within a broad contribution between 1 and 4 eV which is composed of two resonances observed for example for (CO?)(4)(-) at 2.2 eV and 3.1 eV (peak maxima). While the complexes observed around 3.1 eV are generated via the (2)Π(u) resonance as precursor with subsequent intracluster relaxation, the contribution around 2.2 eV can be associated with a resonant scattering feature, recently discovered in single CO? in the selective excitation of the higher energy member of the well known Fermi dyad [M. Allan, Phys. Rev. Lett., 2001, 87, 0332012]. Formation of (CO?)(n)(-) in the threshold region involves vibrational Feshbach resonances (VFRs) as previously discovered via an ultrahigh resolution (1 meV) laser photoelectron attachment method [E. Leber, S. Barsotti, I. I. Fabrikant, J. M. Weber, M.-W. Ruf and H. Hotop, Eur. Phys. J. D, 2000, 12, 125]. The complexes (CO?)(n)O(-) clearly arise from DEA at an individual molecule within the cluster involving both the (2)Π(u) and the core excited resonance.     

Cross sections for electron impact excitation of the vibrationally resolved A 1Pi electronic state of carbon monoxide

The authors report new differential cross section measurements for electron impact excitation of the A (1)Pi(v(')) states of carbon monoxide. The energy range is 20-200 eV. They also reanalyze the A (1)Pi(v(')) manifold cross sections of Middleton et al. [J. Phys. B 26, 1743 (1993)] in order to provide a basis for comparison with our new vibrationally resolved differential cross sections. Excellent agreement is found between the two sets of measurements at all common energies. From 20 to 200 eV the present differential cross sections are extrapolated and integrated, and the corresponding integral excitation cross sections determined. New scaled Born integral cross sections, calculated as a part of the present study, are compared against these experimental integral cross sections, with excellent agreement being found for all the A (1)Pi(v(')=0-7)<--X (1)Sigma(g) (+)(v(")=0) transitions. In addition our scaled Born integral cross sections are found to be in excellent agreement between 300 and 1500 eV with those derived from the previous experiments of Lassettre and Skerbele [J. Chem. Phys. 54, 1597 (1971)] and of Zhong et al. [Phys. Rev. A 55, 1799 (1997)] and from near threshold to 15 eV with those derived from Zobel et al. [J. Phys. B 29, 813 (1996)] and Zetner et al. (J. Phys. B 31, 2395 (1998)].     

Quantum wave packet dynamics of N(2D)+H2 reaction

The quantum wave packet dynamics of the title reaction within the coupled state approximation is examined here and initial state-selected reaction probabilities, integral reaction cross sections, and thermal rate constants are reported. The ab initio potential energy surface of the electronic ground state (1(2)A(")) of the system recently reported by Ho et al. [J. Chem. Phys., 119, 3063 (2003)] is employed in this investigation. All partial wave contributions up to the total angular momentum J=55 were necessary to obtain converged integral reaction cross sections up to a collision energy of 1.0 eV. Thermal rate constants are calculated from the reaction cross sections and compared with the available theoretical and experimental results. Typical resonances formed during the course of the reaction and elucidating the insertion type mechanism for the product formation are calculated. Vibrational energy levels supported by the deep well (approximately 5.5 eV) of the 1(2)A(") potential energy surface of NH(2) are also calculated for the total angular momentum J=0. A statistical analysis of the spacing between the adjacent levels of this energy spectrum is performed and the extent of irregularity in the spectral sequence is assessed.