Quantum Dynamics in an Explicit Solvent Environment: A Photochemical Bond Cleavage Treated with a Combined QD/MD Approach

Solvent Effects on Ultrafast Charge Transfer Population: Insights from the Quantum Dynamics of Guanine-Cytosine in Chloroform

We study the ultrafast photoactivated dynamics of the hydrogen bonded dimer Guanine-Cytosine in chloroform solution, focusing on the population of the Guanine→Cytosine charge transfer state (GC-CT), an important elementary process for the photophysics and photochemistry of nucleic acids. We integrate a quantum dynamics propagation scheme, based on a linear vibronic model parameterized through time dependent density functional theory calculations, with four different solvation models, either implicit or explicit. On average, after 50 fs, 30∼40 % of the bright excited state population has been transferred to GC-CT. This process is thus fast and effective, especially when transferring from the Guanine bright excited states, in line with the available experimental studies. Independent of the adopted solvation model, the population of GC-CT is however disfavoured in solution with respect to the gas phase. We show that dynamical solvation effects are responsible for this puzzling result and assess the different chemical-physical effects modulating the population of CT states on the ultrafast time-scale. We also propose some simple analyses to predict how solvent can affect the population transfer between bright and CT states, showing that the effect of the solute/solvent electrostatic interactions on the energy of the CT state can provide a rather reliable indication of its possible population.

The Quest to Simulate Excited-State Dynamics of Transition Metal Complexes

This Perspective describes current computational efforts in the field of simulating photodynamics of transition metal complexes. We present the typical workflows and feature the strengths and limitations of the different contemporary approaches. From electronic structure methods suitable to describe transition metal complexes to approaches able to simulate their nuclear dynamics under the effect of light, we give particular attention to build a bridge between theory and experiment by critically discussing the different models commonly adopted in the interpretation of spectroscopic experiments and the simulation of particular observables. Thereby, we review all the studies of excited-state dynamics on transition metal complexes, both in gas phase and in solution from reduced to full dimensionality.

Quantum and semiclassical dynamical studies of nonadiabatic processes in solution: achievements and perspectives

We concisely review the main methodological approaches to model nonadiabatic dynamics in isotropic solutions and their applications. Three general classes of models are identified as the most used to include solvent effects in the simulations. The first model describes the solvent as a set of harmonic collective modes coupled to the solute degrees of freedom, and the second as a continuum, while the third explicitly includes solvent molecules in the calculations. The issues related to the use of these models in semiclassical and quantum dynamical simulations are discussed, as well as the main limitations and perspectives of each approach.

Ultrafast Spectroscopy of Photoactive Molecular Systems from First Principles: Where We Stand Today and Where We Are Going

Computational spectroscopy is becoming a mandatory tool for the interpretation of the complex, and often congested, spectral maps delivered by modern non-linear multi-pulse techniques. The fields of Electronic Structure Methods, Non-Adiabatic Molecular Dynamics, and Theoretical Spectroscopy represent the three pillars of the virtual ultrafast optical spectrometer, able to deliver transient spectra in silico from first principles. A successful simulation strategy requires a synergistic approach that balances between the three fields, each one having its very own challenges and bottlenecks. The aim of this Perspective is to demonstrate that, despite these challenges, an impressive agreement between theory and experiment is achievable now regarding the modeling of ultrafast photoinduced processes in complex molecular architectures. Beyond that, some key recent developments in the three fields are presented that we believe will have major impacts on spectroscopic simulations in the very near future. Potential directions of development, pending challenges, and rising opportunities are illustrated.

Ultrafast Photorelaxation of Uracil Embedded in an RNA Strand

Ultrafast photorelaxation of uracil can be hindered by its natural RNA environment. Multiscale quantum dynamical simulations show that the wave packet can be trapped in the photoexcited electronic state, which could potentially lead to photodamage.

Theoretical-computational modeling of charge transfer and intersystem crossing reactions in complex chemical systems

In this paper we present a theoretical-computational methodology specifically aimed at describing processes involving internal conversion or intersystem crossing, from atomistic (semiclassical) simulations and, hence, very suitable for treating complex atomic-molecular systems. The core of the presented approach is the evaluation of the diabatic perturbed energy surfaces of a portion of the whole system, treated at the quantum level and therefore preventively selected, in semi-classical interaction with the atomic-molecular environment. Subsequently, the estimation of the coupling between the diabatic surfaces and the inclusion of the obtained observables within a properly designed kinetic model allows the reconstruction of the whole phenomenology directly comparable to the experimental (typically kinetic) data. Application to two systems has demonstrated that the proposed approach can represent a valuable tool, somewhat complementary to other methods based on explicit quantum-dynamical approaches, for the theoretical-computational investigations of large and complex atomic-molecular systems.

Ultrafast Reactive Quantum Dynamics Coupled to Classical Solvent Dynamics Using an Ehrenfest Approach

The inclusion of solvent effects in the theoretical analysis of molecular processes becomes increasingly important. Currently, it is not feasible to directly include the solvent on the quantum level. We use an Ehrenfest approach to study the coupled time evolution of quantum dynamically treated solutes and classical solvents system. The classical dynamics of the solvent is coupled to the wavepacket dynamics of the solute and rotational and translational degrees of freedom of the solute are included classically. This allows quantum dynamics simulations for ultrafast processes that are decided by environment interactions without explicit separation of time scales. We show the application to the dissociation of ICN in liquid Ar as a proof of principal system and to the more applied example of uracil in water.

Mixed Quantum/Classical Method for Nonadiabatic Quantum Dynamics in Explicit Solvent Models: The ππ*/nπ* Decay of Thymine in Water as a Test Case

We present a novel mixed quantum classical dynamical method to include solvent effects on internal conversion (IC) processes. All the solute degrees of freedom are represented by a wavepacket moving according to nonadiabatic quantum dynamics, while the motion of an explicit solvent model is described by an ensemble of classical trajectories. The mutual coupling of the solute and solvent dynamics is included within a mean-field framework and the quantum and classical equations of motions are solved simultaneously. As a test case we apply our method to the ultrafast ππ* → nπ* decay of thymine in water. Solvent dynamical response modifies IC yield already on the 50 fs time scale. This effect is due to water librational motions that stabilize the most populated state. Pure static disorder, that is, the existence of different solvent configurations when photoexcitation takes place, also has a remarkable impact on the dynamics.

Molecular features in complex environment: Cooperative team players during excited state bond cleavage

Photoinduced bond cleavage is often employed for the generation of highly reactive carbocations in solution and to study their reactivity. Diphenylmethyl derivatives are prominent precursors in polar and moderately polar solvents like acetonitrile or dichloromethane. Depending on the leaving group, the photoinduced bond cleavage occurs on a femtosecond to picosecond time scale and typically leads to two distinguishable products, the desired diphenylmethyl cations (Ph2CH(+)) and as competing by-product the diphenylmethyl radicals ([Formula: see text]). Conical intersections are the chief suspects for such ultrafast branching processes. We show for two typical examples, the neutral diphenylmethylchloride (Ph2CH-Cl) and the charged diphenylmethyltriphenylphosphonium ions ([Formula: see text]) that the role of the conical intersections depends not only on the molecular features but also on the interplay with the environment. It turns out to differ significantly for both precursors. Our analysis is based on quantum chemical and quantum dynamical calculations. For comparison, we use ultrafast transient absorption measurements. In case of Ph2CH-Cl, we can directly connect the observed signals to two early three-state and two-state conical intersections, both close to the Franck-Condon region. In case of the [Formula: see text], dynamic solvent effects are needed to activate a two-state conical intersection at larger distances along the reaction coordinate.

Molecular Dynamics Simulations of the Initial-State Predict Product Distributions of Dediazoniation of Aryldiazonium in Binary Solvents

The dediazoniation of aryldiazonium salts in mixed solvents proceeds by a borderline SN1 and SN2 pathway, and product distribution should be proportional to the composition of the solvation shell of the carbon attached to the -N2 group (ipso carbon). The rates of dediazoniation of 2,4,6-trimethylbenzenediazonium in water, methanol, ethanol, propanol, and acetonitrile were similar, but measured product distributions were noticeably dependent on the nature of the water/cosolvent mixture. Here we demonstrated that solvent distribution in the first solvation shell of the ipso carbon, calculated from classical molecular dynamics simulations, is equal to the measured product distribution. Furthermore, we showed that regardless of the charge distribution of the initial state, i.e., whether the positive charge is smeared over the molecule or localized on phenyl moiety, the solvent distribution around the reaction center is nearly the same.