Influence of Solvation Environment on Excited State Avoided Crossings and Photodissociation Dynamics

How Solvation Alters the Thermodynamics of Asymmetric Bond-Breaking: Quantum Simulation of NaK+ in Liquid Tetrahydrofuran

Gas-phase potential energy surfaces (PESs) are often used to provide an intuitive understanding of molecular chemical reactivity. Most chemical reactions, however, take place in solution, and it is unclear whether gas-phase PESs accurately represent chemical processes in solvent environments. In this work we use quantum simulations to investigate the dissociation energetics of NaK+ in liquid tetrahydrofuran (THF) to understand the degree to which solvent interactions alter the gas-phase picture. Using umbrella sampling and thermodynamic integration techniques, we construct condensed-phase free energy surfaces of NaK+ on THF in both the ground and electronic excited states. We find that solvation by THF completely alters the nature of the NaK+ bond by reordering the thermodynamic dissociation products. Reaching the thermodynamic dissociation limit in THF also requires a long-range charge transfer process that has no counterpart in the gas phase. Gas-phase PESs, even with perturbations, cannot adequately describe the reactivity of simple asymmetric molecules in solution.

Ab initio symmetric quasi-classical approach to investigate molecular Tully models

We perform on-the-fly non-adiabatic molecular dynamics simulations using the symmetrical quasi-classical (SQC) approach with the recently suggested molecular Tully models: ethylene and fulvene. We attempt to provide benchmarks of the SQC methods using both the square and triangle windowing schemes as well as the recently proposed electronic zero-point-energy correction scheme (the so-called γ correction). We use the quasi-diabatic propagation scheme to directly interface the diabatic SQC methods with adiabatic electronic structure calculations. Our results showcase the drastic improvement of the accuracy by using the trajectory-adjusted γ-corrections, which outperform the widely used trajectory surface hopping method with decoherence corrections. These calculations provide useful and non-trivial tests to systematically investigate the numerical performance of various diabatic quantum dynamics approaches, going beyond simple diabatic model systems that have been used as the major workhorse in the quantum dynamics field. At the same time, these available benchmark studies will also likely foster the development of new quantum dynamics approaches based on these techniques.

Recent progress in approximate quantum dynamics methods for the study of proton-coupled electron transfer reactions

Proton-coupled electron transfer (PCET) reactions are ubiquitous natural processes at the heart of energy conversion reactions in photosynthesis and respiration, DNA repair, and diverse enzymatic reactions. Theoretical formulation and computational method developments have eyed modeling of thermal and photoinduced PCET for the last three decades. The accumulation of these studies, collected in dozens of reviews, accounts, and perspectives, has firmly established the influence of quantum effects, including non-adiabatic electronic transitions, vibrational relaxation, zero-point energy, and proton tunneling, on the rate and mechanism of PCET reactions. Here, we focus on some recently-developed methods, spanning the last eight years, that can quantitatively capture these effects in the PCET context and provide efficient means for their qualitative description in complex systems. The theoretical background of each method and their accuracy with respect to exact results are discussed and the results of relevant PCET simulations based on each method are presented.

Quasi-Diabatic Scheme for Nonadiabatic On-the-Fly Simulations

We use the quasi-diabatic (QD) propagation scheme to perform on-the-fly nonadiabatic simulations of the photodynamics of ethylene. The QD scheme enables a seamless interface between accurate diabatic-based quantum dynamics approaches and adiabatic electronic structure calculations, explicitly avoiding any efforts to construct global diabatic states or reformulate the diabatic dynamics approach to the adiabatic representation. Using the partial linearized path-integral approach and the symmetrical quasi-classical approach as the diabatic dynamics methods, the QD propagation scheme enables direct nonadiabatic simulation with complete active space self-consistent field on-the-fly electronic structure calculations. The population dynamics obtained from both approaches are in a close agreement with the quantum wavepacket-based method and outperform the widely used trajectory surface-hopping approach. Further analysis of the ethylene photodeactivation pathways demonstrates the correct predictions of competing processes of nonradiative relaxation mechanism through various conical intersections. This work provides the foundation of using accurate diabatic dynamics approaches and on-the-fly adiabatic electronic structure information to perform ab initio nonadiabatic simulation.

Ultrafast dynamics of halogens in rare gas solids

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

Nonadiabatic Trajectory Studies of NaI(H2O)n Photodissociation Dynamics

We have investigated the photodissociation dynamics of NaI(H(2)O)(n) [n = 1-4] clusters using the molecular dynamics with quantum transitions method and a quantum mechanics/molecular mechanics description of NaI(H(2)O)(n), which involves a semiempirical valence-bond approach to describe the NaI electronic structure and classical solvent-solvent and solute-solvent interaction potentials. Our simulation results show that the NaI(H(2)O)(n) excited-state population decay upon reaching the NaI curve-crossing region increases with cluster size due to the stabilization of the ionic branch of the NaI excited state by the surrounding water molecules, and the resulting increase in nonadiabatic transition probability. After reaching the curve-crossing region for the first time, however, the excited-state population decay resembles that of bare NaI because of rapid evaporation of 99% and 95% of the water molecules for NaI(H(2)O) and NaI(H(2)O)(n) [n = 2-4], respectively. This extensive evaporation is due to the reversed NaI polarity in the Franck-Condon region of the NaI first excited state, which causes strong repulsive NaI-H(2)O forces and induces rapid nonstatistical water evaporation, where product water molecules are formed more rotationally than translationally hot. A few water molecules (5% or less) remain transiently or permanently bound to NaI, forming long-lived clusters, when NaI remains predominantly ionic, i.e., remains in the excited state, after reaching the curve-crossing region. To connect simulation results with experiment, we have simulated femtosecond probe signals resulting from two-photon and one-photon excitation to the X and I NaI(+) probe states. In agreement with experimental findings, the probe signals resulting from the two-photon probe scheme, where excitation occurs from the covalent branch of the excited state, decay exponentially over the NaI first excited-state vibrational period, with very little evidence of long-time dynamics. The one-photon probe scheme (not used for experimental cluster studies) is shown to be less sensitive to solvation, in that excitation energies will remain similar over a range of cluster sizes, as the ionic branch of the excited state and the NaI(+) probe states are stabilized to the same extent by the presence of water molecules. The resulting probe signals are also more revealing of the NaI(H(2)O)(n) photodissociation dynamics than the two-photon probe signals, as they may allow monitoring of solvation effects on the NaI nonadiabatic dynamics and of successive evaporation of water molecules. Time-resolved photoelectron spectra provide limited additional information regarding the NaI(H(2)O)(n) photodissociation dynamics. A key consequence of the rapid water evaporation demonstrated here is that experimentally observed signals may arise from the photodissociation of much larger NaI(H(2)O)(n) parent clusters.

Effective chromophore potential, dissipative trajectories, and vibrational energy relaxation: Br2 in Ar matrix

The intramolecular wave packet dynamics on the electronic B (3pi0) potential of Br2 in solid argon is induced and interrogated by femtosecond pump-probe spectroscopy. An effective potential of the chromophore in the solid is derived from the wave packet period for different excitation photon energies. Deep in the potential well, it is consistent with vibrational energies from wavelength-resolved spectra. It extends to higher energies, where the vibrational bands merge to a continuum, and even beyond the dissociation limit, thus quantifying the cage effect of the argon matrix. This advantage of pump-probe spectroscopy is related to a reduced contribution of homogeneous and inhomogeneous line broadenings. The vibrational energy relaxation rates are determined by a variation of the probe window spatial position via the probe quantum energy. A very large energy loss in the first excursion of the wave packet is observed near the dissociation limit. This strong interaction with the argon matrix is directly displayed in an experimental trajectory.

LAND-map, a linearized approach to nonadiabatic dynamics using the mapping formalism

We present a new approach for calculating quantum time correlation functions for systems whose dynamics exhibits relevant nonadiabatic effects. The method involves partial linearization of the full quantum path-integral expression for the time correlation function written in the nonadiabatic mapping Hamiltonian formalism. Our analysis gives an algorithm which is both numerically efficient and accurate as we demonstrate in test calculations on the spin-boson model where we find results in good agreement with exact calculations. The accuracy of our new approach is comparable to that of calculations performed using other approximate methods over a relatively broad range of model parameters. However, our method converges relatively quickly when compared with most alternative schemes. These findings are very encouraging in view of the application of the new method for studying realistic nonadiabatic model problems in the condensed phase.

Dynamics and the breaking of a driven cage: I2 in solid Ar

Pump-probe measurements of I2 in solid Ar are reported and analyzed to extract a description of cage response to impulsive excitation, from the gentle kick, up to the breaking point. The most informative data are obtained through wavepacket motion on cage-bound, but otherwise dissociative, potentials where the chromophore acts as a transducer to drive the cage and to report on the local dynamics. This general class of dynamics is identified and analyzed as a function of energy in Ar, Kr, and Xe. The overdriven cage rebounds with a characteristic period of 1.2 ps that shows little dependence on excitation amplitude in all hosts. After rebound, the cage rings as a local resonant mode in Ar, with a period of 1 ps and dephasing time of 3 ps. This mode remains at the Debye edge in Kr and Xe, with periods of 630 and 800 fs, and dephasing times of 8 and 6 ps, respectively. In the bound B-state, the cage fluctuates toward its dilated equilibrium structure on a time scale of 3 ps, which is extracted from the down-chirp in the molecular vibrational frequency. When kicked with excess energy of 4 eV, the Ar cage breaks with 50% probability, and the molecule dissociates. The kinetics of polarization selective, multiphoton dissociation with Gaussian laser intensity profiles is delineated and the ballistics of cage breakout is described: The photodissociation proceeds by destruction of the local lattice, by creating interstitials and vacancies. During large amplitude motion on cage-bound potentials, sudden, nonadiabatic spin-flip transitions can be observed and quantified in space and time. The spin-flip occurs with unit probability in Ar when the I*-I bond is stretched beyond 6 A.

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

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