Absorption of sunlight generates renewable electricity and powers the growth of plants, but also causes severe damage both to synthetic materials and biological tissue. The wildly varying outcomes of these light-induced processes are ultimately determined by much slighter differences in their underlying reaction pathways, induced by the transient properties of short-lived and miniscule molecules; a powerful approach to their detection and characterization is offered by ultrafast x-ray spectroscopy, with identification of spectral fingerprints and further guidance from quantum chemical models.

This thesis contains the computational half of three experimentally joint projects that push the limits for detection of electronic, spin and structural dynamics of small molecular systems in solution. A wide selection of theoretical frameworks are combined to model various aspects of the measurements: from multi-configurational descriptions of non-adiabatic couplings in the photo-dynamics and multi-electron transitions in the x-ray spectroscopy, to affordable simulations of extensive aqueous solutions by density functional theory and classical mechanics.

Applied to experimental data, the presented quantum chemical results allowed in particular to: simultaneously identify molecular forms and electronic states of aqueous 2-thiopyridone, to determine a detailed pathway for its excited-state proton-transfer; characterize the charge-transfer state of aqueous ferricyanide, to extend well-known concepts from steady-state spectroscopy into the ultrafast domain; establish the newly implemented framework of multi-configurational Dyson orbitals, as a powerful tool for simulation of photoelectron spectroscopy.

A number of computational predictions are additionally presented for hitherto-unexplored experimental regions, which may help to guide and optimize future measurements.

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.

Ultrafast non-adiabatic dynamics of the small heteroaromatic compound 2-thiopyridone has been studied with surface hopping simulations based on multi-configurational quantum chemistry. Initial excitation of the bright S₂ (π,π^∗) state is found to promptly relax to S₁ (n,π^∗) through in-plane motion. The subsequent dynamics are oppositely driven by out-of-plane motion, which results both in complex population transfers among all of the available states and intersystem crossing predominantly through the ‘El-Sayed forbidden ’S₁ (n,π^∗) to T₂ (n,π^∗) channel, through significant mixing of electronic excitation characters. Despite this complexity, the femto- to picosecond triplet population, expected from several spectroscopic measurements, is well described as a simple exponential decay of the singlet state manifold. No proton transfer is found in the reported trajectories, but two mechanisms for its possible mediation in previously reported experiments are proposed based on the observed structural dynamics: (i) ultrafast intra-molecular transfer driven by the initially coherent in-plane motion and (ii) inter-molecular solvent-mediated transfer driven by the out-of-plane modes that dominate the later motion.

The emergence of X-ray free electron lasers (X-FELs) has made it possible to probe structural dynamics on the femtosecond timescale. This extension of experimental capabilities also calls for a simultaneous development in theory to help interpret the underlying structure and dynamics encoded within the experimental observable. In the ultrafast regime this often requires a time-dependent theoretical treatment that describes nuclear dynamics beyond the Born–Oppenheimer approximation. In this work, we perform quantum dynamics simulations based upon time-evolving Gaussian basis functions (GBFs) and simulate the ultrafast X-ray Absorption Near-Edge Structure (XANES) spectra of photoexcited pyrazine including two strongly coupled electronically excited states and four normal mode degrees of freedom. Two methods to simulate the excited state XANES spectra are applied, the first is based upon the multi-configurational second order perturbation theory restricted active space (RASPT2) method and the second uses a combination of the maximum overlap method (MOM) and time-dependent density functional theory (TDDFT). We demonstrate that despite the simplicity of the MOM/TDDFT method, it captures several qualitative features of the RASPT2 simulations at much reduced computational effort. However, features such as the conical intersection are a particular exception as they require a multi-configurational treatment. For the nuclear dynamics, we demonstrate that even a small number of GBFs can provide reasonable description of the spectroscopic observable. This work provides perspectives for computationally efficient approaches important for addressing larger systems.

The nitrogen and sulfur K-edge X-ray absorption spectra of aqueous 2-thiopyridone, a model system for excited-state proton transfer in several recent time-resolved measurements, have been simulated from ab initio molecular dynamics. Spectral signatures of the local intra- and inter-molecular structure are identified and rationalized, which facilitates experimental interpretation and optimization. In particular, comparison of aqueous and gas phase spectrum simulations assesses the previously unquantified solvatization effects, where hydrogen bonding is found to yield solvatochromatic shifts up to nearly 1 eV of the main peak positions. Thereby, while each K-edge can still decisively determine the local protonation of its core-excited site, only their combined, complementary fingerprints allow separating all of the three relevant molecular forms, giving a complete picture of the proton transfer.

2-Thiopyridone (2-TP), a common model system for excited-state proton transfer, has been simulated in aqueous solution with ab initio molecular dynamics. The interplay of electronic structure, protonation, and solvatization is investigated by comparison of three differently protonated molecular forms and between the lowest singlet and triplet electronic states. An interdependence clearly manifests in the mixed-character T-1 state for the 2-TP form, systematic structural distortions of the 2-mercaptopyridine (2-MP) form, and photobase protolysis of the 2-TP- form, in the aqueous phase. In comparison, simplified continuum models for the solvatization are found to be significantly inaccurate for several of the species. To facilitate future computational studies, we therefore present a minimal representative solvatization complex for each stable form and electronic state. Our findings demonstrate the importance of explicit solvatization of the compound and sets the studies. stage for including it also in future studies.

In this Article we describe the OpenMolcas environment and invite the computational chemistry community to collaborate. The open-source project already includes a large number of new developments realized during the transition from the commercial MOLCAS product to the open-source platform. The paper initially describes the technical details of the new software development platform. This is followed by brief presentations of many new methods, implementations, and features of the OpenMolcas program suite. These developments include novel wave function methods such as stochastic complete active space self-consistent field, density matrix renormalization group (DMRG) methods, and hybrid multiconfigurational wave function and density functional theory models. Some of these implementations include an array of additional options and functionalities. The paper proceeds and describes developments related to explorations of potential energy surfaces. Here we present methods for the optimization of conical intersections, the simulation of adiabatic and nonadiabatic molecular dynamics, and interfaces to tools for semiclassical and quantum mechanical nuclear dynamics. Furthermore, the Article describes features unique to simulations of spectroscopic and magnetic phenomena such as the exact semiclassical description of the interaction between light and matter, various X-ray processes, magnetic circular dichroism, and properties. Finally, the paper describes a number of built-in and add-on features to support the OpenMolcas platform with postcalculation analysis and visualization, a multiscale simulation option using frozen-density embedding theory, and new electronic and muonic basis sets.

Excited-state proton transfer (ESPT) is a fundamental process in biomolecular photochemistry, but its underlying mediators often evade direct observation. We identify a distinct pathway for ESPT in aqueous 2-thiopyridone, by employing transient N1s X-ray absorption spectroscopy and multi-configurational spectrum simulations. Photoexcitations to the singlet S-2 and S-4 states both relax promptly through intersystem crossing to the triplet T-1 state. The T-1 state, through its rapid population and near nanosecond lifetime, mediates nitrogen site deprotonation by ESPT in a secondary intersystem crossing to the S-0 potential energy surface. This conclusively establishes a dominant ESPT pathway for the system in aqueous solution, which is also compatible with previous measurements in acetonitrile. Thereby, the hitherto open questions of the pathway for ESPT in the compound, including its possible dependence on excitation wavelength and choice of solvent, are resolved.

Soft X-ray spectroscopies are ideal probes of the local valence electronic structure of photocatalytically active metal sites. Here, we apply the selectivity of time resolved resonant inelastic X-ray scattering at the iron L-edge to the transient charge distribution of an optically excited charge-transfer state in aqueous ferricyanide. Through comparison to steady-state spectra and quantum chemical calculations, the coupled effects of valence-shell closing and ligand-hole creation are experimentally and theoretically disentangled and described in terms of orbital occupancy, metal-ligand covalency, and ligand field splitting, thereby extending established steady-state concepts to the excited-state domain. pi-Back-donation is found to be mainly determined by the metal site occupation, whereas the ligand hole instead influences sigma-donation. Our results demonstrate how ultrafast resonant inelastic X-ray scattering can help characterize local charge distributions around catalytic metal centers in short-lived charge-transfer excited states, as a step toward future rationalization and tailoring of photocatalytic capabilities of transition-metal complexes.

We describe how inversion symmetry separation of electronic state manifolds in resonant inelastic soft X-ray scattering (RIXS) can be applied to probe excited-state dynamics with compelling selectivity. In a case study of Fe L-3-edge RIXS in the ferricyanide complex Fe(CN)(6)(3-), we demonstrate with multi-configurational restricted active space spectrum simulations how the information content of RIXS spectral fingerprints can be used to unambiguously separate species of different electronic configurations, spin multiplicities, and structures, with possible involvement in the decay dynamics of photo-excited ligand-to-metal charge-transfer. Specifically, we propose that this could be applied to confirm or reject the presence of a hitherto elusive transient Quartet species. Thus, RIXS offers a particular possibility to settle a recent controversy regarding the decay pathway, and we expect the technique to be similarly applicable in other model systems of photo-induced dynamics.