Real-Time Dynamics of the Formation of Hydrated Electrons upon Irradiation of Water Clusters with Extreme Ultraviolet Light

Two rhombic ice phases from aqueous salt solutions under graphene confinement

Water exhibits rich ice phases depending upon its respective formation conditions, and in particular, the two-dimensional ice with nonhexagonal symmetry adsorbed on solids relates to the exceptional arrangement of water molecules. Despite extensive reporting of two-dimensional ice on various solid surfaces, the geometry and thermodynamics of ice formation from an aqueous salt solution are still unknown. In this Letter, we show the formation of single- and two-phase mixed two-dimensional rhombic ice from aqueous salt solutions with different concentrations under strong compressed confinement of graphene at ambient temperature by using classical molecular dynamics simulations and first-principles calculations. The two rhombic ice phases exhibit identical geometry and thermodynamic properties, but different projections of the oxygen atoms against solid surface symmetry, where they relate to the stable and metastable arrangements of water molecules confined between two graphene layers. A single-phase rhombic ice would grow from the confined saturated aqueous solutions since the previously stable rhombic molecular arrangement becomes an unstable high-energy state by introducing salt ions nearby. Our result reveals different rhombic ice phases growing from pure water and aqueous solutions, highlighting the deciding role of salt ions in the ice formation process due to their common presence in liquids.

Method of kinetic energy reconstruction from time-of-flight mass spectra

We present a method for the reconstruction of ion kinetic energy distributions from ion time-of-flight mass spectra through ion trajectory simulations. In particular, this method is applicable to complicated spectrometer geometries with largely anisotropic ion collection efficiencies. A calibration procedure using a single ion mass peak allows the accurate determination of parameters related to the spectrometer calibration, experimental alignment, and instrument response function, which improves the agreement between simulations and experiment. The calibrated simulation is used to generate a set of basis functions for the time-of-flight spectra, which are then used to transform from time-of-flight to kinetic-energy spectra. We demonstrate this reconstruction method on a recent pump-probe experiment by Asmussen et al. [Asmussen et al., Phys. Chem. Chem. Phys., 23, 15138, (2021)] on helium nanodroplets and retrieve time-resolved kinetic-energy-release spectra for the ions from ion time-of-flight spectra.

Theoretical investigation of Aryl/Alkyl halide reduction with hydrated electrons from energy and AIMD aspects

CONTEXT

In this study, the reactions of hydrated electron (e^-(aq)) with alkyl and aryl halides were simulated with an ab initial molecular dynamics (AIMD) method to reveal the underlying mechanism. An original protocol was developed for preparing the proper initial wavefunction guess of AIMD, in which a single electron was curled in a tetrahedral cavity of four water molecules. Our results show that the stability of e^-(aq) increases with the hydrogen bond grid integrity. The organic halides prefer to react with e^-(aq) in neutral or alkaline environment, while they are more likely to react with hydrogen radical (the product of e^-(aq) and proton) under acidic conditions. The reaction between fluorobenzene/fluoromethane and hydrogen radical is considered as the least favorable reaction due to the highest reaction barriers. The bond dissociation energy (BDE) suggested that the cleavage of the carbon-halogen bond of their anion radical might be a thermodynamically favorable reaction. AIMD results indicated that the LUMO or higher orbitals were the e^-(aq) migration destination. The transplanted electron enhanced carbon-halogen bond vibration intensively, leading to bond cleavage. The solvation process of the departing halogen anions was observed in both fluorobenzene and fluoromethane AIMD simulation, indicating that it might have a significant effect on enthalpy. Side reactions and byproducts obtained during the AIMD simulation suggested the complexity of the e^-(aq) reactions and further investigation was needed to fully understand the reaction mechanisms. This study provided theoretical insight into the pollutant environmental fate and constructed a methodological foundation for AIMD simulation of analogous free radical reactions.

METHODS

The theoretical calculation was conducted on the combination of Gaussian16 and ORCA5.0.3 software packages. The initial geometries, as well as the wavefunction initial guesses, were obtained at PBE0/ma-def2-TZVP/IEFPCM-water level in Gaussian16 unless otherwise stated. AIMD simulations were performed at the same level in ORCA. Wavefunction analysis was carried out with Multiwfn. The details methods were described in the section "Computational details" section.

Scavenging of “dry” electrons prior to hydration by azide ions: effect on the formation of H2 in the radiolysis of water by 60Co γ-rays and tritium β-electrons

In this study, we use Monte Carlo track chemistry simulations to show that “dry” secondary electrons, precursors of the “hydrated” electron (e−aq), can be scavenged on the sub-picosecond time scale prior to hydration, by a high concentration (>0.1–1 M) of azide ions (N3−) in water irradiated with 60Co γ-rays and tritium β-electrons at 25 °C. This is a striking result, as N3− is known to react very slowly with e−aq. These processes tend to significantly reduce the yields of H2 as observed experimentally. For both energetic Compton electrons (“linear energy transfer”, LET ∼ 0.3 keV/µm), which are generated by the cobalt-60 γ-rays, and 3H β-electrons (LET ∼ 6 keV/µm), our H2 yield results confirm previous Monte Carlo simulations, which indicated the necessity of including the capture of the precursors to e−aq. Interestingly, our calculations show no significant changes in the scavenging of “dry” electrons at high azide concentrations in passing from γ-radiolysis to tritium β-radiolysis (i.e., with LET). This led us to the conclusion that the higher H2 yield observed experimentally for 3H β-electrons compared with 60Co γ-rays is mainly explained by the difference in the radiation track structures during the chemical stage (>1 ps). The higher LET of tritium β-electrons leads to more molecular products (H2 in this case) in tritium radiolysis than in γ-radiolysis. Finally, a value of ∼0.5 nm was derived for the reaction distance between N3− and the “dry” electron from the H2 yields observed in 60Co γ-radiolysis at high N3− concentrations.

Structure, energy, and bonding in anionic water tetramers obtained by exhaustive search

An analysis of the structures, some energy related properties, and key aspects of the bonding nature of the microsolvated electron with four water molecules is presented. The study is based on an exhaustive potential energy surface scan of the ground state of (H₂O)₄ ^- at the UCCSD(T)/6-311(3+,4+)G(d,p)//UMP2/6-311(3+,4+)G(d,p) level. A total of 18 structures, most of them not reported before, spanning in an energy range of 8.8 kcal mol^-1 were found. The energetic stability of the clusters is dictated by the effect of the excess electron on their structures, on their partial fragmentation, and on the hydrogen bonds' framework. Vertical detachment energies depend on the number of water molecules holding the excess electron in "direct contact" to their two protons at the same time and, to a lesser extent, also depend on the hydrogen bond sequence in the rest of the structure. In general, hydrogen bonds in (H₂O)₄ ^- are of closed shell character, and there are other less common interactions assisted by the excess electron.

Real-Time Observation of Solvation Dynamics of Electron in Actinide Extraction Binary Solutions of Water and n-Tributyl Phosphate

The excess electron in solution is a highly reactive radical involved in various radiation-induced reactions. Its solvation state critically determines the subsequent pathway and rate of transfer. For instance, water plays a dominating role in the electron-induced dealkylation of n-tributyl phosphate in actinide extraction processing. However, the underlying electron solvation processes in such systems are lacking. Herein, we directly observed the solvation dynamics of electrons in H-bonded water and n-tributyl phosphate (TBP) binary solutions with a mole fraction of water (X_w) varying from 0.05 to 0.51 under ambient conditions. Following the evolution of the absorption spectrum of trapped electrons (not fully solvated) with picosecond resolution, we show that electrons statistically distributed would undergo preferential solvation within water molecules extracted in TBP. We determine the time scale of excess electron full solvation from the deconvoluted transient absorption-kinetical data. The process of solvent reorganization accelerates by increasing the water molar fraction, and the rate of this process is 2 orders of magnitude slower compared to bulk water. We assigned the solvation process to hydrogen network reorientation induced by a negative charge of the excess electron that strongly depends on the local water environment. Our findings suggest that water significantly stabilizes the electron in a deeper potential than the pure TBP case. In its new state, the electron is likely to inhibit the dealkylation of extractants in actinide separation.

Ultra-high dose-rate (FLASH) radiotherapy: Generation of early, transient, strongly acidic spikes in the irradiated tumor environment

Monte Carlo simulations of γ/fast electron-radiolysis of water show that the in situ formation of H₃O⁺ temporarily renders each "native" isolated spur/track region very acidic. For pulsed (FLASH) irradiation with high dose rate, this early time, transient "acid-spike" response is shown to extend evenly across the entire irradiated volume. Since pH controls many cellular processes, this study highlights the need to consider these spikes of acidity in understanding the fundamental mechanisms underlying FLASH radiotherapy.

Time-resolved formation of excited atomic and molecular states in XUV-induced nanoplasmas in ammonia clusters

High intensity XUV radiation from a free-electron laser (FEL) was used to create a nanoplasma inside ammonia clusters with the intent of studying the resulting electron-ion interactions and their interplay with plasma evolution. In a plasma-like state, electrons with kinetic energy lower than the local collective Coulomb potential of the positive ionic core are trapped in the cluster and take part in secondary processes (e.g. electron-impact excitation/ionization and electron-ion recombination) which lead to subsequent excited and neutral molecular fragmentation. Using a time-delayed UV laser, the dynamics of the excited atomic and molecular states are probed from -0.1 ps to 18 ps. We identify three different phases of molecular fragmentation that are clearly distinguished by the effect of the probe laser on the ionic and electronic yield. We propose a simple model to rationalize our data and further identify two separate channels leading to the formation of excited hydrogen.

Ultrafast relaxation of photoexcited superfluid He nanodroplets

The relaxation of photoexcited nanosystems is a fundamental process of light–matter interaction. Depending on the couplings of the internal degrees of freedom, relaxation can be ultrafast, converting electronic energy in a few fs, or slow, if the energy is trapped in a metastable state that decouples from its environment. Here, we study helium nanodroplets excited resonantly by femtosecond extreme-ultraviolet (XUV) pulses from a seeded free-electron laser. Despite their superfluid nature, we find that helium nanodroplets in the lowest electronically excited states undergo ultrafast relaxation. By comparing experimental photoelectron spectra with time-dependent density functional theory simulations, we unravel the full relaxation pathway: Following an ultrafast interband transition, a void nanometer-sized bubble forms around the localized excitation (He\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${}^{* }$$\end{document}*) within 1 ps. Subsequently, the bubble collapses and releases metastable He\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${}^{* }$$\end{document}* at the droplet surface. This study highlights the high level of detail achievable in probing the photodynamics of nanosystems using tunable XUV pulses.

There is interest in understanding the relaxation mechanisms of photoexcitation in atoms, molecules and other complex systems. Here the authors unravel the photoexcitation and ultrafast relaxation of superfluid helium nanodroplets using a pump-probe experiment with FEL pulses.

Real-time observation of water radiolysis and hydrated electron formation induced by extreme-ultraviolet pulses

The dominant pathway of radiation damage begins with the ionization of water. Thus far, however, the underlying primary processes could not be conclusively elucidated. Here, we directly study the earliest steps of extreme ultraviolet (XUV)-induced water radiolysis through one-photon excitation of large water clusters using time-resolved photoelectron imaging. Results are presented for H2O and D2O clusters using femtosecond pump pulses centered at 133 or 80 nm. In both excitation schemes, hydrogen or proton transfer is observed to yield a prehydrated electron within 30 to 60 fs, followed by its solvation in 0.3 to 1.0 ps and its decay through geminate recombination on a ∼10-ps time scale. These results are interpreted by comparison with detailed multiconfigurational non-adiabatic ab-initio molecular dynamics calculations. Our results provide the first comprehensive picture of the primary steps of radiation chemistry and radiation damage and demonstrate new approaches for their study with unprecedented time resolution.

Relaxation Dynamics and Genuine Properties of the Solvated Electron in Neutral Water Clusters

We have investigated the solvation dynamics and the genuine binding energy and photoemission anisotropy of the solvated electron in neutral water clusters with a combination of time-resolved photoelectron velocity map imaging and electron scattering simulations. The dynamics was probed with a UV probe pulse following above-band-gap excitation by an EUV pump pulse. The solvation dynamics is completed within about 2 ps. Only a single band is observed in the spectra, with no indication for isomers with distinct binding energies. Data analysis with an electron scattering model reveals a genuine binding energy in the range of 3.55-3.85 eV and a genuine anisotropy parameter in the range of 0.51-0.66 for the ground-state hydrated electron. All of these observations coincide with those for liquid bulk, which is rather unexpected for an average cluster size of 300 molecules.

MELEXIR: maximum entropy Legendre expanded image reconstruction. A fast and efficient method for the analysis of velocity map imaging or photoelectron imaging data

The MELEXIR program obtains a Legendre expansion of the 3D velocity distribution from 2D images of ions or photoelectrons. The maximum entropy algorithm avoids inverse Abel transforms, is fast and applicable to low-intensity images.