The femtosecond excited-state dynamics following resonant photoexcitation enable the selective deformation of N-H and N-C chemical bonds in 2-thiopyridone in aqueous solution with optical or X-ray pulses. In combination with multiconfigurational quantum-chemical calculations, the orbital-specific electronic structure and its ultrafast dynamics accessed with resonant inelastic X-ray scattering at the N 1s level using synchrotron radiation and the soft X-ray free-electron laser LCLS provide direct evidence for this controlled photoinduced molecular deformation and its ultrashort time-scale.

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.

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.

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.

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.

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.

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.

A combination of multi-configurational restricted active space calculations with a Dyson orbital formalism has been applied for accurate simulations of 4d photo-electron spectra of the I-3(-) molecular ion. The analysis based on the occupation numbers of natural orbitals allowed to predict and rationalize the spectral fingerprints of solvent-induced nuclear asymmetry. In particular, it demonstrates how the nuclear asymmetry directly causes an increase of shake-up intensity. The relative intensity of shake-up and main features of the I 4d XPS spectrum could therefore serve as a simplified experimental observable of structural asymmetry, complementary to changes in the shape of the main spectral features.