Laser pulse trains for controlling excited state dynamics of adenine in water

We investigate theoretically the control of the ultrafast excited state dynamics of adenine in water by laser pulse trains, with the aim to extend the excited state lifetime and to suppress nonradiative relaxation processes. For this purpose, we introduce the combination of our field-induced surface hopping method (FISH) with the quantum mechanical-molecular mechanical (QM/MM) technique for simulating the laser-driven dynamics in the condensed phase under explicit inclusion of the solvent environment. Moreover, we employ parametric pulse shaping in the frequency domain in order to design simplified laser pulse trains allowing to establish a direct link between the pulse parameters and the controlled dynamics. We construct pulse trains which achieve a high excitation efficiency and at the same time keep a high excited state population for a significantly extended time period compared to the uncontrolled dynamics. The control mechanism involves a sequential cycling of the population between the lowest and higher excited states, thereby utilizing the properties of the corresponding potential energy surfaces to avoid conical intersections and thus to suppress the nonradiative decay to the ground state. Our findings provide a means to increase the fluorescence yield of molecules with an intrinsically very short excited state lifetime, which can lead to novel applications of shaped laser fields in the context of biosensing.     

On-the-fly, electric-field-driven, coupled electron-nuclear dynamics

An on-the-fly, electric field driven, coupled electron-nuclear dynamics approach is developed and applied to model the photodissociation of water in the A((1)B1) excited state. In this method, a quantum propagator evolves the photon-induced electronic dynamics in the ultrafast time scale, and a quasi-classical surface hopping approach describes the nuclear dynamics in the slower time scale. In addition, strong system-field interactions are explicitly included in the electronic propagator. This theoretical development enables us to study rapid photon-induced bond dissociation dynamics and demonstrates the partial breakdown of electronic coherence as well as electronic population trapping in the excited state when the molecular vibrations detune the system with respect to the applied field. The method offers a practical way to use on-the-fly dynamics for modeling light-molecule interactions that lead to interesting photochemical events.     

Laser-operated chiral molecular switch: quantum simulations for the controlled transformation between achiral and chiral atropisomers

We report quantum dynamical simulations for the laser controlled isomerization of 1-(2-cis-fluoroethenyl)-2-fluorobenzene based on one-dimensional electronic ground and excited state potentials obtained from (TD)DFT calculations. 1-(2-cis-fluoroethenyl)-2-fluorobenzene supports two chiral and one achiral atropisomers, the latter being the most stable isomer at room temperature. Using a linearly polarized IR laser pulse the molecule is excited to an internal rotation around its chiral axis, i.e. around the C-C single bond between phenyl ring and ethenyl group, changing the molecular chirality. A second linearly polarized laser pulse stops the torsion to prepare the desired enantiomeric form of the molecule. This laser control allows the selective switching between the achiral and either the left- or right-handed form of the molecule. Once the chirality is "switched on" linearly polarized UV laser pulses allow the selective change of the chirality using the electronic excited state as intermediate state.     

A comparison between different semiclassical approximations for optical response functions in nonpolar liquid solutions

The temporal behavior of optical response functions (ORFs) reflects the quantum dynamics of an electronic superposition state, and as such lacks a well-defined classical limit. In this paper, we consider the importance of accounting for the quantum nature of the dynamics when calculating ORFs of different types. To this end, we calculated the ORFs associated with the linear absorption spectrum and the nonlinear two-pulse photon-echo experiment, via the following approaches: (1) the semiclassical forward-backward approach; (2) an approach based on linearizing the path-integral forward-backward action in terms of the difference between the forward and backward paths; (3) an approach based on ground state nuclear dynamics. The calculations were performed on a model that consists of a two-state chromophore solvated in a nonpolar liquid. The different methods were found to yield very similar results for the absorption spectrum and "diagonal" two-pulse photon echo (i.e., the homodyne-detected signal at time t=t(0) after the second pulse, where t(0) is the time interval between the two pulses). The different approximations yielded somewhat different results in the case of the time-integrated photon-echo signal. The reasons for the similarity between the predictions of different approximations are also discussed.     

Analysis of adiabatic passage by light-induced potentials with chirped laser pulses in three- and four-level diatomic systems

This paper describes an investigation into the process of adiabatic passage by light-induced potentials (APLIP), which was previously suggested as a method for employing two strong picosecond laser pulses to transfer the population between two electronic states. We have extended earlier numerical studies in order to assess the feasibility of an experimental implementation of the APLIP concept. APLIP has been modeled in a three-level model system based on Na2 with chirped pulses, using laser parameters available from a typical chirped pulse amplified Ti:sapphire laser. The model showed that the APLIP process remains essentially unchanged for chirped pulses of equal magnitude and the opposite, or equal and positive sign of chirp as compared to the transform-limited case. We also examined the case of additional electronic states by introduction of a fourth state that lies close to the "target," i.e., final, state. The investigation showed that there are circumstances in which a significant fraction of the population gets transferred to this state which will disrupt the APLIP process. However, by switching to this fourth state as the target state in an experiment, good transfer efficiency is recovered. The results of the extension of the original APLIP modeling to chirped pulses and additional electronic states indicate that an APLIP experimental realization should be feasible in Na2.     

New Energy-Based Decoherence Correction Approaches for Trajectory Surface Hopping

Inspired by the branching corrected surface hopping (BCSH) method [J. Xu and L. Wang, J. Chem. Phys. 150 , 164101 (2019)], we present two new decoherence time formulas for trajectory surface hopping. Both the proposed linear and exponential formulas characterize the decoherence time as functions of the energy difference between adiabatic states and correctly capture the decoherence effect due to wave packet reflection as predicted by BCSH. The relevant parameters are trained in a series of 200 diverse models with different initial nuclear momenta, and the exact quantum solutions are utilized as references. As demonstrated in the three standard Tully models, the two new approaches exhibit significantly higher reliability than the widely used counterpart algorithm while holding the appealing efficiency, thus promising for nonadiabatic dynamics simulations of general systems.     

Singlet versus triplet dynamics of β-carotene studied by quantum control spectroscopy

Biomolecules very often present complex energy deactivation networks with overlapping electronic absorption bands, making their study a difficult task. This can be especially true in transient absorption spectroscopy when signals from bleach, excited state absorption and stimulated emission contribute to the signal. However, quantum control spectroscopy can be used to discriminate specific electronic states of interest by applying specifically designed laser pulses. Recently, we have shown the control of energy flow in bacterial light-harvesting using shaped pump pulses in the visible and the selective population of pathways in carotenoids using an additional depletion pulse in the transient absorption technique. Here, we apply a closed-loop optimization approach to β-carotene using a spatial light modulator to decipher the energy flow network after a multiphoton excitation with a shaped ultrashort pulse in the near-IR. After excitation, two overlapping bands were detected and identified as the S₁ state and the first triplet state T₁. Using the transient absorption signal at a specific probe delay as feedback, the triplet signal could be optimized over the singlet contribution.     

Mixed quantum-classical dynamics in the adiabatic representation to simulate molecules driven by strong laser pulses

The dynamics of molecules under strong laser pulses is characterized by large Stark effects that modify and reshape the electronic potentials, known as laser-induced potentials (LIPs). If the time scale of the interaction is slow enough that the nuclear positions can adapt to these externally driven changes, the dynamics proceeds by adiabatic following, where the nuclei gain very little kinetic energy during the process. In this regime we show that the molecular dynamics can be simulated quite accurately by a semiclassical surface-hopping scheme formulated in the adiabatic representation. The nuclear motion is then influenced by the gradients of the laser-modified potentials, and nonadiabatic couplings are seen as transitions between the LIPs. As an example, we simulate the process of adiabatic passage by light induced potentials in Na(2) using the surface-hopping technique both in the diabatic representation based on molecular potentials and in the adiabatic representation based on LIPs, showing how the choice of the representation is crucial in reproducing the results obtained by exact quantum dynamical calculations.     

Theoretical description of femtosecond fluorescence depletion spectrum of molecules in solution

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

Photo-induced desorption of NO from NiO(100): calculation of the four-dimensional potential energy surfaces and systematic wave packet studies

The velocity distributions of the laser-induced desorption of NO molecules from an epitaxially grown film of NiO(100) on Ni(100) have been studied [Mull et al., J. Chem. Phys., 1992, 96, 7108]. A pronounced bimodality of velocity distributions has been found, where the NO molecules desorbing with higher velocities exhibit a coupling to the rotational quantum states J. In this article we present simulations of state resolved velocity distributions on a full ab initio level. As a basis for this quantum mechanical treatment a 4D potential energy surface (PES) was constructed for the electronic ground and a representative excited state, using a NiO5Mg(18+)13 cluster. The PESs of the electronic ground and an excited state were calculated at the CASPT2 and the configuration interaction (CI) level of theory, respectively. Multi-dimensional quantum wave packet simulations on these two surfaces were performed for different sets of degrees of freedom. Our key finding is that at least a 3D wave packet simulation, in which the desorption coordinate Z, polar angle theta and lateral coordinate X are included, is necessary to allow the simulation of experimental velocity distributions. Analysis of the wave packet dynamics demonstrates that essentially the lateral coordinate, which was neglected in previous studies [Klüner et al., Phys. Rev. Lett. 1998, 80, 5208], is responsible for the experimentally observed bimodality. An extensive analysis shows that the bimodality is due to a bifurcation of the wave packet on the excited state PES, where the motion of the molecule parallel to the surface plays a decisive role.     

Design of UV laser pulses for the preparation of matrix isolated homonuclear diatomic molecules in selective vibrational superposition states

The preparation of matrix isolated homonuclear diatomic molecules in a vibrational superposition state c0Phie=1,v=0+cjPhie=1,v=j, with large (|c0|2 approximately 1) plus small contributions (|cj|2<1) of the ground v=0 and specific v=j low excited vibrational eigenstates, respectively, in the electronic ground (e=1) state, and without any net population transfer to electronic excited (e>1) states, is an important challenge; it serves as a prerequisite for coherent spin control. For this purpose, the authors investigate two scenarios of laser pulse control, involving sequential or intrapulse pump- and dump-type transitions via excited vibronic states Phiex,k with a dominant singlet or triplet character. The mechanisms are demonstrated by means of quantum simulations for representative nuclear wave packets on coupled potential energy surfaces, using as an example a one-dimensional model for Cl2 in an Ar matrix. A simple three-state model (including Phi1,0, Phi1,j and Phiex,k) allows illuminating analyses and efficient determinations of the parameters of the laser pulses based on the values of the transition energies and dipole couplings of the transient state which are derived from the absorption spectra.     

Ultrafast control of the internuclear distance with parabolic chirped pulses

Recently, control over the bond length of a diatomic molecule with the use of parabolic chirped pulses was predicted on the basis of numerical calculations [Chang; et al. Phys. Rev. A 2010, 82, 063414]. To achieve the required bond elongation, a laser scheme was proposed that implies population inversion and vibrational trapping in a dissociative state. In this work we identify two regimes where the scheme works, called the strong and the weak adiabatic regimes. We define appropriate parameters to identify the thresholds where the different regimes operate. The strong adiabatic regime is characterized by a quasi-static process that requires longer pulses. The molecule is stabilized at a bond distance and at a time directly controlled by the pulse in a time-symmetrical way. In this work we analyze the degree of control over the period and elongation of the bond as a function of the pulse bandwidth. The weak adiabatic regime implies dynamic deformation of the bond, which allows for larger bond stretch and the use of shorter pulses. The dynamics is anharmonic and not time-symmetrical and the final state is a wave packet in the ground potential. We show how the vibrational energy of the wave packet can be controlled by changing the pulse duration.     

Coupled-channels quantum theory of electronic flux density in electronically adiabatic processes: application to the hydrogen molecule ion

This article presents the results of the first quantum simulations of the electronic flux density (j(e)) by the "coupled-channels" (CC) theory, the fundamentals of which are presented in the previous article [Diestler, D. J. J. Phys. Chem. A 2012, DOI: 10.1021/jp207843z]. The principal advantage of the CC scheme is that it employs exclusively standard methods of quantum chemistry and quantum dynamics within the framework of the Born-Oppenheimer approximation (BOA). The CC theory goes beyond the BOA in that it yields a nonzero j(e) for electronically adiabatic processes, in contradistinction to the BOA itself, which always gives j(e) = 0. The CC is applied to oriented H(2)(+) vibrating in the electronic ground state ((2)Σ(g)(+)), for which the nuclear and electronic flux densities evolve on a common time scale of about 22 fs per vibrational period. The system is chosen as a touchstone for the CC theory, because it is the only one for which highly accurate flux densities have been calculated numerically without invoking the BOA [Barth et al, Chem. Phys. Lett. 2009, 481, 118]. Good agreement between CC and accurate results supports the CC approach, another advantage of which is that it allows a transparent interpretation of the temporal and spatial properties of j(e).     

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.     

Dynamics of electronic states and spin-flip for photodissociation of dihalogens in matrices: experiment and semiclassical surface-hopping and quantum model simulations for F2 and ClF in solid Ar

Three approaches are combined to study the electronic states' dynamics in the photodissociation of F(2) and ClF in solid argon. These include (a) semiclassical surface-hopping simulations of the nonadiabatic processes involved. These simulations are carried out for the F(2) molecule in a slab of 255 argon atoms with periodic boundary conditions at the ends. The full manifold of 36 electronic states relevant to the process is included. (b) The second approach involves quantum mechanical reduced-dimensionality models for the initial processes induced by a pump laser pulse, which involve wavepacket propagation for the preoriented ClF in the frozen argon lattice and incorporate the important electronic states. The focus is on the study of quantum coherence effects. (c) The final approach is femtosecond laser pump-probe experiments for ClF in Ar. The combined results for the different systems shed light on general properties of the nonadiabatic processes involved, including the singlet to triplet and intertriplet transition dynamics. The main findings are (1) that the system remains in the initially excited-state only for a very brief, subpicosecond, time period. Thereafter, most of the population is transferred by nonadiabatic transitions to other states, with different time constants depending on the systems. (2) Another finding is that the dynamics is selective with regard to the electronic quantum numbers, including the Lambda and Omega quantum numbers, and the spin of the states. (3) The semiclassical simulations show that prior to the first "collision" of the photodissociated F atom with an Ar atom, the argon atoms can be held frozen, without affecting the process. This justifies the rigid-lattice reduced-dimensionality quantum model for a brief initial time interval. (4) Finally, degeneracies between triplets and singlets are fairly localized, but intertriplet degeneracies and near degeneracies can span an extensive range. The importance of quantum effects in photochemistry of matrix-isolated molecules is discussed in light of the results.     

Overlapping resonances in the control of intramolecular vibrational redistribution

Coherent control of bound state processes via the interfering overlapping resonance scenario [Christopher et al., J. Chem. Phys. 123, 064313 (2006)] is developed to control intramolecular vibrational redistribution (IVR). The approach is applied to the flow of population between bonds in a model of chaotic OCS vibrational dynamics, showing the ability to significantly alter the extent and rate of IVR by varying quantum interference contributions.