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.

Retinal photoisomerization versus counterion protonation in light and dark-adapted bacteriorhodopsin and its primary photoproduct.

Discovered over 50 years ago, bacteriorhodopsin is the first recognized and most widely studied microbial retinal protein. Serving as a light-activated proton pump, it represents the archetypal ion-pumping system. Here we compare the photochemical dynamics of bacteriorhodopsin light and dark-adapted forms with that of the first metastable photocycle intermediate known as "K". We observe that following thermal double isomerization of retinal in the dark from bio-active all-trans 15-anti to 13-cis, 15-syn, photochemistry proceeds even faster than the ~0.5 ps decay of the former, exhibiting ballistic wave packet curve crossing to the ground state. In contrast, photoexcitation of K containing a 13-cis, 15-anti chromophore leads to markedly multi-exponential excited state decay including much slower stages. QM/MM calculations, aimed to interpret these results, highlight the crucial role of protonation, showing that the classic quadrupole counterion model poorly reproduces spectral data and dynamics. Single protonation of ASP212 rectifies discrepancies and predicts triple ground state structural heterogeneity aligning with experimental observations. These findings prompt a reevaluation of counter ion protonation in bacteriorhodopsin and contribute to the broader understanding of its photochemical dynamics.

Reconciling TD-DFT and CASPT2 electronic structure methods for describing the photophysics of DNA.

Time-dependent density functional theory (TD-DFT) and multiconfigurational second-order perturbation theory (CASPT2) are two of the most widely used methods to investigate photoinduced dynamics in DNA-based systems. These methods sometimes give diverse dynamics in physiological environments usually modeled by quantum mechanics/molecular mechanics (QM/MM) protocol. In this work, we demonstrate for the uridine test case that the underlying topology of the potential energy surfaces of electronic states involved in photoinduced relaxation is similar in both electronic structure methods. This is verified by analyzing surface-hopping dynamics performed at the QM/MM level on aqueous solvated uridine at TD-DFT and CASPT2 levels. By constraining the dynamics to remain on π π * state we observe similar fluctuations in energy and relaxation lifetimes in surface-hopping dynamics in both TD-DFT and experimentally validated CASPT2 methods. This finding calls for a systematic comparison of the ES potential energy surfaces of DNA and RNA nucleosides at the single- and multi-reference levels of theory. The anomalous long excited state lifetime at the TD-DFT level is explained by n π * trapping due to the tendency of TD-DFT in QM/MM schemes with electrostatic embedding to underestimate the energy of the π π * state leading to a wrong π π * / n π * energetic order. A study of the FC energetics suggests that improving the description of the surrounding environment through polarizable embedding or by the expansion of QM layer with hydrogen-bonded waters helps restore the correct state order at TD-DFT level. Thus by combining TDDFT with an accurate modeling of the environment, TD-DFT is positioned as the standout protocol to model photoinduced dynamics in DNA-based aggregates and multimers.

Modeling pH-Dependent Biomolecular Photochemistry.

The tuning mechanism of pH can be extremely challenging to model computationally in complex biological systems, especially with respect to the photochemical properties. This article reports a protocol aimed at modeling pH-dependent photodynamics using a combination of constant-pH molecular dynamics and semiclassical nonadiabatic molecular dynamics simulations. With retinal photoisomerization in Anabaena sensory rhodopsin (ASR) as a testbed, we show that our protocol produces pH-dependent photochemical properties, such as the isomerization quantum yield or decay rates. We decompose our results into single-titrated residue contributions, identifying some key tuning amino acids. Additionally, we assess the validity of the single protonation state picture to represent the system at a given pH and propose the most populated protein charge state as a compromise between cost and accuracy.

Time-Resolved X-ray Absorption Spectroscopy: An MCTDH Quantum Dynamics Protocol.

Expressions for linear and nonlinear spectroscopy simulation in the X-ray window in which the time evolution of a photoexcited molecular system is treated via quantum dynamics are derived. By leveraging on the peculiar properties of core-excited/ionized states, first- and third-order response functions are recast in the limit of time-scale separation between the extremely short core-state lifetime and the (comparably longer) electronic-state transfer and nuclear vibrational motion. This work is a natural extension of Segatta et al. (J. Chem. Theory Comput. 2023, 19, 2075-2091), in which some of the present authors coupled MCTDH quantum dynamics to spectroscopy simulation at different levels of sophistication. Full quantum dynamics and approximate expressions are compared by simulating X-ray transient absorption spectroscopy at the carbon K-edge in the pyrene molecule.

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.

The carbonyl-lock mechanism underlying non-aromatic fluorescence in biological matter.

Challenging the basis of our chemical intuition, recent experimental evidence reveals the presence of a new type of intrinsic fluorescence in biomolecules that exists even in the absence of aromatic or electronically conjugated chemical compounds. The origin of this phenomenon has remained elusive so far. In the present study, we identify a mechanism underlying this new type of fluorescence in different biological aggregates. By employing non-adiabatic ab initio molecular dynamics simulations combined with a data-driven approach, we characterize the typical ultrafast non-radiative relaxation pathways active in non-fluorescent peptides. We show that the key vibrational mode for the non-radiative decay towards the ground state is the carbonyl elongation. Non-aromatic fluorescence appears to emerge from blocking this mode with strong local interactions such as hydrogen bonds. While we cannot rule out the existence of alternative non-aromatic fluorescence mechanisms in other systems, we demonstrate that this carbonyl-lock mechanism for trapping the excited state leads to the fluorescence yield increase observed experimentally, and set the stage for design principles to realize novel non-invasive biocompatible probes with applications in bioimaging, sensing, and biophotonics.

Characterizing Conical Intersections in DNA/RNA Nucleobases with Multiconfigurational Wave Functions of Varying Active Space Size.

We characterize the photochemically relevant conical intersections between the lowest-lying accessible electronic excited states of the different DNA/RNA nucleobases using Cholesky decomposition-based complete active space self-consistent field (CASSCF) algorithms. We benchmark two different basis set contractions and several active spaces for each nucleobase and conical intersection type, measuring for the first time how active space size affects conical intersection topographies in these systems and the potential implications these may have toward their description of photoinduced phenomena. Our results show that conical intersection topographies are highly sensitive to the electron correlation included in the model: by changing the amount (and type) of correlated orbitals, conical intersection topographies vastly change, and the changes observed do not follow any converging pattern toward the topographies obtained with the largest and most correlated active spaces. Comparison across systems shows analogous topographies for almost all intersections mediating population transfer to the dark 1nO/Nπ* states, while no similarities are observed for the "ethylene-like" conical intersection ascribed to mediate the ultrafast decay component to the ground state in all DNA/RNA nucleobases. Basis set size seems to have a minor effect, appearing to be relevant only for purine-based derivatives. We rule out structural changes as a key factor in classifying the different conical intersections, which display almost identical geometries across active space and basis set change, and we highlight instead the importance of correctly describing the electronic states involved at these crossing points. Our work shows that careful active space selection is essential to accurately describe conical intersection topographies and therefore to adequately account for their active role in molecular photochemistry.

Quantum-Classical Protocol for Efficient Characterization of Absorption Lineshape and Fluorescence Quenching upon Aggregation: The Case of Zinc Phthalocyanine Dyes.

A quantum-classical protocol that incorporates Jahn-Teller vibronic coupling effects and cluster analysis of molecular dynamics simulations is reported, providing a tool for simulations of absorption spectra and ultrafast nonadiabatic dynamics in large molecular photosystems undergoing aggregation in solution. Employing zinc phthalocyanine dyes as target systems, we demonstrated that the proposed protocol provided fundamental information on vibronic, electronic couplings and thermal dynamical effects that mostly contribute to the absorption spectra lineshape and the fluorescence quenching processes upon dye aggregation. Decomposing the various effects arising upon dimer formation, the structure-property relations associated with their optical responses have been deciphered at atomistic resolution.

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.

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.

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.

The OpenMolcas Web: A Community-Driven Approach to Advancing Computational Chemistry.

The developments of the open-source OpenMolcas chemistry software environment since spring 2020 are described, with a focus on novel functionalities accessible in the stable branch of the package or via interfaces with other packages. These developments span a wide range of topics in computational chemistry and are presented in thematic sections: electronic structure theory, electronic spectroscopy simulations, analytic gradients and molecular structure optimizations, ab initio molecular dynamics, and other new features. This report offers an overview of the chemical phenomena and processes OpenMolcas can address, while showing that OpenMolcas is an attractive platform for state-of-the-art atomistic computer simulations.

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.

Monitoring vibronic coherences and molecular aromaticity in photoexcited cyclooctatetraene with an X-ray probe: a simulation study.

Understanding conical intersection (CI) dynamics and subsequent conformational changes is key for exploring and controlling photo-reactions in aromatic molecules. Monitoring of their time-resolved dynamics remains a formidable experimental challenge. In this study, we simulate the photoinduced S3 to S1 non-adiabatic dynamics of cyclooctatetraene (COT), involving multiple CIs with relaxation times in good agreement with experiment. We further investigate the possibility to directly probe the CI passages in COT by off-resonant X-ray Raman spectroscopy (TRUECARS) and time-resolved X-ray diffraction (TRXD). We find that these signals sensitively monitor key chemical features during the ultrafast dynamics. First, we distinguish two CIs by using TRUECARS signals with their appearances at different Raman shifts. Second, we demonstrate that TRXD, where X-ray photons scatter off electron densities, can resolve ultrafast changes in the aromaticity of COT. It can further distinguish between planar and non-planar geometries explored during the dynamics, as e.g. two different tetraradical-type CIs. The knowledge gained from these measurements can give unique insight into fundamental chemical properties that dynamically change during non-adiabatic passages.

Nonlinear Molecular Electronic Spectroscopy via MCTDH Quantum Dynamics: From Exact to Approximate Expressions.

We present an accurate and efficient approach to computing the linear and nonlinear optical spectroscopy of a closed quantum system subject to impulsive interactions with an incident electromagnetic field. It incorporates the effect of ultrafast nonadiabatic dynamics by means of explicit numerical propagation of the nuclear wave packet. The fundamental expressions for the evaluation of first- and higher-order response functions are recast in a general form that can be used with any quantum dynamics code capable of computing the overlap of nuclear wave packets evolving in different states. Here we present the evaluation of these expressions with the multiconfiguration time-dependent Hartree (MCTDH) method. Application is made to pyrene, excited to its lowest bright excited state S2 which exhibits a sub-100-fs nonadiabatic decay to a dark state S1. The system is described by a linear vibronic coupling Hamiltonian, parametrized with multiconfiguration electronic structure methods. We show that the ultrafast nonadiabatic dynamics can have a remarkable effect on the spectral line shapes that goes beyond simple lifetime broadening. Furthermore, a widely employed approximate expression based on the time scale separation of dephasing and population relaxation is recast in the same theoretical framework. Application to pyrene shows the range of validity of such approximations.

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.

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.

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.