Glycoboehmite (GB) materials are synthesized by a solvothermal reaction to form layered aluminum oxyhydroxide (boehmite) modified by intercalated butanediol molecules. These hybrid materials offer a platform to design materials with potentially novel sorption, wetting, and catalytic properties. Several synthetic methods have been used, resulting in different structural and spectroscopic properties, but atomistic detail is needed to determine the interlayer structure to explore the synthetic control of GB materials. Here, we use classical molecular dynamics (MD) simulations to compare the structural properties of GB interlayers containing chemisorbed butanediol molecules as a function of diol loading. Accompanying quantum (density functional theory, DFT) static calculations and MD simulations are used to validate the classical model and compute the infrared spectra of various models. Classical MD results reveal the existence of two unique interlayer environments at higher butanediol loading, corresponding to smaller (cross-linked) and expanded interlayers. DFT-computed infrared spectra reveal the sensitivity of the aluminol O-H stretch frequencies to the interlayer environment, consistent with the spectrum of the synthesized material. Insight from these simulations will aid in the characterization of the newly synthesized GB materials.
Jove Colon, Carlos F.; Ho, Tuan A.; Lopez, Carlos M.; Rutqvist, Jonny; Guglielmi, Yves; Hu, Mengsu; Sasaki, Tsubasa; Yoon, Sangcheol; Steefel, Carl I.; Tournassat, Christophe; Mital, Utkarsh; Luu, Keurfon; Sauer, Kirsten B.; Caporuscio, Florie A.; Rock, Marlena J.; Zandanel, Amber E.; Zavarin, Mavrik; Wolery, Thomas J.; Chang, Elliot; Han, Sol-Chan; Wainwright, Haruko; Greathouse, Jeffery A.
This report represents the milestone deliverable M2SF-23SN010301072 “Evaluation of Nuclear Spent Fuel Disposal in Clay-Bearing Rock - Process Model Development and Experimental Studies” The report provides a status update of FY23 activities for the work package Argillite Disposal work packages for the DOE-NE Spent Fuel Waste Form Science and Technology (SFWST) Program. Clay-rich geological media (often referred as shale or argillite) are among the most abundant type of sedimentary rock near the Earth’s surface. Argillaceous rock formations have the following advantageous attributes for deep geological nuclear waste disposal: widespread geologic occurrence, found in stable geologic settings, low permeability, self-sealing properties, low effective diffusion coefficient, high sorption capacity, and have the appropriate depth and thickness to host nuclear waste repository concepts. The DOE R&D program under the Spent Fuel Waste Science Technology (SFWST) campaign has made key progress (through experiment, modeling, and testing) in the study of chemical and physical phenomena that could impact the long-term safety assessment of heat-generating nuclear waste disposition in clay/shale/argillaceous rock. International collaboration activities comprising field-scale heater tests, field data monitoring, and laboratory-scale experiments provide key information on changes to the engineered barrier system (EBS) material exposed high thermal loads. Moreover, consideration of direct disposal of large capacity dual-purpose canisters (DPCs) as part of the back-end SNF waste disposition strategy has generated interest in improving our understanding of the effects of elevated temperatures on the engineered barrier system (EBS) design concepts. Chemical and structural analyses of sampled bentonite material from laboratory tests at elevated temperatures are key to the characterization of thermal effects affecting bentonite clay barrier performance. The knowledge provided by these experiments is crucial to constrain the extent of sacrificial zones in the EBS design during the thermal period. Thermal, hydrologic, mechanical, and chemical (THMC) data collected from heater tests and laboratory experiments have been used in the development, validation, and calibration of THMC simulators to model near-field coupled processes. This information leads to the development of simulation approaches to assess issues on coupled processes involving porous media flow, transport, geomechanical phenomena, chemical interactions with barrier/geologic materials, and the development of EBS concepts. These lines of knowledge are central to the design of deep geological backfilled repository concepts where temperature plays a key role in the EBS behavior, potential interactions with host rock, and long-term performance in the safety assessment.
Senanayake, Hasini S.; Wimalasiri, Pubudu N.; Godahewa, Sahan M.; Thompson, Ward H.; Greathouse, Jeffery A.
Here, we present a classical interatomic force field, silica-DDEC, to describe the interactions of amorphous and crystalline silica surfaces, parametrized using density functional theory-based charges. Charge schemes for silica surfaces were developed using the density-derived electrostatic and chemical (DDEC) method, which reproduces atomic charges of the periodic models as well as the electrostatic potential away from the atom sites. Lennard–Jones parameters were determined by requiring the correct description of (i) the amorphous silica density, coordination defects, and local coordination geometry, relative to experimental measurements, and (ii) water-silica interatomic distances compared with ab initio results. Deprotonated surface silanol sites are also described within the model based on DDEC charges. The result is a general electronic structure-derived model for describing fully flexible amorphous and crystalline silica surfaces and interactions of liquids with silica surfaces of varying structure and protonation state.
This report describes research and development (R&D) activities conducted during Fiscal Year 2023 (FY23) in the Advanced Fuels and Advanced Reactor Waste Streams Strategies work package in the Spent Fuel Waste Science and Technology (SFWST) Campaign supported by the United States (U.S.) Department of Energy (DOE). This report is focused on evaluating and cataloguing Advanced Reactor Spent Nuclear Fuel (AR SNF) and Advanced Reactor Waste Streams (ARWS) and creating Back-end Nuclear Fuel Cycle (BENFC) strategies for their disposition. The R&D team for this report is comprised of researchers from Sandia National Laboratories and Enviro Nuclear Services, LLC.
Diffusion properties of bulk fluids have been predicted using empirical expressions and machine learning (ML) models, suggesting that predictions of diffusion also should be possible for fluids in confined environments. The ability to quickly and accurately predict diffusion in porous materials would enable new discoveries and spur development in relevant technologies such as separations, catalysis, batteries, and subsurface applications. Here in this work, we apply artificial neural network (ANN) models to predict the simulated self-diffusion coefficients of real liquids in both bulk and pore environments. The training data sets were generated from molecular dynamics (MD) simulations of Lennard-Jones particles representing a diverse set of 14 molecules ranging from ammonia to dodecane over a range of liquid pressures and temperatures. Planar, cylindrical, and hexagonal pore models consisted of walls composed of carbon atoms. Our simple model for these liquids was primarily used to generate ANN training data, but the simulated self-diffusion coefficients of bulk liquids show excellent agreement with experimental diffusion coefficients. ANN models based on simple descriptors accurately reproduced the MD diffusion data for both bulk and confined liquids, including the trend of increased mobility in large pores relative to the corresponding bulk liquid.
Strong gas-mineral interactions or slow adsorption kinetics require a molecular-level understanding of both adsorption and diffusion for these interactions to be properly described in transport models. In this combined molecular simulation and experimental study, noble gas adsorption and mobility is investigated in two naturally abundant zeolites whose pores are similar in size (clinoptilolite) and greater than (mordenite) the gas diameters. Simulated adsorption isotherms obtained from grand canonical Monte Carlo simulations indicate that both zeolites can accommodate even the largest gas (Rn). However, gas mobility in clinoptilolite is significantly hindered at pore-limiting window sites, as seen from molecular dynamics simulations in both bulk and slab zeolite models. Experimental gas adsorption isotherms for clinoptilolite confirm the presence of a kinetic barrier to Xe uptake, resulting in the unusual property of reverse Kr/Xe selectivity. Finally, a kinetic model is used to fit the simulated gas loading profiles, allowing a comparison of trends in gas diffusivity in the zeolite pores.
Understanding the adsorption of isolated metal cations from water on to mineral surfaces is critical for toxic waste retention and cleanup in the environment. Heterogeneous nucleation of metal oxyhydroxides and other minerals on material surfaces is key to crystal growth and dissolution. The link connecting these two areas, namely cation dimerization and polymerization, is far less understood. In this work we apply ab initio molecular dynamics calculations to examine the coordination structure of hydroxide-bridged Cu(II) dimers, and the free energy changes associated with Cu(II) dimerization on silica surfaces. The dimer dissociation pathway involves sequential breaking of two Cu2+-OH− bonds, yielding three local minima in the free energy profiles associated with 0-2 OH− bridges between the metal cations, and requires the design of a (to our knowledge) novel reaction coordinate for the simulations. Cu(II) adsorbed on silica surfaces are found to exhibit stronger tendency towards dimerization than when residing in water. Cluster-plus-implicit-solvent methods yield incorrect trends if OH− hydration is not correctly depicted. The predicted free energy landscapes are consistent with fast equilibrium times (seconds) among adsorbed structures, and favor Cu2+ dimer formation on silica surfaces over monomer adsorption.
Tracer gases, whether they are chemical or isotopic in nature, are useful tools in examining the flow and transport of gaseous or volatile species in the underground. One application is using detection of short-lived argon and xenon radionuclides to monitor for underground nuclear explosions. However, even chemically inert species, such as the noble gases, have bene observed to exhibit non-conservative behavior when flowing through porous media containing certain materials, such as zeolites, due to gas adsorption processes. This report details the model developed, implemented, and tested in the open source and massively parallel subsurface flow and transport simulator PFLOTRAN for future use in modeling the transport of adsorbing tracer gases.
The structural and dynamical properties of nanoconfined solutions can differ dramatically from those of the corresponding bulk systems. Understanding the changes induced by confinement is central to controlling the behavior of synthetic nanostructured materials and predicting the characteristics of biological and geochemical systems. A key outstanding issue is how the molecular-level behavior of nanoconfined electrolyte solutions is reflected in different experimental, particularly spectroscopic, measurements. This is addressed here through molecular dynamics simulations of the OH stretching infrared (IR) spectroscopy of NaCl, NaBr, and NaI solutions in isotopically dilute HOD/D2O confined in hydroxylated amorphous silica slit pores of width 1-6 nm and pH ∼2. In addition, the water reorientation dynamics and spectral diffusion, accessible by pump-probe anisotropy and two-dimensional IR measurements, are investigated. The aim is to elucidate the effect of salt identity, confinement, and salt concentration on the vibrational spectra. It is found that the IR spectra of the electrolyte solutions are only modestly blue-shifted upon confinement in amorphous silica slit pores, with both the size of the shift and linewidth increasing with the halide size, but these effects are suppressed as the salt concentration is increased. This indicates the limitations of linear IR spectroscopy as a probe of confined water. However, the OH reorientational and spectral diffusion dynamics are significantly slowed by confinement even at the lowest concentrations. The retardation of the dynamics eases with increasing salt concentration and pore width, but it exhibits a more complex behavior as a function of halide.
Structural properties of the anionic surfactant dioctyl sodium sulfosuccinate (AOT or Aerosol-OT) adsorbed on the mica surface were investigated by molecular dynamics simulation, including the effect of surface loading in the presence of monovalent and divalent cations. The simulations confirmed recent neutron reflectivity experiments that revealed the binding of anionic surfactant to the negatively charged surface via adsorbed cations. At low loading, cylindrical micelles formed on the surface, with sulfate head groups bound to the surface by water molecules or adsorbed cations. Cation bridging was observed in the presence of weakly hydrating monovalent cations, while sulfate groups interacted with strongly hydrating divalent cations through water bridges. The adsorbed micelle structure was confirmed experimentally with cryogenic electronic microscopy, which revealed micelles approximately 2 nm in diameter at the basal surface. At higher AOT loading, the simulations reveal adsorbed bilayers with similar surface binding mechanisms. Adsorbed micelles were slightly thicker (2.2–3.0 nm) than the corresponding bilayers (2.0–2.4 nm). Upon heating the low loading systems from 300 K to 350 K, the adsorbed micelles transformed to a more planar configuration resembling bilayers. The driving force for this transition is an increase in the number of sulfate head groups interacting directly with adsorbed cations.
The ability to predict transport properties of liquids quickly and accurately will greatly improve our understanding of fluid properties both in bulk and complex mixtures, as well as in confined environments. Such information could then be used in the design of materials and processes for applications ranging from energy production and storage to manufacturing processes. As a first step, we consider the use of machine learning (ML) methods to predict the diffusion properties of pure liquids. Recent results have shown that Artificial Neural Networks (ANNs) can effectively predict the diffusion of pure compounds based on the use of experimental properties as the model inputs. In the current study, a similar ANN approach is applied to modeling diffusion of pure liquids using fluid properties obtained exclusively from molecular simulations. A diverse set of 102 pure liquids is considered, ranging from small polar molecules (e.g., water) to large nonpolar molecules (e.g., octane). Self-diffusion coefficients were obtained from classical molecular dynamics (MD) simulations. Since nearly all the molecules are organic compounds, a general set of force field parameters for organic molecules was used. The MD methods are validated by comparing physical and thermodynamic properties with experiment. Computational input features for the ANN include physical properties obtained from the MD simulations as well as molecular properties from quantum calculations of individual molecules. Fluid properties describing the local liquid structure were obtained from center of mass radial distribution functions (COM-RDFs). Feature sensitivity analysis revealed that isothermal compressibility, heat of vaporization, and the thermal expansion coefficient were the most impactful properties used as input for the ANN model to predict the MD simulated self-diffusion coefficients. The MD-based ANN successfully predicts the MD self-diffusion coefficients with only a subset (2 to 3) of the available computationally determined input features required. A separate ANN model was developed using literature experimental self-diffusion coefficients as model targets. Although this second ML model was not as successful due to a limited number of data points, a good correlation is still observed between experimental and ML predicted self-diffusion coefficients.
Predicting the diffusion coefficient of fluids under nanoconfinement is important for many applications including the extraction of shale gas from kerogen and product turnover in porous catalysts. Due to the large number of important variables, including pore shape and size, fluid temperature and density, and the fluid-wall interaction strength, simulating diffusion coefficients using molecular dynamics (MD) in a systematic study could prove to be prohibitively expensive. Here, we use machine learning models trained on a subset of MD data to predict the self-diffusion coefficients of Lennard-Jones fluids in pores. Our MD data set contains 2280 simulations of ideal slit pore, cylindrical pore, and hexagonal pore geometries. We use the forward feature selection method to determine the most useful features (i.e., descriptors) for developing an artificial neutral network (ANN) model with an emphasis on easily acquired features. Our model shows good predictive ability with a coefficient of determination (i.e., R2) of ∼0.99 and a mean squared error of ∼2.9 × 10-5. Finally, we propose an alteration to our feature set that will allow the ANN model to be applied to nonideal pore geometries.
The objective of this project was to eliminate and/or render bulk agent unusable by a threat entity via neutralization and/or polymerization of the bulk agent using minimal quantities of additives. We proposed the in situ neutralization and polymerization of bulk chemical agents (CAs) by performing reactions in the existing CA storage container via wet chemical approaches using minimal quantities of chemical based materials. This approach does not require sophisticated equipment, fuel to power generators, electricity to power equipment, or large quantities of decontaminating materials. By utilizing the CA storage container as the batch reactor, the amount of logistical resources can be significantly reduced. Fewer personnel are required since no sophisticated equipment needs to be set up, configured, or operated. Employing the CA storage container as the batch reactor enables the capability to add materials to multiple containers in a short period of time as opposed to processing one container at a time for typical batch reactor approaches. In scenarios where a quick response is required, the material can be added to all the CA containers and left to react on its own without intervention. Any attempt to filter the CA plus material solution will increase the rate of reaction due to increased agitation of the solution.
As a general-purpose force field for molecular simulations of layered materials and their fluid interfaces, Clayff continues to see broad usage in atomistic computational modeling for numerous geoscience and materials science applications due to its (1) success in predicting properties of bulk nanoporous materials and their interfaces, (2) transferability to a range of layered and nanoporous materials, and (3) simple functional form which facilitates incorporation into a variety of simulation codes. Here, we review applications of Clayff to model bulk phases and interfaces not included in the original parameter set and recent modifications for modeling surface terminations such as hydroxylated nanoparticle edges. We conclude with a discussion of expectations for future developments.
We present a combined molecular dynamics (MD) simulation and X-ray absorption fine structure (XAFS) spectroscopic investigation of aqueous iron adsorption on nanoconfined amorphous silica surfaces. The simulation models examine the effects of pore size, pH (surface charge), iron valency, and counter-ion (chloride or hydroxide). The simulation methods were validated by comparing the coordination environment of adsorbed iron with coordination numbers and bond lengths derived from XAFS. In the MD models, nanoconfinement effects on local iron coordination were investigated by comparing results for unconfined silica surfaces and in confined domains within 2 nm, 4 nm, and 8 nm pores. Experimentally, coordination environments of iron adsorbed onto mesoporous silica with 4 nm and 8 nm pores at pH 7.5 were investigated. The effect of pH in the MD models was included by simulating Fe(ii) adsorption onto negatively charged SiO2surfaces and Fe(iii) adsorption on neutral surfaces. The simulation results show that iron adsorption depends significantly on silica surface charge, as expected based on electrostatic interactions. Adsorption on a negatively charged surface is an order of magnitude greater than on the neutral surface, and simulated surface coverages are consistent with experimental results. Pore size effects from the MD simulations were most notable in the adsorption of Fe(ii) at deprotonated surface sites (SiO−), but adsorption trends varied with concentration and aqueous Fe speciation. The coordination environment of adsorbed iron varied significantly with the type of anion. Considerable ion pairing with hydroxide anions led to the formation of oligomeric surface complexes and aqueous species, resulting in larger iron hydroxide clusters at higher surface loadings.
Water in nano-scale confining environments is a key element in many biological, material, and geological systems. The structure and dynamics of the liquid can be dramatically modified under these conditions. Probing these changes can be challenging, but vibrational spectroscopy has emerged as a powerful tool for investigating their behavior. A critical, evolving component of this approachis a detailed understanding of the connection between spectroscopic features and molecular-level details. In this paper, this issue is addressed by using molecular dynamics simulations to simulate the linear infrared (IR) and Raman spectra for isotopically dilute HOD in D2O confined inhydroxylated amorphous silica slit pores. The effect of slit-pore width and hydroxyl density on thesilica surface on the vibrational spectra is also investigated. The primary effect of confinement is a blueshift in the frequency of OH groups donating a hydrogen bond to the silica surface. Thisappears as a slight shift in the total (measurable) spectra but is clearly seen in the distance-based IR and Raman spectra. Analysis indicates that these changes upon confinement are associated withtheweaker hydrogen-bond accepting properties of silica oxygens compared to water molecules.
This report summarizes the results obtained during the LDRD project entitled "Partitioning of Complex Fluids at Mineral Interfaces." This research addressed fundamental aspects of such interfaces, which are relevant to energy-water applications in the subsurface, including fossil energy extraction and carbon sequestration. This project directly addresses the problem of selectivity of complex fluid components at mineral-fluid interfaces, where complex fluids are defined as a mixture of hydrophobic and hydrophilic components: e.g., water, aqueous ions, polar/nonpolar organic compounds. Specifically, this project investigates how adsorption selectivity varies with surface properties and fluid composition. Both experimental and molecular modeling techniques were used to better understand trends in surface wettability on mineral surfaces. The experimental techniques spanned the macroscale (contact angle measurements) to the nanoscale (cryogenic electronic microscopy and vibrational spectroscopy). We focused on an anionic surfactant and a well-characterized mineral phase representative of clay phases present in oil- and gas-producing shale deposits. Collectively, the results consistently demonstrate that the presence of surfactant in the aqueous fluid significantly affects the mineral-fluid interfacial structure. Experimental and molecular modeling results reveal details of the surfactant structure at the interface, and how this structure varies with surfactant coverage and fluid composition.
Gonzalez, Rafael I.; Rojas-Nunez, Javier; Valencia, Felipe J.; Munoz, Francisco; Baltazar, Samuel E.; Allende, Sebastian; Rogan, Jose; Valdivia, Juan A.; Kiwi, Miguel; Ramírez, Ricardo; Greathouse, Jeffery A.
Imogolite is a fascinating inorganic nanotube that is found in nature or synthesized in a laboratory. The synthesis process is carried out in liquid media, and leads to the formation of almost monodisperse diameter nanotubes. Here we investigate, employing classical molecular dynamics simulations, the interaction of water and imogolite for nanotubes of several radii. We established that water penetrates the pores of N = 9 and larger nanotubes, and adopts a coaxial arrangement in it. Also, while water molecules can diffuse along the center of the nanotube, the molecules next to the inner imogolite walls have very low mobility. At the outer nanotube wall, an increase of water density is observed, this effect extends up to 1 nm, beyond which water properties are bulk-like. Both phenomena are affected by the imogolite curvature.
Classical molecular dynamics (MD) simulations were performed to provide a conceptual understanding of the amorphous-crystalline interface for a candidate negative thermal expansion (NTE) material, ZrW2O8. Simulations of pressure-induced amorphization at 300 K indicate that an amorphous phase forms at pressures of 10 GPa and greater, and this phase persists when the pressure is subsequently decreased to 1 bar. However, the crystalline phase is recovered when the slightly distorted 5 GPa phase is relaxed to 1 bar. Simulations were also performed on a two-phase model consisting of the high-pressure amorphous phase in direct contact with the crystalline phase. Upon equilibration at 300 K and 1 bar, the crystalline phase remains unchanged beyond a thin layer of disrupted structure at the crystalline-amorphous interface. Differences in local atomic structure at the interface are quantified from the simulation trajectories.
Here we report molecular level details regarding the adsorption of sarin (GB) gas in a prototypical zirconium-based metal-organic framework (MOF, UiO-66). By combining predictive modeling and experimental spectroscopic techniques, we unambiguously identify several unique bindings sites within the MOF, using the P=O stretch frequency of GB as a probe. Remarkable agreement between predicted and experimental IR spectrum is demonstrated. As previously hypothesized, the undercoordinated Lewis acid metal site is the most favorable binding site. Yet multiple sites participate in the adsorption process; specifically, the Zr-chelated hydroxyl groups form hydrogen bonds with the GB molecule, and GB weakly interacts with fully coordinated metals. Importantly, this work highlights that subtle orientational effects of bound GB are observable via shifts in characteristic vibrational modes; this finding has large implications for degradation rates and opens a new route for future materials design.
This report describes the potential of a novel class of materials—α-ZrW2O8, Zr2WP2O12, and related compounds that contract upon amorphization as possible radionuclide waste-forms. The proposed ceramic waste-forms would consist of zoned grains, or sintered ceramics with center- loaded radionuclides and barren shells. Radiation-induced amorphization would result in core shrinkage but would not fracture the shells or overgrowths, maintaining isolation of the radionuclide. In this report, we have described synthesis techniques to produce phase-pure forms of the materials, and how to fully densify those materials. Structural models for the materials were developed and validated using DFPT approaches, and radionuclide substitution was evaluated; U(IV), Pu(IV), Tc(IV) and Tc(VII) all readily substitute into the material structures. MD modeling indicated that strain associated with radiation-induced amorphization would not affect the integrity of surrounding crystalline materials, and these results were validated via ion beam experimental studies. Finally, we have evaluated the leach rates of the barren materials, as determined by batch and flow-through reactor experiments. ZrW2O8 leaches rapidly, releasing tungstate while Zr is retained as a solid oxide or hydroxide. Tungsten release rates remain elevated over time and are highly sensitive to contact times, suggesting that this material will not be an effective waste-form. Conversely, tungsten releases rates from Zr2WP2O12 rapidly drop, show little dependence on short-term changes in fluid contact time, and in over time, become tied to P release rates. The results presented here suggest that this material may be a viable waste-form for some hard-to-handle radionuclides such as Pu and Tc.
A molecular-scale understanding of the transition between hydration states in clay minerals remains a challenging problem because of the very fast stepwise swelling process observed from X-ray diffraction (XRD) experiments. XRD profile modeling assumes the coexistence of multiple hydration states in a clay sample to fit the experimental XRD pattern obtained under humid conditions. While XRD profile modeling provides a macroscopic understanding of the heterogeneous hydration structure of clay minerals, a microscopic model of the transition between hydration states is still missing. Here, for the first time, we use molecular dynamics simulation to investigate the transition states between a dry interlayer, one-layer hydrate, and two-layer hydrate. We find that the hydrogen bonds that form across the interlayer at the clay particle edge make an important contribution to the energy barrier to interlayer hydration, especially for initial hydration.
The phonon, infrared, and Raman spectroscopic properties of zirconium tungsten phosphate, Zr2(WO4)(PO4)2 (space group Pbcn, IT No. 60; Z = 4), have been extensively investigated using density functional perturbation theory (DFPT) calculations with the Perdew, Burke, and Ernzerhof exchange-correlation functional revised for solids (PBEsol) and validated by experimental characterization of Zr2(WO4)(PO4)2 prepared by hydrothermal synthesis. Using DFPT-simulated infrared, Raman, and phonon density-of-state spectra combined with Fourier transform infrared and Raman measurements, new comprehensive and extensive assignments have been made for the spectra of Zr2(WO4)(PO4)2, resulting in the characterization of its 29 and 34 most intense IR- and Raman-active modes, respectively. DFPT results also reveal that ν1(PO4) symmetric stretching and ν3(PO4) antisymmetric stretching bands have been interchanged in previous Raman experimental assignments. Negative thermal expansion in Zr2(WO4)(PO4)2 appears to have very limited impact on the spectral properties of this compound. This work shows the high accuracy of the PBEsol exchange-correlation functional for studying the spectroscopic properties of crystalline materials using first-principles methods.
We have investigated cubic zirconium tungstate (ZrW2O8) using density functional perturbation theory (DFPT), along with experimental characterization to assess and validate computational results. Cubic zirconium tungstate is among the few known materials exhibiting isotropic negative thermal expansion (NTE) over a broad temperature range, including room temperature where it occurs metastably. Isotropic NTE materials are important for technological applications requiring thermal-expansion compensators in composites designed to have overall zero or adjustable thermal expansion. While cubic zirconium tungstate has attracted considerable attention experimentally, a very few computational studies have been dedicated to this well-known NTE material. Therefore, spectroscopic, mechanical and thermodynamic properties have been derived from DFPT calculations. A systematic comparison of the calculated infrared, Raman, and phonon density-of-state spectra has been made with Fourier transform far-/mid-infrared and Raman data collected in this study, as well as with available inelastic neutron scattering measurements. The thermal evolution of the lattice parameter computed within the quasi-harmonic approximation exhibits negative values below the Debye temperature, consistent with the observed negative thermal expansion characteristics of cubic zirconium tungstate, α-ZrW2O8. These results show that this DFPT approach can be used for studying the spectroscopic, mechanical and thermodynamic properties of prospective NTE ceramic waste forms for encapsulation of radionuclides produced during the nuclear fuel cycle.
We use molecular simulations to provide a conceptual understanding of a crystalline-amorphous interface for a candidate negative thermal expansion (NTE) material. Specifically, classical molecular dynamics (MD) simulations were used to investigate the temperature and pressure dependence on structural properties of ZrW2O8. Polarizability of oxygen atoms was included to better account for the electronic charge distribution within the lattice. Constant-pressure simulations of cubic crystalline ZrW2O8 at ambient pressure reveal a slight NTE behavior, characterized by a small structural rearrangement resulting in oxygen sharing between adjacent WO4 tetrahedra. Periodic quantum calculations confirm that the MD-optimized structure is lower in energy than the idealized structure obtained from neutron diffraction experiments. Additionally, simulations of pressure-induced amorphization of ZrW2O8 at 300 K indicate that an amorphous phase forms at pressures greater than 10 GPa, and this phase persists when the pressure is decreased to 1 bar. Simulations were performed on a hybrid model consisting of amorphous ZrW2O8 in direct contact with the cubic crystalline phase. Upon equilibration at 300 K and 1 bar, the crystalline phase remains unchanged beyond a thin layer of disrupted structure at the amorphous interface. Detailed analysis reveals the transition in metal coordination at the interface.
Observation of vibrational properties of phyllosilicate edges via a combined molecular modeling and experimental approach was performed. Deuterium exchange was utilized to isolate edge vibrational modes from their internal counterparts. The appearance of a specific peak within the broader D2O band indicates the presence of deuteration on the edge surface, and this peak is confirmed with the simulated spectra. These results are the first to unambiguously identify spectroscopic features of phyllosilicate edge sites.
Appropriate waste-forms for radioactive materials must isolate the radionuclides from the environment for long time periods. To accomplish this typically requires low waste-form solubility, to minimize radionuclide release to the environment. However, radiation eventually damages most waste-forms, leading to expansion, crumbling, increased exposed surface area, and faster dissolution. We have evaluated the use of a novel class of materials-ZrW2O8, Zr2P2WO12 and related compounds-that contract upon amorphization. The proposed ceramic waste-forms would consist of zoned grains, or sintered ceramics with center-loaded radionuclides and barren shells. Radiation-induced amorphization would result in core shrinkage but would not fracture the shells or overgrowths, maintaining isolation of the radionuclide. We have synthesized these phases and have evaluated their leach rates. Tungsten forms stable aqueous species at neutral to basic conditions, making it a reliable indicator of phase dissolution. ZrW2O8 leaches rapidly, releasing tungstate while Zr is retained as a solid oxide or hydroxide. Tungsten release rates remain elevated over time and are highly sensitive to contact times, suggesting that this material will not be an effective waste-form. Conversely, tungsten release rates from Zr2P2WO12 rapidly drop and are tied to P release rates; we speculate that a low-solubility protective Zr-phosphate leach layer forms, slowing further dissolution.
The adsorption of chemical warfare agents and their simulants by Zr (UiO-66) and rare-earth (Y, UiO-66-DOBDC analog)-based metal-organic frameworks (MOFs) is explored here using density functional theory. In particular, we investigate the role of linker functional group (OH, H) and metal atom identity on the binding energies of organophosphorous compounds. Commonly used cluster approximations for MOF secondary building units and various optimization constraints are compared with three-dimensional periodic results. An in-depth scan of potential binding sites and orientations reveals little effect due to metal identity, whereas the effect of linker functionalization depends on the substrate. This finding strongly suggests that full linkers and functional groups should be included in cluster models. Importantly, defect sites show considerably improved binding of organophosphorous compounds as compared to ideal clusters. Favorable binding is also demonstrated at two additional adsorption sites, ZrOH and μ3-OH, that likely play a role in the initial adsorption process. The results presented here portray the importance of including full three-dimensional pore structures in the adsorption process of organophosphorous compounds in MOFs; a critical first step in the degradation of these harmful chemicals.
The negative thermal expansion (NTE) material Zr2(WO4)(PO4)2 has been investigated for the first time within the framework of the density functional perturbation theory (DFPT). The structural, mechanical, and thermodynamic properties of this material have been predicted using the Perdew, Burke and Ernzerhof for solid (PBEsol) exchange-correlation functional, which showed superior accuracy over standard functionals in previous computational studies of the NTE material α-ZrW2O8. The bulk modulus calculated for Zr2(WO4)(PO4)2 using the Vinet equation of state at room temperature is K0 = 63.6 GPa, which is in close agreement with the experimental estimate of 61.3(8) at T = 296 K. The computed mean linear coefficient of thermal expansion is -3.1 × 10-6 K-1 in the temperature range ∼0-70 K, in line with the X-ray diffraction measurements. The mean Grüneisen parameter controlling the thermal expansion of Zr2(WO4)(PO4)2 is negative below 205 K, with a minimum of -2.1 at 10 K. The calculated standard molar heat capacity and entropy are CP0 = 287.6 and S0 = 321.9 J·mol-1·K-1, respectively. The results reported in this study demonstrate the accuracy of DFPT/PBEsol for assessing or predicting the relationship between structural and thermomechanical properties of NTE materials.
Nanoporous materials such as metal-organic frameworks (MOFs) have attractive properties for selective capture of chemical warfare agents (CWAs). For obvious reasons, most research on adsorption of CWAs is performed with simulant molecules rather than real agents. This paper examines how effectively common CWA simulants mimic the adsorption properties of sarin and soman. To this end, we perform molecular simulations in the dilute adsorption limit for four simulants [dimethyl methylphosphonate (DMMP), diethyl chlorophosphate (DCP), diisopropyl fluorophosphate, and dimethyl p-nitrophenyl phosphate (DMNP)] and sarin and soman in a set of 2969 MOFs with experimentally known crystal structures. To establish the robustness of the conclusions with respect to the force field used in these simulations, each system was examined with two independent force fields, a "generic" force field and a density functional theory (DFT)-derived force field we established based on extensive dispersion-corrected DFT calculations of adsorption in the well-known MOF UiO-66. Our results show that when judging the performance of adsorbents using the heat of adsorption, DCP and DMMP are the best simulants for the adsorption of sarin, while DMNP is the best simulant for soman. The adsorption properties of DCP or DMMP show a strong correlation with sarin over a range of MOFs, but the correlation between DMNP and soman is considerably weaker. Comparisons of results with both force fields indicate that our main conclusions are robust with respect to the force field used to define adsorbate-MOF interactions.
The adsorption equilibrium constants of monovalent and divalent cations to material surfaces in aqueous media are central to many technological, natural, and geochemical processes. Cation adsorption-desorption is often proposed to occur in concert with proton transfer on hydroxyl-covered mineral surfaces, but to date this cooperative effect has been inferred indirectly. This work applies density functional theory-based molecular dynamics simulations of explicit liquid water/mineral interfaces to calculate metal ion desorption free energies. Monodentate adsorption of Na+, Mg2+, and Cu2+ on partially deprotonated silica surfaces are considered. Na+ is predicted to be unbound, while Cu2+ exhibits binding free energies to surface SiO- groups that are larger than those of Mg2+. The predicted trends agree with competitive adsorption measurements on fumed silica surfaces. As desorption proceeds, Cu2+ dissociates one of the H2O molecules in its first solvation shell, turning into Cu2+(OH-)(H2O)3, while Mg remains Mg2+(H2O)6. The protonation state of the SiO- group at the initial binding site does not vary monotonically with cation desorption.
Cubic zirconium tungstate (α-ZrW2O8), a well-known negative thermal expansion material, has been investigated within the framework of density functional perturbation theory (DFPT), combined with experimental characterization to assess and validate computational results. Using combined Fourier transform infrared measurements and DFPT calculations, new and extensive assignments were made for the far-infrared (<400 cm−1) spectrum of α-ZrW2O8. A systematic comparison of DFPT-simulated infrared, Raman, and phonon density-of-state spectra with Fourier transform far-/mid-infrared and Raman data collected in this study, as well as with available inelastic neutron scattering measurements, shows the superior accuracy of the PBEsol exchange-correlation functional over standard PBE calculations for studying the spectroscopic properties of this material.
Classical molecular dynamics simulation was used to study the adsorption of Na+, Ca2+, Ba2+, and Cl- ions on gibbsite edge (1 0 0), basal (0 0 1), and nanoparticle (NP) surfaces. The gibbsite NP consists of both basal and edge surfaces. Simulation results indicate that Na+ and Cl- ions adsorb on both (1 0 0) and (0 0 1) surfaces as inner-sphere species (i.e., no water molecules between an ion and the surface). Outer-sphere Cl- ions (i.e., one water molecule between an ion and the surface) were also found on these surfaces. On the (1 0 0) edge, Ca2+ ions adsorb as inner-sphere and outer-sphere complexes, whereas on the (0 0 1) surface, outer-sphere Ca2+ ions are the dominant species. Ba2+ ions were found as inner-sphere and outer-sphere complexes on both surfaces. Calculated ion surface coverages indicate that, for all ions, surface coverages are always higher on the basal surface compared to those on the edge surface. More importantly, surface coverages for cations on the gibbsite NP are always higher than those calculated for the (1 0 0) and (0 0 1) surfaces. This enhanced ion adsorption behavior for the NP is due to the significant number of inner-sphere cations found at NP corners. Outer-sphere cations do not contribute to the enhanced surface coverage. In addition, there is no ion adsorption enhancement observed for the Cl- ion. Our work provides a molecular-scale understanding of the relative significance of ion adsorption onto gibbsite basal versus edge surfaces and demonstrates the corner effect on ion adsorption on NPs.
Molecular tracers that can be selectively placed underground and uniquely identified at the surface using simple on-site spectroscopic methods would significantly enhance subsurface fluid monitoring capabilities. To ensure their widespread utility, the solubility of these tracers must be easily tuned to oil- or water-wet conditions as well as reducing or eliminating their propensity to adsorb onto subsurface rock and/or mineral phases. In this work, molecular dynamics simulations were used to investigate the relative solubilities and mineral surface adsorption properties of three candidate tracer compounds comprising Mg–salen derivatives of varying degrees of hydrophilic character. Simulations in water–toluene liquid mixtures indicate that the partitioning of each Mg–salen compound relative to the interface is strongly influenced by the degree of hydrophobicity of the compound. Simulations of these complexes in fluid-filled mineral nanopores containing neutral (kaolinite) and negatively charged (montmorillonite) mineral surfaces reveal that adsorption tendencies depend upon a variety of parameters, including tracer chemical properties, mineral surface type, and solvent type (water or toluene). Simulation snapshots and averaged density profiles reveal insight into the solvation and adsorption mechanisms that control the partitioning of these complexes in mixed liquid phases and nanopore environments. As a result, this work demonstrates the utility of molecular simulation in the design and screening of molecular tracers for use in subsurface applications.
Elastic and thermodynamic properties of negative thermal expansion (NTE) αα-ZrW2O8 have been calculated using PBEsol and PBE exchange-correlation functionals within the framework of density functional perturbation theory (DFPT). Measured elastic constants are reproduced within ~2% with PBEsol and 6% with PBE. The thermal evolution of the Grüneisen parameter computed within the quasi-harmonic approximation exhibits negative values below the Debye temperature, consistent with observation. The standard molar heat capacity is predicted to be C $O\atop{P}$=192.2 and 193.8 J mol-1K-1 with PBEsol and PBE, respectively. These results suggest superior accuracy of DFPT/PBEsol for studying the lattice dynamics, elasticity and thermodynamics of NTE materials.