Multistate, Polarizable QM/MM Embedding Scheme Based on the Direct Reaction Field Method: Solvatochromic Shifts, Analytical Gradients and Optimizations of Conical Intersections in Solution

Solvent Effects on Nonadiabatic Dynamics: Ab Initio Multiple Spawning Propagated on CASPT2/xTB Potentials

An approach to simulate nonadiabatic dynamics in solution is introduced, which relies on the propagation of the nuclear wavepacket with the Ab Initio Multiple Spawning (AIMS) method under the effect of potential energy calculated with a hybrid but fully quantum mechanical scheme (QM/QM'). The electronic energies of the excited states of the chromophore are calculated with multireference perturbation theory (CASPT2), and the embedding molecules are described with a tight binding Hamiltonian (GFN2-xTB). This implementation is fully open source and relies on the combination of PySpawn, OpenMolcas, and xTB. Additionally, ORCA is used to properly generate the initial conditions in solution, showing how the combination of cutting-edge implementations in several commonly used software can push the state of the art of nonadiabatic dynamics in solution toward a new high standard of accuracy. The dynamics of ethylene in vacuum, in acetone, and in chloroform is reported as a test case, with a detailed analysis of the AIMS runs that shows important geometrical and electronic effects of the solvents on the decay mechanism of the chromophore.

Analytic First-Order Derivatives of CASPT2 Combined with the Polarizable Continuum Model

The complete active space second-order perturbation theory (CASPT2) is valuable for accurately predicting electronic structures and transition energies. However, optimizing molecular geometries in the solution phase has proven challenging. In this study, we develop analytic first-order derivatives of CASPT2 using an implicit solvation model, specifically the polarizable continuum model, within the open-source package OpenMolcas. Analytic gradients and nonadiabatic coupling vectors are computed by solving a modified Z-vector equation. Comparisons with existing theoretical and experimental results demonstrate that the solvent effects can be qualitatively captured using the developed method.

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. In this paper we show how the DRF approach can be combined with the electrostatic potential fitted (ESPF) multipole operator description of the QM region charge density, which significantly improves the efficiency of the method, particularly for large MM systems, and for typical calculations effectively eliminates the dependence on MM system size. We also show how the DRF approach can be combined with fluctuating charge descriptions of the polarizable environment, as well as previously used atom-centered dipole-polarizability based models. We further show that the ESPF-DRF method provides an accurate description of molecular interactions in both ground and excited electronic states of the QM system and apply it to predict the gas to aqueous solution solvatochromic shifts in the UV/visible absorption spectrum of acrolein.

Fully Polarizable Multiconfigurational Self-Consistent Field/Fluctuating Charges Approach

A multiscale model based on the coupling of the multiconfigurational self-consistent field (MCSCF) method and the classical atomistic polarizable fluctuating charges (FQ) force field is presented. The resulting MCSCF/FQ approach is validated by exploiting the CASSCF scheme through application to compute vertical excitation energies of formaldehyde and para-nitroaniline in aqueous solution. The procedure is integrated with molecular dynamics simulations to capture the solute's conformational changes and the dynamic aspects of solvation. Comparative analysis with alternative solvent models, gas-phase calculations, and experimental data provides insights into the model's accuracy in reproducing solute-solvent molecular interactions and spectral signals.

Perspective on Theoretical and Experimental Advances in Atmospheric Photochemistry

Research that explores the chemistry of Earth's atmosphere is central to the current understanding of global challenges such as climate change, stratospheric ozone depletion, and poor air quality in urban areas. This research is a synergistic combination of three established domains: earth observation, for example, using satellites, and in situ field measurements; computer modeling of the atmosphere and its chemistry; and laboratory measurements of the properties and reactivity of gas-phase molecules and aerosol particles. The complexity of the interconnected chemical and photochemical reactions which determine the composition of the atmosphere challenges the capacity of laboratory studies to provide the spectroscopic, photochemical, and kinetic data required for computer models. Here, we consider whether predictions from computational chemistry using modern electronic structure theory and nonadiabatic dynamics simulations are becoming sufficiently accurate to supplement quantitative laboratory data for wavelength-dependent absorption cross-sections, photochemical quantum yields, and reaction rate coefficients. Drawing on presentations and discussions from the CECAM workshop on Theoretical and Experimental Advances in Atmospheric Photochemistry held in March 2024, we describe key concepts in the theory of photochemistry, survey the state-of-the-art in computational photochemistry methods, and compare their capabilities with modern experimental laboratory techniques. From such considerations, we offer a perspective on the scope of computational (photo)chemistry methods based on rigorous electronic structure theory to become a fourth core domain of research in atmospheric chemistry.

PyDFT-QMMM: A modular, extensible software framework for DFT-based QM/MM molecular dynamics

PyDFT-QMMM is a Python-based package for performing hybrid quantum mechanics/molecular mechanics (QM/MM) simulations at the density functional level of theory. The program is designed to treat short-range and long-range interactions through user-specified combinations of electrostatic and mechanical embedding procedures within periodic simulation domains, providing necessary interfaces to external quantum chemistry and molecular dynamics software. To enable direct embedding of long-range electrostatics in periodic systems, we have derived and implemented force terms for our previously described QM/MM/PME approach [Pederson and McDaniel, J. Chem. Phys. 156, 174105 (2022)]. Communication with external software packages Psi4 and OpenMM is facilitated through Python application programming interfaces (APIs). The core library contains basic utilities for running QM/MM molecular dynamics simulations, and plug-in entry-points are provided for users to implement custom energy/force calculation and integration routines, within an extensible architecture. The user interacts with PyDFT-QMMM primarily through its Python API, allowing for complex workflow development with Python scripting, for example, interfacing with PLUMED for free energy simulations. We provide benchmarks of forces and energy conservation for the QM/MM/PME and alternative QM/MM electrostatic embedding approaches. We further demonstrate a simple example use case for water solute in a water solvent system, for which radial distribution functions are computed from 100 ps QM/MM simulations; in this example, we highlight how the solvation structure is sensitive to different basis-set choices due to under- or over-polarization of the QM water molecule's electron density.

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.

Environmentally Driven Symmetry Breaking Quenches Dual Fluorescence in Proflavine

Nonadiabatic couplings between several electronic excited states are ubiquitous in many organic chromophores and can significantly influence optical properties. A recent experimental study demonstrated that the proflavine molecule exhibits surprising dual fluorescence in the gas phase, which is suppressed in polar solvent environments. Here, we uncover the origin of this phenomenon by parametrizing a linear-vibronic coupling Hamiltonian from spectral densities of system-bath coupling constructed along molecular dynamics trajectories, fully accounting for interactions with the condensed-phase environment. The finite-temperature absorption, steady-state emission, and time-resolved emission spectra are then computed using powerful, numerically exact tensor network approaches. We find that the dual fluorescence in vacuum is driven by a single well-defined coupling mode but is quenched in solution due to dynamic solvent-driven symmetry breaking that mixes the two low-lying electronic states. We expect the computational framework developed here to be widely applicable to the study of non-Condon effects in complex condensed-phase environments.