A nuclear configuration interaction approach to study nuclear spin effects: an application to ortho- and para-3 He2 @C60.

We introduce a non-orthogonal configuration interaction approach to investigate nuclear quantum effects on energies and densities of confined fermionic nuclei. The Hamiltonian employed draws parallels between confined systems and many-electron atoms, where effective non-Coulombic potentials represent the interactions of the trapped particles. One advantage of this method is its generality, as it offers the potential to study the nuclear quantum effects of various confined species affected by effective isotropic or anisotropic potentials. As a first application, we analyze the quantum states of two 3 He atoms encapsulated in C60 . At the Hartree-Fock level, we observe the breaking of spin and spatial symmetries. To ensure wavefunctions with the correct symmetries, we mix the broken-symmetry Hartree-Fock states within the non-orthogonal configuration interaction expansion. Our proposed approach predicts singly and triply degenerate ground states for the singlet (para-3 He2 @C60 ) and triplet (ortho-3 He2 @C60 ) nuclear spin configurations, respectively. The ortho-3 He2 @C60 ground state is 5.69 cm-1 higher in energy than the para-3 He2 @C60 ground state. The nuclear densities obtained for these states exhibit the icosahedral symmetry of the C60 embedding potential. Importantly, our calculated energies for the lowest 85 states are in close agreement with perturbation theory results based on a harmonic oscillator plus rigid rotor model of 3 He2 @C60 .

TDDFT and the x-ray absorption spectrum of liquid water: Finding the "best" functional.

We investigate the performance of time-dependent density functional theory (TDDFT) for reproducing high-level reference x-ray absorption spectra of liquid water and water clusters. For this, we apply the integrated absolute difference (IAD) metric, previously used for x-ray emission spectra of liquid water [T. Fransson and L. G. M. Pettersson, J. Chem. Theory Comput. 19, 7333-7342 (2023)], in order to investigate which exchange-correlation (xc) functionals yield TDDFT spectra in best agreement to reference, as well as to investigate the suitability of IAD for x-ray absorption spectroscopy spectrum calculations. Considering highly asymmetric and symmetric six-molecule clusters, it is seen that long-range corrected xc-functionals are required to yield good agreement with the reference coupled cluster (CC) and algebraic-diagrammatic construction spectra, with 100% asymptotic Hartree-Fock exchange resulting in the lowest IADs. The xc-functionals with best agreement to reference have been adopted for larger water clusters, yielding results in line with recently published CC theory, but which still show some discrepancies in the relative intensity of the features compared to experiment.

Symmetry-resolved CO desorption and oxidation dynamics on O/Ru(0001) probed at the C K-edge by ultrafast x-ray spectroscopy.

We report on carbon monoxide desorption and oxidation induced by 400 nm femtosecond laser excitation on the O/Ru(0001) surface probed by time-resolved x-ray absorption spectroscopy (TR-XAS) at the carbon K-edge. The experiments were performed under constant background pressures of CO (6 × 10-8 Torr) and O2 (3 × 10-8 Torr). Under these conditions, we detect two transient CO species with narrow 2π* peaks, suggesting little 2π* interaction with the surface. Based on polarization measurements, we find that these two species have opposing orientations: (1) CO favoring a more perpendicular orientation and (2) CO favoring a more parallel orientation with respect to the surface. We also directly detect gas-phase CO2 using a mass spectrometer and observe weak signatures of bent adsorbed CO2 at slightly higher x-ray energies than the 2π* region. These results are compared to previously reported TR-XAS results at the O K-edge, where the CO background pressure was three times lower (2 × 10-8 Torr) while maintaining the same O2 pressure. At the lower CO pressure, in the CO 2π* region, we observed adsorbed CO and a distribution of OC-O bond lengths close to the CO oxidation transition state, with little indication of gas-like CO. The shift toward "gas-like" CO species may be explained by the higher CO exposure, which blocks O adsorption, decreasing O coverage and increasing CO coverage. These effects decrease the CO desorption barrier through dipole-dipole interaction while simultaneously increasing the CO oxidation barrier.

Electrocatalytic Glycerol Oxidation with Concurrent Hydrogen Evolution Utilizing an Efficient MoOx /Pt Catalyst.

Glycerol electrolysis affords a green and energetically favorable route for the production of value-added chemicals at the anode and H2 production in parallel at the cathode. Here, a facile method for trapping Pt nanoparticles at oxygen vacancies of molybdenum oxide (MoOx ) nanosheets, yielding a high-performance MoOx /Pt composite electrocatalyst for both the glycerol oxidation reaction (GOR) and the hydrogen evolution reaction (HER) in alkaline electrolytes, is reported. Combined electrochemical experiments and theoretical calculations reveal the important role of MoOx nanosheets for the adsorption of glycerol molecules in GOR and the dissociation of water molecules in HER, as well as the strong electronic interaction with Pt. The MoOx /Pt composite thus significantly enhances the specific mass activity of Pt and the kinetics for both reactions. With MoOx /Pt electrodes serving as both cathode and anode, two-electrode glycerol electrolysis is achieved at a cell voltage of 0.70 V to reach a current density of 10 mA cm-2 , which is 0.90 V less than that required for water electrolysis.

Nanoscale Spatial Distribution of Supported Nanoparticles Controls Activity and Stability in Powder Catalysts for CO Oxidation and Photocatalytic H2 Evolution.

Supported metal nanoparticles are essential components of high-performing catalysts, and their structures are intensely researched. In comparison, nanoparticle spatial distribution in powder catalysts is conventionally not quantified, and the influence of this collective property on catalyst performance remains poorly investigated. Here, we demonstrate a general colloidal self-assembly method to control uniformity of nanoparticle spatial distribution on common industrial powder supports. We quantify distributions on the nanoscale using image statistics and show that the type of nanospatial distribution determines not only the stability, but also the activity of heterogeneous catalysts. Widely investigated systems (Au-TiO2 for CO oxidation thermocatalysis and Pd-TiO2 for H2 evolution photocatalysis) were used to showcase the universal importance of nanoparticle spatial organization. Spatially and temporally resolved microkinetic modeling revealed that nonuniformly distributed Au nanoparticles suffer from local depletion of surface oxygen, and therefore lower CO oxidation activity, as compared to uniformly distributed nanoparticles. Nanoparticle spatial distribution also determines the stability of Pd-TiO2 photocatalysts, because nonuniformly distributed nanoparticles sinter while uniformly distributed nanoparticles do not. This work introduces new tools to evaluate and understand catalyst collective (ensemble) properties in powder catalysts, which thereby pave the way to more active and stable heterogeneous catalysts.

Time-resolved observation of transient precursor state of CO on Ru(0001) using carbon K-edge spectroscopy.

The transient dynamics of carbon monoxide (CO) molecules on a Ru(0001) surface following femtosecond optical laser pump excitation has been studied by monitoring changes in the unoccupied electronic structure using an ultrafast X-ray free-electron laser (FEL) probe. The particular symmetry of perpendicularly chemisorbed CO on the surface is exploited to investigate how the molecular orientation changes with time by varying the polarization of the FEL pulses. The time evolution of spectral features corresponding to the desorption precursor state was well distinguished due to the narrow line-width of the C K-edge in the X-ray absorption (XA) spectrum, illustrating that CO molecules in the precursor state rotated freely and resided on the surface for several picoseconds. Most of the CO molecules trapped in the precursor state ultimately cooled back down to the chemisorbed state, while we estimate that ∼14.5 ± 4.9% of the molecules in the precursor state desorbed into the gas phase. It was also observed that chemisorbed CO molecules diffused over the metal surface from on-top sites toward highly coordinated sites. In addition, a new "vibrationally hot precursor" state was identified in the polarization-dependent XA spectra.

Diffusive dynamics during the high-to-low density transition in amorphous ice.

Water exists in high- and low-density amorphous ice forms (HDA and LDA), which could correspond to the glassy states of high- (HDL) and low-density liquid (LDL) in the metastable part of the phase diagram. However, the nature of both the glass transition and the high-to-low-density transition are debated and new experimental evidence is needed. Here we combine wide-angle X-ray scattering (WAXS) with X-ray photon-correlation spectroscopy (XPCS) in the small-angle X-ray scattering (SAXS) geometry to probe both the structural and dynamical properties during the high-to-low-density transition in amorphous ice at 1 bar. By analyzing the structure factor and the radial distribution function, the coexistence of two structurally distinct domains is observed at T = 125 K. XPCS probes the dynamics in momentum space, which in the SAXS geometry reflects structural relaxation on the nanometer length scale. The dynamics of HDA are characterized by a slow component with a large time constant, arising from viscoelastic relaxation and stress release from nanometer-sized heterogeneities. Above 110 K a faster, strongly temperature-dependent component appears, with momentum transfer dependence pointing toward nanoscale diffusion. This dynamical component slows down after transition into the low-density form at 130 K, but remains diffusive. The diffusive character of both the high- and low-density forms is discussed among different interpretations and the results are most consistent with the hypothesis of a liquid-liquid transition in the ultraviscous regime.

X-ray absorption spectrum simulations of hexagonal ice.

We calibrate basis sets and performance of two theoretical approaches to compute X-ray absorption spectra (XAS) of condensed water by comparison to experiments on hexagonal ice Ih. We apply both the transition-potential half-core-hole approach and the complex polarization propagator using four different models of the crystal with increasing oxygen and proton disorder but find poor agreement with experiments. We note that there are large variations in experimental spectra depending on detection mode and how the ice samples were prepared, which leads us to critically investigate what structures were actually prepared and measured in each case. This is done by using a Monte Carlo-based fitting technique which fits the spectra based on a library of precomputed spectra and assigns weights to contributions from different model structures. These are then used to generate O-O and O-H radial distribution functions and tetrahedrality parameters associated with each of the measured spectra. We find that all spectra are associated with sharp peaks at the oxygen positions in the perfect lattice, but with significant disorder around these positions. We suggest that presently available XAS of hexagonal ice are not fully representative of the perfect crystalline lattice, but contain varying amounts of defects and possible contributions from low-density amorphous ice.

Translational and rotational dynamics of high and low density TIP4P/2005 water.

We use molecular dynamics simulations using TIP4P/2005 to investigate the self- and distinct-van Hove functions for different local environments of water, classified using the local structure index as an order parameter. The orientational dynamics were studied through the calculation of the time-correlation functions of different-order Legendre polynomials in the OH-bond unit vector. We found that the translational and orientational dynamics are slower for molecules in a low-density local environment and correspondingly the mobility is enhanced upon increasing the local density, consistent with some previous works, but opposite to a recent study on the van Hove function. From the analysis of the distinct dynamics, we find that the second and fourth peaks of the radial distribution function, previously identified as low density-like arrangements, show long persistence in time. The analysis of the time-dependent interparticle distance between the central molecule and the first coordination shell shows that particle identity persists longer than distinct van Hove correlations. The motion of two first-nearest-neighbor molecules thus remains coupled even when this correlation function has been completely decayed. With respect to the orientational dynamics, we show that correlation functions of molecules in a low-density environment decay exponentially, while molecules in a local high-density environment exhibit bi-exponential decay, indicating that dynamic heterogeneity of water is associated with the heterogeneity among high-density and between high-density and low-density species. This bi-exponential behavior is associated with the existence of interstitial waters and the collapse of the second coordination sphere in high-density arrangements, but not with H-bond strength.

Radial distribution functions of water: Models vs experiments.

We study the temperature behavior of the first four peaks of the oxygen-oxygen radial distribution function of water, simulated by the TIP4P/2005, MB-pol, TIP5P, and SPC/E models and compare to experimental X-ray diffraction data, including a new measurement which extends down to 235 K [H. Pathak et al., J. Chem. Phys. 150, 224506 (2019)]. We find the overall best agreement using the MB-pol and TIP4P/2005 models. We observe, upon cooling, a minimum in the position of the second shell simulated with TIP4P/2005 and SPC/E potentials, located close to the temperature of maximum density. We also calculated the two-body entropy and the contributions coming from the first, second, and outer shells to this quantity. We show that, even if the main contribution comes from the first shell, the contribution of the second shell can become important at low temperature. While real water appears to be less ordered at short distance than obtained by any of the potentials, the different water potentials show more or less order compared to the experiments depending on the considered length-scale.

A proposal for the structure of high- and low-density fluctuations in liquid water.

Based on recent experimental data that can be interpreted as indicating the presence of specific structures in liquid water, we build and optimize two structural models which we compare with the available experimental data. To represent the proposed high-density liquid structures, we use a model consisting of chains of water molecules, and for low-density liquid, we investigate fused dodecahedra as templates for tetrahedral fluctuations. The computed infrared spectra of the models are in very good agreement with the extracted experimental spectra for the two components, while the extracted structures from molecular dynamics (MD) simulations give spectra that are intermediate between the experimentally derived spectra. Computed x-ray absorption and emission spectra as well as the O-O radial distribution functions of the proposed structures are not contradicted by experiment. The stability of the proposed dodecahedral template structures is investigated in MD simulations by seeding the starting structure, and remnants found to persist on an ∼30 ps time scale. We discuss the possible significance of such seeds in simulations and whether they can be viable candidates as templates for structural fluctuations below the compressibility minimum of liquid water.

A Chemical View on X-ray Photoelectron Spectroscopy: the ESCA Molecule and Surface-to-Bulk XPS Shifts.

In this paper we remind the reader of a simple, intuitive picture of chemical shifts in X-ray photoelectron spectroscopy (XPS) as the difference in chemical bonding between the probed atom and its neighbor to the right in the periodic table, the so called Z+1 approximation. We use the classical ESCA molecule, ethyl trifluoroacetate, and 4d-transition metals to explicitly demonstrate agreement between core-level shifts computed as differences between final core-hole states and the approach where each core-ionized atom is replaced by a Z+1 atom. In this final state, or total energy picture, the XPS shift arises due to the more or less unfavorable chemical bonding of the effective nitrogen in the carbon geometry for the ESCA molecule. Surface core level shifts in metals are determined by whether the Z+1 atom as an alloy segregates to the surface or is more soluble in the bulk. As further illustration of this more chemical picture, we compare the geometry of C 1s and O 1s core-ionized CO with that of, respectively, NO+ and CF+ . The scope is not to propose a new method to compute XPS shifts but rather to stress the validity of this simple interpretation.

X-ray and Electron Spectroscopy of Water.

Here we present an overview of recent developments of X-ray and electron spectroscopy to probe water at different temperatures. Photon-induced ionization followed by detection of electrons from either the O 1s level or the valence band is the basis of photoelectron spectroscopy. Excitation between the O 1s and the unoccupied states or occupied states is utilized in X-ray absorption and X-ray emission spectroscopies. These techniques probe the electronic structure of the liquid phase and show sensitivity to the local hydrogen-bonding structure. Both experimental aspects related to the measurements and theoretical simulations to assist in the interpretation are discussed in detail. Different model systems are presented such as the different bulk phases of ice and various adsorbed monolayer structures on metal surfaces.

Relationship between x-ray emission and absorption spectroscopy and the local H-bond environment in water.

The connection between specific features in the water X-ray absorption spectrum and X-ray emission spectrum (XES) and the local H-bond coordination is studied based on structures obtained from path-integral molecular dynamics simulations using either the opt-PBE-vdW density functional or the MB-pol force field. Computing the XES spectrum using all molecules in a snapshot results in only one peak in the lone-pair (1b1) region, while the experiment shows two peaks separated by 0.8-0.9 eV. Different H-bond configurations were classified based on the local structure index (LSI) and a geometrical H-bond cone criterion. We find that tetrahedrally coordinated molecules characterized by high LSI values and two strong donated and two strong accepted H-bonds contribute to the low energy 1b1 emission peak and to the post-edge region in absorption. Molecules with the asymmetric H-bond environment with one strong accepted H-bond and one strong donated H-bond and low LSI values give rise to the high energy 1b1 peak in the emission spectrum and mainly contribute to the pre-edge and main-edge in the absorption spectrum. The 1b1 peak splitting can be increased to 0.62 eV by imposing constraints on the H-bond length, i.e., for very tetrahedral structures short H-bonds (less than 2.68 Å) and for very asymmetric structures elongated H-bonds (longer than 2.8 Å). Such structures are present, but underrepresented, in the simulations which give more of an average of the two extremes.

Atom-specific activation in CO oxidation.

We report on atom-specific activation of CO oxidation on Ru(0001) via resonant X-ray excitation. We show that resonant 1s core-level excitation of atomically adsorbed oxygen in the co-adsorbed phase of CO and oxygen directly drives CO oxidation. We separate this direct resonant channel from indirectly driven oxidation via X-ray induced substrate heating. Based on density functional theory calculations, we identify the valence-excited state created by the Auger decay as the driving electronic state for direct CO oxidation. We utilized the fresh-slice multi-pulse mode at the Linac Coherent Light Source that provided time-overlapped and 30 fs delayed pairs of soft X-ray pulses and discuss the prospects of femtosecond X-ray pump X-ray spectroscopy probe, as well as X-ray two-pulse correlation measurements for fundamental investigations of chemical reactions via selective X-ray excitation.

Temperature-Independent Nuclear Quantum Effects on the Structure of Water.

Nuclear quantum effects (NQEs) have a significant influence on the hydrogen bonds in water and aqueous solutions and have thus been the topic of extensive studies. However, the microscopic origin and the corresponding temperature dependence of NQEs have been elusive and still remain the subject of ongoing discussion. Previous x-ray scattering investigations indicate that NQEs on the structure of water exhibit significant temperature dependence [Phys. Rev. Lett. 94, 047801 (2005)PRLTAO0031-900710.1103/PhysRevLett.94.047801]. Here, by performing wide-angle x-ray scattering of H_{2}O and D_{2}O droplets at temperatures from 275 K down to 240 K, we determine the temperature dependence of NQEs on the structure of water down to the deeply supercooled regime. The data reveal that the magnitude of NQEs on the structure of water is temperature independent, as the structure factor of D_{2}O is similar to H_{2}O if the temperature is shifted by a constant 5 K, valid from ambient conditions to the deeply supercooled regime. Analysis of the accelerated growth of tetrahedral structures in supercooled H_{2}O and D_{2}O also shows similar behavior with a clear 5 K shift. The results indicate a constant compensation between NQEs delocalizing the proton in the librational motion away from the bond and in the OH stretch vibrational modes along the bond. This is consistent with the fact that only the vibrational ground state is populated at ambient and supercooled conditions.

MP4 Study of the Anharmonic Coupling of the Shared Proton Stretching Vibration of the Protonated Water Dimer in Equilibrium and Transition States.

The structure and harmonic and anharmonic IR spectra of the protonated water dimer (PWD) were calculated in C1, C2, and Cs symmetry at the MP4/acc-pVTZ level of theory. We found that structure and IR spectra are practically identical in C2 and C1 symmetry, demonstrating that an equilibrium C1 configuration of the PWD is not realized. Anharmonic coupling of the shared proton stretching vibration with all other modes in the PWD in C2 and Cs symmetry was the focus of this investigation. For this purpose, 28 two-dimensional potential energy surfaces (2D PES) were built at the MP4/acc-pVTZ level of theory and the corresponding vibrational Schrödinger equations were solved using the DVR method. Differences in the coupling of the investigated mode with other modes in the C2 and Cs configurations, along with some factors that determine the red- or blue-shift of the stretching vibration frequency, were analyzed. We obtained a rather reasonable value of the stretching frequency of the bridging proton (1058.4 cm-1) unperturbed by Fermi resonance. The Fermi resonance between the fundamental vibration ν7 and the combined vibration ν2 + ν6 of the same symmetry was analyzed through anharmonic second-order perturbation theory calculations, as well as by 3D PES constructed using Q2, Q6, and Q7 as normal coordinates. A significant (up to 50%) transfer of intensity from the fundamental vibration to the combined one was found. We have estimated the frequency of the bridging proton stretching vibration in the Cs configuration of the PWD based on calculations of the intrinsic anharmonicity and anharmonic double modes interactions at the MP4/acc-pVTZ level of theory (1261 cm-1).

Requirements of first-principles calculations of X-ray absorption spectra of liquid water.

A computational benchmark study on X-ray absorption spectra of water has been performed by means of transition-potential density functional theory (TP-DFT), damped time-dependent density functional theory (TDDFT), and damped coupled cluster (CC) linear response theory. For liquid water, using TDDFT with a tailored CAM-B3LYP functional and a polarizable embedding, we find that an embedding with over 2000 water molecules is required to fully converge spectral features for individual molecules, but a substantially smaller embedding can be used within averaging schemes. TP-DFT and TDDFT calculations on 100 MD structures demonstrate that TDDFT produces a spectrum with spectral features in good agreement with experiment, while it is more difficult to fully resolve the spectral features in the TP-DFT spectrum. Similar trends were also observed for calculations of bulk ice. In order to further establish the performance of these methods, small water clusters have been considered also at the CC2 and CCSD levels of theory. Issues regarding the basis set requirements for spectrum simulations of liquid water and the determination of gas-phase ionization potentials are also discussed.

Evaporative cooling of microscopic water droplets in vacuo: Molecular dynamics simulations and kinetic gas theory.

In the present study, we investigate the process of evaporative cooling of nanometer-sized droplets in vacuum using molecular dynamics simulations with the TIP4P/2005 water model. The results are compared to the temperature evolution calculated from the Knudsen theory of evaporation which is derived from kinetic gas theory. The calculated and simulation results are found to be in very good agreement for an evaporation coefficient equal to unity. Our results are of interest to experiments utilizing droplet dispensers as well as to cloud micro-physics.