Measures of stochasticity for nonstationary initial preparations of coupled oscillators

The contributions of different eigenstates of a nonlinearly coupled oscillator system to the expansion of a local wave packet are analyzed from an information theoretical point of view. Such a wave packet can be considered as a nonstationary vibrational state of an electronically excited manifold of a molecule after Franck–Condon type initial preparation. The distributions of these contributions are compared to their individual stochastic ideals using Ruch's concept of “mixing distances.” The stochastic ideals are constructed via a probability diffusion process between neighboring states of the original distributions, representing an initial preparation corresponding to a Hamiltonian with only irregular eigenstates. Gaussian minimum uncertainty wave packets as initial states in a two-dimensional nonlinear oscillator system with classically regular and chaotic energy ranges are studied numerically. It is found that distance measures, partly reflecting the “mixing distance” of a distribution from its stochastic ideal, show a large fluctuation in the classical regular energy range and a small fluctuation in the range where most of the classical trajectories move chaotically. This indicates that for this type of initial preparation process the actual location of the initial state in space plays the dominant role for the dynamics in the low-energy range while for wave packets starting near the dissociation energy of the model system this location becomes unimportant.     

Excitation, dynamics, and control of rotationally autoionizing Rydberg states of H2

The dynamics of rotationally autoionizing Rydberg states of molecular hydrogen is investigated using a time-dependent extension of multichannel quantum defect theory, in which the time-dependent wave packets are constructed using first-order perturbation theory. An analytical expression for the complex excitation function for a sequence of Gaussian excitation pulses is derived and then employed to investigate the influence of pairs of pulses with well-defined phase differences on the decay dynamics and final-state composition.     

Molecular Dynamics of Excited State Intramolecular Charge Transfer in Solution by Coherent Nuclear Wave Packets

Revealing a proper reaction coordinate in a chemical reaction is the key step towards elucidation of the molecular reaction dynamics. In this report, we investigated the dynamics of intramolecular charge transfer (ICT) of 8-aminopyrene-1,3,6-trisulfonic acid (APTS) occurring in the excited state by time-resolved fluorescence (TF) and TF spectra. Accurate reaction rates and rate-dependent nuclear wave packets in the product state allow detailed investigation of the molecular reaction dynamics. The ICT rate is solvent dependent: (34 fs)⁻¹, (87 fs)⁻¹, and (∞)⁻¹ in water, formamide, and dimethylformamide, respectively. By recording spectra of the nuclear wave packets for different reaction rates, chemical species responsible for the emission spectra can be positively identified. The origin of the wave packets can be deduced from the amplitude change of the wave packets at different reaction rates, and the vibrational modes that are associated with the reaction coordinate could be identified. Theoretical calculations of the vibrational reorganization energies reproduce the experimental spectrum of the nuclear wave packets and corroborate the conclusions.     

Electronic coherences and vibrational wave-packets in single molecules studied with femtosecond phase-controlled spectroscopy

Employing femtosecond pulse-shaping techniques we investigate ultrafast, coherent and incoherent dynamics in single molecules at room temperature. In first experiments single molecules are excited into their purely electronic 0-0 transition by phase-locked double-pulse sequences with pulse durations of 75 fs and 20 nm spectral band width. Their femtosecond kinetics can then be understood in terms of a 2-level system and modelled with the optical Bloch equations. We find that we observe the coherence decay in single molecules, and the purely electronic dephasing times can be retrieved directly in the time domain. In addition, the Rabi-frequencies and thus the transition dipole moments of single molecules are determined from these data. Upon excitation of single molecules into a vibrational level of the electronically excited state also incoherent intra-molecular vibrational relaxation is recorded. Increasing the spectral band width of the excitation pulses to up to 120 nm (resulting in a transform-limited pulse width of 15 fs) coherent superpositions of excited state vibrational modes, i.e. vibrational wave packets, are excited. The wave-packet oscillations in the excited state potential energy surface are followed in time by a phase-controlled pump-probe scheme, which permits to record wave packet interference, and to determine the energies of vibrational modes and their coupling strengths to the electronic transition.     

Coherent polyatomic dynamics studied by femtosecond time-resolved photoelectron spectroscopy: dissociation of vibrationally excited CS2 in the 6s and 4d Rydberg states

The dissociation dynamics of the 6s and 4d Rydberg states of carbon disulfide (CS(2)*) are studied by time-resolved photoelectron spectroscopy. The CS(2) is excited by two photons of 267 nm (pump) to the 6s and 4d Rydberg states and probed by ionization with either 800 or 400 nm. The experiments can distinguish and successfully track the time dynamics of both spin [1/2] (upper) and [3/2] (lower) cores of the excited Rydberg states, which are split by 60 meV, by measuring the outgoing electron kinetic energies. Multiple mode vibrational wave packets are created within the Rydberg states and observed through recurrence interferences in the final ion state. Fourier transformation of the temporal response directly reveals the coherent population of several electronic states and vibrational modes. The composition of the wave packet is varied experimentally by tuning the excitation frequency to particular resonances between 264 and 270 nm. The work presented here shows that the decay time of the spin components exhibits sensitivity to the electronic and vibrational states accessed in the pump step. Population of the bending mode results in an excited state lifetime of as little as 530 fs, as compared to a several picosecond lifetime observed for the electronic origin bands. Experiments that probe the neutral state dynamics with 400 nm reveal a possible vibrationally mediated evolution of the wave packet to a different Franck-Condon window as a consequence of Renner-Teller splitting. Upon bending, symmetry lowering from D(infinityh) to C(2v) enables ionization to the CS(2) (+) (B (2)Pi(u)) final state. The dissociation dynamics observed are highly mode specific, as revealed by the frequency and temporal domain analysis of the photoelectron spectra.     

Preparation and decay of unstable molecular states undergoing memory effects: Excitation by coherent electromagnetic pulses of arbitrary strength

It is shown that the preparation and decay of non-markovian dissipative systems can be handled without any restriction on both excitation intensity and memory strength. Any perturbative approximation in the treatment of non-markovian interaction is avoided. The theory is explicity related to physical models relevant to the problem of molecular radiationless decay. In this respect, it is shown that a quantitative evaluation of the “radiative state” time evolution along the lines of the recent paper by Rhodes can easily be performed whatever the strength of external interaction. It is emphasized that no quantum beat phenomenon may be interpreted in an unambiguous way in the absence of any information on the dissipation coupling. As a new effect of coherent excitation on the nature of the excited state, it is shown that increasing light intensity may allow the molecule to be prepared in the “radiative state” without using 4 excitations. It is stressed, furthermore, that the excitation of a discrete set of states by long duration pulses has the effect of invalidating the replacement of a complex molecular system by a more simple one, which in turn is to be regarded as the basic idea of the theory. The above physical circumstance seems to be the unique one under which the present approach is no longer applicable. Evidence for this statement is provided by a quantitative simulation of hypothetical experiments of the same kind as the recent ones by Zewail et al.     

Auger decay calculations with core-hole excited-state molecular-dynamics simulations of water

We report a new theoretical procedure for calculating Auger decay transition rates including effects of core-hole excited-state dynamics. Our procedure was applied to the normal and first resonant Auger processes of gas-phase water and compared to high-resolution experiments. In the normal Auger decay, calculated Auger spectra were found to be insensitive to the dynamics, while the repulsive character of the first resonant core-excited state makes the first resonantly excited Auger decay spectra depend strongly on the dynamics. The ultrafast dissociation of water upon O(1s)-->4a(1) excitation was analyzed and found to be very sensitive to initial vibrational distortions in the ground state which furthermore affect the excitation energy. Our calculated spectra reproduce the experimental Auger spectra except for the Franck-Condon vibrational structure which is not included in the procedure. We found that the Auger decay of OH and O fragments contributes to the total intensity, and that the contribution from these fragments increases with increasing excitation energy.     

Two-dimensional infrared spectral signature and hydration of the oxalate dianion

Ultrafast vibrational spectra of the aqueous oxalate ion in the region of its carboxylate asymmetric stretch modes show novel relaxation processes. Two-dimensional infrared vibrational echo spectra and the vibrational dynamics obtained from them along with measurements of the anisotropy decay provide a picture in which the localization of the oxalate vibrational excitation onto the carboxylate groups occurs in ~450 fs. Molecular dynamics simulations are used to characterize the vibrational dynamics in terms of dihedral angle motion between the two carboxylate planes and solvation dynamics. The localization of the oxalate vibrational excitation onto the carboxylates is induced by the fluctuations in the carboxylate vibrational frequencies which are shown by theory and experiment to have a similar correlation time as the anisotropy decay.     

Laser adiabatic manipulation of the bond length of diatomic molecules with a single chirped pulse

We propose and test numerically a scheme for controlling the bond distance in a diatomic molecule that requires the use of a single chirped pulse. The laser prepares a superposition state of both nuclear and electronic degrees of freedom, where the main character of the electronic wave function is that of an excited dissociative state. The main limitation of the scheme is the need of ultra broadband pulses, where the bandwidth must be of the order of the dissociation energy to achieve large bond elongations. The scheme can be used to deform the bond during the laser excitation to an arbitrary large and constant value, or to allow slow time-dependent bond elongations. Additionally, the scheme can be used to prepare highly excited vibrational wave packets in the ground potential after the pulse is switched off, at the expense of losing some population that dissociates. These wave packets are initially localized at the outer well of the potential, at energies controllable by the excitation process.     

Slow periodic oscillations in time domain dynamics of NO2

The authors investigated the time domain nonadiabatic dynamics of NO2 on the coupled X 2A1 and A 2B2 electronic states by launching wave packets on the excited electronic state and focused on the evolution at long times (t>200 fs), which has received little attention up to now. The authors showed that the initial fast spreading of the wave packets is followed at all energies by slow periodic intramolecular vibronic energy redistribution (IVER) with periods in the range of 0.3 to several tens of picoseconds. These energy transfers lead to oscillations with the same periods in the population of each electronic state. Propagation of wave packets indicates that IVER frequencies also dominate the fluctuations of the squared modulus of the autocorrelation function |A(t)|2 at energies not too high above the bottom of the conical intersection, but this is no longer the case at higher energies. For example, for initial wave packets prepared by almost vertical excitation of the vibrational ground state of the ground electronic surface, the oscillations of |A(t)|2 essentially reflect the detuning from 1:2 resonance between the frequency of the bend and that of the symmetric stretch in the excited electronic state. These theoretical results were used to discuss the possible origin of the low-frequency oscillations which were recently observed in time domain experimental spectra of NO2.     

Nonadiabatic effects on peptide vibrational dynamics induced by conformational changes

Quantum dynamical simulations of vibrational spectroscopy have been carried out for glycine dipeptide (CH(3)-CO-NH-CH(2)-CO-NH-CH(3)). Conformational structure and dynamics are modeled in terms of the two Ramachandran dihedral angles of the molecular backbone. Potential energy surfaces and harmonic frequencies are obtained from electronic structure calculations at the density functional theory (DFT) [B3LYP/6-31+G(d)] level. The ordering of the energetically most stable isomers (C(7) and C(5)) is reversed upon inclusion of the quantum mechanical zero point vibrational energy. Vibrational spectra of various isomers show distinct differences, mainly in the region of the amide modes, thereby relating conformational structures and vibrational spectra. Conformational dynamics is modeled by propagation of quantum mechanical wave packets. Assuming a directed energy transfer to the torsional degrees of freedom, transitions between the C(7) and C(5) minimum energy structures occur on a sub-picosecond time scale (700...800 fs). Vibrationally nonadiabatic effects are investigated for the case of the coupled, fundamentally excited amide I states. Using a two state-two mode model, the resulting wave packet dynamics is found to be strongly nonadiabatic due to the presence of a seam of the two potential energy surfaces. Initially prepared adiabatic vibrational states decay upon conformational change on a time scale of 200...500 fs with population transfer of more than 50% between the coupled amide I states. Also the vibrational energy transport between localized (excitonic) amide I vibrational states is strongly influenced by torsional dynamics of the molecular backbone where both enhanced and reduced decay rates are found. All these observations should allow the detection of conformational changes by means of time-dependent vibrational spectroscopy.     

Theoretical analysis of photoinduced H-atom elimination in thiophenol

The photoinduced hydrogen elimination reaction in thiophenol via the conical intersections of the dissociative (1)πσ? excited state with the bound (1)ππ? excited state and the electronic ground state has been investigated with ab initio electronic-structure calculations and time-dependent quantum wave-packet calculations. A screening of the coupling constants of the symmetry-allowed coupling modes at the (1)ππ?-(1)πσ? and (1)πσ?-S(0) conical intersection shows that the SH torsional mode is by far the most important coupling mode at both conical intersections. A model including three intersecting potential-energy surfaces (S(0), (1)ππ?, (1)πσ?) and two nuclear degrees of freedom (SH stretch and SH torsion) has been constructed on the basis of ab initio complete-active-space self-consistent field and multireference second-order perturbation theory calculations. The nonadiabatic quantum wave-packet dynamics initiated by optical excitation of the (1)ππ? and (1)πσ? states has been explored for this three-state two-coordinate model. The photodissociation dynamics is characterized in terms of snapshots of time-dependent wave packets, time-dependent electronic population probabilities, and the branching ratio of the (2)σ/(2)π electronic states of the thiophenoxyl radical. The dependence of the timescale of the photodissociation process and the branching ratio on the initial excitation of the SH stretching and SH torsional vibrations has been analyzed. It is shown that the node structure, which is imposed on the nuclear wave packets by the initial vibrational preparation as well as by the transitions through the conical intersections, has a profound effect on the photodissociation dynamics. The effect of additional weak coupling modes of CC twist (ν(16a)) and ring-distortion (ν(16b)) character has been investigated with three-dimensional and four-dimensional time-dependent wave-packet calculations, and has been found to be minor.     

A transition state view on reactive scattering: initial state-selected reaction probabilities for the H + CH4 → H2 + CH3 reaction studied in full dimensionality

Initial state-selected reaction probabilities for the H+CH(4)→H(2)+CH(3) reaction are computed for vanishing total angular momentum by full-dimensional calculations employing the multiconfigurational time-dependent Hartree approach. An ensemble of wave packets completely describing reactivity for total energies up to 0.58 eV is constructed in the transition state region by diagonalization of the thermal flux operator. These wave packets are then propagated into the reactant asymptotic region to obtain the initial state-selected reaction probabilities. Reaction probabilities for reactants in all rotational states of the vibrational 1A(1), 1F(2), and 1E levels of methane are presented. Vibrational excitation is found to decrease reactivity when reaction probabilities at equivalent total energies are compared but to increase reaction probabilities when the comparison is done at the basis of equivalent collision energies. Only a fraction of the initial vibrational energy can be utilized to promote the reaction. The effect of rotational excitation on the reactivity differs depending on the initial vibrational state of methane. For the 1A(1) and 1F(2) vibrational states of methane, rotational excitation decreases the reaction probability even when comparing reaction probabilities at equivalent collision energies. In contrast, rotational energy is even more efficient than translational energy in increasing the reaction probability when the reaction starts from the 1E vibrational state of methane. All findings can be explained employing a transition state based interpretation of the reaction process.     

Femtosecond dynamics of isolated phenylcarbenes

Understanding the primary photophysical processes in molecules is essential for interpreting their photochemistry, because molecules rarely react from the initially excited electronic state. In this study the ultrafast excited-state dynamics of chlorophenylcarbene (CPC) and trifluoromethylphenylcarbene (TFPC), two species that are considered as models for carbene dynamics, were investigated by femtosecond time-resolved pump probe spectroscopy in the gas phase. Their dynamics was followed in real time by time-resolved photoionization and photoelectron imaging. CPC was excited at 265 nm into the 3 1A' state, corresponding to excitation from a pi-orbital of the aromatic ring into the LUMO. The LUMO contains a contribution of the p-orbital at the carbene center. Three time constants are apparent in the photoelectron images: A fast decay process with tau1 approximately 40 fs, a second time constant of tau2 approximatley 350 fs, and an additional time constant of tau3 approximately 1 ps. The third time constant is only visible in the time-dependence of low kinetic energy electrons. Due to the dense manifold of excited states between 3.9 and 5 eV, known from ab initio calculations, the recorded time-resolved electron images show broad and unstructured bands. A clear population transfer between the states thus can not directly be observed. The fast deactivation process is linked to either a population transfer between the strongly coupled excited states between 3.9 and 5.0 eV or the movement of the produced wave packet out of the Franck-Condon region. Since the third long time constant is only visible for photoelectrons at low kinetic energy, evidence is given that this time constant corresponds to the lifetime of the lowest excited A 1A' state. The remaining time constant reflects a deactivation of the manifold of states in the range 3.9-5.0 eV down to the A 1A' state.     

Efficient and robust strong-field control of population transfer in sensitizer dyes with designed femtosecond laser pulses

We demonstrate control of electronic population transfer in molecules with the help of appropriately shaped femtosecond laser pulses. To this end we investigate two photosensitizer dyes in solution being prepared in the triplet ground state. Excitation within the triplet system is followed by intersystem crossing and the corresponding singlet fluorescence is monitored as a measure of population transfer in the triplet system. We record control landscapes with respect to the fluorescence intensity on both dyes by a systematic variation of laser pulse shapes combining second order and third order dispersion. In the strong-field regime we find highly structured topologies with large areas of maximum or minimum population transfer being insensitive over a certain range of applied laser intensities thus demonstrating robustness. We then compare our experimental results with simulations on generic molecular potentials by solving the time-dependent Schr?dinger equation for excitation with shaped pulses. Control landscapes with respect to population transfer confirm the general trends from experiments. An analysis of regions with maximum or minimum population transfer indicates that coherent processes are responsible for the outcome of our excitation process. The physical mechanisms of joint motion of ground and excited state wave packets or population of a vibrational eigenstate in the excited state permit us to discuss the molecular dynamics in an atom-like picture.     

Theoretical study of solvent influence on the electronic absorption and emission spectra of kynurenine

Neutral/zwitterionic form equilibrium, excited state wave functions, absorption and emission spectra of kynurenine (KN) in various solvents (water, methanol, ethanol, and dimethylsulfoxide) have been studied theoretically. The ground electronic state geometries have been optimized by density functional theory methods; the geometries of the first two singlets excited electronic states have been optimized using the CASSCF technique. The influence of the solvent was taken into account by the calculation of the solvation free energies using the Polarizable Continuum Model (PCM). The spectra of electronic absorption and fluorescence emission have been calculated by the CS‐INDO S‐CI and SDT‐CI methods [Momicchioli, Baraldi, and Bruni, Chem Phys, 1983, 82, 229]. The calculated data reproduce the experimental positions of maxima and the solvent‐induced shifts of the absorption and emission bands well. The energy gap between the two lowest excited states of KN increases from aprotic to protic solvents. This fact suggests that the “proximity effect” cannot be responsible for the ultrafast decay of KN fluorescence in protic solvents. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Rate kernel theory for pseudo-first-order kinetics of diffusion-influenced reactions and application to fluorescence quenching kinetics

Theoretical foundation of rate kernel equation approaches for diffusion-influenced chemical reactions is presented and applied to explain the kinetics of fluorescence quenching reactions. A many-body master equation is constructed by introducing stochastic terms, which characterize the rates of chemical reactions, into the many-body Smoluchowski equation. A Langevin-type of memory equation for the density fields of reactants evolving under the influence of time-independent perturbation is derived. This equation should be useful in predicting the time evolution of reactant concentrations approaching the steady state attained by the perturbation as well as the steady-state concentrations. The dynamics of fluctuation occurring in equilibrium state can be predicted by the memory equation by turning the perturbation off and consequently may be useful in obtaining the linear response to a time-dependent perturbation. It is found that unimolecular decay processes including the time-independent perturbation can be incorporated into bimolecular reaction kinetics as a Laplace transform variable. As a result, a theory for bimolecular reactions along with the unimolecular process turned off is sufficient to predict overall reaction kinetics including the effects of unimolecular reactions and perturbation. As the present formulation is applied to steady-state kinetics of fluorescence quenching reactions, the exact relation between fluorophore concentrations and the intensity of excitation light is derived.     

The influence of ultrashort excitations on the vibrational dynamics in a diatomic molecule

The problem of vibrational wave packet dynamics in the system of two electronic states of a diatomic molecule, where the states are coupled by infinitely short light pulses, is solved. The electronic states were modeled by shifted harmonic oscillators with different frequencies. Exact expressions for the probability densities of the wave packets in the ground and excited states were derived. The spatial, spectral, and temporal characteristics of the wave packets, namely, the range of motion, spatial width, mean energy, spectral width (the mean number of vibrational states in a wave packet), and the autocorrelation function, were calculated as functions of the molecular parameters (the frequency ratio and the distance between the potential minima) and of the delay time between the light pulses. The possibility of controlling the mean energy and spectral width of the wave packets in the ground electronic state by varying the delay time is considered. It was shown that "squeezed" wave packets can be prepared in the ground electronic state if the upper electronic state is shallow.     

Calculation of infrared absorption spectra using Gaussian variational wave packets

Formulas are developed for calculation of infrared absorption spectra with the help of Gaussian variational wave packets. The excitation is viewed as occurring to a linear combination of basis states corresponding to a Gaussian multiplied by Hermite polynomials of coordinates. Application of the formulas is examplified on HDO isolated in a D2O matrix. A simple approximation is proposed for the fundamental transitions.     

Initial dynamics of the Norrish Type I reaction in acetone: probing wave packet motion

The Norrish Type I reaction in the S(1) (nπ*) state of acetone is a prototype case of ketone photochemistry. On the basis of results from time-resolved mass spectrometry (TRMS) and photoelectron spectroscopy (TRPES) experiments, it was recently suggested that after excitation the wave packet travels toward the S(1) minimum in less than 30 fs and stays there for more than 100 picoseconds [Chem. Phys. Lett.2008, 461, 193]. In this work we present simulated TRMS and TRPES signals based on ab initio multiple spawning simulations of the dynamics during the first 200 fs after excitation, getting quite good agreement with the experimental signals. We can explain the ultrafast decay of the experimental signals in the following manner: the wave packet simply travels, mainly along the deplanarization coordinate, out of the detection window of the ionizing probe. This window is so narrow that subsequent revival of the signal due to the coherent deplanarization vibration is not observed, meaning that from the point of view of the experiment the wave packets travels directly to the S(1) minimum. This result stresses the importance of pursuing a closer link to the experimental signal when using molecular dynamics simulations in interpreting experimental results.