Ultrafast excited-state dynamics of tetraphenylethylene studied by semiclassical simulation

Detailed simulation study is reported for the excited-state dynamics of photoisomerization of cis-tetraphenylethylene (TPE) following excitation by a femtosecond laser pulse. The technique for this investigation is semiclassical dynamics simulation, which is described briefly in the paper. Upon photoexcitation by a femtosecond laser pulse, the stretching motion of the ethylenic bond of TPE is initially excited, leading to a significant lengthening of ethylenic bond in 300 fs. Twisting motion about the ethylenic bond is activated by the energy released from the relaxation of the stretching mode. The 90 degrees twisting about the ethylenic bond from an approximately planar geometry to nearly a perpendicular conformation in the electronically excited state is completed in 600 fs. The torsional dynamics of phenyl rings which is temporally lagging behind occurs at about 5 ps. Finally, the twisted TPE reverts to the initial conformation along the twisting coordinate through nonadiabatic transitions. The simulation results provide a basis for understanding several spectroscopic observations at molecular levels, including ultrafast dynamic Stokes shift, multicomponent fluorescence, viscosity dependence of the fluorescence lifetime, and radiationless decay from electronically excited state to the ground state along the isomerization coordinate.     

Picosecond IR-UV pump-probe spectroscopic study on the intramolecular vibrational energy redistribution of NH2 and CH stretching vibrations of jet-cooled aniline

Intramolecular vibrational energy redistribution (IVR) of the NH2 symmetric and asymmetric stretching vibrations of jet-cooled aniline has been investigated by picosecond time-resolved IR-UV pump-probe spectroscopy. A picosecond IR laser pulse excited the NH2 symmetric or asymmetric stretching vibration of aniline in the electronic ground state and the subsequent time evolutions of the excited level as well as redistributed levels were observed by a picosecond UV pulse. The IVR lifetimes for symmetric and asymmetric stretches were obtained to be 18 and 34 ps, respectively. In addition, we obtained the direct evidence that IVR proceeds via two-step bath states; that is, the NH2 stretch energy first flows into the doorway state and the energy is further dissipated into dense bath states. The rate constants of the second step were estimated to be comparable to or slower than those of the first step IVR. The relaxation behavior was compared with that of IVR of the OH stretching vibration of phenol [Y. Yamada, T. Ebata, M. Kayano, and M. Mikami J. Chem. Phys. 120, 7400 (2004)]. We found that the second step IVR process of aniline is much slower than that of phenol, suggesting a large difference of the "doorway state increasing the dense bath states" anharmonic coupling strength between the two molecules. We also observed IVR of the CH stretching vibrations, which showed much faster IVR behavior than that of the NH2 stretches. The fast relaxation is described by the interference effect, which is caused by the coherent excitation of the quasistationary states.     

Picosecond IR-UV pump-probe spectroscopic study of the dynamics of the vibrational relaxation of jet-cooled phenol. I. Intramolecular vibrational energy redistribution of the OH and CH stretching vibrations of bare phenol

The intramolecular vibrational energy redistribution (IVR) of the OH stretching vibration of jet-cooled phenol-h6 (C6H8OH) and phenol-d8 (C6D8OH) in the electronic ground state has been investigated by picosecond time-resolved IR-UV pump-probe spectroscopy. The OH stretching vibration of phenol was excited with a picosecond IR laser pulse, and the subsequent temporal evolutions of the initially excited level and the redistributed ones due to the IVR were observed by multiphoton ionization detection with a picosecond UV pulse. The IVR lifetime for the OH stretch vibration of phenol-h6 was determined to be 14 ps, while that of the OH stretch for phenol-d8 was found to be 80 ps. This remarkable change of the IVR rate constant upon the dueteration of the CH groups strongly suggests that the "doorway states" for the IVR from the OH level would be the vibrational states involving the CH stretching modes. We also investigated the IVR rate of the CH stretching vibration for phenol-h6. It was found that the IVR lifetime of the CH stretch is less than 5 ps. The fast IVR is described by the strong anharmonic resonance of the CH stretch with many other combinations or overtone bands.     

Modeling the photophysics and photochromic potential of 1,2-dihydronaphthalene (DHN): a combined CASPT2//CASSCF-topological and MMVB-dynamical investigation

The photochemical ring opening of 1,2-dihydronaphthalene (DHN) was investigated using two complementary computational approaches. CASPT2//CASSCF minimum energy paths were characterized for reaction channels on the three lowest-energy singlet excited states, describing initial evolution of the spectroscopic bright (ionic) state and its subsequent decay to dark (covalent) states of benzene-like and hexatriene-like character. Although the benzene-like state is unreactive and can radiate, the hexatriene-like state has indirect access to a low-energy conical intersection seam, at which radiationless decay to the ground state and subsequent product formation can take place. An MMVB molecular dynamics simulation was carried out on the reactive hexatriene-like excited state, suggesting that intramolecular vibrational energy redistribution (IVR) controls the radiationless decay and the photoproduct distribution (which is qualitatively reproduced).     

Ultrafast decay of electronically excited singlet cytosine via a pi,pi* to n(O),pi* state switch

Singlet fluorescence lifetimes of adenosine, cytidine, guanosine, and thymidine, determined by femtosecond pump-probe spectroscopy (Pecourt, J.-M. L.; Peon, J.; Kohler, B. J. Am. Chem. Soc. 2000, 122, 9348. Pecourt, J.-M. L.; Peon, J.; Kohler, B. J. Am. Chem. Soc. 2001, 123, 10370), show that the excited states produced by 263 nm light in these nucleosides decay in the subpicosecond range (290-720 fs). Ultrafast radiationless decay to the ground state greatly reduces the probability of photochemical damage. In this work we present a theoretical study of isolated cytosine, the chromophore of cytidine. The experimental lifetime of 720 fs indicates that there must be an ultrafast decay channel for this species. We have documented the possible decay channels and approximate energetics, using a valence-bond derived analysis to rationalize the structural details of the paths. The mechanism favored by our calculations and the experimental data involves (1) a two-mode decay coordinate composed of initial bond length inversion followed by internal vibrational energy redistribution (IVR) to populate a carbon pyramidalization mode, (2) a state switch between the pi,pi* and nO,pi* (excitation from oxygen lone pair) excited states, and (3) decay to the ground state through a conical intersention. A second decay path through the nN,pi* state (excitation from the nitrogen lone pair), with a higher barrier, involves out-of-plane bending of the amino substituent.     

Direct dynamics study of ultrafast vibrational energy relaxation in ice Ih

Car-Parrinello molecular dynamics (CPMD) and a previously developed wave packet model are used to study ultrafast relaxation in water clusters. Water clusters of 15 water molecules are used to represent ice Ih. The relaxation is studied by exciting a symmetric or an asymmetric stretch mode of the central water molecule. The CPMD results suggest that relaxation occurs within 100 fs. This is in agreement with experimental work by Woutersen and Bakker and the earlier wave packet calculations. The CPMD results further indicate that the excitation energy is transferred both intramolecularly and intermolecularly on roughly the same time scale. The intramolecular energy transfer occurs predominantly between the symmetric and asymmetric modes while the bend mode is largely left unexcited on the short time scale studied here.     

Vibrationally controlled chemistry: mode- and bond-selected reaction of CH3D with Cl

Selective vibrational excitation controls the competition between C-H and C-D bond cleavage in the reaction of CH(3)D with Cl, which forms either HCl + CH(2)D or DCl + CH(3). The reaction of CH(3)D molecules with the first overtone of the C-D stretch (2nu(2)) excited selectively breaks the C-D bond, producing CH(3) exclusively. In contrast, excitation of either the symmetric C-H stretch (nu(1)), the antisymmetric C-H stretch (nu(4)), or a combination of antisymmetric stretch and CH(3) umbrella bend (nu(4) + nu(3)) causes the reaction to cleave only a C-H bond to produce solely CH(2)D. Initial preparation of C-H stretching vibrations with different couplings to the reaction coordinate changes the rate of the H-atom abstraction reaction. Excitation of the symmetric C-H stretch (nu(1)) of CH(3)D accelerates the H-atom abstraction reaction 7 times more than excitation of the antisymmetric C-H stretch (nu(4)) even though the two lie within 80 cm(-1) of the same energy. Ab initio calculations and a simple theoretical model help identify the dynamics behind the observed mode selectivity.     

Mechanistic Origin of the Vibrational Coherence Accompanying the Photoreaction of Biomimetic Molecular Switches

The coherent photoisomerization of a chromophore in condensed phase is a rare process in which light energy is funneled into specific molecular vibrations during electronic relaxation from the excited to the ground state. In this work, we employed ultrafast spectroscopy and computational methods to investigate the molecular origin of the coherent motion accompanying the photoisomerization of indanylidene–pyrroline (IP) molecular switches. UV/Vis femtosecond transient absorption gave evidence for an excited‐ and ground‐state vibrational wave packet, which appears as a general feature of the IP compounds investigated. In close resemblance to the coherent photoisomerization of rhodopsin, the sudden onset of a far‐red‐detuned and rapidly blue‐shifting photoproduct signature indicated that the population arriving on the electronic ground state after nonadiabatic decay through the conical intersection (CI) is still very focused in the form of a vibrational wave packet. Semiclassical trajectories were employed to investigate the reaction mechanism. Their analysis showed that coupled double‐bond twisting and ring inversions, already populated during the excited‐state reactive motion, induced periodic changes in π‐conjugation that modulate the ground‐state absorption after the non‐adiabatic decay. This prediction further supports that the observed ground‐state oscillation results from the reactive motion, which is in line with a biomimetic, coherent photoisomerization scenario. The IP compounds thus appear as a model system to investigate the mechanism of mode‐selective photomechanical energy transduction. The presented mechanism opens new perspectives for energy transduction at the molecular level, with applications to the design of efficient molecular devices.     

Coherent low-frequency motions of hydrogen bonded acetic acid dimers in the liquid phase

Ultrafast vibrational dynamics of cyclic hydrogen bonded dimers and the underlying microscopic interactions are studied in temporally and spectrally resolved pump-probe experiments with 100 fs time resolution. Femtosecond excitation of the O-H and/or O-D stretching mode gives rise to pronounced changes of the O-H/O-D stretching absorption displaying both rate-like kinetic and oscillatory components. A lifetime of 200 fs is measured for the v=1 state of the O-H stretching oscillator. The strong oscillatory absorption changes are due to impulsively driven coherent wave packet motions along several low-frequency modes of the dimer between 50 and 170 cm(-1). Such wave packets generated via coherent excitation of the high-frequency O-H/O-D stretching oscillators represent a clear manifestation of the anharmonic coupling of low- and high-frequency modes. The underdamped low-frequency motions dephase on a time scale of 1-2 ps. Calculations of the vibrational potential energy surface based on density functional theory give the frequencies, anharmonic couplings, and microscopic elongations of the low-frequency modes, among them intermolecular hydrogen bond vibrations. Oscillations due to the excitonic coupling between the two O-H or O-D stretching oscillators are absent as is independently confirmed by experiments on mixed dimers with uncoupled O-H and O-D stretching oscillators.     

Picosecond IR-UV pump-probe spectroscopic study on the vibrational energy flow in isolated molecules and clusters

Intramolecular vibrational energy redistribution (IVR) and vibrational predissociation (VP) from the XH stretching vibrations, where X refers to O or C atom, of aromatic molecules and their hydrogen(H)-bonded clusters are investigated by picosecond time-resolved IR-UV pump probe spectroscopy in a supersonic beam. For bare molecules, we mainly focus on IVR of the OH stretch of phenol. We describe the IVR of the OH stretch by a two-step tier model and examine the effect of the anharmonic coupling strength and the density of states on IVR rate and mechanism by using isotope substitution. In the H-bonded clusters of phenol, we show that the relaxation of the OH stretching vibration can be described by a stepwise process and then discuss which process is sensitive to the H-bonding strength. We discuss the difference/similarity of IVR/VP between the "donor" and the "acceptor" sites in phenol-ethylene cluster by exciting the CH stretch vibrations. Finally, we study the vibrational energy transfer in the isolated molecules having the alkyl chain, namely phenylalcanol (PA). In this system, we measure the rate constant of the vibrational energy transfer between the OH stretch and the vibrations of benzene ring which are connected at the both ends of the alkyl chain. This energy transfer can be called "through-bond IVR". We investigate the three factors which are thought to control the energy transfer rate; (1) "OH <--> next CH(2)" coupling, (2) chain length and (3) conformation. We discuss the energy transfer mechanism in PAs by examining these factors.     

Photoionization-induced dynamics of ammonia: ab initio potential energy surfaces and time-dependent wave packet calculations for the ammonia cation

An analytical anharmonic six-dimensional three-sheeted potential energy surface for the ground and first excited states of the ammonia cation has been developed which is tailored to model the ultrafast photoinduced dynamics. Selected ab initio cuts, obtained by multireference configuration interaction calculations, have been used to determine the parameters of a diabatic representation for this Jahn-Teller and pseudo-Jahn-Teller system. The model includes higher-order coupling terms both for the Jahn-Teller and for the pseudo-Jahn-Teller matrix elements. The relaxation to the ground state is possible via dynamical pseudo-Jahn-Teller couplings involving the asymmetric bending and stretching coordinates. The photoelectron spectrum of NH3 and the internal conversion dynamics of NH3+ have been determined by wave packet propagation calculations employing the multiconfigurational time-dependent Hartree method. Three different time scales are found in the dynamics calculations for the second absorption band. The ultrafast Jahn-Teller dynamics of the two excited states occurs on a 5 fs time scale. The major part of the internal conversion to the ground state takes place within a short time scale of 20 fs. This fast internal conversion is, however, incomplete and the remaining excited state population does not decay completely even within 100 fs.     

Excited state intramolecular proton transfer in Schiff bases. Decay of the locally excited enol state observed by femtosecond resolved fluorescence

Although the late (t>1 ps) photoisomerization steps in Schiff bases have been described in good detail, some aspects of the ultrafast (sub-100 fs) proton transfer process, including the possible existence of an energy barrier, still require experimental assessment. In this contribution we present femtosecond fluorescence up-conversion studies to characterize the excited state enol to cis-keto tautomerization through measurements of the transient molecular emission. Salicylideneaniline and salicylidene-1-naphthylamine were examined in acetonitrile solutions. We have resolved sub-100 fs and sub-0.5 ps emission components which are attributed to the decay of the locally excited enol form and to vibrationally excited states as they transit to the relaxed cis-keto species in the first electronically excited state. From the early spectral evolution, the lack of a deuterium isotope effect, and the kinetics measured with different amounts of excess vibrational energy, it is concluded that the intramolecular proton transfer in the S1 surface occurs as a barrierless process where the initial wave packet evolves in a repulsive potential toward the cis-keto form in a time scale of about 50 fs. The absence of an energy barrier suggests the participation of normal modes which modulate the donor to acceptor distance, thus reducing the potential energy during the intramolecular proton transfer.     

Multiple steps and multiple excitations in photoisomerization of azobenzene

Detailed understanding of femtosecond-scale photochemical processes requires dynamical simulations that are complementary to interpretations based on transitions between energy surfaces. For cis to trans photoisomerization of azobenzene following a 100 fs laser pulse, we find that the mechanism is rotation about the central NN bond, and the process is complete in less than 1 ps. The initial excitation and subsequent de-excitation are each achieved via multiple steps, with the molecule always in a superposition of electronic states, as the 3N nuclear degrees of freedom are excited by the laser pulse.     

Time-dependent treatment of intramolecular energy transfer for HeICl complex

The time-dependent golden wave packet method has been used for calculating the decay widths of vibrational predissociation for HeICl complex in the B state with total angular momentum J=0. This is a good example of intramolecular energy transfer. We examine the dependence of the final rotational distribution (partial decay width) of ICl fragment on the stretching excitation. It is found that computed final rotational distributions are weakly dependent on the vibrational level being excited. Unlike the smoothly varying rotational distribution for lower initial vibrational levels, for higher initial vibrational levels the rotational distribution indicates the very pronounced oscillatory structure. The analysis of the rotational distribution as a function of propagation time reveals the predominant role of the final states interaction in determining the final rotational distribution.     

Vibrational overtone initiated unimolecular dissociation of HOCH2OOH and HOCD2OOH: evidence for mode selective behavior

The vibrational overtone induced unimolecular dissociation of HMHP (HOCH(2)OOH) and HMHP-d(2) (HOCD(2)OOH) into OH and HOCH(2)O (HOCD(2)O) fragments is investigated in the region of the 4nu(OH) and 5nu(OH) bands. The unimolecular dissociation rates in the threshold region, corresponding to the 4nu(OH) band, exhibit measurable differences associated with excitation of the OH stretch of the alcohol versus the peroxide functional group, with the higher energy alcohol OH stretching state exhibiting a slower dissociation rate compared to the lower energy peroxide OH stretch in both HMHP and HMHP-d(2). Predictions using the Rice-Ramsperger-Kassel-Marcus theory give rates that are in reasonably good agreement with the measured dissociation rate for the alcohol OH stretch but considerably differ from the measured rates for the peroxide OH stretch in both isotopomers. The present results are interpreted as suggesting that the extent of intramolecular vibrational energy redistribution (IVR) is different for the two OH stretching states associated with the two functional groups in HMHP, with IVR being substantially less complete for the peroxide OH stretch. Analysis of the OH fragment product state distributions in conjunction with phase-space theory simulation gives a D(0) value of 38+/-0.7 kcal/mole for breaking the peroxide bond in HMHP.     

Intramolecular vibrational energy redistribution in bridged azulene-anthracene compounds: ballistic energy transport through molecular chains

Intramolecular vibrational energy flow in excited bridged azulene-anthracene compounds is investigated by time-resolved pump-probe laser spectroscopy. The bridges consist of molecular chains and are of the type (CH(2))(m) with m up to 6 as well as (CH(2)OCH(2))(n) (n=1,2) and CH(2)SCH(2). After light absorption into the azulene S(1) band and subsequent fast internal conversion, excited molecules are formed where the vibrational energy is localized at the azulene side. The vibrational energy transfer through the molecular bridge to the anthracene side and, finally, to the surrounding medium is followed by probing the red edge of the azulene S(3) absorption band at 300 nm and/or the anthracene S(1) absorption band at 400 nm. In order to separate the time scales for intramolecular and intermolecular energy transfer, most of the experiments were performed in supercritical xenon where vibrational energy transfer to the bath is comparably slow. The intramolecular equilibration proceeds in two steps. About 15%-20% of the excitation energy leaves the azulene side within a short period of 300 fs. This component accompanies the intramolecular vibrational energy redistribution (IVR) within the azulene chromophore and it is caused by dephasing of normal modes contributing to the initial local excitation of the azulene side and extending over large parts of the molecule. Later, IVR in the whole molecule takes place transferring vibrational energy from the azulene through the bridge to the anthracene side and thereby leading to microcanonical equilibrium. The corresponding time constants tau(IVR) for short bridges increase with the chain length. For longer bridges consisting of more than three elements, however, tau(IVR) is constant at around 4-5 ps. Comparison with molecular dynamics simulations suggests that the coupling of these chains to the two chromophores limits the rate of intramolecular vibrational energy transfer. Inside the bridges the energy transport is essentially ballistic and, therefore, tau(IVR) is independent on the length.     

Libration induced stretching mode excitation for pump-probe spectroscopy in pure liquid water

We developed an experimental approach to study pure liquid water in the infrared and avoid thermal effects. This technique is based on libration induced stretching excitation of water molecules. A direct correspondence between frequencies within the libration and OH stretching bands is demonstrated. Energy diffusion is studied in pure liquid water by measuring wave packet dynamics of OH stretching vibrator with infrared femtosecond spectroscopy. Wave packet dynamics reveals ultrafast energy dynamics and reflects 130 fs intermolecular energy transfer between water vibrators. Energy diffusion is almost two orders of magnitude faster than self diffusion in water.     

Vibrational relaxation of CN stretch of pseudo-halide anions (OCN-, SCN-, and SeCN-) in polar solvents

The vibrational relaxation dynamics of pseudo-halide anions XCN- (X = O, S, Se) in polar solvents were studied to understand the effect of charge on solute-to-solvent intermolecular energy transfer (IET) and solvent assisted intramolecular vibrational relaxation (IVR) pathways. The T1 relaxation times of the CN stretch in these anions were measured by IR pump/IR probe spectroscopy, in which the 0-1 transition was excited, and the 0-1 and 1-2 transitions were monitored to follow the recovery of the ground state and decay of the excited state. For these anions in five solvents, H2O, D2O, CH3OH, CH3CN, and (CH3)2SO, relaxation rates followed the trend of OCN- > SCN- > SeCN-. For these anions and isotopes of SCN-, the relaxation rate was a factor of a few (2.5-10) higher in H2O than in D2O. To further probe the solvent isotope effect, the relaxation rates of S12C14N-, S13C14N-, and S12C15N- in deuterated methanols (CH3OH, CH3OD, CH3OH, CD3OD) were compared. Relaxation rate was found to be affected by the change of solvent vibrational band at the CN- stretching mode (CD3 symmetric stretch) and lower frequency regions, suggesting the presence of both direct IET and solvent assisted IVR relaxation pathways. The possible relaxation pathways and mechanisms for the observed trends in solute and solvent dependence were discussed.     

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).