Different chemical dynamics for different conformers of biological molecules: photoionization of glycine

Single-photon ionization dynamics of two conformers of glycine is studied by classical trajectory simulations using the semiempirical PM3 potential surface in "on the fly" calculations. Initial conditions for the trajectories are weighted according to the Wigner distribution function computed for the initial vibrational ground state. Vertical ionization in the spirit of the classical Franck-Condon principle is assumed. The dynamics of the two conformers are compared during the first 10 ps. The comparison shows very different dynamical behavior for the two conformers. In particular, the chemical fragmentation pathways differ in part. Also, one of the conformers gives much higher rates of conformational transitions, while the other conformer gives larger chemical fragmentation yields. The example shows significantly different chemical dynamics for two conformers close in energy and separated by a low barrier.     

Representing and selecting vibrational angular momentum states for quasiclassical trajectory chemical dynamics simulations

Linear molecules with degenerate bending modes have states, which may be represented by the quantum numbers N and L. The former gives the total energy for these modes and the latter identifies their vibrational angular momentum jz. In this work, the classical mechanical analog of the N,L-quantum states is reviewed, and an algorithm is presented for selecting initial conditions for these states in quasiclassical trajectory chemical dynamics simulations. The algorithm is illustrated by choosing initial conditions for the N = 3 and L = 3 and 1 states of CO2. Applications of this algorithm are considered for initial conditions without and with zero-point energy (zpe) included in the vibrational angular momentum states and the C-O stretching modes. The O-atom motions in the x,y-plane are determined for these states from classical trajectories in Cartesian coordinates and are compared with the motion predicted by the normal-mode model. They are only in agreement for the N = L = 3 state without vibrational angular momentum zpe. For the remaining states, the Cartesian O-atom motions are considerably different from the elliptical motion predicted by the normal-mode model. This arises from bend-stretch coupling, including centrifugal distortion, in the Cartesian trajectories, which results in tubular instead of elliptical motion. Including zpe in the C-O stretch modes introduces considerable complexity into the O-atom motions for the vibrational angular momentum states. The short-time O-atom motions for these trajectories are highly irregular and do not appear to have any identifiable characteristics. However, the O-atom motions for trajectories integrated for substantially longer period of times acquire unique properties. With C-O stretch zpe included, the long-time O-atom motion becomes tubular for trajectories integrated to approximately 14 ps for the L = 3 states and to approximately 44 ps for the L = 1 states.     

Vibrational spectroscopy and dynamics in the CH-stretch region of fluorene by IVR-assisted, ionization-gain stimulated Raman spectroscopy

We report stimulated Raman spectra at 0.2 and 0.03 cm(-1) resolution in the CH-stretching region of jet-cooled fluorene. The results were obtained by a version of ionization-gain stimulated Raman spectroscopy in which resonant two-photon ionization probing of the state-population changes arising from stimulated Raman transitions is assisted by the process of intramolecular vibrational redistribution (IVR) in the Raman-excited molecule. The fluorene spectra reveal extensive vibrational coupling interactions involving both the aliphatic and aromatic CH-stretching first excited states with nearby background states. Results pertaining to the symmetric aliphatic CH-stretching fundamental are consistent with a tier model of IVR and point to vibrational energy flow out of the CH stretch on a approximately 1 ps time scale with subsequent redistribution on a approximately 5 ps time scale.     

Excited state dynamics and fragmentation channels of the protonated dipeptide H2N-Leu-Trp-COOH

The excited state dynamics of the isolated and protonated peptide H(2)N-Leu-Trp-COOH are analyzed by fs pump-probe spectroscopy. The peptides are brought into the gas phase by electrospray ionization, and fs pump-probe excitation is detected by fragment ion formation. The pump laser addressed the excited pipi* state of the indole chromophore of the amino acid tryptophan. The subsequent excited state dynamics agreed with a biexponential decay with time constants of 500 fs and 10 ps. This is considerably shorter than the lifetime of neutral tryptophan in solution and in proteins, but similar to isolated, protonated tryptophan. Several models are discussed to explain the experimental results but the detailed quenching mechanism remains unresolved.     

Femtosecond dynamics of Cu(CD3OD)

We report the femtosecond nuclear dynamics of Cu(CD3OD) van der Waals clusters, investigated using photodetachment-photoionization spectroscopy. Photodetachment of an electron from Cu-(CD3OD) with a 150 fs, 398 nm laser pulse produces a vibrationally excited neutral complex that undergoes ligand reorientation and dissociation. The dynamics of Cu(CD3OD) on the neutral surface is interrogated by delayed femtosecond resonant two-photon ionization. Analysis of the resulting time-dependent signals indicates that the nascent Cu(CD3OD) complex dissociates on two distinct time scales of 3 and 30 ps. To understand the origins of the observed time scales, complimentary studies were performed. These included measurement of the photoelectron spectrum of Cu-(CD3OD) as well as a series of calculations of the structure and the electronic and vibrational energies of the anion and neutral complexes. Based on the comparisons of the experimental and calculated results for Cu(CD3OD) with those obtained from earlier studies of Cu(H2O), we conclude that the 3 ps time scale reflects the energy transfer from the rotation of CD3OD in the complex to the dissociation coordinate, while the 30 ps time scale reflects the energy transfer from the excited methyl torsion states to the dissociation coordinate.     

Femtosecond time-resolved geometry relaxation and ultrafast intramolecular energy redistribution in Ag2Au.

The ultrafast dynamics of the bimetallic cluster Ag2Au is investigated by pump-probe negative ion-to-neutral-to-positive ion (NeNePo) spectroscopy. Preparation of the neutral cluster in a highly nonequilibrium state by electron detachment from the mass-selected anion, and subsequent probing of the neutral nuclear dynamics through two-photon ionization to the cationic state, leads to strongly probe-energy-dependent transient cation-abundance signals. The origin of this pronounced time and wavelength dependence of the ionization probability on the femtosecond scale is revealed by ab initio theoretical simulations of the transient spectra. Based on the analysis of underlying dynamics, two fundamental processes involving geometry relaxation from linear to triangular structure followed by ultrafast intramolecular vibrational energy redistribution (IVR) have been identified and for the first time experimentally observed in the frame of NeNePo spectroscopy under conditions close to zero electron kinetic energy.     

吡啶气相分子的共振多光子电离谱

本工作观测了双光子能量在35180-36002cm^-^1范围内吡啶分子的共振多光子电离谱。对电子跃迁S~1(^1B~1)→S~0(^1A~l)的振动带做了归属, 首次在^1B~1态观察到7个新的振动模。讨论了N原子取代对分子结构及振动模的影响。     

Conformational transitions of glycine induced by vibrational excitation of the O-H stretch

Vibrational energy flow and conformational transitions following excitation of the OH stretching mode of the most stable conformer of glycine are studied by classical trajectories. "On the fly" simulations with the PM3 semiempirical electronic structure method for the potential surface are used. Initial conditions are selected to correspond to the ν=1 excitation of the OH stretch. The main findings are: (1) An an equilibrium-like ratio is established between the populations of the 3 lowest-lying conformers after about 10 picoseconds. (2) There is a high probability throughout the 150 ps of the simulations for finding the molecule in geometries far from the equilibrium structures of the lowest-energy conformers. (3) Energy from the initial excited OH (ν=1) stretch flows preferentially to 5 other vibrational modes, including the bending motion of the H atom. (4) RRK theory yields conformational transition rates that deviate substantially from the classical trajectory results. Possible implication of these results for vibrational energy flow and conformational transitions in small biological molecules are discussed.     

Post-transition state dynamics for propene ozonolysis: Intramolecular and unimolecular dynamics of molozonide

A direct chemical dynamics simulation, at the B3LYP6-31G(d) level of theory, was used to study the post-transition state intramolecular and unimolecular dynamics for the O3 + propene reaction. Comparisons of B3LYP6-31G(d) with CCSD(T)/cc-pVTZ and other levels of theory show that the former gives accurate structures and energies for the reaction's stationary points. The direct dynamics simulations are initiated at the anti and syn O3 + propene transition states (TSs) and the TS symmetries are preserved in forming the molozonide intermediates. Anti<-->syn molozonide isomerization has a very low barrier of 2-3 kcalmol and its Rice-Ramsperger-Kassel-Marcus (RRKM) lifetime is 0.3 ps. However, the trajectory isomerization is slower and it is unclear whether this anti<-->syn equilibration is complete when the trajectories are terminated at 1.6 ps. The syn (anti) molozonides dissociate to CH3CHO + H2COO and H2CO + syn (anti) CH3CHOO. The kinetics for the latter reactions are in overall good agreement with RRKM theory, but there is a symmetry preserving non-RRKM dynamical constraint for the former. Dissociation of anti molozonide to CH3CHO + H2COO is enhanced and suppressed, respectively, for the trajectory ensembles initiated at the anti and syn O3 + propene TSs. The dissociation of syn molozonide to CH3CHO + H2COO may also be enhanced for trajectories initiated at the syn O3 + propene TS. At the time the trajectories are terminated at 1.6 ps, the ratio of the trajectory and RRKM values of the CH3CHO + H2COO product yield is 1.6 if the symmetries of the initiation and dissociation TSs are the same and 0.6 if their symmetries are different. There are coherences in the intramolecular energy flow, which depend on molozonide's symmetry (i.e., anti or syn). This symmetry related dynamics is not completely understood, but it is clearly related to the non-RRKM dynamics for anti<-->syn isomerization and anti molozonide dissociation to CH3CHO + H2COO. Correlations are found between the stretching motions of molozonide, indicative of nonchaotic and non-RRKM dynamics. The non-RRKM dynamics of molozonide dissociation partitions vibration energy to H2COO that is larger than statistical partitioning. Though the direct dynamics simulations are classical, better agreement is obtained using quantum instead of classical harmonic RRKM theory. This may result from the neglect of anharmonicity in the RRKM calculations, the non-RRKM dynamics of the classical trajectories, or a combination of these two effects. The trajectories suggest that the equilibrium syn/anti molozonide ratio is approximately 1.1-1.2 times larger than that predicted by the harmonic densities of state, indicating an anharmonic correction.     

DNA excited-state dynamics: ultrafast internal conversion and vibrational cooling in a series of nucleosides

To better understand DNA photodamage, several nucleosides were studied by femtosecond transient absorption spectroscopy. A 263-nm, 150-fs ultraviolet pump pulse excited each nucleoside in aqueous solution, and the subsequent dynamics were followed by transient absorption of a femtosecond continuum pulse at wavelengths between 270 and 700 nm. A transient absorption band with maximum amplitude near 600 nm was detected in protonated guanosine at pH 2. This band decayed in 191 +/- 4 ps in excellent agreement with the known fluorescence lifetime, indicating that it arises from absorption by the lowest excited singlet state. Excited state absorption for guanosine and the other nucleosides at pH 7 was observed in the same spectral region, but decayed on a subpicosecond time scale by internal conversion to the electronic ground state. The cross section for excited state absorption is very weak for all nucleosides studied, making some amount of two-photon ionization of the solvent unavoidable. The excited state lifetimes of Ado, Guo, Cyd, and Thd were determined to be 290, 460, 720, and 540 fs, respectively (uncertainties are +/-40 fs). The decay times are shorter for the purines than for the pyrimidine bases, consistent with their lower propensity for photochemical damage. Following internal conversion, vibrationally highly excited ground state molecules were detected in experiments on Ado and Cyd by hot ground state absorption at ultraviolet wavelengths. The decays are assigned to intermolecular vibrational energy transfer to the solvent. The longest time constant observed for Ado is approximately 2 ps, and we propose that solute-solvent H-bonds are responsible for this fast rate of vibrational cooling. The results show for the first time that excited singlet state dynamics of the DNA bases can be directly studied at room temperature. Like sunscreens that function by light absorption, the bases rapidly convert dangerous electronic energy into heat, and this property is likely to have played a critical role in life's early evolution on earth.     

Ultrafast dynamics of acetylacetone (2,4-pentanedione) in the S2 state

The dynamics of the enolic form of acetylacetone (E-AcAc) was investigated using a femtosecond pump-probe experiment. The pump at 266 nm excited E-AcAc in the first bright state, S2(pi pi*). The resulting dynamics was probed by multiphoton ionization at 800 nm. It was investigated for 80 ps on the S2(pi pi*) and S1(n pi*) potential energy surfaces. An important step is the transfer from S2 to S1 that occurs with a time constant of 1.4 +/- 0.2 ps. Before, the system had left the excitation region in 70 +/- 10 fs. An intermediate step was identified when E-AcAc traveled on the S2 surface. Likely, it corresponds to an accidental resonance in the detection scheme that is met along this path. More importantly, some clues are given that an intramolecular vibrational energy relaxation is observed, which transfers excess vibrational energy from the enolic group O-H to the other modes of the molecule. The present multistep evolution of excited E-AcAc probably also describes, at least qualitatively, the dynamics of other electronically excited beta-diketones.     

Ultrafast dynamics in thiophene investigated by femtosecond pump probe photoelectron spectroscopy and theory

A hybrid of a time-of-flight mass spectrometer and a time-of-flight "magnetic-bottle type" photoelectron (PE) spectrometer is used for fs pump-probe investigations of the excited state dynamics of thiophene. A resonant two-photon ionization spectrum of the onset of the excited states has been recorded with a tunable UV laser of 190 fs pulse width. With the pump laser set to the first intense transition we find by UV probe ionization first a small time shift of the maxima in the PE spectrum and then a fast decay to a low constant intensity level. The fitted time constants are 80+/-10 fs, and 25+/-10 fs, respectively. Theoretical calculations show that upon geometry relaxation the electronic state order changes and conical intersections between excited states exist. We use the vertical state order S1, S2, S3 to define the terms S1, S2, and S3 for the characterization of the electron configuration of these states. On the basis of our theoretical result we discuss the electronic state order in the UV spectra and identify in the photoelectron spectrum the origin of the first cation excited state D1. The fast excited state dynamics agrees best with a vibrational dynamics in the photo-excited S1 (80+/-10 fs) and an ultrafast decay via a conical intersection, presumably a ring opening to the S3 state (25+/-10 fs). The subsequently observed weak constant signal is taken as an indication, that in the gas phase the ring-closure to S0 is slower than 50 ps. An ultrafast equilibrium between S1 and S2 before ring opening is not supported by our data.     

Time scales and pathways of vibrational energy relaxation in liquid CHBr3 and CDBr3

Molecular dynamics simulations are used in conjunction with Landau-Teller, fluctuating Landau-Teller, and time-dependent perturbation theories to investigate energy flow out of various vibrational states of liquid CHBr3 and CDBr3. The CH stretch overtone is found to relax with a time scale of about 1 ps compared to the 50 ps rate for the fundamental. The relaxation pathways and rates for the CD stretch decay in CDBr3 are computed in order to understand the changes arising from deuteration. While the computed relaxation rate agrees well with experiments, the pathway is found to be more complex than anticipated. In addition to the above channels for CH(D) stretch relaxation that involve only the hindered translations and rotations of the solvent, routes involving off-resonant and resonant excitations of solvent vibrational modes are also examined. Finally, the decay of energy from low frequency states to near-lying solute states and solvent vibrations are studied.     

Ultrafast excited state dynamics of modified phthalocyanines: p-HPcZn and p-HPcCo

Two modified metallophthalocyanines (MPcs) containing sulfonic naphthoxy substituents were synthesized. The measurements of transient absorption and time-resolved photoluminescence were used to study the ultrafast response and excited state dynamics of two MPcs in dimethyl sulfoxide (DMSO) solution, which were predominantly in the monomeric form. Under excitation at 400 nm, these molecules experience vibrational relaxation to the bottom of the first excited state and then the excitation rapidly converts to the low-lying charge-transfer (CT) state and finally reaches the triplet states. Under excitation at 800 nm, they show a two-photon absorption character, and their excited state dynamics exhibit strong dependence on the probe wavelength. The main results with 400 nm pumping are similar to the results with 800 nm pumping. For p-HPcZn, weak two-photon photoluminescence was also observed with a lifetime of 52 +/- 2 ps. A four-level model was used to illustrate the excited state dynamics of p-HPcZn, while a five-level model was suggested for p-HPcCo molecule.     

Experimental and Theoretical Investigation of the 3sp(d) Rydberg States of Fenchone by Polarized Laser Resonance-Enhanced-Multiphoton-Ionization and Fourier Transform VUV Absorption Spectroscopy

The VUV absorption spectrum of fenchone is re-examined using synchrotron radiation Fourier transform spectrometry, revealing new vibrational structure. Picosecond laser (2+1) resonance enhanced multiphoton ionization (REMPI) spectroscopy complements this, providing an alternative view of the 3spd Rydberg excitation region. These spectra display broadly similar appearance, with minor differences that are largely explained by referring to calculated one- and two-photon electronic excitation cross-sections. Both show good agreement with Franck-Condon simulations of the relevant vibrational structures. Parent ion REMPI ionization yields with both femtosecond and picosecond excitation laser pulses are studied as a function of laser polarization and intensity, the latter providing insight into the relative two-photon excitation and one-photon ionization rates. The experimental circular-linear dichroism observed in the parent ion yields varies strongly between the 3s and 3p Rydberg states, in good overall agreement with the calculated two-photon excitation circular-linear dichroism, while corroborating other evidence that the 3p_z sub-state plays no more than a very minor role in the (2+1) REMPI spectrum. Vibrationally resolved photoelectron spectra are recorded with picosecond pulse duration (2+1) REMPI at selected intermediate vibrational excitations. The 3s intermediate state displays a very strong Δv=0 propensity on ionization, but the 3p intermediate evidences more complex vibronic dynamics, and we infer some 3p→3s internal conversion prior to ionization.     

A molecular dynamics study of Lys-Trp-Lys: structure and dynamics in solution following photoexcitation

We report studies of the structure and dynamics of a tripeptide Lys-Trp-Lys (KWK) in aqueous solution following photoexcitation by molecular dynamics simulations. For ground-state KWK, we observe three stable conformations with free energy differences of less than 5.2 kJ/mol. Each conformer is stabilized by a pi-cation interaction between one of three protonated amino groups and the indole moiety. For the excited state of tryptophan in KWK, the simulated molecular dynamics of the three isomers are similar, all in good agreement with recent femtosecond experiments (J. Phys. Chem. B 2005, 109, 16901). Specifically, we observe: (1) the fluorescence anisotropy is dominated by a single-exponential component and decays in approximately 130 ps, (2) the total dynamic Stokes shift reaches approximately 2700 cm(-1), and (3) the excited state relaxation dynamics occurs on several time scales ranging from femtoseconds to tens of picoseconds. The relaxation dynamics involve rapid initial response of neighboring water, followed by local motions of flexible peptide chains. These processes drive global restructuring of the tripeptide on a rather flat energy surface, inducing slower dynamics evident in both the water and protein contributions to the stabilization energy of the photoexcited chromophore. The water and protein dynamics are strongly correlated. On a still longer time scale, we observe isomerization of two excited state conformers to the other most stable one, an analogue for evolution of trajectories along the funnel on the rugged free energy landscape to the final "native" state. Our studies suggest new experiments to detect this unique dynamics.     

Structure and vibrational dynamics of model compounds of the [FeFe]-hydrogenase enzyme system via ultrafast two-dimensional infrared spectroscopy

Ultrafast two-dimensional infrared (2D) spectroscopy has been applied to study the structure and vibrational dynamics of (mu-S(CH2)3S)Fe2(CO)6, a model compound of the active site of the [FeFe]-hydrogenase enzyme system. Comparison of 2D-IR spectra of (mu-S(CH2)3S)Fe2(CO)6 with density functional theory calculations has determined that the solution-phase structure of this molecule is similar to that observed in the crystalline phase and in good agreement with gas-phase simulations. In addition, vibrational coupling and rapid (<5 ps) solvent-mediated equilibration of energy between vibrationally excited states of the carbonyl ligands of the di-iron-based active site model are observed prior to slower (approximately 100 ps) relaxation to the ground state. These dynamics are shown to be solvent-dependent and form a basis for the future determination of the vibrational interactions between active site and protein.     

Ion-Velocity Map Imaging Study of Photodissociation Dynamics of Acetaldehyde

The photodissociation dynamics of acetaldehyde in the radical channel CH3＋HCO has been reinvestigated using time-sliced velocity map imaging technique in the photolysis wavelength range of 275-321 nm. The CH3 fragments have been probed via （2＋1） resonance-enhanced multiphoton ionization. Images are measured for CH3 formed in the ground and excited states （v2=0 and 1） of the umbrella vibrational mode. For acetaldehyde dissociation on T1 state after intersystem crossing from S1 state, the products are formed with high translational energy release and low internal excitation. The rotational and vibrational energy of both fragments increases with increasing photodissociation energy. The triplet barrier height is estimated at 3.8814-0.006 eV above the ground state of acetaldehyde.     

Long live vinylidene! A new view of the H(2)C=C: --> HC triple bond CH rearrangement from ab initio molecular dynamics.

We present complete active space self-consistent field (CASSCF) ab initio molecular dynamics (AIMD) simulations of the preparation of the metastable species vinylidene, and its subsequent, highly exothermic isomerization to acetylene, via electron removal from vinylidene anion (D(2)C=C(-) --> D(2)C=C: --> DC triple bond CD). After equilibrating vinylidene anion-d(2) at either 600 +/- 300 K (slightly below the isomerization barrier) or 1440 K +/- 720 K (just above the isomerization barrier), we remove an electron to form a vibrationally excited singlet vinylidene-d(2) and follow its dynamical evolution for 1.0 ps. Remarkably, we find that none of the vinylidenes equilibrated at 600 K and only 20% of the vinylidenes equilibrated at 1440 K isomerized, suggesting average lifetimes >1 ps for vibrationally excited vinylidene-d(2). Since the anion and neutral vinylidene are structurally similar, and yet extremely different geometrically from the isomerization transition state (TS), neutral vinylidene is not formed near the TS so that it must live until it has sufficient instantaneous kinetic energy in the correct vibrational mode(s). The origin of the delay is explained via both orbital rearrangement and intramolecular vibrational energy redistribution (IVR) effects. Unique signatures of the isomerization dynamics are revealed in the anharmonic vibrational frequencies extracted from the AIMD, which should be observable by ultrafast vibrational spectroscopy and in fact are consistent with currently available experimental spectra. Most interestingly, of those trajectories that did isomerize, every one of them violated conventional transition-state theory by recrossing back to vinylidene multiple times, against conventional notions that expect highly exothermic reactions to be irreversible. The dynamical motion responsible for the multiple barrier recrossings involves strong mode-coupling between the vinylidene CD(2) rock and a local acetylene DCC bend mode that has been recently observed experimentally. The multiple barrier recrossings can be used, via a generalized definition of lifetime, to reconcile extremely disparate experimental estimates of vinylidene's lifetime (differing by at least 6 orders of magnitude). Last, a caveat: These results are constrained by the approximations inherent in the simulation (classical nuclear motion, neglect of rotation-vibration coupling, and restriction to C(s) symmetry); refinement of these predictions may be necessary when more exact simulations someday become feasible.     

The laser two-photon photolysis of liquid carbon tetrachloride

The two-photon photolysis of liquid CCl4 with 25 ps pulses of 266 nm light has been studied and compared with similar studies with high energy radiation. Both neutral and ionic species are produced from excited states and ionization. The emphasis of the study is on the ionic processes, while some data related to excited states and free radicals are presented. In both radiolysis and photolysis, a solvent separated charged pair, CCl3+ // Cl-, exhibiting a lambda(max) at 475 nm, is observed that exhibits a total growth over 38 to 100 ps. Solutes with ionization potentials less than that of CCl4 (11.47 eV) reduce the yield of the 475 nm species producing radical cations of the solute. The efficiency of this process is about 10-fold larger in radiolysis compared with photolysis. Analysis of the data suggest that the lower energy of two-photon photolysis produces a charge pair CCl4+ // CCl4-, which decays in about 3 ps to CCl4+ // Cl-. This species then decays to CCl3+ // Cl-. The lifetime of the growth of the 475 nm is measured as 46 ps. These studies clearly show areas where radiolysis and photolysis can be quite similar and also areas where the vast difference in excitation energy introduces stark differences in the observed radiation and photoinduced chemistry.