Modeling solvent effects and convergence of 31P-NMR shielding calculations with COBRAMM.

Solvent effects on 31P-NMR parameters for triphenylphosphine oxide and triphenylphosphine in chloroform have been extensively investigated by testing different solvation models. The solvent is described implicitly, mixed implicitly/explicitly, and using full explicit models. Polarizable continuum model (PCM), molecular dynamic (MD) simulations, and hybrid quantum mechanics/molecular mechanics (QM/MM) calculations are used to disclose the effects of solute/solvent interactions and, more generally, the role of the embedding in NMR simulations. The results show the beneficial effect of carrying out QM/MM optimizations on top of geometries directly extracted from classical MD simulations, used to ensure representative conformational sampling. The nuclear shielding convergence has been tested against a different number of snapshots and with the inclusion of solvent shells into the QM region. An automated MD//QM/MM//GIAO protocol, implemented in the COBRAMM package, is here proposed and tested on trimethyl phosphite showing that our approach boosts the convergence of nuclear shielding satisfactorily. The present work aims to be a stepping-stone to assess proper QM/MM computational strategies in simulating chemical shifts in non-homogeneous systems like supramolecular and biological systems.

Ultrafast X-Ray Probes of Elementary Molecular Events.

Elementary events that determine photochemical outcomes and molecular functionalities happen on the femtosecond and subfemtosecond timescales. Among the most ubiquitous events are the nonadiabatic dynamics taking place at conical intersections. These facilitate ultrafast, nonradiative transitions between electronic states in molecules that can outcompete slower relaxation mechanisms such as fluorescence. The rise of ultrafast X-ray sources, which provide intense light pulses with ever-shorter durations and larger observation bandwidths, has fundamentally revolutionized our spectroscopic capabilities to detect conical intersections. Recent theoretical studies have demonstrated an entirely new signature emerging once a molecule traverses a conical intersection, giving detailed insights into the coupled nuclear and electronic motions that underlie, facilitate, and ultimately determine the ultrafast molecular dynamics. Following a summary of current sources and experiments, we survey these techniques and provide a unified overview of their capabilities. We discuss their potential to dramatically increase our understanding of ultrafast photochemistry.

Photoisomerization transition state manipulation by entangled two-photon absorption.

We demonstrate how two-photon excitation with quantum light can influence elementary photochemical events. The azobenzene trans → cis isomerization following entangled two-photon excitation is simulated using quantum nuclear wave packet dynamics. Photon entanglement modulates the nuclear wave packets by coherently controlling the transition pathways. The photochemical transition state during passage of the reactive conical intersection in azobenzene photoisomerization is strongly affected with a noticeable alteration of the product yield. Quantum entanglement thus provides a novel control knob for photochemical reactions. The distribution of the vibronic coherences during the conical intersection passage strongly depends on the shape of the initial wave packet created upon quantum light excitation. X-ray signals that can experimentally monitor this coherence are simulated.

Imaging conical intersection dynamics during azobenzene photoisomerization by ultrafast X-ray diffraction.

X-ray diffraction is routinely used for structure determination of stationary molecular samples. Modern X-ray photon sources, e.g., from free-electron lasers, enable us to add temporal resolution to these scattering events, thereby providing a movie of atomic motions. We simulate and decipher the various contributions to the X-ray diffraction pattern for the femtosecond isomerization of azobenzene, a textbook photochemical process. A wealth of information is encoded besides real-time monitoring of the molecular charge density for the cis to trans isomerization. In particular, vibronic coherences emerge at the conical intersection, contributing to the total diffraction signal by mixed elastic and inelastic photon scattering. They cause distinct phase modulations in momentum space, which directly reflect the real-space phase modulation of the electronic transition density during the nonadiabatic passage. To overcome the masking by the intense elastic scattering contributions from the electronic populations in the total diffraction signal, we discuss how this information can be retrieved, e.g., by employing very hard X-rays to record large scattering momentum transfers.

Unveiling the spatial distribution of molecular coherences at conical intersections by covariance X-ray diffraction signals.

The outcomes and timescales of molecular nonadiabatic dynamics are decisively impacted by the quantum coherences generated at localized molecular regions. In time-resolved X-ray diffraction imaging, these coherences create distinct signatures via inelastic photon scattering, but they are buried under much stronger background elastic features. Here, we exploit the rich dynamical information encoded in the inelastic patterns, which we reveal by frequency-dispersed covariance ultrafast powder X-ray diffraction of stochastic X-ray free-electron laser pulses. This is demonstrated for the photoisomerization of azobenzene involving the passage through a conical intersection, where the nuclear wave packet branches and explores different quantum pathways. Snapshots of the coherence dynamics are obtained at high frequency shifts, not accessible with conventional diffraction measurements. These provide access to the timing and to the confined spatial distribution of the valence electrons directly involved in the conical intersection passage. This study can be extended to full three-dimensional imaging of conical intersections with ultrafast X-ray and electron diffraction.

Double Thionated Pyrimidine Nucleobases: Molecular Tools with Tunable Photoproperties.

Sulfur-substituted nucleobases are DNA and RNA base derivatives that exhibit extremely efficient photoinduced intersystem crossing (ISC) dynamics into the lowest-energy triplet state. The long-lived and reactive triplet states of sulfur-substituted nucleobases are crucial due to their wide range of potential applications in medicine, structural biology, and the development of organic light-emitting diodes (OLEDs) and other emerging technologies. However, a comprehensive understanding of non-negligible wavelength-dependent changes in the internal conversion (IC) and ISC events is still lacking. Here, we study the underlying mechanism using joint experimental gas-phase time-resolved photoelectron spectroscopy (TRPES) and theoretical quantum chemistry methods. We combine 2,4-dithiouracil (2,4-DTU) TRPES experimental data with computational analysis of the different photodecay processes, which are induced by increasing excitation energies along the entire linear absorption (LA) ultraviolet (UV) spectrum. Our results show how the double-thionated uracil (U), i.e., 2,4-DTU, appears as a versatile photoactivatable instrument. Multiple decay processes can be initiated with different ISC rates or triplet-state lifetimes that resemble the distinctive behavior of the singly substituted 2- or 4-thiouracil (2-TU or 4-TU). We obtained a clear partition of the LA spectrum based on the dominant photoinduced process. Our work clarifies the reasons behind the wavelength-dependent changes in the IC, ISC, and triplet-state lifetimes in doubly thionated U, becoming a biological system of utmost importance for wavelength-controlled applications. These mechanistic details and photoproperties are transferable to closely related molecular systems such as thionated thymines.

Separating Convective from Diffusive Mass Transport Mechanisms in Ionic Liquids by Redox Pro-fluorescence Microscopy.

The study of electrochemical reactivity requires analytical techniques capable of probing the diffusion of reactants and products to and from electrified interfaces. Information on diffusion coefficients is often obtained indirectly by modeling current transients and cyclic voltammetry data, but such measurements lack spatial resolution and are accurate only if mass transport by convection is negligible. Detecting and accounting for adventitious convection in viscous and wet solvents, such as ionic liquids, is technically challenging. We have developed a direct, spatiotemporally resolved optical tracking of diffusion fronts which can detect and resolve convective disturbances to linear diffusion. By tracking the movement of an electrode-generated fluorophore, we demonstrate that parasitic gas evolving reactions lead to 10-fold overestimates of macroscopic diffusion coefficients. A hypothesis is put forward linking large barriers to inner-sphere redox reactions, such as hydrogen gas evolution, to the formation of cation-rich overscreening and crowding double layer structures in imidazolium-based ionic liquids.

The protein environment restricts the intramolecular charge transfer character of the luciferine/luciferase complex.

The electronic characterization of the luciferine/luciferase complex is fundamental to tune its photophysical properties and develop more efficient devices based on this luminiscent system. Here, we apply molecular dynamics simulations, hybrid quantum mechanics/molecular mechanics (QM/MM) calculations and transition density analysis to compute the absorption and emission spectra of luciferine/luciferase and analyze the nature of the relevant electronic state and its behaviour with the intramolecular and intermolecular degrees of freedom. It is found that the torsional motion of the chromophore is hampered by the presence of the enzyme, reducing the intramolecular charge transfer nature of the absorbing and emitting state. In addition, such a reduced charge transfer character does not correlate in a strong way neither with the intramolecular motion of the chromophore nor with the chromophore/amino-acid distances. However, the presence of a polar environment around the oxygen atom of the thiazole ring of the oxyluciferin, coming from both the protein and the solvent, enhances the charge transfer character of the emitting state.

Engineering Azobenzene Derivatives to Control the Photoisomerization Process.

In this work, we show how the structural features of photoactive azobenzene derivatives can influence the photoexcited state behavior and the yield of the trans/cis photoisomerization process. By combining high-resolution transient absorption experiments in the vis-NIR region and quantum chemistry calculations (TDDFT and RASPT2), we address the origin of the transient signals of three poly-substituted push-pull azobenzenes with an increasing strength of the intramolecular interactions stabilizing the planar trans isomer (absence of intramolecular H-bonds, methyl, and traditional H-bond, respectively, for 4-diethyl-4'-nitroazobenzene, Disperse Blue 366, and Disperse Blue 165) and a commercial red dye showing keto-enol tautomerism involving the azo group (Sudan Red G). Our results indicate that the intramolecular H-bonds can act as a "molecular lock" stabilizing the trans isomer and increasing the energy barrier along the photoreactive CNNC torsion coordinate, thus preventing photoisomerization in the Disperse Blue dyes. In contrast, the involvement of the azo group in keto-enol tautomerism can be employed as a strategy to change the nature of the lower excited state and remove the nonproductive symmetric CNN/NNC bending pathway typical of the azo group, thus favoring the productive torsional motion. Taken together, our results can provide guidelines for the structural design of azobenzene-based photoswitches with a tunable excited state behavior.

Modus Operandi of a Pedalo-Type Molecular Switch: Insight from Dynamics and Theoretical Spectroscopy.

Molecular switches which can be triggered by light to interconvert between two or more well-defined conformation differing in their chemical or physical properties are fundamental for the development of materials with on-demand functionalities. Recently, a novel molecular switch based on a the azodicarboxamide core has been reported. It exhibits a volume-conserving conformational change upon excitation, making it a promising candidate for embedding in confined environments. In order to rationally implement and efficiently utilize the azodicarboxamide molecular switch, detailed insight into the coordinates governing the excited-state dynamics is needed. Here, we report a detailed comparative picture of the molecular motion at the atomic level in the presence and absence of explicit solvent. Our hybrid quantum mechanics/molecular mechanics (QM/MM) excited state simulations reveal that, although the energy landscape is slightly modulated by the solvation, the light-induced motion is dominated by a bending-assisted pedalo-type motion independent of the solvation. To support the predicted mechanism, we simulate time-resolved IR spectroscopy from first principles, thereby resolving fingerprints of the light-induced switching process. Our calculated time-resolved data are in good agreement with previously reported measured spectra.

Environment-Driven Coherent Population Transfer Governs the Ultrafast Photophysics of Tryptophan.

By combining UV transient absorption spectroscopy with sub-30-fs temporal resolution and CASPT2/MM calculations, we present a complete description of the primary photoinduced processes in solvated tryptophan. Our results shed new light on the role of the solvent in the relaxation dynamics of tryptophan. We unveil two consecutive coherent population transfer events involving the lowest two singlet excited states: a sub-50-fs nonadiabatic La → Lb transfer through a conical intersection and a subsequent 220 fs reverse Lb → La transfer due to solvent-assisted adiabatic stabilization of the La state. Vibrational fingerprints in the transient absorption spectra provide compelling evidence of a vibronic coherence established between the two excited states from the earliest times after photoexcitation and lasting until the back-transfer to La is complete. The demonstration of response to the environment as a driver of coherent population dynamics among the excited states of tryptophan closes the long debate on its solvent-assisted relaxation mechanisms and extends its application as a local probe of protein dynamics to the ultrafast time scales.

Automatized protocol and interface to simulate QM/MM time-resolved transient absorption at TD-DFT level with COBRAMM.

We present a series of new implementations that we recently introduced in COBRAMM, the open-source academic software developed in our group. The goal of these implementations is to offer an automatized workflow and interface to simulate time-resolved transient absorption (TA) spectra of medium-to-big chromophore embedded in a complex environment. Therefore, the excited states absorption and the stimulated emission are simulated along nonadiabatic dynamics performed with trajectory surface hopping. The possibility of treating systems from medium to big size is given by the use of time-dependent density functional theory (TD-DFT) and the presence of the environment is taken into account employing a hybrid quantum mechanics/molecular mechanics (QM/MM) scheme. The full implementation includes a series of auxiliary scripts to properly setup the QM/MM system, the calculation of the wavefunction overlap along the dynamics for the propagation, the evaluation of the transition dipole moment at linear response TD-DFT level, and scripts to setup, run and analyze the TA from an ensemble of trajectories. Altogether, we believe that our implementation will open the door to the easily simulate the time-resolved TA of systems so far computationally inaccessible.

Modelling quenching mechanisms of disordered molecular systems in the presence of molecular aggregates.

Exciton density dynamics recorded in time-resolved spectroscopic measurements is a useful tool to recover information on energy transfer (ET) processes that can occur at different timescales, up to the ultrafast regime. Macroscopic models of exciton density decays, involving both direct Förster-like ET and diffusion mechanisms for exciton-exciton annihilation, are largely used to fit time-resolved experimental data but generally neglect contributions from molecular aggregates that can work as quenching species. In this work, we introduce a macroscopic model that includes contributions from molecular aggregate quenchers in a disordered molecular system. As an exemplifying case, we considered a homogenous distribution of rhodamine B dyes embedded in organic nanoparticles to set the initial parameters of the proposed model. The influence of such model parameters is systematically analysed, showing that the presence of molecular aggregate quenchers can be monitored by evaluating the exciton density long time decays. We showed that the proposed model can be applied to molecular systems with ultrafast decays, and we anticipated that it could be used in future studies for global fitting of experimental data with potential support from first-principles simulations.

Coherent vibrational modes promote the ultrafast internal conversion and intersystem crossing in thiobases.

Thionated nucleobases are obtained by replacing oxygen with sulphur atoms in the canonical nucleobases. They absorb light efficiently in the near-ultraviolet, populating singlet states which undergo intersystem crossing to the triplet manifold on an ultrashort time scale with a high quantum yield. Nonetheless there are still important open questions about the primary mechanisms responsible for this ultrafast transition. Here we track both the electronic and the vibrational ultrafast excited-state dynamics towards the triplet state for solvated 4-thiothymidine (4TT) and 4-thiouracil (4TU) with sub-30 fs broadband transient absorption spectroscopy in the ultraviolet. A global and target analysis allows us to simultaneously resolve the contributions of the different electronically and vibrationally excited states to the whole data set. Our experimental results, combined with state-of-the-art quantum mechanics/molecular mechanics simulations and Damped Oscillation Associated Spectra (DOAS) and target analysis, support that the relaxation to the triplet state is mediated by conical intersections promoted by vibrational coherences through the population of an intermediate singlet state. In addition, the analysis of the coherent vibrational dynamics reveals that, despite sharing the same relaxation mechanism and similar chemical structures, 4TT and 4TU exhibit rather different geometrical deformations, characterized by the conservation of planarity in 4TU and its partial rupture in 4TT.

Luciferase-free Luciferin Electrochemiluminescence.

Luciferin is one of Nature's most widespread luminophores, and enzymes that catalyze luciferin luminescence are the basis of successful commercial "glow" assays for gene expression and metabolic ATP formation. Herein we report an electrochemical method to promote firefly's luciferin luminescence in the absence of its natural biocatalyst-luciferase. We have gained experimental and computational insights on the mechanism of the enzyme-free luciferin electrochemiluminescence, demonstrated its spectral tuning from green to red by means of electrolyte engineering, proven that the colour change does not require, as still debated, a keto/enol isomerization of the light emitter, and gained evidence of the electrostatic-assisted stabilization of the charge-transfer excited state by double layer electric fields. Luciferin's electrochemiluminescence, as well as the in situ generation of fluorescent oxyluciferin, are applied towards an optical measurement of diffusion coefficients.

Diffractive Imaging of Conical Intersections Amplified by Resonant Infrared Fields.

The fate of virtually all photochemical reactions is determined by conical intersections. These are energetically degenerate regions of molecular potential energy surfaces that strongly couple electronic states, thereby enabling fast relaxation channels. Their direct spectroscopic detection relies on weak features that are often buried beneath stronger, less interesting contributions. For azobenzene photoisomerization, a textbook photochemical reaction, we demonstrate how a resonant infrared field can be employed during the conical intersection passage to significantly enhance its coherence signatures in time-resolved X-ray diffraction while leaving the product yield intact. This transition-state amplification holds promise to bring signals of conical intersections above the detection threshold.

iSPECTRON: A simulation interface for linear and nonlinear spectra with ab-initio quantum chemistry software.

We introduce iSPECTRON, a program that parses data from common quantum chemistry software (NWChem, OpenMolcas, Gaussian, Cobramm, etc.), produces the input files for the simulation of linear and nonlinear spectroscopy of molecules with the Spectron code, and analyzes the spectra with a broad range of tools. Vibronic spectra are expressed in term of the electronic eigenstates, obtained from quantum chemistry computations, and vibrational/bath effects are incorporated in the framework of the displaced harmonic oscillator model, where all required quantities are computed at the Franck-Condon point. The program capabilities are illustrated by simulating linear absorption, transient absorption and two dimensional electronic spectra of the pyrene molecule. Calculations at two levels of electronic structure theory, time-dependent density functional theory (with NWChem) and RASSCF/RASPT2 (with OpenMolcas) are presented and compared where possible. The iSPECTRON program is available online at https://github.com/ispectrongit/iSPECTRON/ and distributed open source under the terms of the Educational Community License version 2.0 (ECL 2.0).

Manipulating Core Excitations in Molecules by X-Ray Cavities.

Core excitations on different atoms are highly localized and therefore decoupled. By placing molecules in an x-ray cavity the core transitions become coupled via the exchange of cavity photons and form delocalized hybrid light-matter excitations known as core polaritons. We demonstrate these effects for the two inequivalent carbon atoms in 1,1-difluoroethylene. Polariton signatures in the x-ray absorption, two-photon absorption, and multidimensional four-wave mixing signals are predicted.

Control of Protonated Schiff Base Excited State Decay within Visual Protein Mimics: A Unified Model for Retinal Chromophores.

Artificial biomimetic chromophore-protein complexes inspired by natural visual pigments can feature color tunability across the full visible spectrum. However, control of excited state dynamics of the retinal chromophore, which is of paramount importance for technological applications, is lacking due to its complex and subtle photophysics/photochemistry. Here, ultrafast transient absorption spectroscopy and quantum mechanics/molecular mechanics simulations are combined for the study of highly tunable rhodopsin mimics, as compared to retinal chromophores in solution. Conical intersections and transient fluorescent intermediates are identified with atomistic resolution, providing unambiguous assignment of their ultrafast excited state absorption features. The results point out that the electrostatic environment of the chromophore, modified by protein point mutations, affects its excited state properties allowing control of its photophysics with same power of chemical modifications of the chromophore. The complex nature of such fine control is a fundamental knowledge for the design of bio-mimetic opto-electronic and photonic devices.

Quantum Chemical Modeling of the Photoinduced Activity of Multichromophoric Biosystems.

Multichromophoric biosystems represent a broad family with very diverse members, ranging from light-harvesting pigment-protein complexes to nucleic acids. The former are designed to capture, harvest, efficiently transport, and transform energy from sunlight for photosynthesis, while the latter should dissipate the absorbed radiation as quickly as possible to prevent photodamages and corruption of the carried genetic information. Because of the unique electronic and structural characteristics, the modeling of their photoinduced activity is a real challenge. Numerous approaches have been devised building on the theoretical development achieved for single chromophores and on model Hamiltonians that capture the essential features of the system. Still, a question remains: is a general strategy for the accurate modeling of multichromophoric systems possible? By using a quantum chemical point of view, here we review the advancements developed so far highlighting differences and similarities with the single chromophore treatment. Finally, we outline the important limitations and challenges that still need to be tackled to reach a complete and accurate picture of their photoinduced properties and dynamics.