New mechanistic insight into electronically excited CO-NiO(100): a quantum dynamical analysis

In a recent study Redlich et al. [Redlich et al., Chem. Phys. Lett. 2006, 420, 110] measured the velocity distribution of CO molecules desorbing from a NiO(100) surface after irradiation with an ultraviolet (UV) laser pulse. Due to the complexity of the involved processes no experimental evidence on the excitation and desorption mechanism could be obtained. In recent ab initio studies Mehdaoui et al. [Mehdaoui et al., Phys. Rev. Lett. 2007, 98, 037601] have shown that a 5sigma --> 2pi* (a (3)Pi) like transition within the CO adsorbate is most likely the crucial excitation step in the CO-NiO(100) system. At first sight this seems unlikely, since the interaction of CO molecules with the NiO(100) surface is very weak (-0.30 eV) and the corresponding CO gas phase transition energy is about 1.5 eV higher than the laser pulse energy of 4.66 eV used in the experiment. In this work we give further insight into relevant electronically excited states and identify the desorption mechanism by analysing the dynamical processes after laser excitation by quantum dynamical wave packet simulations on the basis of three-dimensional (3D) ab initio potential energy surfaces. The results corroborate the so far discussed excitation mechanism, which proposes the formation of a genuine C-Ni bond as the driving force for photodesorption, as the crucial excitation step.     

Photodesorption of diatomic molecules from surfaces: A theoretical approach based on first principles

Photodesorption of small molecules from surfaces is one of the most fundamental processes in surface photochemistry. Despite its apparent simplicity, a microscopic understanding beyond a qualitative picture still poses a true challenge for theory. While the dynamics of nuclear motion can be treated on various levels of sophistication, all approaches suffer from the lack of sufficiently accurate potential energy surfaces, in particular for electronically excited states involved in the desorption scenario.In the last decade, a systematic and accurate methodology has been developed which allows a reliable calculation of accurate ground and excited state potential energy surfaces (PES) for different adsorbate–substrate systems. These potential energy surfaces serve as a prerequisite for subsequent quantum dynamical wave packet calculations, which allow for a direct simulation of experimentally observable quantities such as quantum state resolved velocity distributions.In the first part of this review, we will focus on scalar properties of desorbing diatomic molecules from insulating surfaces, where we also present a recently developed strategy of obtaining accurate potential energy surfaces using quantum chemical approaches. In general, diatomic molecules on large band gap materials such as oxide surfaces are studied which allows the use of sufficiently large cluster models and accurate ab initio methods beyond density functional theory (DFT). In the second part, we will focus on the vectorial aspects of the dynamics of nuclear motion and present simulations of experimentally accessible observables such as velocity distributions, Doppler profiles and alignment parameters. For each system, the microscopic mechanism of photodesorption is elucidated. We will demonstrate that the driving force of surface photochemistry is strongly dependent on details of the electronic structure of the adsorbate–substrate systems. This implies that great caution is advisable if experimental results are interpreted using empirical or semi-empirical models.     

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.     

Quantum mechanical investigations of the N(4S)+O2(X 3Sigmag -)-->NO(X 2Pi)+O(3P) reaction

The reaction between energetic nitrogen atoms and oxygen molecules has received important attention in connection with nitric oxide chemistry in the lower thermosphere. We report time-independent quantum mechanical calculations of the N(4S)+O2-->NO+O reaction employing the X 2A' and a 4A' electronic potential energy surfaces of Sayos et al. [J. Chem. Phys. 117, 670 (2002)]. We confirm the production of highly vibrationally excited NO molecules, consistent with previous semiclassical and more recent time-dependent quantum wave packet studies. Calculations are carried out for total angular momentum quantum number J=0 and cross sections and rate coefficients are extracted using the J-shifting approximation. The results are in good agreement with available experimental and theoretical data.     

Mechanistic aspects of proton chain transfer in the green fluorescent protein. Part II. A comparison of minimal quantum chemical models

In this paper we report the results of extensive quantum chemical reaction pathway calculations for the electronic ground state of several different cluster models that mimic the proton chain transfer path within the green fluorescent protein (GFP). Our principal objective is to establish the robustness with respect to variations in the model of our recent mechanistic inferences for the ground state proton chain transfer [S. Wang and S. C. Smith, J. Phys. Chem. B, 2006, 110, 5084]. Additionally, comparison of our ground state results with the excited state proton transfer (ESPT) study by Vendrell et al. [O. Vendrell, R. Gelabert, M. Moreno and J. M. Lluch, J. Am. Chem. Soc., 2006, 128, 3564] leads to the conclusion that the mechanism of proton chain transfer may be expected to be analogous in ground and excited states, principally because in both cases the loss of the chromophore's phenolic proton contributes strongly to the reaction coordinate only late in the reaction path.     

Wave packet simulation of nonadiabatic dynamics in highly excited 1,3-dibromopropane

We have conducted wave packet simulations of excited-state dynamics of 1,3-dibromopropane (DBP) with the aim of reproducing the experimental results of the gas-phase pump-probe experiment by Kotting et al. [ Kotting, C. ; Diau, E. W.-G. ; S?lling, T. I. ; Zewail, A. H. J. Phys. Chem. A 2002, 106, 7530 ]. In the experiment, DBP is excited to a Rydberg state 8 eV above the ground state. The interpretation of the results is that a torsional motion of the bromomethylene groups with a vibrational period of 680 fs is activated upon excitation. The Rydberg state decays to a valence state, causing a dissociation of one of the carbon bromine bonds on a time scale of 2.5 ps. Building the theoretical framework for the wave packet propagation around this model of the reaction dynamics, the simulations reproduce, to a good extent, the time scales observed in the experiment. Furthermore, the simulations provide insight into how the torsion motion influences the bond breakage, and we can conclude that the mechanism that delays the dissociation is solely the electronic transition from the Rydberg state to the valence state and does not involve, for example, intramolecular vibrational energy redistribution (IVR).     

Three‐body molecular states of the system in the Born–Oppenheimer approximation

In this work, we present a quantum mechanical treatment of the three‐body molecular system in the Born–Oppenheimer (BO) approximation, were the nuclei dynamics is evaluated over the potential energy surfaces (PES) induced by the electronic states. The PES corresponding to the two lowest electronic levels are the ones described by Martinazzo et al. (Chem. Phys. 2003, 287, 335), and are used to write the three‐body Schrdinger equation of the three atomic system. We use the generalized Sturmian functions method to expand the wave functions in each (distinguishable) pair of relative coordinates or Jacobi pairs, and analyze the convergence differences between the series. A partial‐wave decomposition of the potential is proposed to simplify the Hamiltonian matrix element calculation. Bound states are considered for the ground and first excited electronic PES, the spreading of energies after sudden electronic transitions studied, and the break‐up probability induced by the sudden change of the PES evidenced.     

Preparation and resolution of molecular states by coherent sequences of phase-locked ultrashort laser pulses

We study the application of nonlinear wave packet interferometry to the preparation and resolution of the overlaps of nonstationary nuclear wave functions evolving in an excited electronic state of a diatomic molecule. It is shown that possible experiments with two phase-locked ultrashort pulsepairs can be used to determine a specific vibrational wave packet state in terms of coherent states of the ground electronic state. We apply this scheme to an idealized molecule with harmonic potential energy surfaces and to the X <-- B transition states of the iodine molecule. Our results indicate that this scheme is very promising as a potential tool to quantum control.     

Simulation of vibrational dephasing of I(2) in solid Kr using the semiclassical Liouville method

In this paper, we present simulations of the decay of quantum coherence between vibrational states of I(2) in its ground (X) electronic state embedded in a cryogenic Kr matrix. We employ a numerical method based on the semiclassical limit of the quantum Liouville equation, which allows the simulation of the evolution and decay of quantum vibrational coherence using classical trajectories and ensemble averaging. The vibrational level-dependent interaction of the I(2)(X) oscillator with the rare-gas environment is modeled using a recently developed method for constructing state-dependent many-body potentials for quantum vibrations in a many-body classical environment [J. M. Riga, E. Fredj, and C. C. Martens, J. Chem. Phys. 122, 174107 (2005)]. The vibrational dephasing rates gamma(0n) for coherences prepared between the ground vibrational state mid R:0 and excited vibrational state mid R:n are calculated as a function of n and lattice temperature T. Excellent agreement with recent experiments performed by Karavitis et al. [Phys. Chem. Chem. Phys. 7, 791 (2005)] is obtained.     

Theoretical study of the electronic nonadiabatic transitions in the photoelectron spectroscopy of F2O

The photoelectron spectrum of F2O pertaining to ionizations to the ground (X2B1) and low-lying excited electronic states (A2B2, B2A1, and C2A2) of F2O+ is investigated theoretically. The near equilibrium potential energy surfaces of the ground electronic state (X2B1) of F2O and the mentioned ground and excited electronic states of F2O+ reported by Wang et al. ( J. Chem. Phys. 2001, 114, 10682) for the C2v configuration are extended for the Cs geometry assuming a harmonic vibration along the asymmetric stretching mode. The vibronic interactions between the A2B2 and B2A1 electronic states of F2O+ are treated within a linear coupling approach, and the strength of the vibronic coupling parameter is calculated by an ab initio method. The nuclear dynamics is simulated by both time-independent quantum mechanical and time-dependent wave packet approaches. Although the first photoelectron band exhibits resolved vibrational progression along the symmetric stretching mode, the second one is highly overlapping. The latter is attributed to the nonadiabatic interactions among the energetically close A2B2, B2A1, and C2A2 electronic states of F2O+. The theoretical findings are in good accord with the available experimental results.     

Two-dimensional surrogate Hamiltonian investigation of laser-induced desorption of NO/NiO(100)

The photodesorption of NO from NiO(100) is studied from first principles, with electronic relaxation treated by the use of the surrogate Hamiltonian approach. Two nuclear degrees of freedom of the adsorbate-substrate system are taken into account. To perform the quantum dynamical wave-packet calculations, a massively parallel implementation with a one-dimensional data decomposition had to be introduced. The calculated desorption probabilities and velocity distributions are in qualitative agreement with experimental data. The results are compared to those of stochastic wave-packet calculations where a sufficiently large number of quantum trajectories is propagated within a jumping wave-packet scenario.     

The structure of 5-cyanoindole in the ground and the lowest electronically excited singlet states, deduced from rotationally resolved electronic spectroscopy and ab initio theory

The structure and electronic properties of the electronic ground and the lowest excited singlet states of 5-cyanoindole (5CI) were determined using rotationally resolved spectroscopy of the vibrationless electronic origin of 5CI. In contrast to most other indole derivatives, the lowest excited state of 5CI is determined to be of L(a) character. The conventional approximate coupled cluster singles and doubles model (CC2) fails to describe the geometry of the excited state correctly. Nevertheless, scaling the spin components of equal and opposite spins within the CC2 model as proposed by Hellweg et al. (Phys. Chem. Chem. Phys., 2008, 10, 1159) resulted in very good geometry parameters for the excited state.     

The role of the excited electronic states in the C+ + H2O reaction

The electronic excited states of the [COH2]+ system have been studied in order to establish their role in the dynamics of the C+ + H2O-->[COH]+ +H reaction, which is a prototypical ion-molecule reaction. The most relevant minima and saddle points of the lowest excited state have been determined and energy profiles for the lowest excited doublet and quartet electronic states have been computed along the fragmentation and isomerization coordinates. Also, nonadiabatic coupling strengths between the ground and the first excited state have been computed where they can be large. Our analysis suggests that the first excited state could play an important role in the generation of the formyl isomer, which has been detected in crossed beam experiments [D. M. Sonnenfroh et al., J. Chem. Phys. 83, 3985 (1985)], but could not be explained in quasiclassical trajectory computations [Y. Ishikawa et al., Chem. Phys. Lett. 370, 490 (2003); J. R. Flores, J. Chem. Phys. 125, 164309 (2006)].     

A quantum wave packet dynamical study of the electronic and spin-orbit coupling effects on the resonances in Cl(2P) + H2 scattering

Dynamical resonances in Cl(2P) + H2 scattering are investigated with the aid of a time-dependent wave packet approach using the Capecchi-Werner coupled ab initio potential energy surfaces [Phys. Chem. Chem. Phys. 2004, 6, 4975]. The resonances arising from the prereactive van der Waals well (approximately 0.5 kcal/mol) and the transition-state (TS) region of the 2Sigma(1/2) ground spin-orbit (SO) state of the Cl(2P) + H2 system are calculated and assigned by computing their eigenfunctions and lifetimes. The excitation of even quanta along the bending coordinate of the resonances is observed. The resonances exhibit an extended van der Waals progression, which can be attributed to the dissociative states of ClH2. Excitation of H2 vibration is also identified in the high-energy resonances. The effect of the excited 2P(1/2) SO state of Cl on these resonances is examined by considering the electronic and SO coupling in the dynamical simulations. While the electronic coupling has only a minor impact on the resonance structures, the SO coupling has significant effect on them. The nonadiabatic effect due to the SO coupling is stronger, and as a result, the spectrum becomes broad and diffuse particularly at high energies. We also report the photodetachment spectrum of ClD2- and compare the theoretical findings with the available experimental results.     

Experimental and theoretical differential cross sections for the N(2D) + H2 reaction

In this paper, we report a combined experimental and theoretical study on the dynamics of the N(2D) + H2 insertion reaction at a collision energy of 15.9 kJ mol(-1). Product angular and velocity distributions have been obtained in crossed beam experiments and simulated by using the results of quantum mechanical (QM) scattering calculations on the accurate ab initio potential energy surface (PES) of Pederson et al. (J. Chem. Phys. 1999, 110, 9091). Since the QM calculations indicate that there is a significant coupling between the product angular and translational energy distributions, such a coupling has been explicitly included in the simulation of the experimental results. The very good agreement between experiment and QM calculations sustains the accuracy of the NH2 ab initio ground state PES. We also take the opportunity to compare the accurate QM differential cross sections with those obtained by two approximate methods, namely, the widely used quasiclassical trajectory calculations and a rigorous statistical method based on the coupled-channel theory.     

Theoretical investigation of rotationally inelastic collisions of CH2(X) with helium

Following our earlier work on collisions of He with the methylene radical in its excited ?(1)A(1) state [L. Ma, M. H. Alexander, and P. J. Dagdigian, J. Chem. Phys. 134, 154307 (2011)], we investigate here the analogous relaxation of CH(2) in its ground X(3)B(1) electronic state. The molecule is treated as semi-rigid, with fixed bond lengths but a varying bond angle. We use an ab initio potential energy surface (PES) which is averaged over the CH(2) bending angle weighted by the square of the bending wave function. The PES for the interaction of He with CH(2) in the X state is considerably less anisotropic than for interaction with the ? state since the two 2p electrons on the C atom are evenly distributed among the bonding and non-bonding molecular orbitals. We report quantum scattering calculations of state-to-state and total removal cross sections as well as total removal rate constants at room temperature. Because of the less pronounced anisotropy, these cross sections and rate constants are considerably smaller than for collisions of CH(2)(?) with He. Finally, we investigate the dependence of rotational inelasticity on the bending vibrational quantum number.     

Modeling of HeN+ clusters. II. Calculation of He3+ vibrational spectrum

We have computed the vibrational spectrum of the helium ionized trimer He(3)(+) using three different potential energy surfaces [D. T. Chang and G. L. Gellene, J. Chem. Phys. 119, 4694 (2003); E. Scifoni et al., ibid. 125, 164304 (2006); I. Paidarova et al., Chem. Phys. 342, 64 (2007)]. Differences in the details of these potential energy surfaces induce discrepancies between bound state energies of the order of 0.01 eV. The effects of the geometric phase induced by the conical intersection between the ground electronic potential energy surface and the first excited one are studied by computing vibrational spectra with and without this phase. The six lowest vibrational bound states are negligibly affected by the geometric phase. Indeed, they correspond to wavefunctions localized in the vicinity of the linear symmetric configurations and can be assigned well defined vibrational quantum numbers. On the other hand, higher excited states are delocalized, cannot be assigned definite vibrational quantum numbers, and the geometric phase shifts their energies by approximately 0.005 eV.     

Pump-dump iterative squeezing of vibrational wave packets

The free motion of a nonstationary vibrational wave packet in an electronic potential is a source of interesting quantum properties. In this work we propose an iterative scheme that allows continuous stretching and squeezing of a wave packet in the ground or in an excited electronic state, by switching the wave function between both potentials with pi pulses at certain times. Using a simple model of displaced harmonic oscillators and delta pulses, we derive the analytical solution and the conditions for its possible implementation and optimization in different molecules and electronic states. We show that the main constraining parameter is the pulse bandwidth. Although in principle the degree of squeezing (or stretching) is not bounded, the physical resources increase quadratically with the number of iterations, while the achieved squeezing only increases linearly.     

Theoretical description of femtosecond fluorescence depletion spectrum of molecules in solution

A theoretical model used for calculating the fluorescence depletion spectrum (FDS) of molecules in liquids induced by femtosecond pump-probe laser pulses is proposed based on the reduced density matrix theory. The FDS intensity is obtained by calculating the stimulated emission of the excited electronic state. As an application of the theoretical model, the FDS of oxazine 750 (OX-750) molecule in acetone solution is calculated. The simulated FDS agrees with the experimental result of Liu et al. [J. Y. Liu et al., J. Phys. Chem. A 107, 10857 (2003)]. The calculated vibrational relaxation rate is 2.5 ps(-1) for the OX-750 molecule. Vibrational population dynamics and wave packet evolution in the excited state are described in detail. The effect of the probe pulse parameter on the FDS is also discussed.