Efficient Polarizable QM/MM Using the Direct Reaction Field Hamiltonian with Electrostatic Potential Fitted Multipole Operators.

Electronic polarization and dispersion are decisive actors in determining interaction energies between molecules. These interactions have a particularly profound effect on excitation energies of molecules in complex environments, especially when the excitation involves a significant degree of charge reorganization. The direct reaction field (DRF) approach, which has seen a recent revival of interest, provides a powerful framework for describing these interactions in quantum mechanics/molecular mechanics (QM/MM) models of systems, where a small subsystem of interest is described using quantum chemical methods and the remainder is treated with a simple MM force field.

Investigating the Origin of Automatic Rhodopsin Modeling Outliers Using the Microbial Rhodopsin as Testbed.

The automatic rhodopsin modeling (ARM) approach is a computational workflow devised for the automatic buildup of hybrid quantum mechanics/molecular mechanics (QM/MM) models of wild-type rhodopsins and mutants, with the purpose of establishing trends in their photophysical and photochemical properties. Despite the success of ARM in accurately describing the visible light absorption maxima of many rhodopsins, for a few cases, called outliers, it might lead to large deviations with respect to experiments. Applying ARM to rhodopsin (GR), a microbial rhodopsin with important applications in optogenetics, we analyze the origin of such outliers in the absorption energies obtained for GR wild-type and mutants at neutral pH, with a total root-mean-square deviation (RMSD) of 0.

Temperature Controlled Decay and Pendulum Dynamics of Green Fluorescent Protein (GFP) Chromophore.

The excited-state dynamics of the GFP chromophore, HBDI (anionic -hydroxybenzylidene-2,3-dimethylimidazolinone), were investigated through a combination of theoretical nonadiabatic molecular dynamics (NAMD) simulations and femtosecond transient absorption spectroscopy (fs-TA). The NAMD simulations revealed that the primary dynamics in excited states involve the formation of a P-twisted intermediate (S), which undergoes pendulum-like oscillations with respect to ϕ = 90°. This motion serves as a reservoir for the excited-state population and the primary source of fluorescence.

α-Trifluoromethylated Quinolines as Safe and Storable PET-Donor for Radical Polymerizations.

In the quest for powerful, safe, and storable photoinduced-electron transfer (PET) donors, the attention is turned to the α-trihalomethylated amine moiety that is not studied in the context of PET-reductants. The thermal and photophysical properties of α-trifluoromethylated quinolines are thus studied and their reducing abilities evaluated as initiators of polymerization reactions. Polymers of high molecular weights are obtained through a radical polymerization process and the PET-donor can be stored within the monomer for several months without losing its efficiency.

OpenQP: A Quantum Chemical Platform Featuring MRSF-TDDFT with an Emphasis on Open-Source Ecosystem.

The (Open Quantum Platform) is a new open-source quantum chemistry library developed to tackle sustainability and interoperability challenges in the field of computational chemistry. provides various popular quantum chemical theories as autonomous such as energy and gradient calculations of HF, DFT, TDDFT, SF-TDDFT, and MRSF-TDDFT, thereby allowing easy interconnection with third-party software. A scientifically notable feature is the innovative mixed-reference spin-flip time-dependent density functional theory (MRSF-TDDFT) and its customized exchange-correlation functionals such as the DTCAM series of VAEE, XI, XIV, AEE, and VEE, which significantly expand the applicability scope of DFT and TDDFT.

Toward Consistent Predictions of Core/Valence Ionization Potentials and Valence Excitation Energies by MRSF-TDDFT.

Optimizing exchange-correlation functionals for both core/valence ionization potentials (cIPs/vIPs) and valence excitation energies (VEEs) at the same time in the framework of MRSF-TDDFT is self-contradictory. To overcome the challenge, within the previous "adaptive exact exchange" or double-tuning strategy on Coulomb-attenuating XC functionals (CAM), a new XC functional specifically for cIPs and vIPs was first developed by enhancing exact exchange to both short- and long-range regions. The resulting DTCAM-XI functional achieved remarkably high accuracy in its predictions with errors of less than half eV.

Analytic Gradients for the Electrostatic Embedding QM/MM Model in Periodic Boundary Conditions Using Particle-Mesh Ewald Sums and Electrostatic Potential Fitted Charge Operators.

Long-range electrostatic effects are fundamental for describing chemical reactivity in the condensed phase. Here, we present the methodology of an efficient quantum mechanical/molecular mechanical (QM/MM) model in periodic boundary conditions (PBC) compatible with QM/MM boundaries at chemical bonds. The method combines electrostatic potential fitted charge operators and electrostatic potentials derived from the smooth particle-mesh Ewald (PME) sum approach.

Influence of pump laser fluence on ultrafast myoglobin structural dynamics.

High-intensity femtosecond pulses from an X-ray free-electron laser enable pump-probe experiments for the investigation of electronic and nuclear changes during light-induced reactions. On timescales ranging from femtoseconds to milliseconds and for a variety of biological systems, time-resolved serial femtosecond crystallography (TR-SFX) has provided detailed structural data for light-induced isomerization, breakage or formation of chemical bonds and electron transfer. However, all ultrafast TR-SFX studies to date have employed such high pump laser energies that nominally several photons were absorbed per chromophore.

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.

Solvent Effects and pH Dependence of the X-ray Absorption Spectra of Proline from Electrostatic Embedding Quantum Mechanics/Molecular Mechanics and Mixed-Reference Spin-Flip Time-dependent Density-Functional Theory.

The accurate description of solvent effects on X-ray absorption spectra (XAS) is fundamental for comparing the simulated spectra with experiments in solution. Currently, few protocols exist that can efficiently reproduce the effects of the solute/solvent interactions on XAS. Here, we develop an efficient and accurate theoretical protocol for simulating the solvent effects on XAS.

Assessment of the Electron Correlation Treatment on the Quantum-Classical Dynamics of Retinal Protonated Schiff Base Models: XMS-CASPT2, RMS-CASPT2, and REKS Methods.

We compare the performance of three different multiconfigurational wave function-based electronic structure methods and two implementations of the spin-restricted ensemble-referenced Kohn-Sham (REKS) method. The study is characterized by three features: (i) it uses a small set of quantum-classical trajectories rather than potential energy surface mapping, (ii) it focuses, exclusively, on the photoisomerization of retinal protonated Schiff base models, and (iii) it probes the effect of both methyl substitution and the increase in length of the conjugate π-system. For each tested method, the corresponding analytical gradients are used to drive the quantum-classical (Tully's FSSH method) trajectory propagation, including the recent multistate XMS-CASPT2 and RMS-CASPT2 gradients.

Automatic Rhodopsin Modeling with Multiple Protonation Microstates.

Automatic Rhodopsin Modeling (ARM) is a simulation protocol providing QM/MM models of rhodopsins capable of reproducing experimental electronic absorption and emission trends. Currently, ARM is restricted to a single protonation microstate for each rhodopsin model. Herein, we incorporate an extension of the minimal electrostatic model (MEM) into the ARM protocol to account for all relevant protonation microstates at a given pH.

Doubly Tuned Exchange-Correlation Functionals for Mixed-Reference Spin-Flip Time-Dependent Density Functional Theory.

It is demonstrated that significant accuracy improvements in MRSF-TDDFT can be achieved by introducing two different exchange-correlation (XC) functionals for the reference Kohn-Sham DFT and the response part of the calculations, respectively. Accordingly, two new XC functionals of doubly tuned Coulomb attenuated method-vertical excitation energy (DTCAM-VEE) and DTCAM-AEE were developed on the basis of the "adaptive exact exchange (AEE)" concept in the framework of the Coulomb-attenuating XC functionals. The values by DTCAM-VEE are in excellent agreement with those of Thiel's set [mean absolute errors (MAEs) and the interquartile range (IQR) values of 0.

Ultrafast Spin Crossover Photochemical Mechanism in [Fe(2,2'-bipyridine)]] Revealed by Quantum Dynamics.

Photoexcitation of [Fe(2,2'-bipyridine)] induces a subpicosecond spin crossover transformation from a low-spin singlet to a high-spin quintet state. The mechanism involves metal-centered (MC) and metal-ligand charge transfer (MLCT) triplet intermediates, but their individual contributions to this efficient intersystem crossing have been object of debate. Employing quantum wavepacket dynamics, we show that MC triplets are catalyzing the transfer to the high-spin state.

The OpenMolcas : 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.

An efficient electrostatic embedding QM/MM method using periodic boundary conditions based on particle-mesh Ewald sums and electrostatic potential fitted charge operators.

Hybrid quantum mechanics/molecular mechanics (QM/MM) models are successful at describing the properties and reactivity of biological macromolecules. Combining ab initio QM/MM methods and periodic boundary conditions (PBC) is currently the optimal approach for modeling chemical processes in an infinite environment, but frequently, these models are too time-consuming for general applicability to biological systems in a solution. Here, we define a simple and efficient electrostatic embedding QM/MM model in PBC, combining the benefits of electrostatic potential fitted atomic charges and particle-mesh Ewald sums, which can efficiently treat systems of an arbitrary size at a reasonable computational cost.

Photochemistry of Thymine in Solution and DNA Revealed by an Electrostatic Embedding QM/MM Combined with Mixed-Reference Spin-Flip TDDFT.

The photochemistry of nucleobases, important for their role as building blocks of DNA, is largely affected by the electrostatic environment in which they are soaked. For example, despite the numerous studies of thymine in solution and DNA, there is still a debate on the photochemical deactivation pathways after UV absorption. Many theoretical models are oversimplified due to the lack of computationally accurate and efficient electronic structure methodologies that capture excited state electron correlation effects when nucleobases are embedded in large electrostatic media.

Towards the engineering of a photon-only two-stroke rotary molecular motor.

The rational engineering of photoresponsive materials, e.g., light-driven molecular motors, is a challenging task.

On the fluorescence enhancement of arch neuronal optogenetic reporters.

The lack of a theory capable of connecting the amino acid sequence of a light-absorbing protein with its fluorescence brightness is hampering the development of tools for understanding neuronal communications. Here we demonstrate that a theory can be established by constructing quantum chemical models of a set of Archaerhodopsin reporters in their electronically excited state. We found that the experimentally observed increase in fluorescence quantum yield is proportional to the computed decrease in energy difference between the fluorescent state and a nearby photoisomerization channel leading to an exotic diradical of the protein chromophore.

Stressing the differences in alizarin and purpurin dyes through UV-visible light absorption and H-NMR spectroscopies.

Three anthraquinone-based chromophores (9,10-anthraquinone, alizarin, purpurin) are compared from the point of view of their experimental and computed NMR and UV-visible light absorption spectra. Using a hybrid (explicit/implicit) solvent model, each proton chemical shift can be reproduced with an error of less than 7%, even when such protons are engaged in inter-molecular hydrogen bonds with the solvent or when the analyzed sample contains a significant amount of impurities, for instance, 9,10-anthraquinone in purpurin. All the steady-state UV-visible absorption spectra feature a significant vibrational progression in the first absorption band.

Quantum dynamics simulations of the thermal and light-induced high-spin to low-spin relaxation in Fe(bpy) and Fe(mtz).

First row transition metal complexes with d to d electronic configurations exhibit spin-crossover (SCO), which can be induced by external stimuli, such as temperature, pressure and light. The low-spin to high-spin transition has been widely studied, but very little is known about the reverse process. Here, we present a theoretical study of thermal and light-induced high-to-low spin crossover in prototypical Fe(II) complexes.

Ultrafast Intersystem Crossing in Xanthone from Wavepacket Dynamics.

Most aromatic ketones containing first-row elements undergo unexpectedly fast intersystem crossing in a few tens of picoseconds and a quantum yield close to unity. Among them, xanthone (9-xanthen-9-one) possesses one of the fastest singlet-triplet rates of only ∼1.5 ps.

Photochemistry of Thymine in Protic Polar Nanomeric Droplets Using Electrostatic Embeding TD-DFT/MM.

Thymine photochemistry is important for understanding DNA photodamage. In the gas phase, thymine undergoes a fast non-radiative decay from S2 to S1. In the S1 state, it gets trapped for several picoseconds until returning to the ground-state S0.