Multiconfigurational Quantum Chemistry Determinations of Absorption Cross Sections (σ) in the Gas Phase and Molar Extinction Coefficients (ε) in Aqueous Solution and Air-Water Interface

Photoreduction of mercuric bromides in polar ice

In the polar regions, which are vulnerable receptors of mercury pollution, atmospheric mercury depletion events (AMDEs) efficiently convert elemental mercury (Hg(0)) into oxidized mercury (Hg(II)) via bromine oxidation. Hg(II) subsequently deposits onto snow and sea ice. While field observations have shown that a large percentage of deposited mercury is re-emitted from the ice to the atmosphere by a photoinduced process, the fundamental photochemistry that drives the re-emission process remains unknown. Here, using multiconfigurational quantum chemistry, we find that the photoreduction of HgBr2, HgBr3-, and HgBr42- in ice is more efficient than in the gas phase. This results from the influence of water molecules on the molecular geometry and electronic structure of mercuric bromides in ice, which enhances the absorption intensities at wavelengths relevant in the troposphere (λ > 290 nm), as compared to gas phase. Kinetic modeling shows that ~30 to 60% of deposited mercury in AMDEs can be reemitted due to the photoreduction of mercuric bromides in ice, in agreement with field observations. Our results reveal a photoreduction mechanism of sunlight-induced excited state chemistry of mercuric bromides on ice. These findings strongly suggest that this chemistry should be incorporated into atmospheric models to account for ice-atmosphere mercury cycling in the polar environments, currently not considered.

easyPARM: Automated, Versatile, and Reliable Force Field Parameters for Metal-Containing Molecules with Unique Labeling of Coordinating Atoms

The dynamics of metal centers are challenging to describe due to the vast variety of ligands, metals, and coordination spheres, hampering the existence of general databases of transferable force field parameters for classical molecular dynamics simulations. Here, we present easyPARM, a Python-based tool that can calculate force field parameters for a wide range of metal complexes from routine frequency calculations with electronic structure methods. The approach is based on a unique labeling strategy, in which each ligand atom that coordinates the metal receives a unique atom type. This design prevents parameter shortage, labeling duplication, and the necessity to post-process output files, even for very complicated coordination spheres, whose parametrization process remain automatic. The program requires the Cartesian Hessian matrix, the geometry xyz file, and the atomic charges to provide reliable force-field parameters extensively benchmarked against density functional theory dynamics in both the gas and condensed phases. The procedure allows the classical description of metal complexes at a low computational cost with an accuracy as good as the quality of the Hessian matrix obtained by quantum chemistry methods. easyPARM v2.00 reads vibrational frequencies and charges in Gaussian (version 09 or 16) or ORCA (version 5 or 6) format and provides refined force-field parameters in Amber format. These can be directly used in Amber and NAMD molecular dynamics engines or converted to other formats. The tool is available free of charge in the GitHub platform (https://github.com/Abdelazim-Abdelgawwad/easyPARM.git).

Absorption Spectrum of Hydroperoxymethyl Thioformate: A Computational Chemistry Study

Hydroperoxymethyl thioformate (or HPMTF) is a compound relevant to the chemistry of sulfur in the marine atmosphere. The chemical cycling of this molecule in the atmosphere is still uncertain due in part to the lack of accurate knowledge of its photolytic behavior. Only approximations based on the properties of its chromophores are used in previous studies. In this work, we calculated the absorption spectra of the molecule in gas and aqueous phases using the Nuclear Ensemble Approach (NEA) and the CASPT2 method. Furthermore, we used such information to obtain relative photolysis rates. We found that the chromophore approximation overestimates the photolysis rates in the gas phase by twice the value obtained with the NEA-CASPT2 protocol. Furthermore, for the aqueous phase, we predict a lower role of photolysis as compared to the gas phase.

Nuclear Quantum Effects Have a Significant Impact on UV/Vis Absorption Spectra of Chromophores in Water

Despite the broadly acknowledged importance of solvation effects on measured UV/Vis spectra in the context of solvatochromism or chemical reactions in solution, it is still an open challenge to calculate UV/Vis spectra with predictive accuracy. This is particularly true when it comes to the impact of nuclear quantum effects on these experimental observables. In the present work, we calculate the UV/Vis absorption spectrum of indole in aqueous solution with a combination of a correlated wavefunction method for computing electronic excitation energies and enhanced path integral simulations for rigorous sampling of nuclear configurations including the quantum effects in solution. After validating our approach based on gas-phase benchmarking, we demonstrate that the lineshape of the spectrum measured in aqueous solution is quantitatively recovered, without the application of any shifting, scaling, or broadening, only after including nuclear quantum effects in addition to thermal fluctuations and solvation at ambient conditions. Our findings demonstrate that nuclear quantum effects are "visible" in UV/Vis spectra of chromophores measured in solution even at room temperature and, therefore, that they must be considered computationally to achieve predictive accuracy.

2-in-1 Phase Space Sampling for Calculating the Absorption Spectrum of the Hydrated Electron

The investigation of vibrational effects on absorption spectrum calculations often employs Wigner sampling or thermal sampling. While Wigner sampling incorporates zero-point energy, it may not be suitable for flexible systems. Thermal sampling is applicable to anharmonic systems yet treats nuclei classically. The application of generalized smoothed trajectory analysis (GSTA) as a postprocessing method allows for the incorporation of nuclear quantum effects (NQEs), combining the advantages of both sampling methods. We demonstrate this approach in computing the absorption spectrum of a hydrated electron. Theoretical exploration of the hydrated electron and its embryonic forms, such as water cluster anions, poses a significant challenge due to the diffusivity of the excess electron and the continuous motion of water molecules. In many previous studies, the wave nature of atomic nuclei is often neglected, despite the substantial impact of NQEs on thermodynamic and spectroscopic properties, particularly for hydrogen atoms. In our studies, we examine these NQEs for the excess electrons in various water systems. We obtained structures from mixed classical-quantum simulations for water cluster anions and the hydrated electron by incorporating the quantum effects of atomic nuclei with the filtration of the classical trajectories. Absorption spectra were determined at different theoretical levels. Our results indicate significant NQEs, red shift, and broadening of the spectra for hydrated electron systems. This study demonstrates the applicability of GSTA to complex systems, providing insights into NQEs on energetic and structural properties.

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.

Non-Phenomenological Description of the Time-Resolved Emission in Solution with Quantum-Classical Vibronic Approaches-Application to Coumarin C153 in Methanol

We report a joint experimental and theoretical work on the steady-state spectroscopy and time-resolved emission of the coumarin C153 dye in methanol. The lowest energy excited state of this molecule is characterized by an intramolecular charge transfer thus leading to remarkable shifts of the time-resolved emission spectra, dictated by the methanol reorganization dynamics. We selected this system as a prototypical test case for the first application of a novel computational protocol aimed at the prediction of transient emission spectral shapes, including both vibronic and solvent effects, without applying any phenomenological broadening. It combines a recently developed quantum-classical approach, the adiabatic molecular dynamics generalized vertical Hessian method (Ad-MD|gVH), with nonequilibrium molecular dynamics simulations. For the steady-state spectra we show that the Ad-MD|gVH approach is able to reproduce quite accurately the spectral shapes and the Stokes shift, while a ∼0.15 eV error is found on the prediction of the solvent shift going from gas phase to methanol. The spectral shape of the time-resolved emission signals is, overall, well reproduced, although the simulated spectra are slightly too broad and asymmetric at low energies with respect to experiments. As far as the spectral shift is concerned, the calculated spectra from 4 ps to 100 ps are in excellent agreement with experiments, correctly predicting the end of the solvent reorganization after about 20 ps. On the other hand, before 4 ps solvent dynamics is predicted to be too fast in the simulations and, in the sub-ps timescale, the uncertainty due to the experimental time resolution (300 fs) makes the comparison less straightforward. Finally, analysis of the reorganization of the first solvation shell surrounding the excited solute, based on atomic radial distribution functions and orientational correlations, indicates a fast solvent response (≈100 fs) characterized by the strengthening of the carbonyl-methanol hydrogen bond interactions, followed by the solvent reorientation, occurring on the ps timescale, to maximize local dipolar interactions.

Unraveling the contributions to the spectral shape of flexible dyes in solution: insights on the absorption spectrum of an oxyluciferin analogue

We present a computational investigation of the absorption spectrum in water of 5,5-spirocyclopropyl-oxyluciferin (5,5-CprOxyLH), an analogue of the emitter compound responsible for the bioluminescence in fireflies. Several factors participate in determining the 5,5-CprOxyLH's spectral shape: (i) the contribution of the four close-energy excited states, which show significant non-adiabatic couplings, (ii) the flexible molecular structure and (iii) the specific interactions established with the surrounding environment, which strongly couple the protic solvent dynamics with the dye's spectral response. To tackle the challenge to capture and dissect the role of all these effects we preliminarily investigate the role of non-adiabatic couplings with quantum dynamics simulations and a linear vibronic coupling model in the gas phase. Then, we account for both the molecular flexibility and solvent interactions by resorting to a mixed quantum classical protocol, named Adiabatic Molecular Dynamics generalized Vertical Gradient (Ad-MD|gVG), which is built on a method recently proposed by some of us. It is rooted in the partition between stiff degrees of freedom of the dye, accounted for at the vibronic level within the harmonic approximation, and flexible degrees of freedom of the solute (and of the solvent), described classically through a sampling based on Molecular Dynamics (MD). Ad-MD|gVG avoids spurious effects arising in the excited state Hessians due to non-adiabatic couplings, and can therefore be applied to account for the contributions of the first four excited states to the 5,5-CprOxyLH absorption spectrum. The final simulated spectrum is in very good agreement with the experiment, especially when the MD is driven by a refined quantum-mechanically derived force-field. More importantly, the origin of each separate contribution to the spectral shape is appropriately accounted for, paving the way to future applications of the method to more complex systems or alternative spectroscopies, as emission or circular dichroism.

Spatial and Temporal Resolution of the Oxygen-Independent Photoinduced DNA Interstrand Cross-Linking by a Nitroimidazole Derivative

DNA damage is ubiquitous in nature and is at the basis of emergent treatments such as photodynamic therapy, which is based on the activation of highly oxidative reactive oxygen species by photosensitizing O₂. However, hypoxia observed in solid tumors imposes the necessity to devise oxygen-independent modes of action able to induce DNA damage under a low oxygen concentration. The complexity of these DNA damage mechanisms in realistic environments grows exponentially when taking into account light absorption and subsequent excited-state population, photochemical and (photo)-redox reactions, the multiple species involved in different electronic states, noncovalent interactions, multiple reaction steps, and the large number of DNA reactive sites. This work tackles all the intricate reactivity of a photosensitizer based on a nitroimidazole derivative reacting toward DNA in solution under UV light exposition. This is performed through a combination of ground- and excited-state quantum chemistry, classical molecular dynamics, and hybrid QM/MM simulations to rationalize in detail the formation of DNA interstrand cross-links (ICLs) exerted by the noncanonical noncovalent photosensitizer. Unprecedented spatial and temporal resolution of these phenomena is achieved, revealing that the ICL is sequence-specific and that the fastest reactions take place at AT, GC, and GT steps involving either the opposite nucleobases or adjacent Watson-Crick base pairs. The N7 and O6 positions of guanine, the N7 and N3 sites of adenine, the N4 position of cytosine, and the O2 atom of thymine are deemed as the most nucleophile sites and are positively identified to participate in the ICL productions. This work provides a multiscale computational protocol to study DNA reactivity with noncovalent photosensitizers, and contributes to the understanding of therapies based on photoinduced DNA damage at molecular and electronic levels. In addition, we believe the depth understanding of these processes should assist the design of new photosensitizers considering their molecular size, electronic properties, and the observed regioselectivity toward nucleic acids.

Reconstruction of Nuclear Ensemble Approach Electronic Spectra Using Probabilistic Machine Learning

The theoretical prediction of molecular electronic spectra by means of quantum mechanical (QM) computations is fundamental to gain a deep insight into many photophysical and photochemical processes. A computational strategy that is attracting significant attention is the so-called Nuclear Ensemble Approach (NEA), that relies on generating a representative ensemble of nuclear geometries around the equilibrium structure and computing the vertical excitation energies (ΔE) and oscillator strengths (f) and phenomenologically broadening each transition with a line-shaped function with empirical full-width δ. Frequently, the choice of δ is carried out by visually finding the trade-off between artificial vibronic features (small δ) and over-smoothing of electronic signatures (large δ). Nevertheless, this approach is not satisfactory, as it relies on a subjective perception and may lead to spectral inaccuracies overall when the number of sampled configurations is limited due to an excessive computational burden (high-level QM methods, complex systems, solvent effects, etc.). In this work, we have developed and tested a new approach to reconstruct NEA spectra, dubbed GMM-NEA, based on the use of Gaussian Mixture Models (GMMs), a probabilistic machine learning algorithm, that circumvents the phenomenological broadening assumption and, in turn, the use of δ altogether. We show that GMM-NEA systematically outperforms other data-driven models to automatically select δ overall for small datasets. In addition, we report the use of an algorithm to detect anomalous QM computations (outliers) that can affect the overall shape and uncertainty of the NEA spectra. Finally, we apply GMM-NEA to predict the photolysis rate for HgBrOOH, a compound involved in Earth's atmospheric chemistry.

OS100: A Benchmark Set of 100 Digitized UV-Visible Spectra and Derived Experimental Oscillator Strengths

Excited-state quantum chemical calculations usually report excitation energies and oscillator strengths, f, for each electronic transition. On the other hand, UV-visible spectrophotometric experiments measure energy-dependent molar extinction/attenuation coefficients, ε(v), that give absorption band line shapes when plotted. ε(v) and f are related, but this relation is complicated by broadening and solvation effects. We fitted and integrated 100 experimental UV-visible spectra to obtain 164 f_exp values for absorption bands appearing in these spectra. The 100 UV-visible spectra belong to solvated organic molecules ranging in size from 6-34 atoms. We estimated uncertainties in the fitting to indicate confidence level in the reported f_exp values. The corresponding computed oscillator strengths (f_comp) were obtained with time-dependent density functional theory and a polarizable continuum solvent model. By expressing experimental and computed absorption strengths using a common quantity, we directly compared f_comp and f_exp. Although f_comp and f_exp are well correlated (linear regression R² = 0.921), f_comp in most cases overestimated f_exp (regression slope = 1.34). The agreement between absolute f_comp and f_exp values was substantially improved by accounting for a solvent refractive index factor, as suggested in some derivations in the literature. The 100 digitized UV-visible spectra are included as plain text files in the Supporting Information to aid in benchmarking computational or machine learning methods that aim to simulate realistic UV-visible absorption spectra.

Calculating Photoabsorption Cross-Sections for Atmospheric Volatile Organic Compounds

Characterizing the photochemical reactivity of transient volatile organic compounds (VOCs) in our atmosphere begins with a proper understanding of their abilities to absorb sunlight. Unfortunately, the photoabsorption cross-sections for a large number of transient VOCs remain unavailable experimentally due to their short lifetime or high reactivity. While structure-activity relationships (SARs) have been successfully employed to estimate the unknown photoabsorption cross-sections of VOCs, computational photochemistry offers another promising strategy to predict not only the vertical electronic transitions of a given molecule but also the width and shape of the bands forming its absorption spectrum. In this work, we focus on the use of the nuclear ensemble approach (NEA) to determine the photoabsorption cross-section of four exemplary VOCs, namely, acrolein, methylhydroperoxide, 2-hydroperoxy-propanal, and (microsolvated) pyruvic acid. More specifically, we analyze the influence that different strategies for sampling the ground-state nuclear density-Wigner sampling and ab initio molecular dynamics with a quantum thermostat-can have on the simulated absorption spectra. We highlight the potential shortcomings of using uncoupled harmonic modes within Wigner sampling of nuclear density to describe flexible or microsolvated VOCs and some limitations of SARs for multichromophoric VOCs. Our results suggest that the NEA could constitute a powerful tool for the atmospheric community to predict the photoabsorption cross-section for transient VOCs.

Multiscale Modeling of Electronic Spectra Including Nuclear Quantum Effects

Theoretical prediction of electronic absorption spectra without input from experiments is no easy feat, as it requires addressing all of the factors that affect line shapes. In practice, however, the methodologies are limited to treat these ingredients only to a certain extent. Here, we present a multiscale protocol that addresses the temperature, solvent, and nuclear quantum effects as well as anharmonicity and the reconstruction of the final spectra from individual transitions. First, quantum mechanics/molecular mechanics (QM/MM) molecular dynamics is conducted to obtain trajectories of solute-solvent configurations, from which the corresponding quantum-corrected ensembles are generated through the generalized smoothed trajectory analysis (GSTA). The optical spectra of the ensembles are then produced by calculating vertical transitions using time-dependent density-functional theory (TDDFT) with implicit solvation. To obtain the final spectral shapes, the stick spectra from TDDFT are convoluted with Gaussian kernels where the half-widths are determined by a statistically motivated strategy. We have tested our method by calculating the UV-vis spectra of a recently discovered acridine photocatalyst in two redox states. Vibronic progressions and broadenings due to the finite lifetime of the excited states are not included in the methodology yet. Nuclear quantization affects the relative peak intensities and widths, which is necessary to reproduce the experimental spectrum. We have also found that using only the optimized geometry of each molecule works surprisingly well if a proper empirical broadening factor is applied. This is explained by the rigidity of the conjugated chromophore moieties of the selected molecules, which are mainly responsible for the excitations in the spectra. In contrast, we have also shown that other parts of the molecules are flexible enough to feature anharmonicities that impair the use of other techniques such as Wigner sampling.