Michael Holmboe

Porosity investigation of compacted bentonite using XRD profile modeling

Many countries intend to use compacted bentonite as a barrier in their deep geological repositories for nuclear waste. In order to describe and predict hydraulic conductivity or radionuclide transport through the bentonite barrier, fundamental understanding of the microstructure of compacted bentonite is needed. This study examined the interlayer swelling and overall microstructure of Wyoming Bentonite MX-80 and the corresponding homo-ionic Na(+) and Ca(2+) forms, using XRD with samples saturated under confined swelling conditions and free swelling conditions. For the samples saturated under confined conditions, the interparticle, or so-called free or external porosity was estimated by comparing the experimental interlayer distances obtained from one-dimensional XRD profile fitting against the maximum interlayer distances possible for the corresponding water content. The results showed that interlayer porosity dominated total porosity, irrespective of water content, and that the interparticle porosity was lower than previously reported in the literature. At compactions relevant for the saturated bentonite barrier (1.4-1.8 g/cm(3)), the interparticle porosity was estimated to ≤3%.

MOLECULAR DYNAMICS SIMULATIONS OF ANION EXCLUSION IN CLAY INTERLAYER NANOPORES

The aqueous chemistry of water films confined between clay mineral surfaces remains an important unknown in predictions of radioelement migration from radioactive waste repositories. This issue is particularly important in the case of long-lived anionic radioisotopes (129I-, 99TcO4-, 36Cl-) which interact with clay minerals primarily by anion exclusion. For example, models of ion migration in clayey media do not agree as to whether anions are completely or partially excluded from clay interlayer nanopores. In the present study, this key issue was addressed for Cl- using MD simulations for a range of nanopore widths (6 to 15 Å) overlapping the range of average pore widths that exists in engineered clay barriers. The MD simulation results were compared with the predictions of a thermodynamic model (Donnan Equilibrium model) and two pore-scale models based on the Poisson-Boltzmann equation under the assumption that interlayer water behaves as bulk liquid water. The simulations confirmed that anion exclusion from clay interlayers is greater than predicted by the pore-scale models, particularly at the smallest pore size examined. This greater anion exclusion stems from Cl- being more weakly solvated in nano-confined water than it is in bulk liquid water. Anion exclusion predictions based on the Poisson-Boltzmann equation were consistent with the MD simulation results, however, if the predictions included an ion closest approach distance to the clay mineral surface on the order of 2.0 ± 0.8 Å. These findings suggest that clay interlayers approach a state of complete anion exclusion (hence, ideal semi-permeable membrane properties) at a pore width of 4.2 ± 1.5 Å.

COLLOID DIFFUSION IN COMPACTED BENTONITE: MICROSTRUCTURAL CONSTRAINTS

In Sweden and in many other countries, a bentonite barrier will be used in the repository for spent nuclear fuel. In the event of canister failure, colloidal diffusion is a potential, but scarcely studied mechanism of radionuclide migration through the bentonite barrier. Column and in situ experiments are vital in understanding colloid diffusion and in providing information about the micro structure of compacted bentonite and identifying cut-off limits for colloid filtration. This study examined diffusion of negatively charged 2-, 5-, and 15-nm gold colloids in 4-month diffusion experiments using MX-80 Wyoming bentonite compacted to dry densities of 0.6–2.0 g/cm3. Breakthrough of gold colloids was not observed in any of the three diffusion experiments. In a gold-concentration profile analysis, colloid diffusion was only observed for the smallest gold colloids at the lowest dry density used (estimated apparent diffusivity Da ≈ 5 × 10−13 m2/s). The results from a microstructure investigation using low-angle X-ray diffraction suggest that at the lowest dry density used, interlayer transport of the smallest colloids cannot be ruled out as a potential diffusion pathway, in addition to the expected interparticle transport. In all other cases, with either greater dry densities or larger gold colloids, compacted bentonite will effectively prevent diffusion of negatively charged colloids due to filtration.

Partitioning into Colloidal Structures of Fasted State Intestinal Fluid Studied by Molecular Dynamics Simulations

We performed molecular dynamics (MD) simulations to obtain insights into the structure and molecular interactions of colloidal structures present in fasted state intestinal fluid. Drug partitioning and interaction were studied with a mixed system of the bile salt taurocholate (TCH) and 1,2-dilinoleoyl-sn-glycero-3-phosphocholine (DLiPC). Spontaneous aggregation of TCH and DLiPC from unconstrained MD simulations at the united-atom level using the Berger/Gromos54A7 force fields demonstrated that intermolecular hydrogen bonding between TCH molecules was an important factor in determining the overall TCH and DLiPC configuration. In bilayered systems, these intermolecular hydrogen bonds resulted in embedded transmembrane TCH clusters. Free energy simulations using the umbrella sampling technique revealed that the stability of these transmembrane TCH clusters was superior when they consisted of 3 or 4 TCH per bilayer leaflet. All-atom simulations using the Slipids/GAFF force fields showed that the TCH embedded in the bilayer decreased the energy barrier to penetrate the bilayer (ΔGpen) for water, ethanol, and carbamazepine, but not for the more lipophilic felodipine and danazol. This suggests that diffusion of hydrophilic to moderately lipophilic molecules through the bilayer is facilitated by the embedded TCH molecules. However, the effect of embedded TCH on the overall lipid/water partitioning was significant for danazol, indicating that the incorporation of TCH plays a crucial role for the partitioning of lipophilic solutes into e.g. lipidic vesicles existing in fasted state intestinal fluids. To conclude, the MD simulations revealed important intermolecular interactions in lipidic bilayers, both between the bile components themselves and with the drug molecules.

Atomistic simulations of cation hydration in sodium and calcium montmorillonite nanopores

During the last four decades, numerous studies have been directed to the swelling smectite-rich clays in the context of high-level radioactive waste applications and waste-liners for contaminated sites. The swelling properties of clay mineral particles arise due to hydration of the interlayer cations and the diffuse double layers formed near the negatively charged montmorillonite (MMT) surfaces. To accurately study the cation hydration in the interlayer nanopores of MMT, solvent-solute and solvent-clay surface interactions (i.e., the solvation effects and the shape effects) on the atomic level should be taken into account, in contrast to many recent electric double layer based methodologies using continuum models. Therefore, in this research we employed fully atomistic simulations using classical molecular dynamics (MD) simulations, the software package GROMACS along with the CLAYFF forcefield and the SPC/E water model. We present the ion distributions and the deformation of the hydrated coordination structures, i.e., the hydration shells of Na+ and Ca2+ in the interlayer, respectively, for MMT in the first-layer, the second-layer, the third-layer, the fourth-layer, and the fifth-layer (1W, 2W, 3W, 4W, and 5W) hydrate states. Our MD simulations show that Na+ in Na-MMT nanopores have an affinity to the ditrigonal cavities of the clay layers and form transient inner-sphere complexes at about 3.8 Å from clay midplane at water contents less than the 5W hydration state. However, these phenomena are not observed in Ca-MMT regardless of swelling states. For Na-MMT, each Na+ is coordinated to four water molecules and one oxygen atom of the clay basal-plane in the first hydration shell at the 1W hydration state, and with five to six water molecules in the first hydration shell within a radius of 3.1 Å at all higher water contents. In Ca-MMT, however each Ca2+ is coordinated to approximately seven water molecules in the first hydration shell at the 1W hydration state and about eight water molecules in the first hydration shell within a radius of 3.3 Å at all higher hydration states. Moreover, the MD results show that the complete hydration shells are nearly spherical with an orthogonal coordination sphere. They could only be formed when the basal spacing d001 ≥ 18.7 Å, i.e., approximately, the interlayer separation h ≥ 10 Å. Comparison between DFT and MD simulations shows that DFT failed to reproduce the outer-sphere complexes in the Stern-layer (within ∼5.0 Å from the clay basal-plane), observed in the MD simulations.