Nuclear dynamics in resonant inelastic X-ray scattering and X-ray absorption of methanol

New Implementation of an Equation-of-Motion Coupled-Cluster Damped-Response Framework with Illustrative Applications to Resonant Inelastic X-ray Scattering

We present an implementation of a damped response framework for calculating resonant inelastic X-ray scattering (RIXS) at the equation-of-motion coupled-cluster singles and doubles (CCSD) and second-order approximate coupled-cluster singles and doubles (CC2) levels of theory in the open-source program e^T. This framework lays the foundation for future extension to higher excitation methods (notably, the coupled-cluster singles and doubles with perturbative triples, CC3) and to multilevel approaches. Our implementation adopts a fully relaxed ground state and different variants of the core-valence separation projection technique to address convergence issues. Illustrative results are compared with those obtained within the frozen-core core-valence separated approach, available in Q-Chem, as well as with experiment. The performance of the CC2 method is evaluated in comparison with that of CCSD. It is found that, while the CC2 method is noticeably inferior to CCSD for X-ray absorption spectra, the quality of the CC2 RIXS spectra is often comparable to that of the CCSD level of theory, when the same valence excited states are probed. Finally, we present preliminary RIXS results for a solvated molecule in aqueous solution.

R-Group stabilization in methylated formamides observed by resonant inelastic X-ray scattering

The inherent stability of methylated formamides is traced to a stabilization of the deep-lying σ-framework by resonant inelastic X-ray scattering at the nitrogen K-edge. Charge transfer from the amide nitrogen to the methyl groups underlie this stabilization mechanism that leaves the aldehyde group essentially unaltered and explains the stability of secondary and tertiary amides.

The inherent stability of methylated formamides is traced to a stabilization of the deep-lying σ-framework by resonant inelastic X-ray scattering at the nitrogen K-edge.

Fundamental electronic changes upon intersystem crossing in large aromatic photosensitizers: free base 5,10,15,20-tetrakis(4-carboxylatophenyl)porphyrin

Free base 5,10,15,20-tetrakis(4-carboxylatophenyl)porphyrin stands for the class of powerful porphyrin photosensitizers for singlet oxygen generation and light-harvesting. The atomic level selectivity of dynamic UV pump – N K-edge probe X-ray absorption spectroscopy in combination with time-dependent density functional theory (TD-DFT) gives direct access to the crucial excited molecular states within the unusual relaxation pathway. The efficient intersystem crossing, that is El-Sayed forbidden and not facilitated by a heavy atom is confirmed to be the result of the long singlet excited state lifetime (Q_x 4.9 ns) and thermal effects. Overall, the interplay of stabilization by conservation of angular momenta and vibronic relaxation drive the de-excitation in these chromophores.

The crucial transient states of free-base porphyrins are characterized by time-resolved X-ray absorption spectroscopy unraveling their unusual relaxation pathway.

Targeting Individual Tautomers in Equilibrium by Resonant Inelastic X-ray Scattering

Tautomerism is one of the most important forms of isomerism, owing to the facile interconversion between species and the large differences in chemical properties introduced by the proton transfer connecting the tautomers. Spectroscopic techniques are often used for the characterization of tautomers. In this context, separating the overlapping spectral response of coexisting tautomers is a long-standing challenge in chemistry. Here, we demonstrate that by using resonant inelastic X-ray scattering tuned to the core excited states at the site of proton exchange between tautomers one is able to experimentally disentangle the manifold of valence excited states of each tautomer in a mixture. The technique is applied to the prototypical keto-enol equilibrium of 3-hydroxypyridine in aqueous solution. We detect transitions from the occupied orbitals into the LUMO for each tautomer in solution, which report on intrinsic and hydrogen-bond-induced orbital polarization within the π and σ manifolds at the proton-transfer site.

Interpretation of the X-Ray Emission Spectra of Liquid Water through Temperature and Isotope Dependence

The interpretation of x-ray emission spectroscopy (XES) spectra in terms of their sensitivity to the hydrogen bonding and the consequent microheterogeneity in liquid water has been debated over a decade. To shed a light on this problem, we report the theoretical reproduction of the debated 1b_{1} peaks observed in the XES spectra of liquid water using semiclassical Kramers-Heisenberg formula. The essence of the temperature and isotope dependence of the 1b_{1} double peaks is explained by molecular dynamics simulations including full vibrational (O─H stretching, bending, and) modes, rotational combined with the density functional theory and core-hole induced dynamics. Some inconsistencies exist with the experimental XES profile, which illustrates the need to employ a more precise theoretical calculations for both geometry sampling and electronic structure using a more sophisticated procedure.

Extended quasiparticle approach to non-resonant and resonant X-ray emission spectroscopy

The initial state of X-ray emission spectroscopy (XES) and resonant inelastic X-ray emission spectroscopy (RIXS) is a highly excited eigenstate with a deep core hole after a X-ray photoelectron spectroscopy...

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.

The porphyrin center as a regulator for metal-ligand covalency and π hybridization in the entire molecule

The central moiety of porphyrins is shown to control the charge state of the inner complex and links it by covalent interaction to the peripheral substituents. This link, which enables the versatile functions of porphyrins, is not picked up in the established, reduced four orbital picture [Gouterman, J. Mol. Spectrosc., 1961, 6, 138]. X-ray absorption spectroscopy at the N K-edge with density functional theory approaches gives access to the full electronic structure, in particular the π* manifold beyond the Gouterman orbitals. Systematic variation of the central moiety highlights two linked, governing trends: The ionicity of the porphyrin center increases from the aminic N-H to N-Cu to N-Zn to N-Mg to the iminic N:. At the same time covalency with peripheral substituents increases and compensates the buildup of high charge density at the coordinated nitrogen sites.

Predictions of Pre-edge Features in Time-Resolved Near-Edge X-ray Absorption Fine Structure Spectroscopy from Hole-Hole Tamm-Dancoff-Approximated Density Functional Theory

Time-resolved near-edge X-ray absorption fine structure (TR-NEXAFS) spectroscopy is a powerful technique for studying photochemical reaction dynamics with femtosecond time resolution. In order to avoid ambiguity in TR-NEXAFS spectra from nonadiabatic dynamics simulations, core- and valence-excited states must be evaluated on equal footing and those valence states must also define the potential energy surfaces used in the nonadiabatic dynamics simulation. In this work, we demonstrate that hole-hole Tamm-Dancoff-approximated density functional theory (hh-TDA) is capable of directly simulating TR-NEXAFS spectroscopies. We apply hh-TDA to the excited-state dynamics of acrolein. We identify two pre-edge features in the oxygen K-edge TR-NEXAFS spectrum associated with the S₂ (ππ*) and S₁ (nπ*) excited states. We show that these features can be used to follow the internal conversion dynamics between the lowest three electronic states of acrolein. Due to the low, O(N²) apparent computational complexity of hh-TDA and our GPU-accelerated implementation, this method is promising for the simulation of pre-edge features in TR-NEXAFS spectra of large molecules and molecules in the condensed phase.

Breaking the Symmetry of Pyrimidine: Solvent Effects and Core-Excited State Dynamics

Symmetry and its breaking crucially define the chemical properties of molecules and their functionality. Resonant inelastic X-ray scattering is a local electronic structure probe reporting on molecular symmetry and its dynamical breaking within the femtosecond scattering duration. Here, we study pyrimidine, a system from the C_2v point group, in an aqueous solution environment, using scattering though its 2a₂ resonance. Despite the absence of clean parity selection rules for decay transitions from in-plane orbitals, scattering channels including decay from the 7b₂ and 11a₁ orbitals with nitrogen lone pair character are a direct probe for molecular symmetry. Computed spectra of explicitly solvated molecules sampled from a molecular dynamics simulation are combined with the results of a quantum dynamical description of the X-ray scattering process. We observe dominant signatures of core-excited Jahn-Teller induced symmetry breaking for resonant excitation. Solvent contributions are separable by shortening of the effective scattering duration through excitation energy detuning.

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.

Valence Photoelectron Spectra of Liquid Methanol and Ethanol Measured Using He II Radiation

High-resolution photoelectron (PE) spectra of liquid methanol and ethanol were measured using a liquid microjet and He IIα radiation (40.813 eV). The vertical ionization energy and the ionization threshold were determined as 9.70 ± 0.07 and 8.69 ± 0.07 eV for methanol and 9.52 ± 0.07 and 8.52 ± 0.07 eV for ethanol, respectively. Individual photoemission bands observed for the liquids are well correlated with those in PE spectra of the gaseous samples also measured in the present study, except that the liquid band positions were shifted on average by -1.23 eV for methanol and -1.10 eV for ethanol as compared to the gas. The 5a' and 7a' bands of liquid methanol exhibit specifically larger broadening than other bands, for which we attempted spectral fitting with two components, similarly with the case of the 3a₁ band of liquid water. PE spectra of both liquid and gaseous ethanol are congested partly due to the presence of the trans and gauche isomers; however, the overall band positions are generally in good agreement with predictions based on quantum chemical calculations. Comparison of the measured PE spectra with experimental and simulated X-ray emission spectra indicate that spectral differences in the lowest ionization band of both methanol and ethanol originate from involvement of nuclear dynamics in the X-ray emission process.

How Hydrogen Bonding Amplifies Isomeric Differences in Pyridones toward Strong Changes in Acidity and Tautomerism

Steric hindrance of hydration and hydrogen bond enhancement by localized charges have been identified as key factors for the massive chemical differences between the hydroxypyridine/pyridone isomers in aqueous solution. While all isomers occur mainly in the hydroxypyridine form in the gas phase, they differ by more than 3 orders of magnitude both in their acidity and tautomeric equilibrium constants upon hydration. By monitoring the electronic and solvation structures as a function of the protonation state and the O^– substitution position on the pyridine ring, the amplification of the isomeric differences in aqueous solution has been investigated. Near-edge X-ray absorption fine structure (NEXAFS) measurements at the N K-edge served as the probe of the chemical state. The combination of molecular dynamics simulations, complete active space self-consistent field (CASSCF), and time-dependent density functional theory (TD-DFT) spectral calculations contributes to unraveling the principles of tautomerism and acidity in multiple biochemical systems based on tautomerism.

Vibrational resonant inelastic X-ray scattering in liquid acetic acid: a ruler for molecular chain lengths

Quenching of vibrational excitations in resonant inelastic X-ray scattering (RIXS) spectra of liquid acetic acid is observed. At the oxygen core resonance associated with localized excitations at the O-H bond, the spectra lack the typical progression of vibrational excitations observed in RIXS spectra of comparable systems. We interpret this phenomenon as due to strong rehybridization of the unoccupied molecular orbitals as a result of hydrogen bonding, which however cannot be observed in x-ray absorption but only by means of RIXS. This allows us to address the molecular structure of the liquid, and to determine a lower limit for the average molecular chain length.

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.