Electronic structure, bonding and stability of fumarate, maleate, and succinate dianions from X-ray spectroscopy

The electronic structure of the fumarate, maleate, and succinate dianions in the context of their stability is determined in a joint experimental and computational study with X-ray absorption spectroscopy and resonant inelastic X-ray scattering at the O K-edge. The study reveals differences in the electronic states and molecular orbitals of the three molecules. In particular, maleate has a non-degenerate oxygen core-orbital with an energy difference of approximately 0.15 eV, visible in a two peak structure in XAS. Polarization-dependent RIXS provides information on the orientation of the occupied valence molecular orbitals with respect to the carboxylate group plane and shows a gradually increasing energy gap between the HOMO and excited π* LUMO from fumarate to maleate to succinate. We also demonstrate the energy excitation dependence of the RIXS spectra of maleate, with the total inelastic RIXS profile shifting towards higher energy loss as the detuning is increased from negative to positive values. Our findings show that maleate is less stable than fumarate and succinate due to the presence of electronic density on its HOMO orbital on the CC bond between carboxylate groups, which can lead to weaker bonding of maleate with molecules or ions.

Soft X-ray Spectroscopy Simulations with Multiconfigurational Wave Function Theory: Spectrum Completeness, Sub-eV Accuracy, and Quantitative Reproduction of Line Shapes

Multireference methods are known for their ability to accurately treat states of very different nature in many molecular systems, facilitating high-quality simulations of a large variety of spectroscopic techniques. Here, we couple the multiconfigurational restricted active space self-consistent field RASSCF/RASPT2 method (of the CASSCF/CASPT2 methods family) to the displaced harmonic oscillator (DHO) model, to simulate soft X-ray spectroscopy. We applied such an RASSCF/RASPT2+DHO approach at the K-edges of various second-row elements for a set of small organic molecules that have been recently investigated at other levels of theory. X-ray absorption near-edge structure (XANES) and X-ray photoelectron spectroscopy (XPS) are simulated with a sub-eV accuracy and a correct description of the spectral line shapes. The method is extremely sensitive to the observed spectral shifts on a series of differently fluorinated ethylene systems, provides spectral fingerprints to distinguish between stable conformers of the glycine molecule, and accurately captures the vibrationally resolved carbon K-edge spectrum of formaldehyde. Differences with other theoretical methods are demonstrated, which show the advantages of employing a multireference/multiconfigurational approach. A protocol to systematically increase the number of core-excited states considered while maintaining a contained computational cost is presented. Insight is eventually provided for the effects caused by removing core-electrons from a given atom in terms of bond rearrangement and influence on the resulting spectral shapes within a unitary orbital-based framework for both XPS and XANES spectra.

Stochastic Generalized Active Space Self-Consistent Field: Theory and Application

An algorithm to perform stochastic generalized active space calculations, Stochastic-GAS, is presented, that uses the Slater determinant based FCIQMC algorithm as configuration interaction eigensolver. Stochastic-GAS allows the construction and stochastic optimization of preselected truncated configuration interaction wave functions, either to reduce the computational costs of large active space wave function optimizations, or to probe the role of specific electron correlation pathways. As for the conventional GAS procedure, the preselection of the truncated wave function is based on the selection of multiple active subspaces while imposing restrictions on the interspace excitations. Both local and cumulative minimum and maximum occupation number constraints are supported by Stochastic-GAS. The occupation number constraints are efficiently encoded in precomputed probability distributions, using the precomputed heat bath algorithm, which removes nearly all runtime overhead of GAS. This strategy effectively allows the FCIQMC dynamics to a priori exclude electronic configurations that are not allowed by GAS restrictions. Stochastic-GAS reduced density matrices are stochastically sampled, allowing orbital relaxations via Stochastic-GASSCF, and direct evaluation of properties that can be extracted from density matrices, such as the spin expectation value. Three test case applications have been chosen to demonstrate the flexibility of Stochastic-GAS: (a) the Stochastic-GASSCF [5·(6, 6)] optimization of a stack of five benzene molecules, that shows the applicability of Stochastic-GAS toward fragment-based chemical systems; (b) an uncontracted stochastic MRCISD calculation that correlates 96 electrons and 159 molecular orbitals, and uses a large (32, 34) active space reference wave function for an Fe(II)-porphyrin model system, showing how GAS can be applied to systematically recover dynamic electron correlation, and how in the specific case of the Fe(II)-porphyrin dynamic correlation further differentially stabilizes the ³E_g over the ⁵A_1g spin state; (c) the study of an Fe₄S₄ cluster's spin-ladder energetics via highly truncated stochastic-GAS [4·(5, 5)] wave functions, where we show how GAS can be applied to understand the competing spin-exchange and charge-transfer correlating mechanisms in stabilizing different spin-states.

X-ray emission spectroscopy: a genetic algorithm to disentangle core–hole-induced dynamics

A genetic algorithm (GA) is developed and applied to make proper connections of final-state potential-energy surfaces and X-ray emission (XES) cross sections between steps in the time-propagation of H-bonded systems after a core–hole is created. We show that this modification results in significantly improved resolution of spectral features in XES with the semiclassical Kramers–Heisenberg approach which takes into account important interference effects. We demonstrate the effects on a water pentamer model as well as on two 17-molecules water clusters representing, respectively, tetrahedral (D2A2) and asymmetric (D1A1) H-bonding environments. For D2A2, the applied procedure improves significantly the obtained intensities, whereas for D1A1 the effects are smaller due to milder dynamics during the core–hole life-time as only one hydrogen is involved. We reinvestigate XES for liquid ethanol and, by properly disentangling the relevant states in the dense manifold of states using the GA, now resolve the important 3a′′ state as a peak rather than a shoulder. Furthermore, by applying the SpecSwap-RMC procedure, we reweigh the distribution of structures in the sampling of the liquid to fit to experiment and estimate the ratio between the main anti and gauche conformers in the liquid at room temperature. This combination of techniques will be generally applicable to challenging problems in liquid-phase spectroscopy.

Coupling Methylammonium and Formamidinium Cations with Halide Anions: Hybrid Orbitals, Hydrogen Bonding, and the Role of Dynamics

The electronic structures of four precursors for organic-inorganic hybrid perovskites, namely, methylammonium chloride and iodide, as well as formamidinium bromide and iodide, are investigated by X-ray emission (XE) spectroscopy at the carbon and nitrogen K-edges. The XE spectra are analyzed based on density functional theory calculations. We simulate the XE spectra at the Kohn-Sham level for ground-state geometries and carry out detailed analyses of the molecular orbitals and the electronic density of states to give a thorough understanding of the spectra. Major parts of the spectra can be described by the model of the corresponding isolated organic cation, whereas high-emission energy peaks in the nitrogen K-edge XE spectra arise from electronic transitions involving hybrids of the molecular and atomic orbitals of the cations and halides, respectively. We find that the interaction of the methylammonium cation is stronger with the chlorine than with the iodine anion. Furthermore, our detailed theoretical analysis highlights the strong influence of ultrafast proton dynamics in the core-excited states, which is an intrinsic effect of the XE process. The inclusion of this effect is necessary for an accurate description of the experimental nitrogen K-edge X-ray emission spectra and gives information on the hydrogen-bonding strengths in the different precursor materials.

Analytical gradients for core-excited states in the algebraic diagrammatic construction (ADC) framework

Expressions for analytical molecular gradients of core-excited states have been derived and implemented for the hierarchy of algebraic diagrammatic construction (ADC) methods up to extended second-order within the core-valence separation (CVS) approximation. We illustrate the use of CVS-ADC gradients by determining relaxed core-excited state potential energy surfaces and optimized geometries for water, formic acid, and benzene. For water, our results show that in the dissociative lowest core-excited state, a linear configuration is preferred. For formic acid, we find that the O K-edge lowest core-excited state is non-planar, a fact that is not captured by the equivalent core approximation where the core-excited atom with its hole is replaced by the "Z + 1" neighboring atom in the periodic table. For benzene, the core-excited state gradients are presented along the Jahn-Teller distorted geometry of the 1s → π* excited state. Our development may pave a new path to studying the dynamics of molecules in their core-excited states.

Site-specific generation of excited state wavepackets with high-intensity attosecond x rays

High-intensity attosecond x rays can produce coherent superpositions of valence-excited states through two-photon Raman transitions. The broad-bandwidth, high-field nature of the pulses results in a multitude of accessible excited states. Multiconfigurational quantum chemistry with the time-dependent Schrödinger equation is used to examine population transfer dynamics in stimulated x-ray Raman scattering of the nitric oxide oxygen and nitrogen K-edges. Two pulse schemes initiate wavepackets of different characters and demonstrate how chemical differences between core-excitation pathways affect the dynamics. The population transfer to valence-excited states is found to be sensitive to the electronic structure and pulse conditions, highlighting complexities attributed to the Rabi frequency. The orthogonally polarized two-color-pulse setup has increased selectivity while facilitating longer, less intense pulses than the one-pulse setup. Population transfer in the 1s → Rydberg region is more effective but less selective at the nitrogen K-edge; the selectivity is reduced by double core-excited states. Result interpretation is aided by resonant inelastic x-ray scattering maps.

Hydrogen bond effects in multimode nuclear dynamics of acetic acid observed via resonant x-ray scattering

A theoretical and experimental study of the gas phase and liquid acetic acid based on resonant inelastic x-ray scattering (RIXS) spectroscopy is presented. We combine and compare different levels of theory for an isolated molecule for a comprehensive analysis, including electronic and vibrational degrees of freedom. The excitation energy scan over the oxygen K-edge absorption reveals nuclear dynamic effects in the core-excited and final electronic states. The theoretical simulations for the monomer and two different forms of the dimer are compared against high-resolution experimental data for pure liquid acetic acid. We show that the theoretical model based on a dimer describes the hydrogen bond formation in the liquid phase well and that this bond formation sufficiently alters the RIXS spectra, allowing us to trace these effects directly from the experiment. Multimode vibrational dynamics is accounted for in our simulations by using a hybrid time-dependent stationary approach for the quantum nuclear wave packet simulations, showing the important role it plays in RIXS.

TD-DFT simulations of K-edge resonant inelastic X-ray scattering within the restricted subspace approximation

Truncation of orbital subspaces in TD-DFT yields an accurate description of RIXS spectra for soft X-ray K-edges.

Observation of Double Excitations in the Resonant Inelastic X-ray Scattering of Nitric Oxide

The nitrogen K-edge resonant inelastic X-ray scattering (RIXS) map of nitric oxide (NO) has been measured and simulated to provide a detailed analysis of the observed features. High-resolution experimental RIXS maps were collected using an in situ gas flow cell and a high-transmission soft X-ray spectrometer. Accurate descriptions of the ground, excited, and core-excited states are based upon restricted active space self-consistent-field calculations using second order multiconfigurational perturbation theory. The nitrogen K-edge RIXS map of NO shows a range of features that can be assigned to intermediate states arising from 1s → π* and 1s → Rydberg excitations; additional bands are attributed to doubly excited intermediate states comprising 1s → π* and π → π* excitations. These results provide a detailed picture of RIXS for an open-shell molecule and an extensive description of the core-excited electronic structure of NO, an important molecule in many chemical and biological processes.

A simple molecular orbital picture of RIXS distilled from many-body damped response theory

Ab initio calculations of resonant inelastic x-ray scattering (RIXS) often rely on damped response theory, which prevents the divergence of response solutions in the resonant regime. Within the damped response theory formalism, RIXS moments are expressed as the sum over all electronic states of the system [sum-over-states (SOS) expressions]. By invoking resonance arguments, this expression can be reduced to a few terms, an approximation commonly exploited for the interpretation of computed cross sections. We present an alternative approach: a rigorous formalism for deriving a simple molecular orbital picture of the RIXS process from many-body calculations using the damped response theory. In practical implementations, the SOS expressions of RIXS moments are recast in terms of matrix elements between the zero-order wave functions and first-order frequency-dependent response wave functions of the initial and final states such that the RIXS moments can be evaluated using complex response one-particle transition density matrices (1PTDMs). Visualization of these 1PTDMs connects the RIXS process with the changes in electronic density. We demonstrate that the real and imaginary components of the response 1PTDMs can be interpreted as contributions of the undamped off-resonance and damped near-resonance SOS terms, respectively. By analyzing these 1PTDMs in terms of natural transition orbitals, we derive a rigorous, black-box mapping of the RIXS process into a molecular orbital picture. We illustrate the utility of the new tool by analyzing RIXS transitions in the OH radical, benzene, para-nitroaniline, and 4-amino-4'-nitrostilbene. These examples highlight the significance of both the near-resonance and off-resonance channels.

Exploring the capabilities of optical pump X-ray probe NEXAFS spectroscopy to track photo-induced dynamics mediated by conical intersections

In the present contribution we introduce an accurate theoretical approach for the simulation of NEXAFS spectra of organic molecules, employing azobenzene as a test case.

How to stay out of trouble in RIXS calculations within equation-of-motion coupled-cluster damped response theory? Safe hitchhiking in the excitation manifold by means of core–valence separation

We present a novel approach with robust convergence of the response equations for computing resonant inelastic X-ray scattering (RIXS) cross sections within the equation-of-motion coupled-cluster (EOM-CC) framework.

Probing hydrogen bond strength in liquid water by resonant inelastic X-ray scattering

Local probes of the electronic ground state are essential for understanding hydrogen bonding in aqueous environments. When tuned to the dissociative core-excited state at the O1s pre-edge of water, resonant inelastic X-ray scattering back to the electronic ground state exhibits a long vibrational progression due to ultrafast nuclear dynamics. We show how the coherent evolution of the OH bonds around the core-excited oxygen provides access to high vibrational levels in liquid water. The OH bonds stretch into the long-range part of the potential energy curve, which makes the X-ray probe more sensitive than infra-red spectroscopy to the local environment. We exploit this property to effectively probe hydrogen bond strength via the distribution of intramolecular OH potentials derived from measurements. In contrast, the dynamical splitting in the spectral feature of the lowest valence-excited state arises from the short-range part of the OH potential curve and is rather insensitive to hydrogen bonding.