COBRAMM 2.0 — A software interface for tailoring molecular electronic structure calculations and running nanoscale (QM/MM) simulations

MiMiC: A high-performance framework for multiscale molecular dynamics simulations

MiMiC is a framework for performing multiscale simulations in which loosely coupled external programs describe individual subsystems at different resolutions and levels of theory. To make it highly efficient and flexible, we adopt an interoperable approach based on a multiple-program multiple-data (MPMD) paradigm, serving as an intermediary responsible for fast data exchange and interactions between the subsystems. The main goal of MiMiC is to avoid interfering with the underlying parallelization of the external programs, including the operability on hybrid architectures (e.g., CPU/GPU), and keep their setup and execution as close as possible to the original. At the moment, MiMiC offers an efficient implementation of electrostatic embedding quantum mechanics/molecular mechanics (QM/MM) that has demonstrated unprecedented parallel scaling in simulations of large biomolecules using CPMD and GROMACS as QM and MM engines, respectively. However, as it is designed for high flexibility with general multiscale models in mind, it can be straightforwardly extended beyond QM/MM. In this article, we illustrate the software design and the features of the framework, which make it a compelling choice for multiscale simulations in the upcoming era of exascale high-performance computing.

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.

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.

Sub-100-fs energy transfer in coenzyme NADH is a coherent process assisted by a charge-transfer state

Excitation energy transfer (EET) is a key photoinduced process in biological chromophoric assemblies. Here we investigate the factors which can drive EET into efficient ultrafast sub-ps regimes. We demonstrate how a coherent transport of electronic population could facilitate this in water solvated NADH coenzyme and uncover the role of an intermediate dark charge-transfer state. High temporal resolution ultrafast optical spectroscopy gives a 54±11 fs time constant for the EET process. Nonadiabatic quantum dynamical simulations computed through the time-evolution of multidimensional wavepackets suggest that the population transfer is mediated by photoexcited molecular vibrations due to strong coupling between the electronic states. The polar aqueous solvent environment leads to the active participation of a dark charge transfer state, accelerating the vibronically coherent EET process in favorably stacked conformers and solvent cavities. Our work demonstrates how the interplay of structural and environmental factors leads to diverse pathways for the EET process in flexible heterodimers and provides general insights relevant for coherent EET processes in stacked multichromophoric aggregates like DNA strands.

Multiscale biomolecular simulations in the exascale era

The complexity of biological systems and processes, spanning molecular to macroscopic scales, necessitates the use of multiscale simulations to get a comprehensive understanding. Quantum mechanics/molecular mechanics (QM/MM) molecular dynamics (MD) simulations are crucial for capturing processes beyond the reach of classical MD simulations. The advent of exascale computing offers unprecedented opportunities for scientific exploration, not least within life sciences, where simulations are essential to unravel intricate molecular mechanisms underlying biological processes. However, leveraging the immense computational power of exascale computing requires innovative algorithms and software designs. In this context, we discuss the current status and future prospects of multiscale biomolecular simulations on exascale supercomputers with a focus on QM/MM MD. We highlight our own efforts in developing a versatile and high-performance multiscale simulation framework with the aim of efficient utilization of state-of-the-art supercomputers. We showcase its application in uncovering complex biological mechanisms and its potential for leveraging exascale computing.

Computational Spectroscopy of Aqueous Solutions: The Underlying Role of Conformational Sampling

Fully atomistic multiscale polarizable quantum mechanics (QM)/molecular mechanics (MM) approaches, combined with techniques to sample the solute-solvent phase space, constitute the most accurate method to compute spectral signals in aqueous solution. Conventional sampling strategies, such as classical molecular dynamics (MD), may encounter drawbacks when the conformational space is particularly complex, and transition barriers between conformers are high. This can lead to inaccurate sampling, which can potentially impact the accuracy of spectral calculations. For this reason, in this work, we compare classical MD with enhanced sampling techniques, i.e., replica exchange MD and metadynamics. In particular, we show how the different sampling techniques affect computed UV, electronic circular dichroism, nuclear magnetic resonance shielding, and optical rotatory dispersion of N-acetylproline-amide in aqueous solution. Such a system is a model peptide characterized by complex conformational variability. Calculated values suggest that spectral properties are influenced by solute conformers, relative population, and solvent effects; therefore, particular care needs to be paid for when choosing the sampling technique.

Interoperable workflows by exchanging grid-based data between quantum-chemical program packages

Quantum-chemical subsystem and embedding methods require complex workflows that may involve multiple quantum-chemical program packages. Moreover, such workflows require the exchange of voluminous data that go beyond simple quantities, such as molecular structures and energies. Here, we describe our approach for addressing this interoperability challenge by exchanging electron densities and embedding potentials as grid-based data. We describe the approach that we have implemented to this end in a dedicated code, PyEmbed, currently part of a Python scripting framework. We discuss how it has facilitated the development of quantum-chemical subsystem and embedding methods and highlight several applications that have been enabled by PyEmbed, including wave-function theory (WFT) in density-functional theory (DFT) embedding schemes mixing non-relativistic and relativistic electronic structure methods, real-time time-dependent DFT-in-DFT approaches, the density-based many-body expansion, and workflows including real-space data analysis and visualization. Our approach demonstrates, in particular, the merits of exchanging (complex) grid-based data and, in general, the potential of modular software development in quantum chemistry, which hinges upon libraries that facilitate interoperability.

Ultrafast photochemistry and electron-diffraction spectra in n → (3s) Rydberg excited cyclobutanone resolved at the multireference perturbative level

We study the ultrafast time evolution of cyclobutanone excited to the singlet n → Rydberg state through non-adiabatic surface-hopping simulationsperformed at extended multi-state complete active space second-order perturbation (XMS-CASPT2) level of theory. These dynamics predict relaxation to the ground-state with a timescale of 822 ± 45 fs with minimal involvement of the triplets. The major relaxation path to the ground-state involves a three-state degeneracy region and leads to a variety of fragmented photoproducts. We simulate the resulting time-resolved electron-diffraction spectra, which track the relaxation of the excited state and the formation of various photoproducts in the ground state.

Multiexcitonic and optically bright states in subunits of pentacene crystals: A hybrid DFT/MRCI and molecular mechanics study

A hybrid quantum mechanics/molecular mechanics setup was used to model electronically excited pentacene in the crystal phase. Particularly interesting in the context of singlet fission (SF) is the energetic location of the antiferromagnetically coupled multiexcitonic singlet state, 1(TT), and the ferromagnetically coupled analog in relation to the optically bright singlet state. To provide photophysical properties of the accessible spin manifold, combined density functional theory and multi-reference configuration interaction calculations were performed on pentacene dimers and a trimer, electrostatically embedded in the crystal. The likelihood of a quintet intermediate in the SF process was estimated by computing singlet-quintet electron spin-spin couplings employing the Breit-Pauli Hamiltonian. The performance of the applied methods was assessed on the pentacene monomer. The character of the optically bright state and the energetic location of the 1(TT) state depend strongly on the relative orientation of the pentacene units. In the V-shaped dimers and in the trimer, the optically bright state is dominated by local and charge transfer (CT) excitations, with admixtures of doubly excited configurations. The CT excitations gain weight upon geometry relaxation, thus supporting a CT-mediated SF mechanism as the primary step of the SF process. For the slip-stacked dimer, the energetic order of the bright and the 1(TT) states swaps upon geometry relaxation, indicating strong nonadiabatic coupling close to the Franck-Condon region-a prerequisite for a coherent SF process. The multiexcitonic singlet, triplet, and quintet states are energetically too far apart and their spin-spin couplings are too small to bring about a noteworthy multiplicity mixing.

Potential energy interpolation with target-customized weighting coordinates: application to excited-state dynamics of photoactive yellow protein chromophore in water

Interpolation of potential energy surfaces (PESs) can provide a practical route to performing molecular dynamics simulations with a reliability matching a high-level quantum chemical calculation. An obstacle to its widespread use is perhaps the lack of general and optimal interpolation settings that can be applied in a black-box manner for any given molecular system. How to set up the weights for interpolation is one such task, and we still need to diversify the approaches in order to treat various systems. Here, we develop a new interpolation weighting scheme, which allows us to choose the weighting coordinates in a system-specific manner, by amplifying the contribution from specific internal coordinates. The new weighting scheme with an appropriate selection of coordinates is proved to be effective in reducing the interpolation error along the reaction pathway. As a demonstration, we consider the photoactive yellow protein chromophore system, as it constitutes itself as an interesting target that bears long-standing questions related to excited-state dynamics inside protein environments. We build its two-state diabatic interpolated PES with the new weighting scheme. We indeed see the utility of our scheme by conducting nonadiabatic molecular dynamics simulations with the required semi-global PES based on a limited number of data points.

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 byn π * 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.

Galván I, Rivalta I, Aquilante F, Garavelli M, Lindh R, Segarra-Martí J. 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.

Galván I, Alavi A, Aleotti F, Aquilante F, Autschbach J, Avagliano D, Baiardi A, Bao JJ, Battaglia S, Birnoschi L, Blanco-González A, Bokarev SI, Broer R, Cacciari R, Calio PB, Carlson RK, Carvalho Couto R, Cerdán L, Chibotaru LF, Chilton NF, Church JR, Conti I, Coriani S, Cuéllar-Zuquin J, Daoud RE, Dattani N, Decleva P, de Graaf C, Delcey M, De Vico L, Dobrautz W, Dong SS, Feng R, Ferré N, Filatov(Gulak) M, Gagliardi L, Garavelli M, González L, Guan Y, Guo M, Hennefarth MR, Hermes MR, Hoyer CE, Huix-Rotllant M, Jaiswal VK, Kaiser A, Kaliakin DS, Khamesian M, King DS, Kochetov V, Krośnicki M, Kumaar AA, Larsson ED, Lehtola S, Lepetit MB, Lischka H, López Ríos P, Lundberg M, Ma D, Mai S, Marquetand P, Merritt ICD, Montorsi F, Mörchen M, Nenov A, Nguyen VHA, Nishimoto Y, Oakley MS, Olivucci M, Oppel M, Padula D, Pandharkar R, Phung QM, Plasser F, Raggi G, Rebolini E, Reiher M, Rivalta I, Roca-Sanjuán D, Romig T, Safari AA, Sánchez-Mansilla A, Sand AM, Schapiro I, Scott TR, Segarra-Martí J, Segatta F, Sergentu DC, Sharma P, Shepard R, Shu Y, Staab JK, Straatsma TP, Sørensen LK, Tenorio BNC, Truhlar DG, Ungur L, Vacher M, Veryazov V, Voß TA, Weser O, Wu D, Yang X, Yarkony D, Zhou C, Zobel JP, Lindh R. 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.

NEXMD v2.0 Software Package for Nonadiabatic Excited State Molecular Dynamics Simulations

We present NEXMD version 2.0, the second release of the NEXMD (Nonadiabatic EXcited-state Molecular Dynamics) software package. Across a variety of new features, NEXMD v2.0 incorporates new implementations of two hybrid quantum-classical dynamics methods, namely, Ehrenfest dynamics (EHR) and the Ab-Initio Multiple Cloning sampling technique for Multiconfigurational Ehrenfest quantum dynamics (MCE-AIMC or simply AIMC), which are alternative options to the previously implemented trajectory surface hopping (TSH) method. To illustrate these methodologies, we outline a direct comparison of these three hybrid quantum-classical dynamics methods as implemented in the same NEXMD framework, discussing their weaknesses and strengths, using the modeled photodynamics of a polyphenylene ethylene dendrimer building block as a representative example. We also describe the expanded normal-mode analysis and constraints for both the ground and excited states, newly implemented in the NEXMD v2.0 framework, which allow for a deeper analysis of the main vibrational motions involved in vibronic dynamics. Overall, NEXMD v2.0 expands the range of applications of NEXMD to a larger variety of multichromophore organic molecules and photophysical processes involving quantum coherences and persistent couplings between electronic excited states and nuclear velocity.

Enhanced Solid-State Fluorescence of Flavin Derivatives by Incorporation in the Metal-Organic Frameworks MIL-53(Al) and MOF-5

The flavin derivatives 10-methyl-isoalloxazine (MIA) and 6-fluoro-10-methyl-isoalloxazine (6F-MIA) were incorporated in two alternative metal-organic frameworks, (MOFs) MIL-53(Al) and MOF-5. We used a post-synthetic, diffusion-based incorporation into microcrystalline MIL-53 powders with one-dimensional (1D) pores and an in-situ approach during the synthesis of MOF-5 with its 3D channel network. The maximum amount of flavin dye incorporation is 3.9 wt% for MIA@MIL-53(Al) and 1.5 wt% for 6F-MIA@MIL-53(Al), 0.85 wt% for MIA@MOF-5 and 5.2 wt% for 6F-MIA@MOF-5. For the high incorporation yields the probability to have more than one dye molecule in a pore volume is significant. As compared to the flavins in solution, the fluorescence spectrum of these flavin@MOF composites is broadened at the bathocromic side especially for MIA. Time-resolved spectroscopy showed that multi-exponential fluorescence lifetimes were needed to describe the decays. The fluorescence-weighted lifetime of flavin@MOF of 4 ± 1 ns also corresponds to those in solution but is significantly prolonged compared to the solid flavin dyes with less than 1 ns, thereby confirming the concept of "solid solutions" for dye@MOF composites. The fluorescence quantum yield (Φ_F) of the flavin@MOF composites is about half of the solution but is significantly higher compared to the solid flavin dyes. Both the fluorescence lifetime and quantum yield of flavin@MOF decrease with the flavin loading in MIL-53 due to the formation of various J-aggregates. Theoretical calculations using plane-wave and QM/MM methods are in good correspondence with the experimental results and explain the electronic structures as well as the photophysical properties of crystalline MIA and the flavin@MOF composites. In the solid flavins, π-stacking interactions of the molecules lead to a charge transfer state with low oscillator strength resulting in aggregation-caused quenching (ACQ) with low lifetimes and quantum yields. In the MOF pores, single flavin molecules represent a major population and the computed MIA@MOF structures do not find π-stacking interactions with the pore walls but only weak van-der-Waals contacts which reasons the enhanced fluorescence lifetime and quantum yield of the flavins in the composites compared to their neat solid state. To analyze the orientation of flavins in MOFs, we measured fluorescence anisotropy images of single flavin@MOF-5 crystals and a static ensemble flavin@MIL53 microcrystals, respectively. Based on image information, anisotropy distributions and overall curve of the time-resolved anisotropy curves combined with theoretical calculations, we can prove that all fluorescent flavins species have a defined and rather homogeneous orientation in the MOF framework. In MIL-53, the transition dipole moments of flavins are orientated along the 1D channel axis, whereas in MOF-5 we resolved an average orientation that is tilted with respect to the cubic crystal lattice. Notably, the more hydrophobic 6F-MIA exhibits a higher degree order than MIA. The flexible MOF MIL-53(Al) was optimized essentially to the experimental large-pore form in the guest-free state with QuantumEspresso (QE) and with MIA molecules in the pores the structure contracted to close to the experimental narrow-pore form which was also confirmed by PXRD. In summary, the incorporation of flavins in MOFs yields solid-state materials with enhanced rigidity, stabilized conformation, defined orientation and reduced aggregations of the flavins, leading to increased fluorescence lifetime and quantum yield as controllable photo-luminescent and photo-physical properties.

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.

Newton-X Platform: New Software Developments for Surface Hopping and Nuclear Ensembles

Newton-X is an open-source computational platform to perform nonadiabatic molecular dynamics based on surface hopping and spectrum simulations using the nuclear ensemble approach. Both are among the most common methodologies in computational chemistry for photophysical and photochemical investigations. This paper describes the main features of these methods and how they are implemented in Newton-X. It emphasizes the newest developments, including zero-point-energy leakage correction, dynamics on complex-valued potential energy surfaces, dynamics induced by incoherent light, dynamics based on machine-learning potentials, exciton dynamics of multiple chromophores, and supervised and unsupervised machine learning techniques. Newton-X is interfaced with several third-party quantum-chemistry programs, spanning a broad spectrum of electronic structure methods.

INAQS, a Generic Interface for Nonadiabatic QM/MM Dynamics: Design, Implementation, and Validation for GROMACS/Q-CHEM simulations

The accurate description of large molecular systems in complex environments remains an ongoing challenge for the field of computational chemistry. This problem is even more pronounced for photoinduced processes, as multiple excited electronic states and their corresponding nonadiabatic couplings must be taken into account. Multiscale approaches such as hybrid quantum mechanics/molecular mechanics (QM/MM) offer a balanced compromise between accuracy and computational burden. Here, we introduce an open-source software package (INAQS) for nonadiabatic QM/MM simulations that bridges the sampling capabilities of the GROMACS MD package and the excited-state infrastructure of the Q-CHEM electronic structure software. The interface is simple and can be adapted easily to other MD codes. The code supports a variety of different trajectory-based molecular dynamics, ranging from Born-Oppenheimer to surface hopping dynamics. To illustrate the power of this combination, we simulate electronic absorption spectra, free-energy surfaces along a reaction coordinate, and the excited-state dynamics of 1,3-cyclohexadiene in solution.

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.

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 L_a → L_b transfer through a conical intersection and a subsequent 220 fs reverse L_b → L_a transfer due to solvent-assisted adiabatic stabilization of the L_a 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 L_a 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.

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'.

Time-Resolved Optical Pump-Resonant X-ray Probe Spectroscopy of 4-Thiouracil: A Simulation Study

We theoretically monitor the photoinduced ππ* → nπ* internal conversion process in 4-thiouracil (4TU), triggered by an optical pump. The element-sensitive spectroscopic signatures are recorded by a resonant X-ray probe tuned to the sulfur, oxygen, or nitrogen K-edge. We employ high-level electronic structure methods optimized for core-excited electronic structure calculation combined with quantum nuclear wavepacket dynamics computed on two relevant nuclear modes, fully accounting for their quantum nature of nuclear motions. We critically discuss the capabilities and limitations of the resonant technique. For sulfur and nitrogen, we document a pre-edge spectral window free from ground-state background and rich with ππ* and nπ* absorption features. The lowest sulfur K-edge shows strong absorption for both ππ* and nπ*. In the lowest nitrogen K-edge window, we resolve a state-specific fingerprint of the ππ* and an approximate timing of the conical intersection via its depletion. A spectral signature of the nπ* transition, not accessible by UV-vis spectroscopy, is identified. The oxygen K-edge is not sensitive to molecular deformations and gives steady transient absorption features without spectral dynamics. The ππ*/nπ* coherence information is masked by more intense contributions from populations. Altogether, element-specific time-resolved resonant X-ray spectroscopy provides a detailed picture of the electronic excited-state dynamics and therefore a sensitive window into the photophysics of thiobases.

Fewest switches surface hopping with Baeck-An couplings

In the Baeck-An (BA) approximation, first-order nonadiabatic coupling vectors are given in terms of adiabatic energy gaps and the second derivative of the gaps with respect to the coupling coordinate. In this paper, a time-dependent (TD) BA approximation is derived, where the couplings are computed from the energy gaps and their second time-derivatives. TD-BA couplings can be directly used in fewest switches surface hopping, enabling nonadiabatic dynamics with any electronic structure methods able to provide excitation energies and energy gradients. Test results of surface hopping with TD-BA couplings for ethylene and fulvene show that the TD-BA approximation delivers a qualitatively correct picture of the dynamics and a semiquantitative agreement with reference data computed with exact couplings. Nevertheless, TD-BA does not perform well in situations conjugating strong couplings and small velocities. Considered the uncertainties in the method, TD-BA couplings could be a competitive approach for inexpensive, exploratory dynamics with a small trajectories ensemble. We also assessed the potential use of TD-BA couplings for surface hopping dynamics with time-dependent density functional theory (TDDFT), but the results are not encouraging due to singlet instabilities near the crossing seam with the ground state.

Fewest switches surface hopping with Baeck-An couplings

In the Baeck-An (BA) approximation, first-order nonadiabatic coupling vectors are given in terms of adiabatic energy gaps and the second derivative of the gaps with respect to the coupling coordinate. In this paper, a time-dependent (TD) BA approximation is derived, where the couplings are computed from the energy gaps and their second time-derivatives. TD-BA couplings can be directly used in fewest switches surface hopping, enabling nonadiabatic dynamics with any electronic structure methods able to provide excitation energies and energy gradients. Test results of surface hopping with TD-BA couplings for ethylene and fulvene show that the TD-BA approximation delivers a qualitatively correct picture of the dynamics and a semiquantitative agreement with reference data computed with exact couplings. Nevertheless, TD-BA does not perform well in situations conjugating strong couplings and small velocities. Considered the uncertainties in the method, TD-BA couplings could be a competitive approach for inexpensive, exploratory dynamics with a small trajectories ensemble. We also assessed the potential use of TD-BA couplings for surface hopping dynamics with time-dependent density functional theory (TDDFT), but the results are not encouraging due to singlet instabilities near the crossing seam with the ground state.

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...

Tracking excited state decay mechanisms of pyrimidine nucleosides in real time

DNA owes its remarkable photostability to its building blocks-the nucleosides-that efficiently dissipate the energy acquired upon ultraviolet light absorption. The mechanism occurring on a sub-picosecond time scale has been a matter of intense debate. Here we combine sub-30-fs transient absorption spectroscopy experiments with broad spectral coverage and state-of-the-art mixed quantum-classical dynamics with spectral signal simulations to resolve the early steps of the deactivation mechanisms of uridine (Urd) and 5-methyluridine (5mUrd) in aqueous solution. We track the wave packet motion from the Franck-Condon region to the conical intersections (CIs) with the ground state and observe spectral signatures of excited-state vibrational modes. 5mUrd exhibits an order of magnitude longer lifetime with respect to Urd due to the solvent reorganization needed to facilitate bulky methyl group motions leading to the CI. This activates potentially lesion-inducing dynamics such as ring opening. Involvement of the 1nπ* state is found to be negligible.

Unified Description of Ultrafast Excited State Decay Processes in Epigenetic Deoxycytidine Derivatives

Epigenetic DNA modifications play a fundamental role in modulating gene expression and regulating cellular and developmental biological processes, thereby forming a second layer of information in DNA. The epigenetic 2'-deoxycytidine modification 5-methyl-2'-deoxycytidine, together with its enzymatic oxidation products (5-hydroxymethyl-2'-deoxycytidine, 5-formyl-2'-deoxycytidine, and 5-carboxyl-2'-deoxycytidine), are closely related to deactivation and reactivation of DNA transcription. Here, we combine sub-30-fs transient absorption spectroscopy with high-level correlated multiconfigurational CASPT2/MM computational methods, explicitly including the solvent, to obtain a unified picture of the photophysics of deoxycytidine-derived epigenetic DNA nucleosides. We assign all the observed time constants and identify the excited state relaxation pathways, including the competition of intersystem crossing and internal conversion for 5-formyl-2'-deoxycytidine and ballistic decay to the ground state for 5-carboxy-2'-deoxycytidine. Our work contributes to shed light on the role of epigenetic derivatives in DNA photodamage as well as on their possible therapeutic use.

Computational spectroscopy of complex systems

Numerous linear and non-linear spectroscopic techniques have been developed to elucidate structural and functional information of complex systems ranging from natural systems, such as proteins and light-harvesting systems, to synthetic systems, such as solar cell materials and light-emitting diodes. The obtained experimental data can be challenging to interpret due to the complexity and potential overlapping spectral signatures. Therefore, computational spectroscopy plays a crucial role in the interpretation and understanding of spectral observables of complex systems. Computational modeling of various spectroscopic techniques has seen significant developments in the past decade, when it comes to the systems that can be addressed, the size and complexity of the sample types, the accuracy of the methods, and the spectroscopic techniques that can be addressed. In this Perspective, I will review the computational spectroscopy methods that have been developed and applied for infrared and visible spectroscopies in the condensed phase. I will discuss some of the questions that this has allowed answering. Finally, I will discuss current and future challenges and how these may be addressed.

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 S₁ to T₁ in solution is provided.

Velocity Adjustment in Surface Hopping: Ethylene as a Case Study of the Maximum Error Caused by Direction Choice

The most common surface hopping dynamics algorithms require velocity adjustment after hopping to ensure total-energy conservation. Based on the semiclassical analysis, this adjustment must be made parallel to the nonadiabatic coupling vector's direction. Nevertheless, this direction is not always known, and the common practice has been to adjust the velocity in either the linear momentum or velocity directions. This paper benchmarks surface hopping dynamics of photoexcited ethylene with velocity adjustment in several directions, including those of the nonadiabatic coupling vector, the momentum, and the energy gradient difference. It is shown that differences in time constants and structural evolution fall within the statistical uncertainty of the method considering up to 500 trajectories in each dynamics set, rendering the three approaches statistically equivalent. For larger ensembles beyond 1000 trajectories, significant differences between the results arise, limiting the validity of adjustment in alternative directions. Other possible adjustment directions (velocity, single-state gradients, angular momentum) are evaluated as well. Given the small size of ethylene, the results reported in this paper should be considered an upper limit for the error caused by the choice of the velocity-adjustment direction on surface hopping dynamics.

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).

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.

Photoinduced Forward and Backward Pedalo-Type Motion of a Molecular Switch

Photoresponsive molecular switches enable spatial and temporal control of molecular processes and are therefore crucial for the development of smart functional materials. Because the light-induced dynamics of these switching units are at the core of the resulting functionality, a detailed insight into their structural time evolution is fundamental for molecular embedding. Here, we performed a hybrid quantum mechanics (CASPT2 and TDDFT)/molecular mechanics (QM/MM) study to elucidate the photodynamics of an azodicarboxamide-based molecular switch, which is a promising candidate for implementation in highly dense environments such as polymers. In particular, we report a detailed picture of the molecular motion at the atomic level based on a relevant number of excited-state trajectories. We show that the azodicarboxamide-based molecular switch undergoes both a forward and backward pedalo-type motion upon excitation. Trans-cis photoisomerization on the other hand, which is well-known to occur for other azo-based chromophores, is shown to be a negligible pathway. By validating the volume-conserving pedalo-type motion, we provide a rational basis for the design of novel types of photoresponsive functional materials in which the active component must operate in a confined space.

Modern quantum chemistry with [Open]Molcas

MOLCAS/OpenMolcas is an ab initio electronic structure program providing a large set of computational methods from Hartree-Fock and density functional theory to various implementations of multiconfigurational theory. This article provides a comprehensive overview of the main features of the code, specifically reviewing the use of the code in previously reported chemical applications as well as more recent applications including the calculation of magnetic properties from optimized density matrix renormalization group wave functions.

Non-adiabatic Excited-State Molecular Dynamics: Theory and Applications for Modeling Photophysics in Extended Molecular Materials

Optically active molecular materials, such as organic conjugated polymers and biological systems, are characterized by strong coupling between electronic and vibrational degrees of freedom. Typically, simulations must go beyond the Born-Oppenheimer approximation to account for non-adiabatic coupling between excited states. Indeed, non-adiabatic dynamics is commonly associated with exciton dynamics and photophysics involving charge and energy transfer, as well as exciton dissociation and charge recombination. Understanding the photoinduced dynamics in such materials is vital to providing an accurate description of exciton formation, evolution, and decay. This interdisciplinary field has matured significantly over the past decades. Formulation of new theoretical frameworks, development of more efficient and accurate computational algorithms, and evolution of high-performance computer hardware has extended these simulations to very large molecular systems with hundreds of atoms, including numerous studies of organic semiconductors and biomolecules. In this Review, we will describe recent theoretical advances including treatment of electronic decoherence in surface-hopping methods, the role of solvent effects, trivial unavoided crossings, analysis of data based on transition densities, and efficient computational implementations of these numerical methods. We also emphasize newly developed semiclassical approaches, based on the Gaussian approximation, which retain phase and width information to account for significant decoherence and interference effects while maintaining the high efficiency of surface-hopping approaches. The above developments have been employed to successfully describe photophysics in a variety of molecular materials.

First-principles characterization of the singlet excited state manifold in DNA/RNA nucleobases

TD-DFT characterization of the high-energy singlet excited state manifold of the canonical DNA/RNA nucleobasesin vacuumis assessed against RASPT2 reference computations for reliable simulations of linear and non-linear electronic spectra.

Modeling multidimensional spectral lineshapes from first principles: application to water-solvated adenine

We theoretically describe spectral lineshape from first principles, providing insight into solvent–solute interactions in terms of static and dynamic disorder and how these shape experimental signals in linear and non-linear optical spectroscopies.

Intersystem crossing in thiobases proceeds by a dark intermediate state

4-thiouracil (4TU) is studied by transient absorption spectroscopy employing sub-20 fs UV-pulses and hybrid QM(CASPT2) / MM(AMBER) computations (static and dynamic), evidencing that, along the photoexcited relaxation pathway, intersystem crossing originates from a dark intermediate state.

UV-light induced vibrational coherences explain Kasha rule violation in frans-azobenzene

Sub-20-fs transient absorption spectroscopy and simulations show that CNN-bendings dominate the sub-ps dynamics of ππ*-excited trans-azobenzene, thereby driving the system to the ground state through a non-productive decay channel in violation of the Kasha rule.

Two-dimensional UV spectroscopy: a new insight into the structure and dynamics of biomolecules

Two-dimensional ultraviolet spectroscopy has the potential to deliver rich structural and dynamical information on biomolecules such as DNA and proteins.

Impact of fluorination on the photophysics of the flavin chromophore: a quantum chemical perspective

10-Methylisoalloxazine (MIA) and its fluorinated derivatives (6-9F-MIA) were investigated by means of quantum chemistry, looking into the influence of fluorination on fluorescence, absorption and inter-system crossing (ISC) in vacuum and in aqueous solution.

Conical intersection dynamics of pyrimidine nucleosides tracked with sub-20-fs UV pulses

By combining transient absorption spectroscopy with sub-20-fs UV pulses and ab initio numerical simulations we follow the ultrafast dynamics in pyrimidine nucleosides and visualize the passage through conical intersections presiding excited state deactivation.