Search NASA⌕ Search

SEARCH · Search NASA

Results for “Nuclear Configuration Interaction”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Seismic response of vertical dry storage casks under three-dimensional earthquake motions

Ensuring the long-term seismic safety of dry storage casks (DSCs) is becoming increasingly critical as these systems evolve from temporary to de facto permanent repositories for spent nuclear fuels. Traditional seismic soil–structure interaction (SSI) assessment methods use one-dimensional deconvolution or simplified boundary conditions to model incident waves. Although computationally appealing, simplifying assumptions may alter the seismic risk by neglecting the full complexity of three-dimensional (3D) wave propagation effects. To address this challenge, this paper introduces a novel high-fidelity computational framework that leverages the Domain Reduction Method (DRM) with perfectly matched layers (PML) to accurately transfer complex, 3D seismic wavefields from regional-scale fault-rupture simulations into local-scale finite element models of DSCs. Using broadband, physics-based ground motions from a generic M w 7.0 strike-slip event, both single-cask and multi-cask configurations were investigated under near- and far-field conditions. Emphasis is placed on capturing complex SSI, spatial variability in the ground motion, and nonlinear phenomena such as cask rocking and sliding. Numerical results demonstrate that near-field conditions, where forward directivity and fling-step effects dominate, lead to significantly higher DSC rocking and sliding. Far-field cases, by contrast, generally exhibit modest responses. Incorporating SSI tends to amplify or alter DSC response spectra and introduce response variability, which underscores the need for site-specific evaluations and robust modeling approaches to ensure the seismic integrity of DSCs in interim spent fuel storage installations.

Das, Tonmoy↗
DOE OSTI2021

Valence shell photoelectron angular distributions and vibrationally resolved spectra of imidazole: A combined experimental–theoretical study

Linearly polarized synchrotron radiation has been used to record polarization dependent valence shell photoelectron spectra of imidazole in the photon energy range 21–100 eV. These have allowed the photoelectron angular distributions, as characterized by the anisotropy parameter β, and the electronic state intensity branching ratios to be determined. Complementing these experimental data, theoretical photoionization partial cross sections and β-parameters have been calculated for the outer valence shell orbitals. The assignment of the structure appearing in the experimental photoelectron spectra has been guided by vertical ionization energies and spectral intensities calculated by various theoretical methods that incorporate electron correlation and orbital relaxation. Strong orbital relaxation effects have been found for the 15a', nitrogen lone-pair orbital. The calculations also predict that configuration mixing leads to the formation of several low-lying satellite states. The vibrational structure associated with ionization out of a particular orbital has been simulated within the Franck–Condon model using harmonic vibrational modes. The adiabatic approximation appears to be valid for the X 2A" state, with the β-parameter for this state being independent of the level of vibrational excitation. However, for all the other outer valence ionic states, a disparity occurs between the observed and the simulated vibrational structure, and the measured β-parameters are at variance with the behavior expected at the level of the Franck–Condon approximation. Furthermore, these inconsistencies suggest that the excited electronic states may be interacting vibronically such that the nuclear dynamics occur over coupled potential energy surfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

New isomeric transition in Mg 36 : Bridging the N = 20 and N = 28 islands of inversion

Here, we observed a new isomeric γ transition at 168 keV in 36 Mg, with a half-life of T1/2=90⁢($^{+410}_{-50}$) ns. We propose that the observed transition deexcites a new 0 + isomeric state at 833 keV and populates the previously known first 2 + state. The existence of this isomer is consistent with the predictions of the large-scale shell-model calculations of 36 Mg using the sdpf-u-mix interaction. The observed excitation energy of the second 0 + state is caused by the small energy separation between two prolate-deformed configurations where the intruder configuration corresponds to two-neutron excitations from the sd to the pf shell. Within this interpretation, 36 Mg becomes the crossing point between nuclei in which ground state deformed/superdeformed configurations are caused by the dominance of N=20 intruders ( 32,34 Mg) and nuclei where deformed configurations are associated with the breaking of the N=28 closure and a large occupancy of the 1⁢p 3/2 neutron orbit ( 38 Mg and beyond). We found the lack of three-body monopole corrections in other effective interactions results in a predominance of N=20 intruder configurations past 38 Mg incompatible with our observation. We conclude that 36 Mg bridges the N=20 and N=28 islands of inversion, forming the so-called big island of deformation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2015

High-temperature strain cell for tomographic imaging

This disclosure provides systems, methods, and apparatus related to the high temperature mechanical testing of materials. In one aspect, a method includes providing an apparatus. The apparatus may include a chamber. The chamber may comprise a top portion and a bottom portion, with the top portion and the bottom portion each joined to a window material. A first cooled fixture and a second cooled fixture may be mounted to the chamber and configured to hold the sample in the chamber. A plurality of heating lamps may be mounted to the chamber and positioned to heat the sample. The sample may be placed in the first and the second cooled fixtures. The sample may be heated to a specific temperature using the heating lamps. Radiation may be directed though the window material, the radiation thereafter interacting with the sample and exiting the chamber through the window material.

42 ENGINEERING↗
DOE OSTI2016

An exact factorization perspective on quantum interferences in nonadiabatic dynamics

Nonadiabatic quantum interferences emerge whenever nuclear wavefunctions in different electronic states meet and interact in a nonadiabatic region. In this work, we analyze how nonadiabatic quantum interferences translate in the context of the exact factorization of the molecular wavefunction. In particular, we focus our attention on the shape of the time-dependent potential energy surface—the exact surface on which the nuclear dynamics takes place. We use a one-dimensional exactly solvable model to reproduce different conditions for quantum interferences, whose characteristic features already appear in one-dimension. The time-dependent potential energy surface develops complex features when strong interferences are present, in clear contrast to the observed behavior in simple nonadiabatic crossing cases. Nevertheless, independent classical trajectories propagated on the exact time-dependent potential energy surface reasonably conserve a distribution in configuration space that mimics one of the exact nuclear probability densities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

The NuMI Flux Prediction at ICARUS

The Deep Underground Neutrino Experiment (DUNE) is a next-generation long-baseline neutrino oscillation experiment seeking to probe fundamental symmetries within the structure of the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) mixing matrix, and perform precision measurements its parameters including the neutrino mass ordering via the sign of $\Delta m^2_{31}$, and the charge-parity violating phase, $\delta_{CP}$. To make these measurements with high precision, DUNE will require external $\nu$-Ar scattering cross section data as a crucial input to the oscillation fit. Imaging Cosmic And Rare Underground Signals (ICARUS) is a 476 t liquid argon neutrino detector located at Fermi National Accelerator Laboratory (FNAL) where it is serving as the far detector for the Short-Baseline Neutrino (SBN) program along the Booster Neutrino Beam (BNB) axis. ICARUS additionally lies 795 m downstream and 100.1 mrad off-axis of the Neutrinos at the Main Injector (NuMI) neutrino beam. From this position, ICARUS is exposed to a large flux of NuMI (anti-)electron and (anti-)muon neutrinos, and poses a unique opportunity to provide high-statistics measurements of quasi-elastic and single pion-production cross sections for four neutrino flavors ($\nu_{\mu}$, $\nu_{e}$, $\bar{\nu}_{\mu}$, $\bar{\nu}_{e}$). This dissertation is centered around accurately characterizing the models and estimating their precision for use in making these measurements. This includes identifying major sources of uncertainty in the models such that they can be properly propagated to the cross section measurements. Specifically, this work focused on the model of the NuMI beamline and its impact on the neutrino fluxes, but also delved into the detector response model and its impact on reconstructed observables in the detector. Significant efforts were made to improve the characterization to enhance precision and thus reduce the level of propagated uncertainty. In particular, the NuMI flux was determined to be composed of 57% $\nu_{\mu}$, 38% $\bar{\nu}_{\mu}$, 3% $\nu_{e}$, and 2% $\bar{\nu}_{e}$ while the horns are operating in the positive-particle focusing configuration. The total uncertainty on the $\nu_{\mu} + \bar{\nu}_{\mu}$ ($\nu_{e} + \bar{\nu}_{e}$) flux while operating in the forward horn operating mode was determined to be 10.84% (9.04%). Compared to the on-axis flux, mesons that eventually decay to neutrinos more frequently reinteract within the NuMI structure, resulting in elevated uncertainty as these processes are not well-constrained by existing hadron interaction cross section measurements. Covariance matrices were calculated to propagate the flux uncertainty characterization to NuMI analyses.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2018

Bradwell - The first UK Nuclear Site to Enter Care and Maintenance - 18152

Magnox Ltd, a Cavendish Fluor company, is currently decommissioning twelve nuclear sites on behalf of the UK Nuclear Decommissioning Authority (NDA). Ten of these sites are former Power Stations; the other two are closed Research and Development facilities. An overall strategy for decommissioning the Magnox sites has been agreed between Magnox and the NDA. For all but one of these sites the decommissioning strategy involves an extended quiescent period, known as Care and Maintenance (C and M), during which remaining radiological hazards reduce naturally through radioactive decay. After this quiescent period the sites will be restored to an agreed end state through further decommissioning and clearance work. Sites in C and M will not be permanently manned. They will be passively safe such that prompt action will not be required in response to an alarm or fault. Preparing a nuclear site for C and M involves several inter-related strands of work including definition of the intended physical configuration of the site during C and M (the C and M entry state), establishing a safety case to support that configuration, carrying out the physical works needed to deliver the C and M entry state, preparing arrangements and a team to manage the site during C and M and consolidating the knowledge and information needed to support C and M and final clearance of the site. To enter C and M there needs to be a carefully controlled transition of site, personnel and management system into the agreed C and M configuration. Bradwell, a site on the east coast of England, was selected as part of the overall decommissioning strategy for 'early' entry into C and M; Bradwell is scheduled to enter C and M in late 2018, several years in advance of the next scheduled site, with the intent of carrying learning forward from Bradwell to the subsequent sites. Bradwell is now in the final stages of preparation for C and M. Decommissioning commenced on completion of defueling in 2005. A clearly defined C and M entry state was established early in decommissioning and, through interaction with Regulators and internal and external stakeholders, developed as a shared target, thus providing certainty to the scope of the overall decommissioning project. The site has subsequently been systematically decommissioned to meet this target C and M entry state. All legacy waste has been retrieved and packaged into robust disposal containers, the reactor buildings have a new weather envelope, and much of the site is empty space with concrete hardstanding where facilities used to exist and bare rooms where plant used to operate. The remaining radiological and environmental hazards are not significant. During C and M, which is anticipated to last for 70 years, Bradwell will be managed remotely without the requirement for personnel to be present on site, with the possible exception of an ongoing security presence when the site first enters C and M. In terms of UK decommissioning, many of the projects undertaken at Bradwell were first of a kind. The company operates under a programmed organisational and delivery model, which will enable innovation and learning from Bradwell to be carried forward across the whole Magnox fleet. In addition to decommissioning tools and techniques, processes to manage decommissioning, optimise the C and M entry state to reduce the burden on the UK taxpayer and to transition a site into C and M have been developed. This, in addition to taking people from Bradwell and seeding them into Magnox and the broader UK nuclear industry, should help to provide long term cost savings and delivery efficiencies across the entire NDA estate. This paper describes the journey of Bradwell Site through decommissioning and into C and M, highlighting the innovation and challenges that have contributed to the success of delivering the first UK nuclear site into Care and Maintenance. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2016

A kinematically complete, interdisciplinary, and co-institutional measurement of the 19 F(α,n) cross section for nuclear safeguards science

Alpha particles emitted from the decay of uranium in a UF 6 matrix can interact with fluorine and generate neutrons via the 19 F(α,n) 22 Na reaction. These neutrons can be used to determine the uranium content in a UF 6 storage cylinder. The accuracy of this self-interrogating, non-destructive assay (NDA) technique is, however, limited by the uncertainty of the 19 F(α,n) 22 Na cross section. We have performed complementary measurements of the 19 F(α,n) 22 Na reaction with both 4 He and 19 F beams to improve the precision of the 19 F(α,n) 22 Na cross section over the alpha energy range that encompasses common actinide alpha decay needed for NDA studies. We have determined an absolute cross section for the 19 F(α,n) 22 Na reaction to an average precision of 7.6% over the alpha energy range of 3.9 – 6.7 MeV. We utilized this cross section in a simulation of a 100 g spherical UF 6 assembly and obtained a change in neutron emission rate values of approximately 10-12%, and a significant (factor of 3.6) decrease in the neutron emission rate uncertainty (from 50-51% to 13-14%), compared to simulations using the old cross section. Our new absolute cross section enables improved interpretations of NDAs of containers of arbitrary size and configuration.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Ground Motion Inputs for the Seismic Shake Table Test

Currently, spent nuclear fuel (SNF) is stored in on-site independent spent-fuel storage installations (ISFSIs) at seventythree (73) nuclear power plants (NPPs) in the US. Because a site for geologic repository for permanent disposal of SNF has not been constructed, the SNF will remain in dry storage significantly longer than planned. During this time, the ISFSIs, and potentially consolidated storage facilities, will experience earthquakes of different magnitudes. The dry storage systems are designed and licensed to withstand large seismic loads. When dry storage systems experience seismic loads, there are little data on the response of SNF assemblies contained within them. The Spent Fuel Waste Disposition (SFWD) program is planning to conduct a full-scale seismic shake table test to close the gap related to the seismic loads on the fuel assemblies in dry storage systems. This test will allow for quantifying the strains and accelerations on surrogate fuel assembly hardware and cladding during earthquakes of different magnitudes and frequency content. The main component of the test unit will be the full-scale NUHOMS 32 PTH2 dry storage canister. The canister will be loaded with three surrogate fuel assemblies and twenty-nine dummy assemblies. Two dry storage configurations will be tested – horizontal and vertical above-ground concrete overpacks. These configurations cover 91% of the current dry storage configurations. The major input into the shake table test are the seismic excitations or the earthquake ground motions – acceleration time histories in two horizontal and one vertical direction that will be applied to the shake table surface during the tests. The shake table surface represents the top of the concrete pad on which a dry storage system is placed. The goal of the ground motion task is to develop the ground motions that would be representative of the range of seismotectonic and other conditions that any site in the Western US (WUS) or Central Eastern US (CEUS) might entail. This task is challenging because of the large number of the ISFSI sites, variety of seismotectonic and site conditions, and effects that soil amplification, soil-structure interaction, and pad flexibility may have on the ground motions.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2015

Dark matter directional detection in non-relativistic effective theories

We extend the formalism of dark matter directional detection to arbitrary one-body dark matter-nucleon interactions. The new theoretical framework generalizes the one currently used, which is based on 2 types of dark matter-nucleon interaction only. It includes 14 dark matter-nucleon interaction operators, 8 isotope-dependent nuclear response functions, and the Radon transform of the first 2 moments of the dark matter velocity distribution. We calculate the recoil energy spectra at dark matter directional detectors made of CF{sub 4}, CS{sub 2} and {sup 3}He for the 14 dark matter-nucleon interactions, using nuclear response functions recently obtained through numerical nuclear structure calculations. We highlight the new features of the proposed theoretical framework, and present our results for a spherical dark matter halo and for a stream of dark matter particles. This study lays the foundations for model independent analyses of dark matter directional detection experiments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Ion complexation waves emerge at the curved interfaces of layered minerals

Abstract Visualizing hydrated interfaces is of widespread interest across the physical sciences and is a particularly acute need for layered minerals, whose properties are governed by the structure of the electric double layer (EDL) where mineral and solution meet. Here, we show that cryo electron microscopy and tomography enable direct imaging of the EDL at montmorillonite interfaces in monovalent electrolytes with ångstrom resolution over micron length scales. A learning-based multiple-scattering reconstruction method for cryo electron tomography reveals ions bound asymmetrically on opposite sides of curved, exfoliated layers. We observe conserved ion-density asymmetry across stacks of interacting layers in cryo electron microscopy that is associated with configurations of inner- and outer-sphere ion-water-mineral complexes that we term complexation waves. Coherent X-ray scattering confirms that complexation waves propagate at room-temperature via a competition between ion dehydration and charge interactions that are coupled across opposing sides of a layer, driving dynamic transitions between stacked and aggregated states via layer exfoliation.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Close-out report with links to abstracts

This grant provided A/V support for two technical meetings of the Edge Coordinating Committee: (1) Nov 13, 2013 (co-located with the APS-DPP meeting in Denver, CO) https://ecc.mit.edu/fall-2013-technical-meeting#overlay-context=ecc-meetings; (2) April 28-May 1, 2015 (embedded sessions in the Transport Task Force Meeting, Salem, MA) http://www-internal.psfc.mit.edu/TTF2015/index.html. The ultimate goal of the U.S. Transport Task Force is to develop a physics-based understanding of particle, momentum and heat transport in magnetic fusion devices. This understanding should be of sufficient depth that it allows the development of predictive models of plasma transport that can be validated against experiment, and then used to anticipate the future performance of burning plasmas in ITER, as well as to provide guidance for the design of next-step fusion nuclear science facilities. To achieve success in transport science, it is essential to characterize local fluctuations and transport in toroidal plasmas, to understand the basic mechanisms responsible for transport, and ultimately to control these transport processes. These goals must be pursued in multiple areas, and these topics evolve in order to reflect current interests.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

From asymmetric to symmetric fission in the fermium isotopes within the time-dependent generator-coordinate-method formalism

Background: Predicting the properties of neutron-rich nuclei far from the valley of stability is one of the major challenges of modern nuclear theory. In heavy and superheavy nuclei, a difference of only a few neutrons is sufficient to change the dominant fission mode. A theoretical approach capable of predicting such rapid transitions for neutron-rich systems would be a valuable tool to better understand $r$-process nucleosynthesis or the decay of superheavy elements. Purpose: In this work, we investigate for the first time the transition from asymmetric to symmetric fission through the calculation of primary fission yields with the time-dependent generator coordinate method (TDGCM). We choose here the transition in neutron-rich fermium isotopes, which was the first to be observed experimentally in the late 1970s and is often used as a benchmark for theoretical studies. Methods: We compute the primary fission fragment mass and charge yields for 254 Fm, 256 Fm, and 258 Fm from the TDGCM under the Gaussian overlap approximation. The static part of the calculation (generation of a potential energy surface) consists of a series of constrained Hartree-Fock-Bogoliubov calculations based on the D1S, D1M, or D1N parametrization of the Gogny effective interaction in a two-center harmonic oscillator basis. The two-dimensional dynamics in the collective space spanned by the quadrupole and octupole moments ( $\hat{Q}_{20}$, $\hat{Q}_{30}$) is then computed with the finite element solver Felix-2.0. Results: The available experimental data and the TDGCM post-dictions are consistent and agree especially on the position in the fermium isotopic chain at which the transition occurs. In addition, the TDGCM predicts two distinct asymmetric modes for the fission of 254 Fm. Conclusions: Thanks to its intrinsic accounting of shell effects and to its ability to describe the dynamics of the system up to configurations close to scission, the TDGCM is able to describe qualitatively the fission yield transition in the neutron-rich fermium isotopes. This makes it a promising tool to study the evolution of the fission yields far from the valley of stability. The main limitation of the method lies in the presence of discontinuities in the two-dimensional manifold of generator states.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2014

Microstructure from joint analysis of experimental data and ab initio interactions: Hydrogenated amorphous silicon

A study of the formation of voids and molecular hydrogen in hydrogenated amorphous silicon is presented based upon a hybrid approach that involves inversion of experimental nuclear magnetic resonance data in conjunction with ab initio total-energy relaxations in an augmented solution space. The novelty of this approach is that the voids and molecular hydrogen appear naturally in the model networks unlike conventional approaches, where voids are created artificially by removing silicon atoms from the networks. Two representative models with 16 and 18 at. % of hydrogen are studied in this work. The result shows that the microstructure of the a-Si:H network consists of several microvoids and few molecular hydrogen for concentration above 15 at. % H. The microvoids are highly irregular in shape and size, and have a linear dimension of 5–7 Å. The internal surface of a microvoid is found to be decorated with 4–9 hydrogen atoms in the form of monohydride Si–H configurations as observed in nuclear magnetic resonance experiments. The microstructure consists of (0.9–1.4)% hydrogen molecules of total hydrogen in the networks. These observations are consistent with the outcome of infrared spectroscopy, nuclear magnetic resonance, and calorimetry experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Octupole correlations near Te 110

The lifetime of the 2 + and 9 – , 11 – , 13 – , 15 – states in the neutron-deficient 110 Te was measured for the first time using the recoil distance Doppler shift technique. The reported value of the reduced transition probability B(E2; $0$$^{+}_{g.s}$ → 2 + ) = 4.3(8) × 10 3 e 2 fm 4 supports the systematic for even-mass Te isotopes and was interpreted in the framework of the large-scale shell model and cranked shell model calculations. The measured reduced transition probabilities in the negative-parity yrast band revealed the upward trend towards the high spins. Furthermore, the enhanced collectivity is discussed in terms of the tilted axis cranking approach and the symmetry configuration mixing method with the Gogny D1S interaction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Cross-shell excitations in Ca 46 studied with fusion reactions induced by a reaccelerated rare isotope beam

Discovering unexplored high-spin states in neutron-rich nuclei can open up a new direction to study band structure and the associated shell structure in isospin-asymmetric many-body systems. However, experimental reach has so far been limited to neutron-deficient or stable nuclei which are preferentially produced in fusion reactions used in such studies. Here in this paper, we report the first γ-ray spectroscopy with fusion reactions using a reaccelerated rare-isotope beam of 45 K performed at the ReA3 facility of the National Superconducting Cyclotron Laboratory. Using particle and γ-ray coincidence techniques, three new higher-lying states around 6 MeV and five new γ-ray transitions were identified for 46 Ca, suggesting three independent band structures formed from different particle-hole configurations. The rotational-like band built on the 0$^+_2$ state is established up to the tentatively assigned 6$^+_2$ state. New results are compared to large-scale shell model calculations, confirming the validity of the effective interaction describing particle-hole excitations across the Z=20 and N=28 shell gaps in the vicinity of doubly-magic 48 Ca.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2014

Advanced techniques for characterization of ion beam modified materials

Understanding the mechanisms of damage formation in materials irradiated with energetic ions is essential for the field of ion-beam materials modification and engineering. Utilizing incident ions, electrons, photons, and positrons, various analysis techniques, including Rutherford backscattering spectrometry (RBS), electron RBS, Raman spectroscopy, high-resolution X-ray diffraction, small-angle X-ray scattering, and positron annihilation spectroscopy, are routinely used or gaining increasing attention in characterizing ion beam modified materials. The distinctive information, recent developments, and some perspectives in these techniques are reviewed in this paper. Applications of these techniques are discussed to demonstrate their unique ability for studying ion-solid interactions and the corresponding radiation effects in modified depths ranging from a few nm to a few tens of μm, and to provide information on electronic and atomic structure of the materials, defect configuration and concentration, as well as phase stability, amorphization and recrystallization processes. Finally, such knowledge contributes to our fundamental understanding over a wide range of extreme conditions essential for enhancing material performance and also for design and synthesis of new materials to address a broad variety of future energy applications.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Hydration structure and energetics of iron-exchanged montmorillonite

The interaction between iron released from corroded steel canisters and bentonite is a key process influencing the long-term performance of nuclear waste repositories. In particular, the migration of Fe²⁺ into montmorillonite (Mnt) interlayers may alter its hydration, swelling, and ion-transport properties. In the present study, molecular dynamics simulations were performed to investigate the hydration behavior, structural response, and transport properties of Fe-exchanged montmorillonite (Fe-Mnt) under varying hydration states. The simulations focus on short- to intermediate-time-scale Fe 2+ and Fe 3+ interlayer exchange and hydration effects, and do not consider long-term structural substitution, Fe-bearing clay phase stabilization, or secondary iron mineral precipitation. Systems containing Na + -, Fe 2+ -, and Fe 3+ -Mnt were examined using both periodic and edge-exposed configurations to evaluate interlayer structure, ion exchange, and free energy of Fe intercalation. The results show that Fe ions influence the interlayer spacing primarily at low water contents (<1 bilayer), where Fe-Mnt exhibits a d-spacing 1–2 Å larger than Na-Mnt due to stronger hydration. The calculated hydration energies follow the order Fe 2+ $<$Fe 3+ $<$Na + . Both water and ion diffusion coefficients decrease upon Fe ion intercalation, with Fe 2+ ions diffusing an order of magnitude more slowly than those of Na + . Free energy profiles further confirm that Fe 2+ and Fe 3+ ions are thermodynamically favored in the interlayer, with Fe 3+ being the most stable. In conclusion, these findings provide molecular-scale insights into the mechanisms of Fe–Na exchange and their implications for bentonite alteration in repository environments.

clay↗