Photodissociation and recombination of F2 molecule in Ar54 cluster: Nonadiabatic molecular dynamics simulations

My Trajectory in Molecular Reaction Dynamics and Spectroscopy

This is the story of a career in theoretical chemistry during a time of dramatic changes in the field due to phenomenal growth in the availability of computational power. It is likewise the story of the highly gifted graduate students and postdoctoral fellows that I was fortunate to mentor throughout my career. It includes reminiscences of the great mentors that I had and of the exciting collaborations with both experimentalists and theorists on which I built much of my research. This is an account of the developments of exciting scientific disciplines in which I was involved: vibrational spectroscopy, molecular reaction mechanisms and dynamics, e.g., in atmospheric chemistry, and the prediction of new, exotic molecules, in particular noble gas molecules. From my very first project to my current work, my career in science has brought me the excitement and fascination of research. What a wonderful pursuit!

Nonequilibrium Solvent Effects during Photodissociation in Liquids: Dynamical Energy Surfaces, Caging, and Chemical Identity

In the gas phase, potential energy surfaces can be used to provide insight into the details of photochemical reaction dynamics. In solution, however, it is unclear what potential energy surfaces, if any, can be used to describe even simple chemical reactions such as the photodissociation of a diatomic solute. In this paper, we use mixed quantum/classical (MQC) molecular dynamics (MD) to study the photodissociation of Na₂⁺ in both liquid Ar and liquid tetrahydrofuran (THF). We examine both the gas-phase potential surfaces and potentials of mean force (PMF), which assume that the solvent remains at equilibrium with the solute throughout the photodissociation process and show that neither resemble a nonequilibrium dynamical energy surface that is generated by taking the time integral of work. For the photodissociation of Na₂⁺ in liquid Ar, the dynamical energy surface shows clear signatures of solvent caging, and the degree of caging is directly related to the mass of the solvent atoms. For Na₂⁺ in liquid THF, local specific interactions between the solute and solvent lead to changes in chemical identity that create a kinetic trap that effectively prevents the molecule from dissociating. The results show that nonequilibrium effects play an important role even in simple solution-phase reactions, requiring the use of dynamical energy surface to understand such chemical events.

Electron-induced chemistry of methyl chloride caged within amorphous solid water

The interaction of low energy electrons (1.0-25 eV) with methyl-chloride (CD3Cl) molecules, caged within Amorphous Solid Water (ASW) films, 10-120 monolayer (ML) thick, has been studied on top of a Ru(0001) substrate under Ultra High Vacuum (UHV) conditions. While exposing the ASW film to 3 eV electrons a static electric field up to 8 × 10(8) V∕m is developed inside the ASW film due to the accumulation of trapped electrons that produce a plate capacitor voltage of exactly 3 V. At the same time while the electrons continuously strike the ASW surface, they are transmitted through the ASW film at currents of ca. 3 × 10(-7) A. These electrons transiently attach to the caged CD3Cl molecules leading to C-Cl bond scission via Dissociative Electron Attachment (DEA) process. The electron induced dissociation cross sections and product formation rate constants at 3.0 eV incident electrons at ASW film thicknesses of 10 ML and 40 ML were derived from model simulations supported by Thermal Programmed Desorption (TPD) experimental data. For 3.0 eV electrons the CD3Cl dissociation cross section is 3.5 × 10(-16) cm(2), regardless of ASW film thickness. TPD measurements reveal that the primary product is deuterated methane (D3CH) and the minor one is deuterated ethane (C2D6).

Quantum dynamics of solid Ne upon photo-excitation of a NO impurity: a Gaussian wave packet approach

A high-dimensional quantum wave packet approach based on Gaussian wave packets in Cartesian coordinates is presented. In this method, the high-dimensional wave packet is expressed as a product of time-dependent complex Gaussian functions, which describe the motion of individual atoms. It is applied to the ultrafast geometrical rearrangement dynamics of NO doped cryogenic Ne matrices after femtosecond laser pulse excitation. The static deformation of the solid due to the impurity as well as the dynamical response after femtosecond excitation are analyzed and compared to reduced dimensionality studies. The advantages and limitations of this method are analyzed in the perspective of future applications to other quantum solids.

Initial Processes of Laser Induced Diatomic Molecular Photodissociation in Matrices: Quantum Simulations for F2 in Ar in Reduced Dimensionality

Quantum dynamical simulations of ultrafast processes of F2 in an Ar matrix induced by a pump laser pulse show that photoexcitation (i) is followed by (ii) a stretch of the dissociative bond, (iii) collision with, vibrational excitations of, and possible exit through the nearest neighbour atoms of the surrounding cage, and (iv) finally, subsequent rescattering of the dissociated atoms causing vibrations of additional matrix atoms. If the laser pulse ended before the onset of process (iii), the initial steps (i) and (ii) are similar to photodissociations in the gas phase and are, therefore, well described in terms of the wavepacket dynamics of the diatomic molecule in reduced dimensionality. For specific symmetry constraints, the initial processes (i), (ii), and (iii) may be described in terms of just 1- and 2- or 3-dimensional wavepacket dynamics, respectively. At the end of the laser pulse, this low-dimensional wavepacket may serve to prepare the high-dimensional one for the subsequent laser-free dynamics of the total system.

Ultrafast dynamics of halogens in rare gas solids

We perform time resolved pump-probe spectroscopy on small halogen molecules ClF, Cl2, Br2, and I2 embedded in rare gas solids (RGS). We find that dissociation, angular depolarization, and the decoherence of the molecule is strongly influenced by the cage structure. The well ordered crystalline environment facilitates the modelling of the experimental angular distribution of the molecular axis after the collision with the rare gas cage. The observation of many subsequent vibrational wave packet oscillations allows the construction of anharmonic potentials and indicate a long vibrational coherence time. We control the vibrational wave packet revivals, thereby gaining information about the vibrational decoherence. The coherence times are remarkable larger when compared to the liquid or high pressure gas phase. This fact is attributed to the highly symmetric molecular environment of the RGS. The decoherence and energy relaxation data agree well with a perturbative model for moderate vibrational excitation and follow a classical model in the strong excitation limit. Furthermore, a wave packet interferometry scheme is applied to deduce electronic coherence times. The positions of those cage atoms, excited by the molecular electronic transitions are modulated by long living coherent phonons of the RGS, which we can probe via the molecular charge transfer states.

Coherent spin control of matrix isolated molecules by IR+UV laser pulses: quantum simulations for ClF in Ar

Two coherent sequential IR+UV laser pulses may be used to generate two time-dependent nuclear wave functions in electronic excited triplet and singlet states via single (UV) and two photon (IR+UV) excitation pathways, exploiting spin-orbit coupling and vibrational pre-excitation, respectively. These wave functions evolve from different Franck-Condon domains until they overlap in a domain of bond stretching with efficient intersystem crossing. Here, the coherence of the laser pulses is turned into optimal interferences of the wave packets, yielding the total wave packet at the target place, time, and with dominant target spin. The time resolution of spin control is few femtoseconds. The mechanism is demonstrated by means of quantum model simulations for ClF in an Ar matrix.

Army ants algorithm for rare event sampling of delocalized nonadiabatic transitions by trajectory surface hopping and the estimation of sampling errors by the bootstrap method

The most widely used algorithm for Monte Carlo sampling of electronic transitions in trajectory surface hopping (TSH) calculations is the so-called anteater algorithm, which is inefficient for sampling low-probability nonadiabatic events. We present a new sampling scheme (called the army ants algorithm) for carrying out TSH calculations that is applicable to systems with any strength of coupling. The army ants algorithm is a form of rare event sampling whose efficiency is controlled by an input parameter. By choosing a suitable value of the input parameter the army ants algorithm can be reduced to the anteater algorithm (which is efficient for strongly coupled cases), and by optimizing the parameter the army ants algorithm may be efficiently applied to systems with low-probability events. To demonstrate the efficiency of the army ants algorithm, we performed atom-diatom scattering calculations on a model system involving weakly coupled electronic states. Fully converged quantum mechanical calculations were performed, and the probabilities for nonadiabatic reaction and nonreactive deexcitation (quenching) were found to be on the order of 10(-8). For such low-probability events the anteater sampling scheme requires a large number of trajectories ( approximately 10(10)) to obtain good statistics and converged semiclassical results. In contrast by using the new army ants algorithm converged results were obtained by running 10(5) trajectories. Furthermore, the results were found to be in excellent agreement with the quantum mechanical results. Sampling errors were estimated using the bootstrap method, which is validated for use with the army ants algorithm.

Cage motions induced by electronic and vibrational excitations: Cl2 in Ar

Femtosecond dynamics of molecular vibrations as well as cage motions in the B<--X transition of Cl2 in solid Ar have been investigated. We observed molecular vibrational wave-packet motion in experimental pump-probe spectra and an additional oscillation with a 500 fs period which is assigned to the zone-boundary phonon of the Ar crystal. The cage motion is impulsively driven by the B<--X transition due to the expansion of the electronic cloud of the chromophore. To clarify the underlying mechanism, we performed simulations based on the diatomics-in-molecules method which takes into account the different shapes of the Cl2 electronic wave function in the B and X states as well as the anisotropic interaction with the matrix. The simulation results show that Ar atom motion in the (100) plane is initiated by the electronic transition and that only those Ar atoms oscillate coherently with an approximately 500 fs period which are essentially decoupled from the molecular vibration. Their phase and time evolution are in good agreement with the experimentally observed oscillation, supporting the assignment as a displacive excitation of coherent phonons.

Ultrafast solvent-induced spin-flip and nonadiabatic coupling: ClF in argon solids

Femtosecond pump-probe spectra show direct evidence for ultrafast solvent-induced spin flip in photodissociation-recombination events of ClF, a light diatomic molecule, for which the spin-orbit coupling is weak. The bound triplet states ((3)Pi) of ClF are probed and the dynamics for excitation to the singlet state ((1)Pi(1)) is compared with excitation to the triplet state B((3)Pi(0)). The population initially excited to the singlet state (1)Pi(1) is transferred to the bound triplet states (3)Pi within tau(f)=0.5 ps. Oscillations in the spectra indicate wave packet dynamics with the triplet state period of 300 to 400 fs in both cases. According to simulations of F(2)/Ar, most of the initially excited singlet state population is converted to repulsive and weakly bound triplet states within approximately 60 fs. In the first ps, 40% of the triplet population accumulates in the weakly bound (3)Pi states, in good accord with the experiment.