{'fraction_non_alphanumeric': 0.03982017809285035, 'fraction_numerical': 0.02709432004841359, 'mean_word_length': 4.730875941339675, 'pattern_counts': {'":': 0, '<': 2, '<?xml version=': 0, '>': 8, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Uranium (VI) Adsorbate Structures on Portlandite [Ca(OH)2] Type Surfaces Determined by Computational Modelling and X-ray Absorption Spectroscopy

Portlandite [Ca(OH)2] is a potentially dominant solid phase in the high pH fluids expected within the cementitious engineered barriers of Geological Disposal Facilities (GDF). This study combined X-ray Absorption Spectroscopy with computational modelling in order to provide atomicscale data which improves our understanding of how a critically important radionuclide (U) will be adsorbed onto this phase under conditions relevant to a GDF environment. Such data are fundamental for predicting radionuclide mass transfer. Surface coordination chemistry and speciation of uranium with portlandite [Ca(OH)2] under alkaline groundwater conditions (ca. pH 12) were determined by both in situ and ex situ grazing incidence extended X-ray absorption fine structure analysis (EXAFS) and by computational modelling at the atomic level. Free energies of sorption of aqueous uranyl hydroxides, [UO2(OH)n] (n = 0–5) with the (001), (100) and (203) or (101) surfaces of portlandite are predicted from the potential of mean force using classical molecular umbrella sampling simulation methods and the structural interactions are further explored using fully periodic density functional theory computations. Although uranyl is predicted to only weakly adsorb to the (001) and (100) clean surfaces, there should be significantly stronger interactions with the (203/101) surface or at hydroxyl vacancies, both prevalent under groundwater conditions. The uranyl surface complex is typically found to include four equatorially coordinated hydroxyl ligands, forming an inner-sphere sorbate by direct interaction of a uranyl oxygen with surface calcium ions in both the (001) and (203/101) cases. In contrast, on the (100) surface, uranyl is sorbed with its axis more parallel to the surface plane. The EXAFS data are largely consistent with a surface structural layer or film similar to calcium uranate, but also show distinct uranyl characteristics, with the uranyl ion exhibiting the classic dioxygenyl oxygens at 1.8 Å and between four and five equatorial oxygen atoms at distances between 2.28 and 2.35 Å from the central U absorber. These experimental data are wholly consistent with the adsorbate configuration predicted by the computational models. These findings suggest that, under the strongly alkaline conditions of a cementitious backfill engineered barrier, there would be significant uptake of uranyl by portlandite to inhibit the mobility of U(VI) from the near field of a geological disposal facility.


Introduction
The interaction of radionuclides with cementitious engineered barrier systems and the surrounding mineral environment is of fundamental importance for the safe disposal of intermediate-level radioactive wastes [1]. Typically, the most ubiquitous radionuclide by mass in intermediate-level wasteforms to be disposed of will be uranium. In its higher oxidation state, U(VI) can be readily solubilised in aqueous conditions and complexes of uranyl ions (UO 2 2+ ) have the potential to be highly mobile [2]. Additionally, under the high pH and mildly reduced conditions associated with cementitious geological disposal, U(VI) is reportedly stabilised [3].
This study will focus on the higher pH environments anticipated over the longer term in intermediate-level waste geological disposal facilities (GDF) which, in many cases, will rely upon cement-based backfills and wasteforms to prevent or limit the migration of the radionuclides to the wider environment. Typical cement degradation conditions are highly alkaline, ranging from pH > 10 in CSH gels through to pH 13 when KOH/NaOH dominate dissolution. Portlandite [Ca(OH) 2 ] is a potentially significant mineral surface in these materials during the period that intermediate-level wasteforms (ILW) will be emplaced, and there is a paucity of quantitative understanding of how the waste species will interact or be retarded by this mineral phase. The extreme pH conditions expected in the most proximate conditions in a GDF are the most difficult to model or access experimentally, so these conditions are the focus of this study. Thus, a detailed and critical understanding of the uptake of uranyl (VI) ions by portlandite is important in any safety case for a cementitious GDF. The objective of this study is to determine adsorbate structures under extremely high pH conditions for uranium on reactive hydroxylated portlandite-type surfaces.
In order to better understand the nature of the interactions of uranium with portlandite, we have employed both experimental grazing incidence X-ray absorption spectroscopy and computational modelling at the atomic level. Our computational models have been constructed to assess the interactions that might be expected under aqueous conditions with pH circa 12 (the point of zero charge pH~12.7 of portlandite [4]) and where, in the absence of high concentrations of carbonate, the uranyl hydroxides [UO 2 (OH) n ] 2-n (n = 0-5) will be the dominant species.
We shall show that, although we might expect the kinetics and surface chemistry of portlandite to be relatively simple, there are a wide range of specific U(VI) interactions to consider. Specifically, precipitation, recrystallization or incorporation of U(VI) near the surface are likely since sorption of uranyl ions is particularly favourable on roughened surfaces exposed to an aqueous environment.
Galmarini et al. [5] have shown that whilst the (001) surface of portlandite is the most stable in the gas phase, the (101) surface [6] (originally identified as the similar (203) surface), is the most stable in aqueous solutions. This latter surface is of particular interest here, as it implies that surface calcium ions are exposed to solvent, rather than a dense hydroxyl layer, and these may act as possible enhanced binding sites for uranyl hydroxide anionic species.
Potentials of mean force, which allow us to determine the free energy of adsorption of the uranyl species under aqueous conditions, have been computed for the three main surfaces of portlandite: (001), (100) and (101). In addition, adsorption will be considered in the region of hydroxyl vacancy defects of these solvated surfaces, including a step on the (101) surface (i.e., 203 plane), to better understand potentially less ideal conditions. Density functional theory is then applied to some configurations for a more accurate assessment of the inner-sphere interactions of uranyl with the surface-in particular, whether uranyl has a tendency to be retained perpendicular to the surface (interacting via the axial oxygen) or in a more parallel fashion (with the surface substituting for equatorial ligands).
Interactions of uranyl complexes with mineral surfaces: Uranyl has a propensity for uptake to a variety of common mineral surfaces, such as the carbonates including calcite or aragonite [7,8], silicates [9][10][11] and hydroxides such as brucite [12,13]. However, surprisingly, there has been relatively little research on the interactions of U(VI) with portlandite [Ca(OH) 2 ], with previous studies focussing on modelling [14][15][16][17] or determining its structural properties [18][19][20] and crystal morphology [5,6]. This is in large part due to the experimental difficulties associated with the high solubility and rapid dissolution kinetics of both CaO and Ca(OH) 2 . Despite these experimental challenges, information about surface-specific active sites is important in order to make predictions about how uranyl adsorbs and may be incorporated into the crystal structure. Previous experimental studies have examined the uptake of uranyl in batch experiments [21] and have also highlighted the relatively high uptake of U(VI) in degraded portlandite phases [22].
To that end, ex-situ and in-situ grazing incidence EXAFS (reflEXAFS) has been used in this study to constrain the adsorption of uranyl onto single crystal (hydroxylated) CaO(111), which, upon hydroxylation, will serve as a structural and chemical equivalent to the portlandite (001) surface. The advantage of this approach is that the data collected represents uranyl adsorption on one orientation so that the local coordination on a specific plane can be precisely constrained. Another advantage is that adsorption of U(VI) can be monitored in real time, which can give information about adsorption kinetics that otherwise would be lost in batch experiments.
Computational modelling can be a complementary aid to experimental approaches to model U(VI) uptake, and, in particular, atomistic molecular dynamics simulations and density functional theory quantum chemistry computations can provide direct atomic detail relating to the specific processes that take place. It is important to recognise that the properties of natural systems in equilibrium are those of thermodynamic ensembles, which are strongly dependent on a variety of possible atomic configurations over an extended timescale. Computationally, this will necessitate a range of computational models to be considered. The ligand exchange and sorption mechanisms are relatively rare events and are unlikely to occur within the pico-to nano-second timescales employed in molecular dynamics simulations [23]. The potential of mean force approach utilised in this study employs the constraint of a "reaction or collective coordinate" to mitigate some of these timescale limitations; however, although it does allow for some statistical sampling, we will only employ one reaction constraint in each simulation and a range of possible speciation models will need to be considered to describe the broader picture. We shall follow the general classical simulation approach utilised in previous sorption studies of uranyl on calcite [24,25], augmented by more accurate first-principles density functional theory calculations of a selection of inner-sphere complex models to probe potential surface-adsorbed states of uranyl. Essentially, the MD and PMF computations provide constraints on the free energy change associated with bringing the uranyl ion onto the mineral surface, with a range of values computed as a function of distance. This allows for a prediction of which surface statistically provides the lowest energy adsorbate layer. The DFT models then provide detailed information about the lowest energy chemical coordination structure for uranyl on each mineral surface of interest. Combining these two methods gives us complementary information about uranyl surface-attachment onto portlandite. Spectroscopic measurements then provide a test of the computational results and will allow us to confirm whether the computationally predicted configurations are discernible on real surfaces. Although the sorption of uranyl on portlandite has not been previously modelled, uranyl species have been shown to adsorb directly onto a range of related mineral surfaces, with a range of computational studies, including sorption of uranyl to calcite [24,26], montmorillonite [27], quartz [28], feldspar [29] and aluminosilicates [30]. These previous studies will be used to help constrain our experimental and computational observations.

Atomistic Models
All computations are initially based upon the crystal structure of bulk portlandite determined by neutron diffraction [20], which has a rhombohedral structure with lattice parameters a = b = 3.589 Å, c = 4.911 Å, α = β = 90 • and γ = 120 • . In this work, classical simulations of the bulk structure or surfaces of portlandite used the force field of Freeman et al. [31] and previously used by Galminari [5] to study its crystal morphology. Water was modelled using the four site TIP4P/2005 model [32] and uranyl was modelled using the parameters of Greathouse et al. [28], modified to include interactions with TIP4P/2005 water and portlandite. These force-field parameters are given in the supporting information.
The (001), (100) and (203) surfaces of portlandite were cleaved from the bulk unit cell using Accelrys Materials Studio [33] and a 50 Å vacuum gap was created above each surface to remove interactions between periodic images of the slab in the z-direction. A structure of 8 × 8 supercells of these surfaces was then generated, each eight Ca(OH) 2 layers thick (2560 atoms total), and the uranyl species (UO 2 (OH) n 2-n ) n = 0-5 were placed into the vacuum gap. The cells were then solvated with TIP4P/2005 water and the whole system was equilibrated using classical molecular dynamics simulations at 298 K.
We note that the (203) surface in the simulation cell differs from a (101) surface by a single step; however, in all the non-defect surface computations, adsorption was not sampled in the vicinity of the step and all interactions were analogous to simulation on a pure (101) surface. Hence, we shall refer to computations involving this surface as (203/101) rather than (203). (See Figure S1 for a comparison of the (101) and (203) surfaces.) Previous calcium oxide sorption experiments have shown that, particularly at high pH, the mineral surface will be rough and quite dynamic, with dissolution and re-crystallisation occurring at relatively short timescales [21]. In order to assess how uranyl sorption may be affected in this situation, a series of modified surfaces were modelled in which a point defect (vacancy) was created in each of the surfaces, (001) and (100), by removing a single hydroxide ion from the exposed face. In the case of the (203/101) surface, the vacancy was introduced into the step that occurs explicitly in the unit cell of the (203) surface and sorption was explicitly sampled around this site. We shall refer to the 3 defect surfaces as (001-OH), (100-OH) and (203-OH), respectively. Although we shall refer to the hydroxyl vacancies as defects, it should be borne in mind that the surfaces in our models are fully hydrated and any vacancy can be bound explicitly by water during equilibration, analogous to lower pH conditions where protonation of some of the hydroxyl surface is expected. Similarly, the (203/101) surfaces will be hydrated and water will solvate exposed calcium ions. Thus, an important differentiating feature of the (001), (100) and (203/101) model surfaces will be the density of hydroxyl groups: these are 0.0896 Å −2 , 0.0567 Å −2 and 0.1038 Å −2 , respectively.

Classical MD and PMF Simulations
Potential of mean force (PMF) calculations have been performed to predict the work (free energy) required to bring the 6 uranyl species, [UO 2 (OH) n ] 2-n (n = 0-5), from being fully solvated in solution towards and into contact with the three (001), (100) and (203/101) surfaces. For the defect surfaces, only the UO 2 (OH) 4 2 complex was explicitly considered, although no specific restraints were included to maintain this coordination environment.
All classical molecular dynamics (MD) simulations were performed using DL_POLY Classic [34] with the PLUMED [35] plugin used to include the umbrella sampling potential. A timestep of 0.5 fs was used in all simulations. The simulations were all performed at a constant temperature of 298 K within the NVT ensemble using the Nose-Hoover [36] thermostat with a 0.1 ps relaxation time. Each simulation consisted of 100 ps for equilibration and 1 ns for production, during which statistics were collected.
A series of umbrella sampling molecular dynamics simulations were performed to compute a potential of mean force (PMF) to quantify the free energy of taking each uranyl complex from an aqueous solution to a surface. Our sampling protocols follow previous work using these techniques [24,37], where a harmonic potential of the form: was used to constrain the collective variable or reaction coordinate. For each PMF calculation, approximately 16 simulation windows were used, covering a coordinate that brings the uranyl complex from the solution to the surface of the slab, with additional windows added to improve the sampling in regions of high force. The coordinate was defined as the distance between the uranium atom and the centre of mass of the calcium atoms in the slab. The surface of the slab is defined as the z-coordinate of the hydrogen atoms at the exposed face of the slab in gas phase. The number of windows, the interval between each one and the value of the constraint parameter k varied such as to include a complete sampling of the reaction coordinate, with a value of k ranging from 0.5 to 3.0 eV. The weighted histogram analysis method [38,39] (WHAM) was used to combine and analyse the window distributions obtained from each simulation window and to calculate the free energy profile. Additional PMF simulations were performed to assess the hydroxyl-water ligand exchange energy (see Supplementary Materials for further details).
Geometry optimisations were performed for models of the UO 2 (OH) 4 2− complex on clean (001, 100 and 203/101) and defect (001, 100 and 203) portlandite surfaces, with the uranyl axis either parallel or perpendicular to the surface (based on observed orientations from the classical trajectories). In each of these calculations, the surface slabs contained 8 Ca(OH) 2 layers that were analogous to the classical simulations, the absorbed complex was surrounded by thirty-eight explicit water molecules and periodic boundary conditions were applied.
In Table 1, we summarise a selection of gas-phase PBE functional structural optimisations with different functionals of UO 2 (OH) 4 2− in a typical conformation, very close in energy to the global minimum, observed in the simulations with 3 of the hydroxide H-atoms oriented in the same direction. EXAFS data from a previous study are included for comparison. These DFT calculations were performed using Gaussian 09 [46] with the pure PBE, hybrid PBE0, B3LYP and M06 functionals [47], a large-core SDD basis set and pseudopotential for uranium [48] and the 6-31G(d,p) basis set for oxygen, hydrogen and subsequently for Ca. We note, for later comparison to our EXAFS uranyl-portlandite distributions, that the PBE functional is likely to overestimate the U-O ax bond lengths and U-O eq bonds by around 0.03 and 0.02 Å, respectively, in agreement with previous work [49]. The r(U-O ax ) and r(U-O eq ) bond lengths of UO 2 (OH) 4 2− with PBE in the gas phase, or when water is implicitly included via the polarizable continuum model (PBE/PCM), are respectively 1.86 Å and 1.85 Å, or slightly longer than was observed experimentally in solution at 1.83 Å. We also note that hybrid functionals do result in shorter axial lengths and the equatorial bond lengths are slightly longer; unfortunately, the additional computational demands of these functionals precluded their use with the surface slab together with explicit water.
Two further gas-phase complexes were also considered (see Table 1) in order to give an estimate of the Ca-U and U-U interaction distances when Ca 2+ directly interacts with UO 2 (OH) 4 2− or if UO 2 (OH) 4 dimerises to (UO 2 ) 2 (OH) 6 2− , in which there are 4 equatorial OH and 2 bridging OH groups. Dimerisation results in the U-O ax and U-O eq distances being shortened, with respect to the monomer, by 0.03 Å, and a U-U distance of 3.98 Å was predicted. The effect of the Ca 2+ counter ion, binding with a Ca-U distance of 3.14 Å, is to lengthen the U-O ax to which it directly interacts, and to shorten the other; the U-O eq bonds to which the Ca 2+ interacts are both considerably lengthened to about 2.50 Å and the remaining bonds slightly shortened.

Experimental
Surface coordination chemistry and speciation of uranium were determined by both in-situ and ex-situ reflection geometry extended X-ray absorption fine structure analysis (reflEXAFS) on oriented hydroxylated single crystals of CaO (111), equivalent to the Ca(OH) 2 (001) surface. In order to obtain a surface similar to that of portlandite (001), commercially bought oriented single crystals of CaO (111) were used in this study and then hydroxylated. Oriented single crystals (10 × 10 × 1 mm) were produced, cut and polished by SurfaceNet GmbH Germany. The surface root mean square (RMS) starting roughness was between 20 and 30 Å.

Ex Situ Adsorption Experiments
A stock solution of 25,000 ppm natural UO 2 Cl 2 was prepared. Two sets of experiments were completed (i) with 2 mM CO 3 2− and (ii) without CO 3 2− . In the first experiment, 1 mL (25 ppm) U(VI) of stock was spiked in a 2 mM K 2 CO 3 solution. The pH was adjusted dropwise using HCl to pH 10.6 (±0.1) to avoid the possible precipitation of compreignacite (K 2 (UO 2 ) 6 O 4 (OH) 6 ·7H 2 O). The dominant uranium species in this solution would be UO 2 (CO 3 ) 3 4− (>99%, PHREEQC, SIT database [51]). The U(VI) solution was transferred into a centrifuge tube along with the crystal. The total contact time of the crystal solution was 48 h (±2 h), after which the crystal was removed using Teflon-coated tweezers and wicked dry (without touching the crystal surface) to be transferred onto a domed sample holder suitable for reflEXAFS. The sample holder was then sealed using Kapton tape and bagged.
The second experiment involved spiking 0.4 mL (10 ppm) of U(VI) stock solution into a decarbonized pH 12 (0.01 M KOH) solution. Decarbonized water was obtained by boiling MilliQ water for one hour while purging it with N 2 gas. The dominant solute species were calculated to be~51% UO 2 (OH) 3 − and~49% UO 2 (OH) 4 2− . During the reaction, the solution remained clear and no precipitate was formed. The total contact time of the crystal solution was 48 h (±2 h), and it was then prepared as for experiment (i).

In Situ Adsorption Experiments
In line with the ex situ batch experiments, parallel solution chemistry was used here. In addition, an extra experiment with 10 ppm U(VI) at pH 12 in 0.01 M LiOH was used to further enhance the solubility of U(VI) at pH 12. These results were compared with the other experiments in order to account for any possible effects of compreignacite precipitation.
These experiments were performed at beam line I18 at Diamond Light Source Ltd., Didcot, UK. Real-time adsorption measurements were performed using a micro-fluid cell with a domed X-ray transparent window designed for actinide-bearing solutions. The sample cell was purged prior to reaction with N 2 gas. A syringe pump was connected to the cell using Teflon tubing and a reactant solution was input at a constant fluid flow rate of 1.5 mL/h for two hours. Uranium fluorescence-intensity data were collected during reaction in order to monitor surface uptake. After reaction, the surface was rinsed with 0.4 mL of CO 2 free MilliQ water at a flow rate of 1.5 mL/h and then purged with N 2 gas. This removed solute species of uranyl from the cell and left behind a thin film containing adsorbed uranyl.

reflEXAFS Measurements
Beam line I18 was configured to perform reflEXAFS. The beam line uses a Si (111) monochromator to select the energy of the beam, which was calibrated against the first derivative of yttrium (Y), defined at 17,038 eV. Fluorescence-yield data were collected using a four element Si drift Vortex detector, which was positioned vertically and directly above the crystal surface. Upstream vertical slits of 50 µm were used to define the incident beam.
Crystal alignment was achieved by first completing a rocking curve measurement so that the crystal surface was positioned flat and cutting half of the incident beam intensity. The angle of the crystal surface was then oriented to an incident angle of~30 millidegrees, but in all cases kept below the critical angle of total external reflection so that the fluorescence signal remained optimized.
The EXAFS, χ(k), were refined from the raw data by standard procedures, including background subtraction, data normalization and fitting. Data analysis was performed using the software packages (d)Athena and (d)Artemis 0.9.18 [52].

Bulk, Surface, and Uranium-Incorporated Portlandite
In order to assess the variance between the classical force-field and quantum DFT models, the bulk structure of portlandite was optimised using the PBE functional (Table S1). The force-field approach predicts lattice parameters of a = b = 3.662 Å and c = 4.740 Å, while DFT predicts a = b = 3.644 Å and c = 4.993 Å, both within a few percent of the literature values [19]. Surface energies of the (001), (100) and (203) surfaces (gas phase, i.e., with a vacuum gap) of portlandite were calculated using both the classical force-field and DFT (PBE) methods. The surface energy is calculated using where E slab is the energy of the slab, E bulk is the energy of the bulk unit cell, n is the number of bulk units in the slab and S is the surface area of the slab defined as AB sin γ, where A and B are the lengths of the A and B lattice vectors and γ is the angle between them. The surface energies (Tables S1-S3) for the two methods show the same trend with the (001) surface being the most stable in the gas phase and the classical values are in agreement with previous values [5]. An important consideration for the interaction of uranyl with portlandite is whether U(VI) can or will be incorporated within the bulk or surface layers of the mineral. Table 1 and Table S1 include the optimised results from the DFT (PBE) computations on two possible incorporated structures: (a) U-portlandite, where a single calcium(II) in a 5 × 5 × 5 supercell of portlandite was substituted with a U(VI) ion and (b) UO 2 (OH) 4 2--portlandite, where a single Ca(OH) 6 unit was substituted by UO 2 (OH) 4 2− . For the U(VI)-substituted portlandite model, the overall structure remains similar to pure portlandite with supercell lattice parameters a = b = 17.945 and c = 24.555 (Table S1) and U-O distances of 2.197 Å and 2.577 Å (Table 1), compared to the corresponding Ca-O distance of 2.363 Å in bulk portlandite [17]. For the UO 2 substituted-portlandite system, the optimized supercell lattice parameters are a = 18.189 Å and c = 24.700 Å (Table S1), again almost identical to those predicted for bulk portlandite. However, importantly in this case, without the two protons to form hydroxyl, the complex retains its uranyl character but is distorted with U-O ax and U-OH distances of 1.87 Å and 2.26 Å, respectively, while the U-Ca distances are 3.751 Å, 3.755 Å, 3.739 Å and 3.749 Å in plane (Table 1).

Speciation in Aqueous Solution
Before modelling the interactions of uranyl with the surface of portlandite, it is necessary to first decide which aqueous species will dominate the solute inventory. Speciation models, which include equilibrium data for UO 2 (OH) 3 − and UO 2 (OH) 4 2− [53], and in which CO 2 /carbonate is limited, indicate that, under highly alkaline conditions (approaching the point of zero charge of portlandite at around pH 12.7 [4]), UO 2 (OH) 4 −2 is likely to be the most prevalent species in solution (see the Pourbaix diagram given in Figure S2a) [54]. These predictions are in line with previous studies [55,56], although at slightly lower pH, UO 2 (OH) 3 − will become prevalent. EXAFS and NMR studies, as well as quantum chemical calculations, have previously studied the [UO 2 (OH) 4 ] 2− species and have shown it to be a stable ion [49,[57][58][59][60][61][62]. Nevertheless, we shall also consider the more unusual UO 2 (OH) 5 3− ion, although this may only exist at extremely high pH [63].
In order to validate that our classical potential is able to predict the likely coordination of hydroxide ligands to uranyl, since these ligands would have the opportunity to exchange with water-or surface-bound portlandite hydroxides when forming innersphere adsorbed species, we have considered the following equilibria in aqueous conditions: (a) UO 2 (OH) 5 3 As noted in the previous section, a PMF has been computed for each of these three reactions using a coordinate corresponding to the dissociation of a hydroxyl ligand from the uranium centre. In each case, as with the sorption mechanism to portlandite, water molecules from solution are permitted to coordinate with the uranyl at every stage of the reaction but the classical potentials used do not permit protonation of any species; this is a major limitation of the model, such that the energy barriers to hydroxyl removal will likely be significantly overestimated in comparison with a real mechanism in which proton exchange would be likely. However, since the relative stability of the various hydroxyl species does not depend on the mechanistic path between them, the difference in the free energy between the bound and unbound forms in each PMF can be combined to give the relative stability of these species. The resulting PMFs for the first three equilibria are summarised in Figure 1 and clearly show that the most stable species is correctly predicted to be UO 2 (OH) 4 2− , in agreement with available experimental equilibrium data.
Relative to UO 2 (OH) 4 2− , species UO 2 (OH) 3 − , UO 2 (OH) 2 and UO 2 (OH) 5 3− are predicted to be <56 kJ·mol −1 ,~174 kJ·mol −1 and >25 kJ·mol −1 less stable, respectively. (More details on the aqueous complex calculations are given in the Supplementary Materials.) Taken overall, our aqueous complex calculations indicate that it is relatively easy and energetically favourable to remove the first hydroxyl group from a UO 2 (OH) 5 3− complex, but subsequent hydroxyl groups will less readily exchange for water; we therefore expect surface-adsorbed uranyl species to preferentially coordinate equatorially with four hydroxyl ligands. The high computational barriers for hydroxyl ligand exchange, either with water or surface hydroxyls, should be kept in mind in the surface-adsorption simulations since the surface species will not spontaneously or readily form inner-sphere complexes. Hence, standard molecular dynamics methods will be inappropriate to predict these species reliably, but the addition of the umbrella potential in our PMD simulations, constraining and forcibly reducing the surface to uranium distance, will encourage some ligand exchange when favourable. On the other hand, exchange of water ligands for hydroxyl is more energetically feasible and will be better described by the classical potentials, with surface hydroxyls more readily able to displace water.

Sorption of Uranyl Hydroxides on Portlandite
In this section, PMF simulations of the free energy change of bringing six different uranyl species into close proximity to the portlandite surface will be discussed. As noted in previous sections, the barriers for hydroxyl-to-water exchanges may preclude significant exchanges in the timescales of the MD simulations and, in reality, the mechanism for sorption on portlandite and the number of OH − /H 2 O ligands involved will be highly pH-dependent. In order to capture the essence of the likely thermodynamic outcome for the uptake of uranyl, we separately considered PMF simulations involving six uranyl hydroxide species, including the thermodynamically less-favoured hydroxylated species containing water ligands. Rather than just focussing on the sorption of UO 2 (OH) 4 2− , we have obtained surface-sorbed states for each of the species, thus accounting for species which may be produced via protonation of an OH − group and subsequent H 2 O-hydroxyl exchange near the surface. All of the bound surface species are, in principle, thermodynamically accessible, although the formation of the less-hydroxylated species would be unfavourable in a high pH environment, particularly from a kinetic perspective; therefore, those species are included simply for comparison.
To assess the importance of each surface species, the free energy of adsorption relative to the fully solvated aqueous species (for example, at a separation from the surface of >6 Å) has been obtained from the PMFs and is presented for each initial aqueous species in figure format, with the lowest energy configurations for each attached species depicted as insets on a set of free energy diagrams. In each figure, a free-energy profile is shown for adsorption to one of the three (001), (100) or (203/101) surfaces; however, the zero distance for each surface will depend on the surface topology and thus this is not directly comparable. A typical structure for each minimum energy region of the profile is shown on each figure; however, it should be noted that these are not stationary point structures and a more appropriate description in terms of the coordination numbers and positions of the peaks in the radial distribution functions are presented in Table 2 to understand the nature of the adsorbed species. In particular, the internuclear distances from the central uranium atom to three types of oxygen are shown: (a) U-OH eq (UO 2 ), where the OH − ligand can be attributed to the aqueous uranyl species; (b) U-OH eq (*), where the OH − ligand can be considered to be part of the Ca(OH) 2 surface; and (c) U-H 2 O eq , where a water molecule is coordinated to the uranyl. Figures for the first three reference species are given in the Supplementary Materials (Figures S3-S5), while those for the two species that are pertinent to an alkaline GDF-type fluid on pristine surfaces are given below in the main text (Figures 2 and 3). Also included is a figure detailing the calculations for UO 2 (OH) 4 2− on defect surfaces (Figure 4).
Preliminary PMF calculations detailed the adsorption of UO2 2+ and other species with low amounts of hydroxylation, as would be expected in the near-neutral pH region, and also the UO 2 (OH) 5 3− hich is probably negligible under GDF conditions (see the Preliminary PMF Adsorption Calculations section of the Supplementary Materials, Figure S6). Here, we now present the calculated adsorption behaviour for the species that represent the dominant uranyl complexes to be expected in a GDF scenario. Table 2. Analysis of the radial distribution functions from the classical simulations and PMF shown in Figures 1, 2 and S3-S6. Coordination numbers (CN) and distances (R) refer to the equatorial hydroxyl (or water) ligands of the sorbed uranyl complex and are separated into either those principally coordinated to uranyl (UO 2 ) or to the surface (*). The average coordination distances are shown for the principal RDF peaks obtained from windows representative of the sorbed complex. Where values are shown in parentheses, these correspond to a parallel ( ) uranyl orientation and values not in parentheses correspond to a perpendicular (⊥) oriented uranyl. Figures S7-S78 give details of all RDF calculations.

Ca(OH) 2 Surface
Aqueous Species U-OH eq (UO 2 ) 1 U-OH eq (*) U-H 2 O eq U-Ca    For the case of an alkaline pH up to approximately pH 12, the species of UO 2 (OH) 3 − is expected to dominate uranyl speciation. The free energy of adsorption onto the (100) surface is most favourable for this species at −45 kJ mol −1 (Figure 2). UO 2 (OH) 3 − adsorbs by dissociating a water ligand from the uranium atom and coordinating to an exposed hydroxyl group in the surface, while another hydroxyl group from the complex coordinates to two exposed calcium atoms, effectively forming a bound UO 2 (OH) 4 2− complex. On the (203/101) surface, a less energetically favourable (−16 kJ mol −1 ) but similar process occurs, whereby the species dissociates a water ligand and coordinates to an exposed hydroxyl group in the surface, while one of the complex hydroxyl groups coordinates to one exposed calcium ion in the surface, thereby bridging the calcium and uranium atoms. UO 2 (OH) 3 − complexes show favourable adsorption to the (001) surface at approximately −7.5 kJ mol −1 . The most favourable configuration has uranyl binding in a perpendicular orientation with a direct interaction of the uranyl oxygen to the surface.
Above pH 12, and probably most pertinent to the hyperalkaline conditions expected in a GDF, the UO 2 (OH) 4 2− species needs to be considered. Figure 3 presents the PMF calculations for this critical species. Again, the (203/101) shows a distinct minimum. Here, the [UO 2 (OH) 4 ] 2− complex adsorbs by coordinating one of its hydroxyl groups to an exposed calcium atom in the slab, bridging it and the uranium atom, while the complex orients to form hydrogen bonds with surface hydroxyl groups. This produces the most favourable free energy change out of all three surfaces at −22 kJ mol −1 . Adsorption on the (100) is less favourable, with a weak minimum at~4.7 Å distance, with adsorption proceeding by replacing a hydroxyl group from the surface (which then migrates and coordinates with two exposed calcium atoms in the surface). Adsorption onto the pristine (001) is unfavourable for this species.
Below, we summarize the general behaviour of these uranyl species on each crystallographic surface.
(001) surface: The PMFs for the sorption of UO 2 2+ , UO 2 (OH) + , UO 2 (OH) 2 and UO 2 (OH) 3 − complexes all show favourable adsorption to the surface at a PMF (collective) coordinate of~3 Å with Ca-U distances of~3.7 to 4.2 Å and have uranyl binding with the axial oxygens in a perpendicular orientation (with a direct interaction of a uranyl oxygen to the surface). All of the perpendicular complexes adsorb by hydrogen bonding to the surface hydroxyls and, with the exception of UO 2 (OH) 2 (by approximately 3 kJ mol −1 ), this is slightly more favourable than forming an inner-sphere complex. At a shorter PMF coordinate of~1.3 Å with Ca-U distances < 5.5 Å, UO 2 2+ , UO 2 (OH) + and UO 2 (OH) 2 are also able to form stable bridged inner-sphere complexes with the O=U=O parallel to the surface by exchange of equatorial water with two surface hydroxyl groups. Since the classical energies of the perpendicular and parallel sorbed complexes are very similar, they will be explored in more detail using DFT; however, the weak binding of UO 2 2+ and UO 2 (OH) + suggests that they would easily desorb from the (001) surface. For UO 2 (OH) 4 2− and UO 2 (OH) 5 3− , perpendicular adsorption to the (001) surface is found to be unfavourable in these calculations. UO 2 (OH) 3 − , UO 2 (OH) 4 2− and UO 2 (OH) 5 3− are furthermore predicted to be unable to form stable parallel inner-sphere complexes.
(203/101) surface: All considered complexes favourably absorb to the (203/101). The first three, UO 2 2+ , UO 2 (OH) + and UO 2 (OH) 2 , all have the ability to burrow into the surface to form a stable adsorbate. The more hydroxylated species, UO 2 (OH) 3 − and UO 2 (OH) 4 −2 , do not display burrowing but are still able to form stable surface species which display relatively large negative free energy changes-especially in the case of UO 2 (OH) 4 −2 , which is the most probable species to be present in a hyperalkaline fluid. Adsorption of UO 2 (OH) 5 3− on this surface is not favourable. Although sorption onto the fully hydroxylated (100) surface is relatively weak for UO 2 (OH) 4 −2 and UO 2 (OH) 5 3− , sorption of these two important species is certainly thermodynamically favourable onto the pristine (203/101) surface, where the lower density of hydroxyl ligands permits bridging and ligand exchange with the uranyl. A five-hydroxyl uranyl species is clearly less favoured, and one of its hydroxyl groups will readily dissociate if there is an available site on the surface for it to bind to. Overall, it is clear that no matter which aqueous species is considered, a four-hydroxide coordinated inner-sphere adsorbed structure would be ultimately preferred, and at high pH, the free energy of sorption to (clean) portlandite may be in the region of approximately −20 kJ mol −1 , if dominated by [UO 2 (OH) 4 ] −2 onto (203/101) faces, as expected.

Sorption on Defect Portlandite Surfaces
Although the clean surfaces provide a useful baseline for sorption, the surface of portlandite will be quite dynamic in aqueous solutions [64]. PMF profiles for sorption of the UO 2 (OH) 4 2− complex to the (001-OH), (100-OH) and (203-OH) surfaces, where a surface OH − has been protonated to H 2 O, are shown in Figure 4. In each case, there are distinct minima corresponding to inner-sphere complexes bound by −183 kJ mol −1 , −40 kJ mol −1 and −221 kJ mol −1 , respectively, considerably stronger than to the clean surfaces. The radial distribution functions, summarised in Table 3, indicate that uranyl maintains a coordination number of 4 OH − around uranium at a distance of 2.28 Å for all three surfaces and no ligand exchange. Each complex binds strongly through a direct interaction of an equatorial uranyl hydroxyl directly to one or more exposed surface calcium ions, displacing the surface water at the vacancy site. In the majority of surface complexes, the uranyl axis is oriented relatively parallel to each surface and the overall structures are similar to UO 2 (OH) 3 − complexes adsorbed to the clean portlandite surfaces. Table 3. Analysis of the radial distribution functions from the classical simulations and PMF of the OH-vacancy defect surfaces shown in Figure 6. Coordination numbers (CN) and distances (R) refer to the equatorial hydroxyl (or water) ligands of the sorbed uranyl complex and are separated into either those principally coordinated to uranyl (≈) or to the surface (>). The average coordination distances are shown for the principal RDF peaks obtained from windows representative of the sorbed complex.

DFT Study of the Surface Complexes
While the PMF computations provide an important quantitative assessment of the free energy change associated with bringing the various uranyl species onto the portlandite surfaces and thus give a useful overview of potential adsorbate-surface affinities, DFT computations allow us to optimize the surface geometries of possible adsorbate structures in a more robust manner. Here, we focus on the dominant aqueous species and provide details of sorbate configuration on pristine and defect (001), (100), and (203/101) surfaces. In Table 4, we present key structural details for the 12 periodic cluster models corresponding to the sorption of UO 2 (OH) 4 2− with the uranyl axis either perpendicular or parallel to the three clean ( Figure 5) or three defect ( Figure 6) surfaces, as observed in the classical studies. The schematic figures show the optimised complexes but omit water and the majority of the surface for clarity.   All of the optimised structures retained a uranyl equatorial coordination number of four for hydroxyl and there were no water molecules in the inner shell. For comparison with the classical PMF results, the equatorially coordinated hydroxyl groups have been partitioned into those bridging to the surface U-OH eq (>) and those exposed to solvent, U-OH eq (≈). The optimised conformations (I) to (XII) can be summarised as follows (the lower-energy configuration between each pair of possible orientations is indicated with an asterisk): 1.
(001, ⊥)*-uranyl oxygen coordinates to the hydrogen atom of a surface hydroxyl group at a distance of 1.80 Å; 2.
(100, ⊥)-uranyl oxygen coordinates to the three exposed calcium ions with O- (203/101, ⊥)-uranyl oxygen atom coordinates to the hydrogen atom of a surface hydroxyl group at a distance of 1.75 Å; 6.
(203/101, )-one uranyl hydroxide group coordinates to an exposed calcium ion and is hydrogen bonded to a surface hydroxyl group with distances of 2.41 Å and 1.83 Å, respectively; a second uranyl hydroxide group also hydrogen bonds to another surface hydroxyl group at a distance of 1.78 Å; 7.
(001-OH, )-a hydroxyl group from the uranyl complex substitutes for the missing surface hydroxyl group with O-Ca distances of 2.46 Å and 2.60 Å; 9.
(100-OH, ⊥)-uranyl oxygen coordinates to an exposed calcium atom at a distance of 2.40 Å; 10. (100-OH, )*-two uranyl hydroxide groups coordinate to two exposed calcium ions at distances of 2.40 Å and 2.60 Å; 11. (203-OH, ⊥)*-uranyl oxygen atom coordinates to two exposed calcium ions at distances of 2.49 Å and 2.65 Å, replacing the missing hydroxyl group; 12. (203-OH, )-a hydroxyl group from the uranyl complex substitutes for the missing surface hydroxyl group with O-Ca distances of 2.56 Å and 2.66 Å, while another hydroxyl group from the uranyl complex hydrogen bonds to a surface hydroxyl group at a distance of 1.75 Å.
The perpendicular-oriented uranyl structures to the (001) and (100) clean surfaces (>-OH) appear to be weakly interacting outer-sphere complexes with >-OH − . . . O=U=O separated by 1.75-1.80 Å. For the other surfaces, including the defect surfaces, the uranyl oxygen can clearly interact more directly with surface calcium ions to form more innersphere-like complexes in these cases. When oriented parallel to the surface, the uranyl more clearly forms inner-sphere complexes with at least one bridging OH − at the vacancy site. The U-OH bond lengths were found to average between 2.24 Å and 2.31 Å, typical of r(U-OH) in isolated UO 2 (OH) 4 2− , using the PBE functional (see Table 1). The axial U-O bond lengths were all found to average in the range of 1.84 Å to 1.89 Å. Although longer than typical U-O ax lengths experimentally observed in uranyl phases (1.72-1.82 Å) and more consistent with uranate phases (1.86-1.97 Å), here, the computational value with the PBE functional is more typical of uranyl in UO 2 (OH) 4 2− . Overall, it is interesting to note the long Ca-U distances in all complexes: at the lower end, in the defect (100) and (203) cases, they are shorter at 3.51 to 3.74 Å, and at the upper end, 5.06 to 5.80 Å, when perpendicular to clean (100) and (001) surfaces. The other cases fall between these extremes at 3.96 and 4.20 Å. It should be noted that the optimised structures are assumed to be local minima and it is not meaningful to consider these structures as absolute models of sorption since more extensive sampling of alternative structures was not possible at this level.
The relative potential energies of the two orientations of uranyl for each surface are compared in Table 5; however, due to the different surface compositions and the large energetic dependence on the configurations of the water molecules, it is not possible to directly compare energies between different surfaces-the values are intended to compare only the two different orientations on the same surface. In order to assess the relative interaction energy of the perpendicular (⊥) and parallel ( ) orientations on the same surface, we have re-computed single-point energies of the surface and uranyl complex at the DFToptimised structure, but have omitted all of the water molecules. Thus, the energetic differences are considerably larger in magnitude than those obtained as free energies by the classical simulations since they do not include any stabilisation due to solvation. The results are quite conclusive in the prediction of favourable (gas phase) interactions of uranyl in a perpendicular relative to a parallel orientation with both the (001) and (203/101) surfaces by −33 to −39 kJ mol −1 and −82 to −93 kJ mol −1 , respectively. For the (100) surface configurations, the situation is reversed and a parallel uranyl orientation is preferred over the perpendicular by 78 to 144 kJ mol −1 . These findings are broadly in line with the results of the classical simulations for the full sorption of the complexes, where the (001) and (203/101) structures favour direct interaction of the surface with the uranyl oxygen. However, the DFT results predict slightly different bonding environments, with (001, ⊥), (100, ), and (203/101, ⊥) whereas the PMF computations predict (001, ), (100, ) and (203/101, sub-parallel). Table 5. Relative potential energies of the surface complexes excluding water molecules at the DFT(PBE)-optimised structures. The energy difference (E(⊥)-E( ) is for the perpendicular orientation relative to the parallel orientation.

Surface
Energy Difference/(kJ mol −1 ) (001) −38.7  Table 6 shows the results from fitting both the ex-situ and in-situ EXAFS data for U(VI) with hydroxylated CaO. Four sets of reaction conditions were fully analysed:

Experimental EXAFS Analysis
13. pH 12 in KOH, 48 h reaction time before ex situ data acquisition; 14. pH 11 with 2 mM CO 3 2− , pre-prepared ex situ 48 h reaction time before data acquisition; 15. pH 12 in KOH, reacted in situ, sorption for approximately 2 h before data acquisition; 16. pH 12 in LiOH, reacted in situ, sorption for approximately 2 h before data acquisition. For the in situ experiments, both XANES spectra of the U(VI)-reacted CaO (111) surfaces showed signs of broadening of the U(VI) absorption peak at the white line (Figure 7, feature A), very similar to XANES reference spectra of Ca-uranate [9] and spectra of uranates published by Bots et al. [65]. This is indicated on Figure 7 by the increase in intensity at feature 'B'. As a result, the feature in the shoulder region (feature C in Figure 7) that is typical of an undistorted U-O axial of the uranyl moiety disappears in both in situ experiments. This implies that the U-O axial bond lengths are elongated after reaction and that the U(VI) bonding environment is distorted upon reaction with the mineral surface. Undistorted uranyl has short U=O bond lengths (1.77-1.81 Å) [66], which produces the diagnostic shoulder feature clearly visible in the uranyl reference spectrum (bottom curve on Figure 7).
In both in situ experiments, the EXAFS fits (see Figure 8) confirmed a significant increase in the U-O axial bond lengths to between 1.86 and 1.93 Å. In addition, four to five equatorial oxygen atoms at a distance of 2.23-2.35 Å could be fitted, which supports the idea that UO 2 (OH) 4 2− is the dominant species involved in adsorption to the mineral surface. Despite the fact that all solutions were undersaturated with compreignacite (K 2 (UO 2 ) 6 O 4 (OH) 6 ·7H 2 O), a further~four Ca atoms and 1-2 U atoms could also be fitted within the spectra. This could indicate some surface co-precipitation of U(VI). The measured CN and radial distances would be in agreement with the formation of calcium uranate phases previously published [10,65,67]. These data will be further discussed in the next section.  In contrast, the ex situ spectra retain feature C (Figure 7), typical of the U=O axial of uranyl, especially with carbonate present. Unlike the in situ experiments, where the surface was freshly reacted and carefully rinsed, the ex situ crystal surfaces roughened significantly during reaction and were not rinsed with MilliQ water. This difference in U(VI) speciation between in situ and ex situ experiments is likely due to differences in sample preparation and surface reactivity with bulk experiments presumably retaining undistorted bulk uranyl, whilst the in situ experiments allow for dynamic reaction of the uranate species to be identified. Linear combination fits of the XANES confirmed a mixture of uranyl and uranate with uranyl being predominantly present in the ex situ samples. This is why the XANES spectra are different; therefore, we conclude that the ex situ measurements provide less reliable configurations of adsorbed species in the presence of an aqueous fluid. We conclude that the in situ EXAFS results under equivalent conditions are more likely to give a reliable understanding of the adsorption of the uranyl-hydroxide complexes, especially with as reactive a surface as CaO(111). As shown by fitting the in situ EXAFS data ( Figure 8, Table 6), adsorption of UO 2 (OH) 4 2− onto a surface equivalent to the portlandite (001) results in an adsorbate complex with elongated axial oxygen distances, discernible Ca and U backscatterers consistent with surface attachment and fourfold equatorial oxygen coordination with distances of 2.35 to 2.28 Å. All of these spectral details are consistent with the lowest energy configuration calculated by DFT on the portlandite (001) surface.

Discussion
The high pH experimental EXAFS data clearly fit to a coordination environment containing two shorter axial U=O bonds and four longer equatorial U-O bonds, which, as stated, is consistent with the computational predictions of both aqueous and surface complexes favouring coordination to four hydroxyl ligands. The average DFT bond lengths (1.84-1.89 Å) for the U=O bonds are close to those measured by the experimental analyses (1.86-1.93 Å), with both indicating a tendency to be longer than aqueous [UO 2 (OH) 4 ] 2− , at 1.83 and~1.85 Å for the experiment and the PBE functional, respectively. The calculated equatorial bonds are in very good agreement with DFT predictions in the range of 2.18-2.38 Å, again matching those from the experimental EXAFS data at 2.28-2.35 Å.
Perhaps the most important aspect of the EXAFS analyses is that the data are not consistent with the U-O environment, which would result from precipitation-deposited crystalline calcium uranate (CaUO 4 ). This is because, although calcium uranate has two distinct U-O coordination distances in its crystalline form, at 1.966 and 2.282 Å [68,69], these are distinct from the distances in uranyl. Typically, U−O axial bond lengths in uranate phases are within the range (1.86−1.97 Å) [10,67,70], longer than the U−O axial bond lengths in uranyl phases (1.70-1.82 Å) [68,69], and the U-O eq distances for several layered metal uranate phases have been reported within the range (2.15-2.30 Å), shorter in general than in uranyl phases (2.27-2.49 Å).
Although the in situ pH 12 (LiOH) EXAFS gives U-Ca scattering (at 3.59-3.78 Å) close to the calcium uranate system (3.653 Å), these data are also similar to the inner-sphere models predicted for uranyl sorption. U-Ca distances consistent with the experimental data are consistently predicted for a wide range of the sorption complexes observed in the classical simulations (3.68 to 4.33 Å) and in the more likely inner-sphere complexes studied by DFT (3.51 to 4.20 Å). Furthermore, the longer U-Ca value, which is observed in the other pH 12 in situ (KOH) experiment (4.08 Å), is not consistent with uranate but does match with inner-sphere DFT models. The shorter ex situ U-Ca distance (3.08 Å) is not observed in any of our sorption models but may be a result of scattering from Ca 2+ coordinated to uranyl anionic species, such as CaUO 2 (OH) 4 . The DFT structure of CaUO 2 (OH) 4 illustrates such a short U-Ca distance (3.14 Å).
The in situ EXAFS U-U distances at pH 12 (3.83-3.84 Å) are similar to uranate (3.850 Å) [68] and do not match uranyl minerals such as metastudtite (UO 4 .2H 2 O) [71], which has a longer U-U distance (4.21 Å). It is interesting to compare the experimental U-U data with our computational predictions where U(VI) or uranyl is incorporated within a crystalline portlandite environment. When incorporated, we predict U-Ca distances between 3.74 and 3.76 Å, at the longer end of the range seen in the in situ pH 12 (LiOH) experiment and notably longer than typical Ca-Ca distance in portlandite (3.6 Å). Although, in the incorporated models, we have not embedded two uranium species, DFT computations on the molecular uranyl dimer complex (UO 2 ) 2 (OH) 6 2− resulted in a U-U distance of 3.98 Å, which suggests the longer experimental values are unlikely to be due to dimeric sorption complexes. Taking a holistic view, the U-O axial experimental observations strongly suggest that U(VI) associated with the portlandite (001) surface is consistent with a mixture of uranyl and uranate-like coordination environments. The computational studies predict a range of possibly stable inner-sphere sorbed uranyl-hydroxyl complexes, particularly to (203/101) or defect surfaces, which are consistent with the experimental fit to 4 longer equatorial U-O bonds. It is also quite clear that the surface complexation at a high pH cannot be described by a single model and a range of adsorbed uranyl, incorporated uranyl and uranate precipitate are all likely to contribute to the uptake of uranium at the surface of portlandite under these experimental conditions.

Conclusions
Density functional theory and classical molecular dynamics calculations have been performed on various uranyl hydroxide species in the presence of water and portlandite. Uranyl has been shown to adopt a four-coordinate structure, with respect to hydroxyl groups, both in solution and when adsorbed to portlandite surfaces, consistent with experimental EXAFS data at high pH.
Although uranyl is found to only weakly adsorb to the (001) and (100) clean surfaces, significantly stronger interactions are predicted with the (203/101) surfaces prevalent under groundwater conditions, at OH-vacancy sites or on a rough surface. The uranyl ions are typically found to form sorbed species with four equatorially coordinated hydroxyl ligands with direct interaction of a uranyl oxygen with surface calcium ions in the (001) and (203/101) cases, but uranyl is likely to be sorbed more parallel to the (100) surface.
The experimental EXAFS data looking at the CaO (111) analogue of the portlandite (001) surface are largely consistent with a six-coordinate structural layer or a deposit similar to calcium uranate, although there is also strong evidence for uranyl-type coordination, as would be expected, to originate from the adsorbates predicted by the computational models. It is quite clear from the X-ray absorption experiments that the short preparation timescales of the in situ and ex situ experiments, of 2 to 48 h, respectively, indicate that the uptake of uranyl by portlandite is extremely rapid, particularly in the context of GDF timescales.

Funding:
The authors gratefully acknowledge financial support for this research from the UK Nuclear Decommissioning Agency (RWMD) via EPSRC grant EP/1036389/1.

Data Availability Statement:
Data used in this study may be made available by request to the corresponding authors.