Recombination of the Hydrated Electron at High Temperature and Pressure in Hydrogenated Alkaline Water

GEANT4-DNA simulation of temperature-dependent and pH-dependent yields of chemical radiolytic species

Objective.GEANT4-DNA can simulate radiation chemical yield (G-value) for radiolytic species such as the hydrated electron (eaq-) with the independent reaction times (IRT) method, however, only at room temperature and neutral pH. This work aims to modify the GEANT4-DNA source code to enable the calculation ofG-values for radiolytic species at different temperatures and pH values.Approach.In the GEANT4-DNA source code, values of chemical parameters such as reaction rate constant, diffusion coefficient, Onsager radius, and water density were replaced by corresponding temperature-dependent polynomials. The initial concentration of hydrogen ion (H+)/hydronium ion (H3O+) was scaled for a desired pH using the relationship pH = -log10[H+]. To validate our modifications, two sets of simulations were performed. (A) A water cube with 1.0 km sides and a pH of 7 was irradiated with an isotropic electron source of 1 MeV. The end time was 1μs. The temperatures varied from 25 °C to 150 °C. (B) The same setup as (A) was used, however, the temperature was set to 25 °C while the pH varied from 5 to 9. The results were compared with published experimental and simulated work.Main results.The IRT method in GEANT4-DNA was successfully modified to simulateG-values for radiolytic species at different temperatures and pH values. Our temperature-dependent results agreed with experimental data within 0.64%-9.79%, and with simulated data within 3.52%-12.47%. The pH-dependent results agreed well with experimental data within 0.52% to 3.19% except at a pH of 5 (15.99%) and with simulated data within 4.40%-5.53%. The uncertainties were below ±0.20%. Overall our results agreed better with experimental than simulation data.Significance.Modifications in the GEANT4-DNA code enabled the calculation ofG-values for radiolytic species at different temperatures and pH values.

Molecular Dynamics Characterization of Dielectron Hydration in Liquid Water with Unique Double Proton Transfers

Radiation chemistry of water and aqueous solutions has always been an interesting scientific issue owing to involving electronic excitations, ionization of solvated species, and formation of radiolytic species and many elementary reactions, but the underlying mechanisms are still poorly understood. Here, we for the first time molecular dynamics characterize the hydration dynamics of two correlated electrons and their triggered unique phenomena in liquid water associated with radiolysis of water using the combined hybrid functional and nonlocal dispersion functional. Hydration of two electrons may experience two distinctly different mechanisms, one forming a spin-paired closed-shell unicaged dielectron hydrate (e₂^2-_aq) and the other forming a spin-paired metastable open-shell bicaged hydrated electron pair (e^-_aq···e^-_aq) which exhibits intriguing antiferromagnetic spin coupling dynamics (in a range of -40 cm^-1 to -500 cm^-1). e^-_aq···e^-_aq can recombine to e₂^2-_aq through a unique solvent fluctuation-controlled gradual-flowing mechanism, and enlarging fluctuation can promote the conversion. Interestingly, we directly observe that e₂^2-_aq as the precursor can trigger hydrogen evolution via unique continuous spontaneous double proton transfer to the dielectron with a short-lived H^-_aq intermediate, but e^-_aq···e^-_aq does not directly. This is the first direct observation for the connection between e₂^2-_aq and spontaneous hydrogen evolution including participation of H^-_aq in aqueous solution, bridging relevant experimental phenomena. This work also evidences an unnoticed process, the double proton transfer mediated charge separation, and presents the first detailed analysis regarding the evolution dynamics of e₂^2-_aq for the understanding of the radiolysis reactions in aqueous solutions.

Partial Molar Volume of the Hydrated Electron

The partial molar volume of the hydrated electron was investigated with pulse radiolysis and transient absorption by measuring the pressure dependence of the equilibrium constant for e^-_aq + NH₄⁺ ⇔ H + NH₃. At 2 kbar pressure, the equilibrium constant decreases relative to 1 bar by only 6%. Using tabulated molar volumes for ammonia and ammonium, we have the result V̅(e^-_aq) - V̅(H) = 11.3 cm³/mol at 25 °C, confirming that V̅(e^-_aq) is positive and even larger than the hydrophobic H atom. Assuming on the basis of recent molecular dynamics simulations that the molar volume of the H atom is somewhat less than that of H₂, we estimate V̅(e^-_aq) = 26 ± 6 cm³/mol. The positive molar volume is consistent with an electron that exists largely in a small solvent void (cavity), ruling out a recent model ( Larsen , R. E. ; Glover , W. J. ; Schwartz , B. J. Science 2010 , 329 , 65 - 69 ) that suggests a noncavity structure with negative molar volume.

Low linear energy transfer radiolysis of supercritical water at 400 °C: in situ generation of ultrafast, transient, density-dependent "acid spikes"

There is growing interest in the radiation chemistry of supercritical water (SCW), as its use as a coolant in a nuclear reactor (Generation IV) is the logical evolution of the current (Generation III or less) water-cooled reactors. However, current knowledge about the potential effects of water radiolysis in a Gen-IV supercritical water-cooled reactor (SCWR) is incomplete. In this work, Monte Carlo track chemistry simulations of the low linear energy transfer (LET) radiolysis of SCW (H2O) at 400 °C are used in combination with a spherical "spur" model to study the effect of water density on the in situ radiolytic formation of H3O+ ions and the corresponding abrupt, transient, highly acidic pH response ("acid spikes") that is observed immediately after irradiation. The magnitude and duration of this acidic pH effect depend on the water density in the considered range of 0.15-0.6 g cm-3. It is strongest at times less than a few tens of picoseconds with the pH remaining nearly constant at ∼1.6 and 1.9 for the highest ("liquid-like") and lowest ("gas-like") density, respectively. At longer times, the pH gradually increases for all densities and finally reaches a constant value corresponding to the non-radiolytic, pre-irradiation concentration of H3O+, due to the autoprotolysis of water. Our results show that the lower the density of the water, the longer the time required to reach this constant value, ranging from ∼50 ns at 0.6 g cm-3 (pH ∼ 5.6) to ∼1 μs at 0.15 g cm-3 (pH ∼ 8.5). The generation of these highly acidic pH fluctuations around the "native" radiation tracks, though local and transient, raises questions about the potential implications of this effect in proposed Gen-IV SCW-cooled reactors regarding corrosion and degradation of materials.

Tuning the size of gold nanoparticles produced by multiple filamentation of femtosecond laser pulses in aqueous solutions

In the present study, we consider the self-regulated generation of spatially homogeneous low density plasma (LDP) micro-channels as a high intensity ionization source arising from the multi-filamentation of powerful femtosecond (fs) laser pulses in aqueous solutions. We investigate the modulation of the femtosecond laser multiple filamentation for tuning the size of gold nanoparticles (AuNPs) synthesized in an irradiated gold chloride solution. Previous studies on the radiation-induced synthesis of colloidal gold by more conventional ionization sources, such as high energy γ-rays and electron beams, highlighted the dependence of the size distribution of AuNPs on the density of energy deposited per unit of time, i.e. the dose rate. The present method of laser-induced production of AuNPs rests on a similar radiation-assisted process, i.e. the reduction of the solvated trivalent gold ions by the hydrated electrons produced upon ionization of water. We find that trivial optical manipulation varies the rate of deposited energy by laser irradiation, which can be considered equivalent to a variation of the dose rate. We investigate the influence of varying the density of energy deposited on the laser-induced gold cluster size distribution and made a comparison with the high energy radiation-induced synthesis of AuNPs. Here, our results highlight that the present method of laser irradiation, in the regime of LDP generation, mimics the radiolysis of water at an adjustable high dose rate. More generally, these spatially and temporally resolved plasmas could be developed as a tool for the unprecedented control of chemistry under ionizing radiation.

Prediction of rate constants of important chemical reactions in water radiation chemistry in sub and supercritical water – non-equilibrium reactions

The rate constants for reactions involved in the radiolysis of water under relevant thermodynamic conditions in supercritical water-cooled reactors are estimated for inputs in simulations of the radiation chemistry in Generation IV nuclear reactors. We have discussed the mechanism of each chemical reaction with a focus on non-equilibrium reactions. We found most of the reactions are activation controlled above the critical point and that the rate constants are not significantly pressure dependent below 300 °C. This work will aid industry with developing chemical control strategies to suppress the concentration of eroding species.

Radiolysis of Supercritical Water at 400°C: A Sensitivity Study of the Density Dependence of the Yield of Hydrated Electrons on the (eaq−+eaq−) Reaction Rate Constant

The temperature dependence of the rate constant (k) of the bimolecular reaction of two hydrated electrons (eaq−) measured in alkaline water exhibits an abrupt drop between 150°C and 200°C; above 250°C, it is too small to be measured reliably. Although this result is well established, the applicability of this sudden drop in k(eaq−+eaq−)) above ∼150°C to neutral or slightly acidic solution, as recommended by some authors, still remains uncertain. In fact, the recent work suggested that in near-neutral water the abrupt change in k above ∼150°C does not occur and that k should increase, rather than decrease, at temperatures greater than 150°C with roughly the same Arrhenius dependence of the data below 150°C. In view of this uncertainty of k, Monte Carlo simulations were used in this study to examine the sensitivity of the density dependence of the yield of eaq− in the low–linear energy transfer (LET) radiolysis of supercritical water (H2O) at 400°C on variations in the temperature dependence of k. Two different values of the eaq− self-reaction rate constant at 400°C were used: one was based on the temperature dependence of k above 150°C as measured in alkaline water (4.2×108  M−1 s−1), and the other was based on an Arrhenius extrapolation of the values below 150°C (2.5×1011  M−1 s−1). In both cases, the density dependences of our calculated eaq− yields at ∼60  ps and 1 ns were found to compare fairly well with the available picosecond pulse radiolysis experimental data (for D2O) for the entire water density range studied (∼0.15–0.6  g/cm3). Only a small effect of k on the variation of G(eaq−)) as a function of density at 60 ps and 1 ns could be observed. In conclusion, our present calculations did not allow us to unambiguously confirm (or deny) the applicability of the predicted sudden drop of k(eaq−+eaq−) at ∼150°C in near-neutral water.

Yields of H2 and hydrated electrons in low-LET radiolysis of water determined by Monte Carlo track chemistry simulations using phenol/N2O aqueous solutions up to 350 °C

The effect of temperature on the yields of H₂ and hydrated electrons in the low linear energy transfer radiolysis of water has been modeled by Monte Carlo track chemistry simulations using phenol/N₂O aqueous solutions from 25 up to 350 °C.

On the spur lifetime and its temperature dependence in the low linear energy transfer radiolysis of water

In the spirit of the radiation chemical "spur model", the lifetime of a spur (τ(s)) is an important indicator of overlapping spurs and the establishment of homogeneity in the distribution of reactive species created by the action of low linear energy transfer (LET) radiation (such as fast electrons or γ irradiation). In fact, τ(s) gives the time required for the changeover from nonhomogeneous spur kinetics to homogeneous kinetics in the bulk solution, thus defining the so-called primary (or "escape") radical and molecular yields of radiolysis, which are obviously basic to the quantitative understanding of any irradiated chemical system. In this work, τ(s) and its temperature dependence have been determined for the low-LET radiolysis of deaerated 0.4 M aqueous solutions of H(2)SO(4) and pure liquid water up to 350 °C using a simple model of energy deposition initially in spurs, followed by random diffusion of the species of the spur during track expansion until spur overlap is complete. Unlike our previous τ(s) calculations, based on irradiated Fricke dosimeter simulations, the current model is free from any effects due to the presence of oxygen or the use of scavengers. In acidic solutions, the spur lifetime values thus obtained are in very good agreement with our previous calculations (after making appropriate corrections, however, to account for the possibility of competition between oxygen and Fe(2+) ions for H˙ atoms in the Fricke dosimeter, an effect which was not included in our original simulations). In this way, we confirm the validity of our previous approach. As expected, in the case of pure, oxygen-free water, our calculated times required to reach complete spur overlap are essentially the same (within uncertainty limits) as those found in acidic solutions. This explicitly reflects the fact that the diffusion coefficients for the hydrated electron and the H˙ atom that are involved in the overall calculation of the lifetime of spurs in neutral or acidic media, respectively, are of similar magnitude over the 25-350 °C temperature range studied.

Dynamics and reactivity of trapped electrons on supported ice crystallites

The solvation dynamics and reactivity of localized excess electrons in aqueous environments have attracted great attention in many areas of physics, chemistry, and biology. This manifold attraction results from the importance of water as a solvent in nature as well as from the key role of low-energy electrons in many chemical reactions. One prominent example is the electron-induced dissociation of chlorofluorocarbons (CFCs). Low-energy electrons are also critical in the radiation chemistry that occurs in nuclear reactors. Excess electrons in an aqueous environment are localized and stabilized by the local rearrangement of the surrounding water dipoles. Such solvated or hydrated electrons are known to play an important role in systems such as biochemical reactions and atmospheric chemistry. Despite numerous studies over many years, little is known about the microscopic details of these electron-induced chemical processes, and interest in the fundamental processes involved in the reactivity of trapped electrons continues. In this Account, we present a surface science study of the dynamics and reactivity of such localized low-energy electrons at D(2)O crystallites that are supported by a Ru(001) single crystal metal surface. This approach enables us to investigate the generation and relaxation dynamics as well as dissociative electron attachment (DEA) reaction of excess electrons under well-defined conditions. They are generated by photoexcitation in the metal template and transferred to trapping sites at the vacuum interface of crystalline D(2)O islands. In these traps, the electrons are effectively decoupled from the electronic states of the metal template, leading to extraordinarily long excited state lifetimes on the order of minutes. Using these long-lived, low-energy electrons, we study the DEA to CFCl(3) that is coadsorbed at very low concentrations (∼10(12) cm(-2)). Using rate equations and direct measurement of the change of surface dipole moment, we estimated the electron surface density for DEA, yielding cross sections that are orders of magnitude higher than the electron density measured in the gas phase.

Radiolysis of supercritical water at 400 °C and liquid-like densities near 0.5 g/cm3 — A Monte Carlo calculation

Monte Carlo simulations are used to calculate the primary radical yields [Formula: see text], g(•OH), the sum [[Formula: see text] + g(•OH) + g(H•)], and the ratio g(H•)/[Formula: see text] in the low linear energy transfer (LET) radiolysis of supercritical water (SCW) at 400 °C in the high-density, liquid-like region near ∼0.5 g/cm3. Using all the currently available information on the reactivities and diffusion coefficients of the radiation-induced species under these conditions, and assuming the aqueous medium to be a “continuum”, a good accord is found between our calculations and the available experimental data. In particular, our computed [Formula: see text] yields at 60 ps and 1 ns compare very well with recently reported direct time-dependent [Formula: see text] yield measurements in SCW (D2O) at 400 °C and 0.570 g/cm3 using picosecond pulse radiolysis experiments.

Carbonate radical formation in radiolysis of sodium carbonate and bicarbonate solutions up to 250 degrees C and the mechanism of its second order decay

Pulse radiolysis experiments published several years ago (J. Phys. Chem. A, 2002, 106, 2430) raised the possibility that the carbonate radical formed from reaction of *OH radicals with either HCO(3)(-) or CO(3)(2-) might actually exist predominantly as a dimer form, for example, *(CO(3))(2)(3-). In this work we re-examine the data upon which this suggestion was based and find that the original data analysis is flawed. A major omission of the original analysis is the recombination reaction *OH + *CO(3)(-) --> HOOCO(2)(-). Upon reanalysis of the published data for sodium bicarbonate solutions and analysis of new transient absorption data we are able to establish the rate constant for this reaction up to 250 degrees C. The mechanism for the second-order self-recombination of the carbonate radical has never been convincingly demonstrated. From a combination of literature data and new transient absorption experiments in the 1-400 ms regime, we are able to show that the mechanism involves pre-equilibrium formation of a C(2)O(6)(2-) dimer, which dissociates to CO(2) and peroxymonocarbonate anion: *CO3(-)+*CO3(-)<-->C2O6(2-)-->CO2+O2COO(2-) *CO3(-) reacts with the product peroxymonocarbonate anion, producing a peroxymonocarbonate radical *O2COO(-), which can also recombine with the carbonate radical: *CO3(-)+CO4(2-)-->*CO4(-)+CO3(2-) *CO3(-)+CO4(-)-->C2O7(2-).