WFOT: A Wave Function Overlap Tool between Single- and Multi-Reference Electronic Structure Methods for Spectroscopy Simulation.

We report the development of a novel diagnostic tool, named wave function overlap tool (WFOT), designed to evaluate the overlap between wave functions computed at single-reference [i.e., time-dependent density functional theory or configuration interaction singles (CIS)] and multireference (i.e., CASSCF/CASPT2) electronic structure levels of theory. It relies on truncating the single- and multireference WFs to CIS-like expansions spanning the same configurational space and maximizing the molecular orbital overlap by means of a unitary transformation. To demonstrate the functionality of the tool, we calculate the transient spectrum of acetylacetone by evaluating excited state absorption signals with multireference quality on top of single-reference on-the-fly dynamics simulations. Semiautomatic spectra generation is facilitated by interfacing the tool with the COBRAMM package, which also allows one to use WFOT with several quantum chemistry codes such as Gaussian, NWChem, and OpenMolcas. Other exciting possibilities for the utilization of the code beyond the simulation of transient absorption spectroscopy are eventually discussed.

In Silico Chemical Experiments in the Age of AI: From Quantum Chemistry to Machine Learning and Back.

Computational chemistry is an indispensable tool for understanding molecules and predicting chemical properties. However, traditional computational methods face significant challenges due to the difficulty of solving the Schrödinger equations and the increasing computational cost with the size of the molecular system. In response, there has been a surge of interest in leveraging artificial intelligence (AI) and machine learning (ML) techniques to in silico experiments. Integrating AI and ML into computational chemistry increases the scalability and speed of the exploration of chemical space. However, challenges remain, particularly regarding the reproducibility and transferability of ML models. This review highlights the evolution of ML in learning from, complementing, or replacing traditional computational chemistry for energy and property predictions. Starting from models trained entirely on numerical data, a journey set forth toward the ideal model incorporating or learning the physical laws of quantum mechanics. This paper also reviews existing computational methods and ML models and their intertwining, outlines a roadmap for future research, and identifies areas for improvement and innovation. Ultimately, the goal is to develop AI architectures capable of predicting accurate and transferable solutions to the Schrödinger equation, thereby revolutionizing in silico experiments within chemistry and materials science.

DELFI: a computer oracle for recommending density functionals for excited states calculations.

Density functional theory (DFT) is the workhorse of computational quantum chemistry. One of its main limitations is that choosing the right functional is a non-trivial task left for human experts. The choice is particularly hard for excited state calculations when using its time-dependent formulation (TD-DFT). This is due to the approximations of the method, but also because the photophysical properties of a molecule are defined by a manifold of states that all need to be properly described. This includes not only the relative energy of the states, but also capturing the correct character, order, and intensity of the transitions. In this work, we developed a neural network to recommend functionals to be used on molecules for TD-DFT calculations, by simultaneously considering all these properties for a manifold of states. This was possible by developing a scoring system to define the accuracy of an excited state's calculation against a higher-accuracy reference. The scoring system is generalizable to any level of theory; we here applied it to evaluate the performance of common functionals of different rungs against a higher accuracy method on a large set of organic molecules. The results are collected in a database that we released and made open, providing four million data points to the community for future applications. The scoring system assigns a value between zero and one hundred to each functional for each molecule, transforming the complicated task of learning photophysical properties into a simpler regression task. We used the dataset to train a graph attention neural network to predict the scores for unseen molecules. We call this oracle DELFI (Data-driven EvaLuation of Functionals by Inference), which can be used to quickly screen and predict the ranking of functionals to calculate the optical properties of organic molecules. We validated DELFI in two in silico experiments: choosing a common functional for a series of spiropyran-merocyanine isomers and a unique functional to screen a large dataset of over 50 000 organic photovoltaic molecules, for which an extensive benchmark would be unfeasible. A corresponding web application allows DELFI to be easily run and the results to be analyzed, alleviating the hurdle of choosing the right functional for TD-DFT calculations.

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.

Good Vibrations Report on the DNA Quadruplex Binding of an Excited State Amplified Ruthenium Polypyridyl IR Probe.

The nitrile containing Ru(II)polypyridyl complex [Ru(phen)2(11,12-dCN-dppz)]2+ (1) is shown to act as a sensitive infrared probe of G-quadruplex (G4) structures. UV-visible absorption spectroscopy reveals enantiomer sensitive binding for the hybrid htel(K) and antiparallel htel(Na) G4s formed by the human telomer sequence d[AG3(TTAG3)3]. Time-resolved infrared (TRIR) of 1 upon 400 nm excitation indicates dominant interactions with the guanine bases in the case of Λ-1/htel(K), Δ-1/htel(K), and Λ-1/htel(Na) binding, whereas Δ-1/htel(Na) binding is associated with interactions with thymine and adenine bases in the loop. The intense nitrile transient at 2232 cm-1 undergoes a linear shift to lower frequency as the solution hydrogen bonding environment decreases in DMSO/water mixtures. This shift is used as a sensitive reporter of the nitrile environment within the binding pocket. The lifetime of 1 in D2O (ca. 100 ps) is found to increase upon DNA binding, and monitoring of the nitrile and ligand transients as well as the diagnostic DNA bleach bands shows that this increase is related to greater protection from the solvent environment. Molecular dynamics simulations together with binding energy calculations identify the most favorable binding site for each system, which are in excellent agreement with the observed TRIR solution study. This study shows the power of combining the environmental sensitivity of an infrared (IR) probe in its excited state with the TRIR DNA "site effect" to gain important information about the binding site of photoactive agents and points to the potential of such amplified IR probes as sensitive reporters of biological environments.

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.

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.

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.

Sampling effects in quantum mechanical/molecular mechanics trajectory surface hopping non-adiabatic dynamics.

The impact of different initial conditions in non-adiabatic trajectory surface hopping dynamics within a hybrid quantum mechanical/molecular mechanics scheme is investigated. The influence of a quantum sampling, based on a Wigner distribution, a fully thermal sampling, based on classical molecular dynamics, and a quantum sampled system, but thermally equilibrated with the environment, is investigated on the relaxation dynamics of solvated fulvene after light irradiation. We find that the decay from the first singlet excited state to the ground state shows high dependency on the initial condition and simulation parameters. The three sampling methods lead to different distributions of initial geometries and momenta, which then affect the fate of the excited state dynamics. We evaluated both the effect of sampling geometries and momenta, analysing how the ultrafast decay of fulvene changes accordingly. The results are expected to be of interest to decide how to initialize non-adiabatic dynamics in the presence of the environment. This article is part of the theme issue 'Chemistry without the Born-Oppenheimer approximation'.

QM/MM Nonadiabatic Dynamics: the SHARC/COBRAMM Approach.

We present the SHARC/COBRAMM approach to enable easy and efficient excited-state dynamics simulations at different levels of electronic structure theory in the presence of complex environments using a quantum mechanics/molecular mechanics (QM/MM) setup. SHARC is a trajectory surface-hoping method that can incorporate the simultaneous effects of nonadiabatic and spin-orbit couplings in the excited-state dynamics of molecular systems. COBRAMM allows ground- and excited-state QM/MM calculations using a subtractive scheme, with electrostatic embedding and a hydrogen link-atom approach. The combination of both free and open-source program packages provides a modular and extensive framework to model nonadiabatic processes after light irradiation from the atomistic scale to the nano-scale. As an example, the relaxation of acrolein from S1 to T1 in solution is provided.

Ultrafast Dynamics of the Isolated Adenosine-5'-triphosphate Dianion Probed by Time-Resolved Photoelectron Imaging.

The excited state dynamics of the doubly deprotonated dianion of adenosine-5'-triphosphate, [ATP-H2]2-, has been spectroscopically explored by time-resolved photoelectron spectroscopy following excitation at 4.66 eV. Time-resolved photoelectron spectra show that two competing processes occur for the initially populated 1ππ* state. The first is rapid electron emission by tunneling through a repulsive Coulomb barrier as the 1ππ* state is a resonance. The second is nuclear motion on the 1ππ* state surface leading to an intermediate that no longer tunnels and subsequently decays by internal conversion to the ground electronic state. The spectral signatures of the features are similar to those observed for other adenine-derivatives, suggesting that this nucleobase is quite insensitive to the nearby negative charges localized on the phosphates, except of course for the appearance of the additional electron tunneling channel, which is open in the dianion.

Enhanced Rigidity Changes Ultraviolet Absorption: Effect of a Merocyanine Binder on G-Quadruplex Photophysics.

The urge to discover selective fluorescent binders to G-quadruplexes (G4s) for rapid diagnosis must be linked to understand the effect that those have on the DNA photophysics. Herein, we report on the electronic excited states of a bound merocyanine dye to c-Myc G4 using extensive multiscale quantum mechanics/molecular mechanics calculations. We find that the absorption spectra of c-Myc G4, both without and with the intercalated dye, are mainly composed of exciton states and mixed local/charge-transfer states. The presence of merocyanine hardly affects the energy range of the guanine absorption or the number of guanines excited. However, it triggers a substantial amount (16%) of detrimental pure charge-transfer states involving oxidized guanines. We identify the rigidity introduced by the probe in G4, reducing the overlap among guanines, as the one responsible for the changes in the exciton and charge-transfer states, ultimately leading to a redshift of the absorption maximum.

Site-Specific Photo-oxidation of the Isolated Adenosine-5'-triphosphate Dianion Determined by Photoelectron Imaging.

Photoelectron imaging of the isolated adenosine-5'-triphosphate dianion excited to the 1ππ* states reveals that electron emission is predominantly parallel to the polarization axis of the light and arises from subpicosecond electron tunneling through the repulsive Coulomb barrier (RCB). The computed RCB shows that the most probable electron emission site is on the amino group of adenine. This is consistent with the photoelectron imaging: excitation to the 1ππ* states leads to an aligned ensemble distributed predominantly parallel to the long axis of adenine; the subsequent electron tunneling site is along this axis; and the negatively charged phosphate groups guide the outgoing electron mostly along this axis at long range. Imaging of electron tunneling from polyanions combined with computational chemistry may offer a general route for probing the intrinsic photo-oxidation site and dynamics as well as the overall structure of complex isolated species.

Spiropyran Meets Guanine Quadruplexes: Isomerization Mechanism and DNA Binding Modes of Quinolizidine-Substituted Spiropyran Probes.

The recent delivery of a fluorescent quinolizidine-substituted spiropyran, which is able to switch in vivo and bind to guanine quadruplexes (G4) at physiological pH values, urged us to elucidate its molecular switching and binding mechanism. Combining multiscale dynamical studies and accurate quantum chemical calculations, we show that, both in water and in the G4 environment, the switching of the spiropyran ring is not promoted by an initial protonation event-as expected by the effect of low pH solutions-but that the deprotonated merocyanine form is an intermediate of the reaction leading to the protonated open species. Additionally, we investigate the binding of both deprotonated and protonated open forms of merocyanine to c-MYC G4s. Both species bind to G4s albeit with different hydrogen-bond patterns and provide distinct rotamers around the exocyclic double bond of the merocyanine forms. Altogether, our study sheds light on the pharmacophoric points for the binding of these probes to DNA, and thereby, contributes to future developments of new G4 binders of the remarkable family of quinolizidine-substituted spiropyrans.

Revealing Ultrafast Population Transfer between Nearly Degenerate Electronic States.

The response of a molecule to photoexcitation is governed by the coupling of its electronic states. However, if the energetic spacing between the electronically excited states at the Franck-Condon window becomes sufficiently small, it is infeasible to selectively excite and monitor individual states with conventional time-resolved spectroscopy, preventing insight into the energy transfer and relaxation dynamics of the molecule. Here, we demonstrate how the combination of time-resolved spectroscopy and extensive surface hopping dynamics simulations with a global fit approach on individually excited ensembles overcomes this limitation and resolves the dynamics in the n3p Rydberg states in acetone. Photoelectron transients of the three closely spaced states n3px, n3py, and n3pz are used to validate the theoretical results, which in turn allow retrieving a comprehensive kinetic model describing the mutual interactions of these states for the first time.

Directional and regioselective hole injection of spiropyran photoswitches intercalated into A/T-duplex DNA.

The electron-hole injection from a family of spiropyran photoswitches into A/T-duplex DNA has been investigated at the molecular level for the first time. Multiscale computations coupled with automatized quantitative wavefunction analysis reveal a pronounced directionality and regioselectivity towards the template strand of the duplex DNA. Our findings suggest that this directional and regioselective photoinduced electron-hole transfer could thus be exploited to tailor the charge transport processes in DNA in specific applications.

DNA-binding mechanism of spiropyran photoswitches: the role of electrostatics.

The binding mechanism of the protonated open form of three spiropyran derivatives into a 12-mer (poly-dAT)2 has been unveiled by means of computational methods. It is found that the electrostatic term in the probe:DNA binding energy, modulates the binding mode, providing new guidelines for the design of spiropyran photoswitches with specific binding modes to DNA.