A radiochemical lab-on-a-chip paired with computer vision to unlock the crystallization kinetics of (Ba,Ra)SO4.

(Ra,Ba)SO4 solid solutions are commonly encountered as problematic scales in subsurface energy-related applications, e.g., geothermal systems, hydraulic fracturing, conventional oil and gas, etc. Despite its relevance, its crystallization kinetics were never determined because of radium (226), high radioactivity (3.7 × 1010 Bq g-1), and utilization in contemporary research, therefore constrained to trace amounts (< 10-8 M) with the composition of BaxRa1-xSO4 commonly restricted to x > 0.99. What if lab-on-a-chip technology could create new opportunities, enabling the study of highly radioactive radium beyond traces to access new information? In this work, we developed a lab-on-a-chip experiment paired with computer vision to evaluate the crystal growth rate of (Ba,Ra)SO4 solid solutions. The computer vision algorithm enhances experimental throughput, yielding robust statistical insights and further advancing the efficiency of such experiments. The 3D analysis results of the precipitated crystals using confocal Raman spectroscopy suggested that {210} faces grew twice as fast as {001} faces, mirroring a common observation reported for pure barite. The crystal growth rate of (Ba0.5Ra0.5)SO4 follows a second-order reaction with a kinetic constant equal to (1.23 ± 0.09) × 10-10 mol m-2 s-1.

Microfluidic investigation of pore-size dependency of barite nucleation.

The understanding and prediction of mineral precipitation processes in porous media are relevant for various energy-related subsurface applications. While it is well known that thermodynamic effects can inhibit crystallization in pores with sizes <0.1 µm, the retarded observation of mineral precipitation as function of pore size is less explored. Using barite as an example and based on a series of microfluidic experiments with well-defined pore sizes and shapes, we show that retardation of observation of barite crystallite can already start in pores of 1 µm size, with the probability of nucleation scaling with the pore volume. In general, it can be expected that mineralization occurs preferentially in larger pores in rock matrices, but other parameters such as the exchange of the fluids with respect to reaction time, as well as shape, roughness, and surface functional properties of the pores may affect the crystallization process which can reverse this trend.

Caesium and iodine release from spent mixed oxide fuels under repository relevant conditions: Initial leaching results.

Autoclave leaching experiments are conducted on three well-characterised, irradiated, and cladded mixed oxide fuel-rod segments with burnups ranging from 29 GWd/tHM to 52 GWd/tHM to investigate the instant release fraction of fission gases and long-lived fission products and to assess the long-term fuel matrix corrosion. The segments are exposed to bicarbonate solutions as reference groundwater at neutral pH and a synthetic young cementitious water at pH 13.5 under reducing atmosphere (4 vol% H2 in Ar at 40 bar pressure), since 2018. The initial leaching results for the fission products caesium and iodine as representative elements of the instant release fraction were found to depend on the leachate composition as well as on the fuel burnup.

Monazite-Type SmPO4 as Potential Nuclear Waste Form: Insights into Radiation Effects from Ion-Beam Irradiation and Atomistic Simulations.

Single-phase monazite-type ceramics are considered as potential host matrices for the conditioning of separated plutonium and minor actinides. Sm-orthophosphates were synthesised and their behaviour under irradiation was investigated with respect to their long-term performance in the repository environment. Sintered SmPO4 pellets and thin lamellae were irradiated with 1, 3.5, and 7 MeV Au ions, up to fluences of 5.1 × 1014 ions cm-2 to simulate ballistic effects of recoiling nuclei resulting from α-decay of incorporated actinides. Threshold displacement energies for monazite-type SmPO4 subsequently used in SRIM/TRIM simulations were derived from atomistic simulations. Raman spectra obtained from irradiated lamellae revealed vast amorphisation at the highest fluence used, although local annealing effects were observed. The broadened, but still discernible, band of the symmetrical stretching vibration in SmPO4 and the negligible increase in P-O bond lengths suggest that amorphisation of monazite is mainly due to a breaking of Ln-O bonds. PO4 groups show structural disorder in the local environment but seem to behave as tight units. Annealing effects observed during the irradiation experiment and the distinctively lower dose rates incurred in actinide bearing waste forms and potential α-radiation-induced annealing effects indicate that SmPO4-based waste forms have a high potential for withstanding amorphisation.

Thermodiffusion of ions in nanoconfined aqueous electrolytes.

Understanding of thermal effects on ion transport in porous media is very important for environmental applications. The movement of ions along a temperature gradient is named thermophoresis or thermodiffusion. In nanoporous media, where the interaction of ions with solid-liquid interfaces has a significant influence on their migration, the theoretical understanding of thermodiffusion is still incomplete. Herein, we present experimental results for the thermodiffusion of cations in saturated nanoporous silica by the through-diffusion method. Both the experimental data and theoretical analysis indicate that the temperature-induced polarization of surface charges strongly influences ionic transport. Stated simply, the electric field in a liquid electrolyte confined in nanopores changes when the applied temperature gradients are altered, thereby affecting the motion of the nanoconfined ionic species. By applying an external temperature field, the gradient of the surface charge density leads to the charged aqueous species exhibiting strong temperature gradient-dependent electrophoretic mobility. When the thickness of the electrical double layer is comparable to the size of the nanopores, the theory used herein indicates that this kind of nonisothermal ionic mobility is up to one order of magnitude larger than classical thermophoretic mobility. This study improves the understanding of the underlying mechanisms that govern the transport of ions in nanoporous media, which could set the stage for diffusional metamaterials induced by specific thermal fields.

Thermodynamic and Structural Modelling of Non-Stoichiometric Ln-Doped UO2 Solid Solutions, Ln = {La, Pr, Nd, Gd}.

Available data on the dependence of the equilibrium chemical potential of oxygen on degrees of doping, z, and non-stoichiometry, x, y, in U1-z Ln z O2+0.5(x-y) fluorite solid solutions and data on the dependence of the lattice parameter, a, on the same variables are combined within a unified structural-thermodynamic model. The thermodynamic model fits experimental isotherms of the oxygen potential under the assumptions of a non-ideal mixing of the endmembers, UO2, UO2.5, UO1.5, LnO1.5, and Ln 0.5U0.5O2, and of a significant reduction in the configurational entropy arising from short-range ordering (SRO) within cation-anion distributions. The structural model further investigates the SRO in terms of constraints on admissible values of cation coordination numbers and, building on these constraints, fits the lattice parameter as a function of z, y, and x. Linking together the thermodynamic and structural models allows predicting the lattice parameter as a function of z, T and the oxygen partial pressure. The model elucidates contrasting structural and thermodynamic changes due to the doping with LaO1.5, on the one hand, and with NdO1.5 and GdO1.5, on the other hand. An increased oxidation resistance in the case of Gd and Nd is attributed to strain effects caused by the lattice contraction due to the doping and to an increased thermodynamic cost of a further contraction required by the oxidation.

A lab-on-a-chip approach integrating in-situ characterization and reactive transport modelling diagnostics to unravel (Ba,Sr)SO4 oscillatory zoning.

The co-precipitation of sulphate minerals such as celestine and barite is widely studied because their formation is ubiquitous in natural and anthropogenic systems. Co-precipitation in porous media results in crystallization of solid solutions yielding characteristics such as oscillatory zoning that are rarely observed in bulk solution or in batch experiments. In the past, the precipitation of compositionally-zoned (Ba,Sr)SO4 crystals was observed post-mortem in macroscopic silica gel counter-diffusion experiments. Their formation was originally explained by the difference in the solubility products of the end-members combined with diffusion-limited transport of solutes to the mineral-fluid interface, while a later study favored the idea of kinetically controlled reactions. With recent advances combining in-operando microfluidic experiments and reactive transport modelling, it is now possible to verify hypotheses on the driving forces of transport-coupled geochemical processes. We developed a "lab on a chip" experiment that enabled the systematic study of the nucleation and growth of oscillatory-zoned (Ba,Sr)SO4 crystals in a microfluidic reactor. The compositions of the solid solutions were determined by in-situ Raman spectroscopy. Our investigation shows (1) that the composition of the nucleating phases can be approximated using classical nucleation theory, (2) that the oscillatory zoning is not solely controlled by the limited diffusional transport of solutes, and (3) that nucleation kinetics plays a major role in the switch between different stoichiometric compositions. The zoning phenomena is governed by the complex interplay between the diffusion of reactants and the crystallization kinetics as well as other factors, e.g. surface tension and lattice mismatch.

Insight Into Disorder, Stress and Strain of Radiation Damaged Pyrochlores: A Possible Mechanism for the Appearance of Defect Fluorite.

We have examined the irradiation response of a titanate and zirconate pyrochlore-both of which are well studied in the literature individually-in an attempt to define the appearance of defect fluorite in zirconate pyrochlores. To our knowledge this study is unique in that it attempts to discover the mechanism of formation by a comparison of the different systems exposed to the same conditions and then examined via a range of techniques that cover a wide length scale. The conditions of approximately 1 displacement per atom via He2+ ions were used to simulate long term waste storage conditions as outlined by previous results from Ewing in a large enough sample volume to allow for neutron diffraction, as not attempted previously. The titanate sample, used as a baseline comparison since it readily becomes amorphous under these conditions behaved as expected. In contrast, the zirconate sample accumulates tensile stress in the absence of detectable strain. We propose this is analogous to the lanthanide zirconate pyrochlores examined by Simeone et al. where they reported the appearance of defect fluorite diffraction patterns due to a reduction in grain size. Radiation damage and stress results in the grains breaking into even smaller crystallites, thus creating even smaller coherent diffraction domains. An (ErNd)2(ZrTi)2O7 pyrochlore was synthesized to examine which mechanism might dominate, amorphization or stress/strain build up. Although strain was detected in the pristine sample via Synchrotron X-ray diffraction it was not of sufficient quality to perform a full analysis on.

Chromium Doped UO2-Based Ceramics: Synthesis and Characterization of Model Materials for Modern Nuclear Fuels.

Cr-doped UO2 as a modern nuclear fuel type has been demonstrated to increase the in-reactor fuel performance compared to conventional nuclear fuels. Little is known about the long-term stability of spent Cr-doped UO2 nuclear fuels in a deep geological disposal facility. The investigation of suitable model materials in a step wise bottom-up approach can provide insights into the corrosion behavior of spent Cr-doped nuclear fuels. Here, we present new wet chemical approaches providing the basis for such model systems, namely co-precipitation and wet coating. Both were successfully tested and optimized, based on detailed analyses of all synthesis steps and parameters: Cr-doping method, thermal treatment, reduction of U3O8 to UO2, green body production, and pellet sintering. Both methods enable the production of suitable model systems with a similar microstructure and density as a reference sample from AREVA. In comparison with results from the classical powder route, similar trends upon grain size and lattice parameter were determined. The results of this investigation highlight the significance of subtly different synthesis routes on the properties of Cr-doped UO2 ceramics. They enable a reproducible tailor-made well-defined microstructure, a homogeneous doping, for example, with lanthanides or alpha sources, the introduction of metallic particles, and a dust-free preparation.

Dissolution Kinetics of International Simple Glass and Formation of Secondary Phases at Very High Surface Area to Solution Ratio in Young Cement Water.

Static dissolution experiments were carried out with the reference International Simple Glass under hyperalkaline pH at 70 °C and very high SA/V ratio. Three aspects of glass dissolution behavior were investigated, (1) the rate drop regime and the residual rate (stage II), (2) the formation of secondary phases including thermodynamic aspects, and (3) the microstructure of the interface of altered glass and secondary phases. A very low residual rate of 6 × 10-6 g/m2d was determined based on boron release, which was several orders of magnitude lower than the initial rate established between the start of the experiments and the first sampling on day 59. The presence of a porous layer with a thickness varying between 80 nm and 250 nm and a pore size between 10 nm and 50 nm was observed. CSH phases with a low Ca/Si ratio of 0.3-0.4 and zeolites were also visible at the surface of the altered glass grains, but no glass alteration resumption occurred, probably due to an important pH decrease already at day 59. Thermodynamic calculations assuming congruent glass dissolution and precipitation of the dissolved aqueous species confirmed the precipitation of CSH phases and zeolites.

Microfluidic flow-through reactor and 3D Raman imaging for in situ assessment of mineral reactivity in porous and fractured porous media.

An in-depth understanding of dissolution and precipitation of minerals in porous and fractured porous media and the complex feedback on the transport of fluids is essential for various subsurface applications. In this context, we developed a novel non-destructive "lab-on-chip" approach for quantitative in situ assessments of mineralogical changes in porous media. Our experimental approach involves a microfluidic flow-through reactor of reactive homogeneous and heterogeneous (fractured) porous media coupled with high-resolution imaging. Here, the reactive medium consists of compacted celestine grains seeded in a reservoir within the microfluidic chip. This medium reacts with a barium chloride solution injected into the microreactor at a constant flow rate, leading to the dissolution of celestine and growth of barite. Various seeding processes of the mineral grains allow the creation of homogeneous reactive porous media or the introduction of large heterogeneities such as fractures. Hence, our approach enables high-resolution investigations of reactive transport in fractured porous media. The use of confocal Raman spectroscopic techniques enables the spatio-temporal visualization of the mineral transformation at the pore-scale in two- and three-dimensions. Moreover, advanced pore-scale modelling correlates the hydrological heterogeneities to the geochemical observations in the micro-reactor, which explains the observed discrepancies between homogeneous and heterogeneous reactive media. Eventually, the proposed methodology can be applied to other chemical systems to provide new insights into hydro-geochemical coupling in porous and fractured porous media as well as high-fidelity datasets to benchmark reactive transport codes that are currently under development.

Rare-Earth Orthophosphates From Atomistic Simulations.

Lanthanide phosphates (LnPO 4) are considered as a potential nuclear waste form for immobilization of Pu and minor actinides (Np, Am, and Cm). In that respect, in the recent years we have applied advanced atomistic simulation methods to investigate various properties of these materials on the atomic scale. In particular, we computed several structural, thermochemical, thermodynamic and radiation damage related parameters. From a theoretical point of view, these materials turn out to be excellent systems for testing quantum mechanics-based computational methods for strongly correlated electronic systems. On the other hand, by conducting joint atomistic modeling and experimental research, we have been able to obtain enhanced understanding of the properties of lanthanide phosphates. Here we discuss joint initiatives directed at understanding the thermodynamically driven long-term performance of these materials, including long-term stability of solid solutions with actinides and studies of structural incorporation of f elements into these materials. In particular, we discuss the maximum load of Pu into the lanthanide-phosphate monazites. We also address the importance of our results for applications of lanthanide-phosphates beyond nuclear waste applications, in particular the monazite-xenotime systems in geothermometry. For this we have derived a state-of-the-art model of monazite-xenotime solubilities. Last but not least, we discuss the advantage of usage of atomistic simulations and the modern computational facilities for understanding of behavior of nuclear waste-related materials.

Transport of oxygen into granitic rocks: Role of physical and mineralogical heterogeneity.

The rock matrix of granites is expected to be an important buffer against the dispersion of contaminants, e.g. radionuclides, and against the ingress of oxygenated glacial meltwater. The influence of matrix heterogeneity on O2 diffusive transport is assessed here by means of numerical experiments based on a micro-Discrete Fracture Network (micro-DFN) representation of the diffusion-available pore space along with random realisations of idealized biotite grains, to simulate the heterogeneous nature of granitic rocks. A homogeneous-based analytical solution is also presented and used to assess possible deviations of the numerical experiments from the assumption of homogeneity. The analytical solution is also used to test upscaled values of mineral surface area. The numerical experiments show that the matrix behaves as a composite system, with the coexistence of fast and slow diffusive pathways. This behavior is more evident at low Damköhler numbers. Our interpretation of the numerical experiments points out the importance to properly characterise the heterogeneity of the rock matrix.

Comparison of Uranium(VI) and Thorium(IV) Silicates Synthesized via Mixed Fluxes Techniques.

Two uranium and two thorium silicates were obtained using high temperature mixed fluxes methods. K14(UO2)3Si10O30 crystallizes in the P21/ c space group and contains open-branched sechser (six) single silicate chains, whereas K2(UO2)Si2O6 crystallizes in the C2/ c space group and is built of unbranched achter (eight) silicate chains. The crystals of K14(UO2)3Si10O30 and K2(UO2)Si2O6 are related by increasing U/Si molar ratios, and both structures contain the same secondary building units (SBUs), [USi6] heptamers. The triangle diagram for all known A+-UO22+-SiO44- phases demonstrates the high polymerization level of silicate groups in the system, which was compared with the family of A+-UO22+-BO33-/BO45- compounds. For both thorium silicates, the transformation of K2ThSi2O7 to K2ThSi3O9 was found to be a factor of the reaction time. K2ThSi2O7 crystallizes in the C2/ c space group and belongs to the Na2SiVISi2O7 structure type. Its 3D framework consists of diorthosilicate Si2O7 group and ThO6 octahedra. Noncentrosymmetric K2ThSi3O9 crystallizes in the hexagonal P63 space group and adopts mineral wadeite-type structure based upon triorthosilicate Si3O9 rings and ThO6 octahedra. The coordination environment of thorium for all existing oxo-anion compounds including B, Si/Ge, P/As, Cr/Mo/W, and S/Se/Te are summarized and analyzed. Additionally, spectroscopic properties of all novel materials have been studied.

Hydrothermal Synthesis, Study, and Classification of Microporous Uranium Silicates and Germanates.

Four novel uranyl silicates and germanates with framework structures, K4Na2(UO2)3(Si2O7)2·3H2O, K4Na2(UO2)3(Ge2O7)2·3H2O, H3O(UO2)2(HGe2O7)·2H2O, and Na2(UO2)GeO4, have been synthesized by means of the hydrothermal method. The structures of the title compounds were refined by single-crystal X-ray diffraction and characterized by Raman spectroscopy. We used the method of secondary building units (SBUs) for a crystal chemical analysis of the 3D framework and their topologies. The framework of the K4Na2(UO2)3(T2O7)2·3H2O (T = Si, Ge) series exhibits large 14-membered rings and smaller 8-membered rings which are built upon [UT4] pentamers. The internal size of the largest pores is approximately 12.39 × 3.33 Å2. H3O(UO2)2(HGe2O7)·2H2O is based on 10-membered rings with intermediate sized pores. They are built upon [U2Ge2] tetramers with 7-fold-coordinated U. The internal dimension of the pores in H3O(UO2)2(HGe2O7)·2H2O is smaller compared to the K4Na2(UO2)3(T2O7)2·3H2O (T = Si, Ge) series with ∼5.91 × 5.33 Å2. Its topology is similar to several uranium germanate synthetic phases and silicate minerals, especially α- and ß-uranophane which are constructed from similar building units. A novel 3D framework type of Na2(UO2)GeO4 with 8-membered rings demonstrates the smallest free volume in the family of porous uranium germanates. It crystallizes in tetragonal symmetry and is built upon corner sharing of [UGe4] pentamers. The size of the channels is ∼6.76 × 4.27 Å2. The vibrational bands in Raman spectra were associated with pyro-(Si2O7)6- and -(Ge2O7)6- groups, with the Ge-OH bond and with H3O+ cations, confirming the results of the X-ray crystallographic structural characterization. We systemized existing uranyl silicates and germanates based on their building units and chemical composition. We found a simple structural dependence between synthetic conditions and chemical composition.

Unexpected Behavior of Np in Oxo-selenate/Oxo-selenite Systems.

A study of neptunium (Np) chemistry in the complex oxo-selenium system has been performed. Hereby, two sets of precipitation experiments were conducted, investigating the influence of the initial oxidation state of selenium using SeIVO2 and H2SeVIO4 with NpV in alkali nitrate solution, keeping the ratio of Np/Se constant. Surprising results were observed. Five novel neptunium and selenium bearing compounds have been obtained by slow evaporation from aqueous solution. The novel NpIV phase K4-x[Np(SeO3)4-x(HSeO3)x]·(H2O)1.5 (1) crystallizes in green-colored, plate-shaped crystals and was obtained by adding SeO2 and ANO3 to a NpV stock solution. Single-crystal X-ray diffraction reveals one-dimensional chain structures composed of square antiprismatic NpO8 polyhedra linked via four trigonal pyramidal SeO3 and HSeO3 units. Raman spectral analysis supports the presence of both selenite and hydroselenite due to the presence of corresponding modes within the spectra. The addition of selenic acid to a NpV stock solution resulted in the precipitation of elongated rose prisms of K2[(NpO2)2(SeO4)3(H2O)2]·(H2O)1.5 (2), Rb2[(NpO2)2(SeO4)3(H2O)2]·(H2O)2 (3) and K9[(NpO2)9(SeO4)13.5(H2O)6]·(H2O)12 (4) as well as light red plates of Cs2[(NpO2)2(SeO4)3] (5). To our knowledge, this is the first report of NpVI selenates. All four structures show two-dimensional layered structures with alkali cations acting as charge balancing counter cations. Hereby the layers of compounds 2 and 3 are found to be orientational geometric isomers. Distinctly different phenomena are made responsible for the phase formation within these systems. The kinetically driven process of NpV disproportionation led to the formation of the NpIV selenites in the SeIV-based system, whereas the oxidation of NpV by reduction of nitrate in acidic conditions is responsible for the formation of the NpVI selenates in the SeVI system. The influence of air oxygen is also discussed for the latter reaction.

Microtomography-based Inter-Granular Network for the simulation of radionuclide diffusion and sorption in a granitic rock.

Field investigation studies, conducted in the context of safety analyses of deep geological repositories for nuclear waste, have pointed out that in fractured crystalline rocks sorbing radionuclides can diffuse surprisingly long distances deep into the intact rock matrix; i.e. much longer distances than those predicted by reactive transport models based on a homogeneous description of the properties of the rock matrix. Here, we focus on cesium diffusion and use detailed micro characterisation data, based on micro computed tomography, along with a grain-scale Inter-Granular Network model, to offer a plausible explanation for the anomalously long cesium penetration profiles observed in these in-situ experiments. The sparse distribution of chemically reactive grains (i.e. grains belonging to sorbing mineral phases) is shown to have a strong control on the diffusive patterns of sorbing radionuclides. The computed penetration profiles of cesium agree well with an analytical model based on two parallel diffusive pathways. This agreement, along with visual inspection of the spatial distribution of cesium concentration, indicates that for sorbing radionuclides the medium indeed behaves as a composite system, with most of the mass being retained close to the injection boundary and a non-negligible part diffusing faster along preferential diffusive pathways.

Porous Uranyl Borophosphates with Unique Three-Dimensional Open-Framework Structures.

Two novel alkali-metal uranyl borophosphates have been prepared and characterized for the first time, namely, K5(UO2)2[B2P3O12(OH)]2(OH)(H2O)2 and K2(UO2)12[B(H2PO4)4](PO4)8(OH)(H2O)6 denoted as KUPB1 and KUPB2, respectively. KUPB1 was obtained hydrothermally at 220 °C and crystallizes in a monoclinic structure in the chiral space group P21. The unit cell parameters of KUPB1 are a = 6.7623(2) Å, b = 19.5584(7) Å, c = 11.0110(4) Å, α = γ = 90°, ß = 95.579(3)°, and V = 1449.42(8) Å3. It features a unique three-dimensional (3D) open-framework structure, composed of two corner-sharing linked one-dimensional (1D) anionic borophosphates (BP), [B2P3O13]5-, along the a axis and uranyl phosphate (UP), [(UO2)(PO4)3]7-, chains along the c axis, further bridged by PO4 tetrahedra. Multi-intersectional channels can be observed within the structure, in which the largest 11-ring (11-R) tunnel size is ∼7.0 Å × 8.8 Å. Its simplified framework can be described as a new 4-nodal net topological type with a point symbol of {4.84.10}{42.6}2{43.62.83.102}{82.10}. By modification of the synthetic conditions of KUPB1 through an increase in the amount of H3BO3 as flux 4-fold and a reduction of water as the reaction medium, the novel compound KUPB2 is generated. The unit cell parameters of KUPB2 are a = b = 21.8747(3) Å, c = 7.0652(2) Å, α = ß = γ = 90°, and V = 3380.72(12) Å3. KUPB2 crystallizes in a tetragonal structure in the polar space group I4̅2m, and its structure is based on a highly complex 3D framework, {(UO2)12[B(PO4)4](PO4)8}9-, in which 1D 8-R UP [(UO2)(PO4)]- tubes can be observed along the c axis. The [(UO2)(PO4)]- tubes consist of three uranyl chains along the c axis, which are linked alternately by [PO4]3- tetrahedra. Those isolated 1D [(UO2)(PO4)]- tubes are further bridged through [(UO2)4B(PO4)4]- clusters, forming an exceptional 3D open-framework structure. Its simplified cation network is a new 5-nodal net topological type such as {32.43.5.62.7.8}8{34.45.54.62}8{4.62.83}4{42.6}4{44.62}. Their facile hydrothermal synthetic routes, porous structure topology, thermal stability, and Raman spectroscopy properties are reported and discussed.

Continuum-based DFN-consistent numerical framework for the simulation of oxygen infiltration into fractured crystalline rocks.

We present an enhanced continuum-based approach for the modelling of groundwater flow coupled with reactive transport in crystalline fractured rocks. In the proposed formulation, flow, transport and geochemical parameters are represented onto a numerical grid using Discrete Fracture Network (DFN) derived parameters. The geochemical reactions are further constrained by field observations of mineral distribution. To illustrate how the approach can be used to include physical and geochemical complexities into reactive transport calculations, we have analysed the potential ingress of oxygenated glacial-meltwater in a heterogeneous fractured rock using the Forsmark site (Sweden) as an example. The results of high-performance reactive transport calculations show that, after a quick oxygen penetration, steady state conditions are attained where abiotic reactions (i.e. the dissolution of chlorite and the homogeneous oxidation of aqueous iron(II) ions) counterbalance advective oxygen fluxes. The results show that most of the chlorite becomes depleted in the highly conductive deformation zones where higher mineral surface areas are available for reactions.

Thorium Chemistry in Oxo-Tellurium System under Extreme Conditions.

Through the use of a high-temperature/high-pressure synthesis method, four thorium oxo-tellurium compounds with different tellurium valence states were isolated. The novel inorganic phases illustrate the intrinsic complexity of the actinide tellurium chemistry under extreme conditions of pressure and temperature. Th2Te3O11 is the first instance of a mixed-valent oxo-tellurium compound, and at the same time, Te exhibits three different coordination environments (TeIVO3, TeIVO4, and TeVIO6) within a single structure. These three types of Te polyhedra are further fused together, resulting in a [Te3O11]8- fragment. Na4Th2(TeVI3O15) and K2Th(TeVIO4)3 are the first alkaline thorium tellurates described in the literature. Both compounds are constructed from ThO9 tricapped trigonal prisms and TeVIO6 octahedra. Na4Th2(TeVI3O15) is a three-dimensional framework based on Th2O15 and Te2O10 dimers, while K2Th(TeVIO4)3 contains tungsten oxide bronze like Te layers linked by ThO9 polyhedra. The structure of ß-Th(TeIVO3)(SO4) is built from infinite thorium chains cross-linked by TeIVO32- and SO42- anions. Close structural analysis suggests that ß-Th(TeIVO3)(SO4) is highly related to the structure of α-Th(SeO4)2. Additionally, the Raman spectra are recorded and the characteristic peaks are assigned based on a comparison of reported tellurites or tellurates.