Valence and Core-Level X-ray Photoelectron Spectroscopy of a Liquid Ammonia Microjet

Electron-induced hydroamination of ethane as compared to ethene: implications for the reaction mechanism

The properties of carbonaceous materials with respect to various applications are enhanced by incorporation of nitrogen-containing moieties like, for instance, amino groups. Therefore, processes that allow the introduction of such functional groups into hydrocarbon compounds are of utmost interest. Previous studies have demonstrated that hydroamination reactions which couple amines to unsaturated sites within hydrocarbon molecules do not only proceed in the presence of suitably tailored catalysts but can also be induced and controlled by electron irradiation. However, studies on electron-induced hydroaminations so far were guided by the hypothesis that unsaturated hydrocarbons are required for the reaction while the reaction would be much less efficient in the case of saturated hydrocarbons. The present work evaluates the validity of this hypothesis by post-irradiation thermal desorption experiments that monitor the electron energy-dependent yield of ethylamine after electron irradiation of mixed C2H4:NH3 and C2H6:NH3 ices with the same composition and thickness. The results reveal that, in contrast to the initial assumption, ethylamine is formed with similar efficiency in both mixed ices. From the dependence of the product yields on the electron energy, we conclude that the reaction in both cases is predominantly driven by electron ionization of NH3. Ethylamine is formed via alternative reaction mechanisms by which the resulting NH2˙ radicals add to C2H4 and C2H6, respectively. The similar efficiency of amine formation in unsaturated and saturated hydrocarbons demonstrates that electron irradiation in the presence of NH3 is a more versatile tool for introducing nitrogen into carbonaceous materials than previously anticipated.

Solvated dielectrons from optical excitation: An effective source of low-energy electrons

Low-energy electrons dissolved in liquid ammonia or aqueous media are powerful reducing agents that promote challenging reduction reactions, but can also cause radiation damage to biological tissue. Knowledge of the underlying mechanistic processes remains incomplete, in particular with respect to the details and energetics of the electron transfer steps. Here, we show how ultraviolet (UV) photoexcitation of metal-ammonia clusters could be used to generate tunable low-energy electrons in situ. Specifically, we identified UV light-induced generation of spin-paired solvated dielectrons and their subsequent relaxation by an unconventional electron-transfer-mediated decay as an efficient low-energy electron source. The process is robust and straightforward to induce, with the prospect of improving our understanding of radiation damage and fostering mechanistic studies of solvated electron reduction reactions.

Exploring Gas-Liquid Reactions with Microjets: Lessons We Are Learning

Liquid water is all around us: at the beach, in a cloud, from a faucet, inside a spray tower, covering our lungs. It is fascinating to imagine what happens to a reactive gas molecule as it approaches the surface of water in these examples. Some incoming molecules may first be deflected away after colliding with an evaporating water molecule. Those that do strike surface H₂O or other surface species may bounce directly off or become momentarily trapped through hydrogen bonding or other attractive forces. The adsorbed gas molecule can then desorb immediately or instead dissolve, passing through the interfacial region and into the bulk, perhaps diffusing back to the surface and evaporating before reacting. Alternatively, it may react with solute or water molecules in the interfacial or bulk regions, and a reaction intermediate or the final product may then desorb into the gas phase. Building a "blow by blow" picture of these pathways is challenging for vacuum-based techniques because of the high vapor pressure of water. In particular, collisions within the thick vapor cloud created by evaporating molecules just above the surface scramble the trajectories and internal states of the gaseous target molecules, hindering construction of gas-liquid reaction mechanisms at the atomic scale that we strive to map out.The introduction of the microjet in 1988 by Faubel, Schlemmer, and Toennies opened up entirely new possibilities. Their inspired solution seems so simple: narrow the end of a glass tube to a diameter smaller than the mean free path of the vapor molecules and then push the liquid through the tube at speeds of a car on a highway. The narrow liquid stream creates a sparse vapor cloud, with water molecules spaced far enough apart that they and the reactant gases interact, at most, weakly. Experimentalists, however, confront a host of challenges: nozzle clogging, unstable jetting, searching for vacuum-compatible solutions, measuring low signal levels, and teasing out artifacts because the slender jet is the smallest surface in the vacuum chamber. In this Account, we describe lessons that we are learning as we explore gases (DCl, (HCOOH)₂, N₂O₅) reacting with water molecules and solute ions in the near-interfacial region of these fast-flowing aqueous microjets.

Hydrogen Bonding in Liquid Ammonia

The nature of hydrogen bonding in condensed ammonia phases, liquid and crystalline ammonia has been a topic of much investigation. Here, we use quantum molecular dynamics simulations to investigate hydrogen bond structure and lifetimes in two ammonia phases: liquid ammonia and crystalline ammonia-I. Unlike liquid water, which has two covalently bonded hydrogen and two hydrogen bonds per oxygen atom, each nitrogen atom in liquid ammonia is found to have only one hydrogen bond at 2.24 Å. The computed lifetime of the hydrogen bond is t ≅ 0.1 ps. In contrast to crystalline water-ice, we find that hydrogen bonding is practically nonexistent in crystalline ammonia-I.

Observation of intermolecular Coulombic decay and shake-up satellites in liquid ammonia

We report the first nitrogen 1s Auger-Meitner electron spectrum from a liquid ammonia microjet at a temperature of ∼223 K (-50 °C) and compare it with the simultaneously measured spectrum for gas-phase ammonia. The spectra from both phases are interpreted with the assistance of high-level electronic structure and ab initio molecular dynamics calculations. In addition to the regular Auger-Meitner-electron features, we observe electron emission at kinetic energies of 374-388 eV, above the leading Auger-Meitner peak (3a₁ ²). Based on the electronic structure calculations, we assign this peak to a shake-up satellite in the gas phase, i.e., Auger-Meitner emission from an intermediate state with additional valence excitation present. The high-energy contribution is significantly enhanced in the liquid phase. We consider various mechanisms contributing to this feature. First, in analogy with other hydrogen-bonded liquids (noticeably water), the high-energy signal may be a signature for an ultrafast proton transfer taking place before the electronic decay (proton transfer mediated charge separation). The ab initio dynamical calculations show, however, that such a process is much slower than electronic decay and is, thus, very unlikely. Next, we consider a non-local version of the Auger-Meitner decay, the Intermolecular Coulombic Decay. The electronic structure calculations support an important contribution of this purely electronic mechanism. Finally, we discuss a non-local enhancement of the shake-up processes.

Imaging of Chemical Kinetics at the Water-Water Interface in a Free-Flowing Liquid Flat-Jet

We present chemical kinetics measurements of the luminol oxidation chemiluminescence (CL) reaction at the interface between two aqueous solutions, using liquid jet technology. Free-flowing liquid microjets are a relatively recent development that have found their way into a growing number of applications in spectroscopy and dynamics. A variant thereof, called flat-jet, is obtained when two cylindrical jets of a liquid are crossed, leading to a chain of planar leaf-shaped structures of the flowing liquid. We here show that in the first leaf of this chain, the fluids do not exhibit turbulent mixing, providing a clean interface between the liquids from the impinging jets. We also show, using the example of the luminol CL reaction, how this setup can be used to obtain kinetics information from friction-less flow and by circumventing the requirement for rapid mixing by intentionally suppressing all turbulent mixing and instead relying on diffusion.

Photoelectron circular dichroism in angle-resolved photoemission from liquid fenchone

We present an experimental X-ray photoelectron circular dichroism (PECD) study of liquid fenchone at the C 1s edge. A novel setup to enable PECD measurements on a liquid microjet [Malerz et al., Rev. Sci. Instrum., 2022, 93, 015101] was used. For the C 1s line assigned to fenchone's carbonyl carbon, a non-vanishing asymmetry is found in the intensity of photoelectron spectra acquired under a fixed angle in the backward-scattering plane. This experiment paves the way towards an innovative probe of the chirality of organic/biological molecules in aqueous solution.

We present the first X-ray photoelectron circular dichroism (PECD) study from a liquid phase sample, exemplified for liquid fenchone at the C 1s edge.

Simulating the Feasibility of Using Liquid Micro-Jets for Determining Electron–Liquid Scattering Cross-Sections

The extraction of electron–liquid phase cross-sections (surface and bulk) is proposed through the measurement of (differential) energy loss spectra for electrons scattered from a liquid micro-jet. The signature physical elements of the scattering processes on the energy loss spectra are highlighted using a Monte Carlo simulation technique, originally developed for simulating electron transport in liquids. Machine learning techniques are applied to the simulated electron energy loss spectra, to invert the data and extract the cross-sections. The extraction of the elastic cross-section for neon was determined within 9% accuracy over the energy range 1–100 eV. The extension toward the simultaneous determination of elastic and ionisation cross-sections resulted in a decrease in accuracy, now to within 18% accuracy for elastic scattering and 1% for ionisation. Additional methods are explored to enhance the accuracy of the simultaneous extraction of liquid phase cross-sections.

Real-Time Equation-of-Motion CCSD Cumulant Green's Function

Many-body excitations in X-ray photoemission spectra have been difficult to simulate from first principles. We have recently developed a cumulant-based one-electron Green's function method using the real-time coupled-cluster-singles equation-of-motion approach (RT-EOM-CCS) that provides a general framework for treating these problems. Here we extend this approach to include double excitations in the ground-state energy and in the coupled cluster amplitudes, which have been implemented using subroutines generated by the Tensor Contraction Engine (TCE). As in the case of the singles approximation, RT-EOM-CCSD yields a nonperturbative cumulant form of the Green's function in terms of the time-dependent cluster amplitudes, adding nonlinear corrections to the traditional cumulant forms. The extended approach is applied to the core-hole spectral function for small molecular systems. We find that, when core-optimized basis sets are used, the doubles contributions reduce the mean absolute errors in the core binding energies of the 10e systems from 0.8 to 0.3 eV. They also significantly improve the quasiparticle-satellite gap by reducing its overestimation from about 3-5 to about 0-1 eV in CH₄, NH₃, and H₂O, and also improving the overall shape of the satellite features. Finally, we demonstrate the application of the new implementation to the larger, classical XPS ESCA series of molecules and show that the singles approximation can be paired with a modest basis set to study carbon speciation.

Electronic structure of the solvated benzene radical anion

The benzene radical anion is a molecular ion pertinent to several organic reactions, including the Birch reduction of benzene in liquid ammonia. The species exhibits a dynamic Jahn-Teller effect due to its open-shell nature and undergoes pseudorotation of its geometry. Here, we characterize the complex electronic structure of this condensed-phase system based on ab initio molecular dynamics simulations and GW calculations of the benzene radical anion solvated in liquid ammonia. Using detailed analysis of the molecular and electronic structure, we find that the spatial character of the excess electron of the solvated radical anion follows the underlying Jahn-Teller distortions of the molecular geometry. We decompose the electronic density of states to isolate the contribution of the solute and to examine the response of the solvent to its presence. Our findings show the correspondence between instantaneous molecular structure and spin density; provide important insights into the electronic stability of the species, revealing that it is, indeed, a bound state in the condensed phase; and offer electronic densities of states that aid in the interpretation of experimental photoelectron spectra.

Photoelectron Spectroscopy of Benzene in the Liquid Phase and Dissolved in Liquid Ammonia

We report valence band photoelectron spectroscopy measurements of gas-phase and liquid-phase benzene as well as those of benzene dissolved in liquid ammonia, complemented by electronic structure calculations. The origins of the sizable gas-to-liquid-phase shifts in electron binding energies deduced from the benzene valence band spectral features are quantitatively characterized in terms of the Born-Haber solvation model. This model also allows to rationalize the observation of almost identical shifts in liquid ammonia and benzene despite the fact that the former solvent is polar while the latter is not. For neutral solutes like benzene, it is the electronic polarization response determined by the high frequency dielectric constant of the solvent, which is practically the same in the two liquids, that primarily determines the observed gas-to-liquid shifts.

Equation-of-Motion Coupled-Cluster Cumulant Green's Function for Excited States and X-Ray Spectra

Green’s function methods provide a robust, general framework within many-body theory for treating electron correlation in both excited states and x-ray spectra. Conventional methods using the Dyson equation or the cumulant expansion are typically based on the GW self-energy approximation. In order to extend this approximation in molecular systems, a non-perturbative real-time coupled-cluster cumulant Green’s function approach has been introduced, where the cumulant is obtained as the solution to a system of coupled first order, non-linear differential equations. This approach naturally includes non-linear corrections to conventional cumulant Green’s function techniques where the cumulant is linear in the GW self-energy. The method yields the spectral function for the core Green’s function, which is directly related to the x-ray photoemission spectra (XPS) of molecular systems. The approach also yields very good results for binding energies and satellite excitations. The x-ray absorption spectrum (XAS) is then calculated using a convolution of the core spectral function and an effective, one-body XAS. Here this approach is extended to include the full coupled-cluster-singles (CCS) core Green’s function by including the complete form of the non-linear contributions to the cumulant as well as all single, double, and triple cluster excitations in the CC amplitude equations. This approach naturally builds in orthogonality and shake-up effects analogous to those in the Mahan-Noizeres-de Dominicis edge singularity corrections that enhance the XAS near the edge. The method is illustrated for the XPS and XAS of NH₃.

Benzene Radical Anion Microsolvated in Ammonia Clusters: Modeling the Transition from an Unbound Resonance to a Bound Species

The benzene radical anion, well-known in organic chemistry as the first intermediate in the Birch reduction of benzene in liquid ammonia, exhibits intriguing properties from the point of view of quantum chemistry. Notably, it has the character of a metastable shape resonance in the gas phase, while measurements in solution find it to be experimentally detectable and stable. In this light, our previous calculations performed in bulk liquid ammonia explicitly reveal that solvation leads to stabilization. Here, we focus on the transition of the benzene radical anion from an unstable gas-phase ion to a fully solvated bound species by explicit ionization calculations of the radical anion solvated in molecular clusters of increasing size. The computational cost of the largest systems is mitigated by combining density functional theory with auxiliary methods including effective fragment potentials or approximating the bulk by polarizable continuum models. Using this methodology, we obtain the cluster size dependence of the vertical binding energy of the benzene radical anion converging to the value of -2.3 eV at a modest computational cost.

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.

Core level photoelectron spectroscopy of heterogeneous reactions at liquid-vapor interfaces: Current status, challenges, and prospects

Liquid-vapor interfaces, particularly those between aqueous solutions and air, drive numerous important chemical and physical processes in the atmosphere and in the environment. X-ray photoelectron spectroscopy is an excellent method for the investigation of these interfaces due to its surface sensitivity, elemental and chemical specificity, and the possibility to obtain information on the depth distribution of solute and solvent species in the interfacial region. In this Perspective, we review the progress that was made in this field over the past decades and discuss the challenges that need to be overcome for investigations of heterogeneous reactions at liquid-vapor interfaces under close-to-realistic environmental conditions. We close with an outlook on where some of the most exciting and promising developments might lie in this field.

Real-Time Coupled-Cluster Approach for the Cumulant Green's Function

Green's function methods within many-body perturbation theory provide a general framework for treating electronic correlations in excited states and spectra. Here, we develop the cumulant form of the one-electron Green's function using a real-time coupled-cluster equation-of-motion approach, in an extension of our previous study (Rehr J.; et al. J. Chem. Phys. 2020, 152, 174113). The approach yields a nonperturbative expression for the cumulant in terms of the solution to a set of coupled first-order, nonlinear differential equations. The method thereby adds nonlinear corrections to traditional cumulant methods, which are linear in the self-energy. The approach is applied to the core-hole Green's function and is illustrated for a number of small molecular systems. For these systems, we find that the nonlinear contributions yield significant improvements, both for quasiparticle properties such as core-level binding energies and for inelastic losses that correspond to satellites observed in photoemission spectra.

Benzene Radical Anion in the Context of the Birch Reduction: When Solvation Is the Key

The benzene radical anion is an important intermediate in the Birch reduction of benzene by solvated electrons in liquid ammonia. Beyond organic chemistry, it is an intriguing subject of spectroscopic and theoretical studies due to its rich structural and dynamical behavior. In the gas phase, the species appears as a metastable shape resonance, while in the condensed phase, it remains stable. Here, we approach the system by ab initio molecular dynamics in liquid ammonia and demonstrate that the inclusion of solvent is crucial and indeed leads to stability. Beyond the mere existence of the radical anion species, our simulations explore its condensed-phase behavior at the molecular level and offer new insights into its properties. These include the dynamic Jahn-Teller distortions, vibrational spectra in liquid ammonia, and the structure of the solvent shell, including the motif of a π-hydrogen bond between ammonia molecules and the aromatic ring.

Deeply cooled and temperature controlled microjets: Liquid ammonia solutions released into vacuum for analysis by photoelectron spectroscopy

A versatile, temperature controlled apparatus is presented, which generates deeply cooled liquid microjets of condensed gases, expelling them via a small aperture into vacuum for use in photoelectron spectroscopy (PES). The functionality of the design is demonstrated by temperature- and concentration-dependent PES measurements of liquid ammonia and solutions of KI and NH₄I in liquid ammonia. The experimental setup is not limited to the usage of liquid ammonia solutions solely.

Ionization Energy of Liquid Water Revisited

The ionization energy of liquid water is one of its most fundamental properties, an important benchmark for first-principles electronic-structure calculations and a crucial reference in the growing field of liquid-phase photoelectron spectroscopy. Despite this significance, a consensus on its value appears to be missing in the literature. Therefore, we use a monochromatized high-harmonic light source to perform detailed measurements of the ionization energy of liquid water in the presence of a tunable bias voltage applied to the liquid jet. Our results suggest that this simple method is sufficient to simultaneously compensate the effects of the streaming potential and that of the vacuum-level offset between the liquid and the photoelectron spectrometer. Our measurements yield corrected values of the vertical and adiabatic ionization energies of the 1b₁ band of bulk liquid water of 11.67(15) and 10.12(15) eV, respectively. Our method is broadly applicable and is likely to result in corrections to the measured ionization energies of solvated species as well.

Effective one-particle energies from generalized Kohn-Sham random phase approximation: A direct approach for computing and analyzing core ionization energies

Generalized-Kohn-Sham (GKS) orbital energies obtained self-consistently from the random phase approximation energy functional with a semicanonical projection (spRPA) were recently shown to rival the accuracy of GW quasiparticle energies for valence ionization potentials. Here, we extend the scope of GKS-spRPA correlated one-particle energies from frontier-orbital ionization to core orbital ionization energies, which are notoriously difficult for GW and other response methods due to strong orbital relaxation effects. For a benchmark consisting of 23 1s core electron binding energies (CEBEs) of second-row elements, chemical shifts estimated from GKS-spRPA one-particle energies yield mean absolute deviations from experiment of 0.2 eV, which are significantly more accurate than the standard GW and comparable to Δ self-consistent field theory without semiempirical adjustment of the energy functional. For small ammonia clusters and cytosine tautomers, GKS-spRPA based chemical shifts capture subtle variations in covalent and noncovalent bonding environments; GKS-spRPA 1s CEBEs for these systems agree with equation-of-motion coupled cluster singles and doubles and ADC(4) results within 0.2-0.3 eV. Two perturbative approximations to GKS-spRPA orbital energies, which reduce the scaling from O(N⁶) to O(N⁵) and O(N⁴), are introduced and tested. We illustrate the application of GKS-spRPA orbital energies to larger systems by using oxygen 1s CEBEs to probe solvation and packing effects in condensed phases of water. GKS-spRPA predicts a lowering of the oxygen 1s CEBE of approximately 1.6-1.7 eV in solid and liquid phases, consistent with liquid-jet X-ray photoelectron spectroscopy and gas phase cluster experiments. The results are rationalized by partitioning GKS-spRPA electron binding energies into static, relaxation, and correlation parts.