Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “geologic”

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 1,459 records · Page 81
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Understanding the Hydronium Cation in the Solid State: A Study in Synthetic Hydronium Uranyl Phosphate and Arsenate Mineral Systems and Their Irradiation Stability

Here, we present neutron diffraction, electronic structure calculations, and optical vibrational spectroscopic characterization of hydronium uranyl phosphate (HUP), the synthetic analog of chernikovite ((H 3 O)(UO 2 )(PO 4 )·3H 2 O) to gain insight into structural features of the hydronium cation in the solid state. HUP crystallizes in P2 1 /c instead of previously reported P4/ncc 1 at room temperature and Pccn below 302 K. Lower symmetry was required due to hydronium and interstitial water positions within the structure. Crystallographic positions of hydronium cations determined from powder neutron diffraction data are bolstered by Raman and attenuated total reflectance infrared spectroscopic measurements and are further informed by density functional theory with phonon eigenvector analysis for spectral assignments. Finally, HUP and its arsenate analog (trogerite, (H 3 O)(UO 2 )(AsO 4 )·3H 2 O), (HUAs) were studied using He 2+ irradiation as an analog for α radiolysis to investigate the irradiation stability of these phases and the stability of hydronium cations in the solid state. Dose studies were employed wherein each sample was irradiated to 5, 10, 15, 25, and 50 MGy. Structural insight regarding irradiated materials is gained using optical vibrational spectroscopy and powder X-ray diffraction. Surprising irradiation stability of HUAs was found up to 50 MGy of dose, which may have important implications for understanding and modeling the geologic stability of legacy U.

actinides↗
DOE OSTI2020

NWPEsSe: an Adaptive-Learning Global Optimization Algorithm for Nanosized Cluster Systems

Global optimization constitutes an important and fundamental problem in theoretical studies in many chemical fields, such as catalysis, materials or separations problems. In this paper, a novel algorithm has been developed for the global optimization of large systems including neat and ligated clusters in gas phase, and supported clusters in periodic boundary conditions. The method is based on an updated artificial bee colony (ABC) algorithm method, that allows for adaptive-learning during the search process. The new algorithm is tested against four classes of systems of diverse chemical nature: gas phase Au_55, ligated Au_8^(2+), Au_8 supported on graphene oxide and defected rutile, and a large cluster assembly ?[Co?_6 Te_8 (PEt_3 )_6][C_60 ]_n, with sizes ranging between 1 to 3 nm and containing up to 1300 atoms. Reliable global minima (GMs) are obtained for all cases, either confirming published data or reporting new lower energy structures. The algorithm and interface to other codes in the form of an independent program, Northwest Potential Energy Search Engine (NWPEsSe), is freely available and it provides a powerful and efficient approach for global optimization of nanosized cluster systems. The work described in this publication was performed at Pacific Northwest National Laboratory (PNNL), which is operated by Battelle for the United States Department of Energy (DOE) under Contract DE-AC05-76RL0180. J. Z. and V.-A. G. acknowledge support from DOE, Office of Science, Office of Basic Energy Sci-ences, Chemical, Geological and Biological Sciences Division and computing resources from PNNL’s Research Computing Facility and the National Energy Research Scientific Computing Center.

Zhang, Jun↗
DOE OSTI2021

Density Functional Tight-Binding Simulations Reveal the Presence of Surface Defects on the Quartz (101)–Water Interface

Understanding the structure and reactivity of quartz–water interfaces is critical for numerous applications in the geological, environmental, and biological sciences. However, disagreements on the atomic-level structure of the interfaces between experiments and simulations are hampering our ability to predict the surface reactivity. Here, we used density functional tight-binding (DFTB)-based molecular dynamics to simulate a series of quartz (101) surfaces having different types and densities of surface defects in water and compared them with the structures determined by X-ray reflectivity measurements. The DFTB simulations are able to reproduce previous classical and quantum mechanical predictions of the pristine quartz (101)–water interface that disagree with experimental observations. To remedy this situation, a set of defective quartz surfaces having various surface silicon (Si) vacancies were built as indicated by recent experimental studies. We found that the rotation of surface [SiO 4 ] tetrahedra near Si vacancies can lead to outward displacements of Si atoms similar to those observed in the experiments. The presence of additional surface Si vacancies caused inward relaxations of terminal oxygens through the formation of hydrogen bonds. Finally, the overall results indicate that the quartz (101)–water interface may include a mixture of geminal (≡Si–(OH) 2 )- and vicinal (≡Si–OH)-type silanol groups together with the presence of surface Si vacancies.

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

Molecular Dynamics Simulations of Calcite Fracture in Water

Calcite (CaCO 3 ) is one of the most common minerals in geologic and engineered systems. It is often in contact with aqueous solutions, causing chemically assisted fracture that is critical to understanding the stability of subsurface systems and manmade structures. Calcite fracture was evaluated with reactive molecular dynamics simulations, including the impacts of crack tip geometry (notch), the presence of water, and surface hydroxyl groups. Chemo-mechanical weakening was assessed by comparing the loads where fracture began to propagate. Our analyses show that in the presence of a notch, the load at which crack growth begins is lower, compared to the effect of water or surface hydroxyls. Additionally, the breaking of two adjacent Ca–O bonds is the kinetic limitation for crack initiation, since transiently broken bonds can reform, not resulting in crack growth. In aqueous environments, fresh (not hydroxylated) calcite surfaces exhibited water strengthening. Manual addition of H + and/or OH – species on the (104) calcite surface resulted in chemo-mechanical weakening of calcite by 9%. Achieving full hydroxylation of the calcite surface was thermodynamically and kinetically limited, with only 0.17–0.01 OH/nm 2 surface hydroxylation observed on the (104) surface at the end of the simulations. In conclusion, the limited reactivity of pure water with the calcite surface restricts the chemo-mechanical effects and suggests that reactions between physiosorbed water and localized structural defects may be dominating the chemo-mechanical process in the studies where water weakening has been reported.

Calcite↗
DOE OSTI2020

Electronic and Vibrational Contributions to the Bulk Stabilities of Trivalent 3d Transition Metal Oxyhydroxides from Electronic Structure Calculations

Any understanding of the thermodynamic driving forces active in the global cycles of earth abundant metals requires attendant knowledge of the bulk mineral free energies. This can be particularly challenging for metal (oxyhydr)oxide phases which at low temperature aqueous conditions can be produced at a variety of hydration states and particle sizes. Here we present high quality density functional theory calculations of the relative bulk stabilities of a geologically important family of polymorphic trivalent metal oxyhydroxide crystal structures MOOH (M=Cr, Al, Fe). To improve understanding of their relative stabilities we permute the metal cations into the various possible canonical structure types including goethite (isostructural to diaspore and bracewillite), lepidocrocite (isosctructural to boehmite), guyaniite, and grimaldiite. Electronic, vibrational, and magnetic contributions are considered on a consistent theoretical footing using hybrid functionals and converged plane-wave bases. We find the considered bulk structures to be thermodynamically degenerate within the precision of the available theory, despite the fact that the vibrational and magnetic contributions vary substantially among the cations. A significant conclusion is that Cr3+ supports strong hydrogen bonds, low O-H stretching frequencies and facile H movement relative to the other cations.

Prange, Micah P.↗
DOE OSTI2020

Using Atom Dynamics to Map the Defect Structure Around an Impurity in Nano-Hematite

The geochemical cycling of actinides such as U is strongly controlled by Fe (oxyhydr)oxides. Indeed, the strong affinity of these minerals for U has factored into the design of geologic repositories for the long-term storage of nuclear waste. Many decades of work have focused on detailing the local coordination environment U associated with hematite (a-Fe2O3), a common and stable Fe(III) oxide, in order to better predict the fate and transport of U in the sub-surface. Here, for the first time, molecular dynamic simulations and atomically resolved electron microscopy of uranium atom mobility under the beam were used to map the topology of defect clusters surrounding structurally incorporated U in hematite. The ability to observe vacancies by means of metal atom hopping can be used to directly probe defects in other materials, with promising applications in geochemistry and materials science.

Ilton, Eugene S.↗
DOE OSTI2023

Insights into Waterflooding in Hydrocarbon-Bearing Nanochannels of Varying Cross Sections from Mesoscopic Multiphase Flow Simulations

Waterflooding is one of the geotechnique used to recover fuel sources from nanoporous geological formations. The scientific understanding of the process that involves the multiphase flow of nanoconfined fluids, however, has lagged, mainly due to the complex nanopore geometries and chemical compositions. To enable benchmarked flow of nanoconfined fluids, the architected geomaterials, such as synthesized mesoporous silica with tunable pore shapes and surface chemical properties, are used for designing and conducting experiments and simulations. This work uses a modified many-body dissipative particle dynamics (mDPD) model with accurately calibrated parameters to perform parametric flow simulations for studying the influences of waterflooding driven power, pore shape, surface roughness, and surface wettability on the multiphase flow in heptane-saturated silica nanochannels. Remarkably, up to 80\% reduction in the effective permeability is found for water-driven heptane flow in a baseline 4.5 nm-wide slit channel, when compared with the Hagen–Poiseuille equation. In the 4.5 nm-wide channels with architected surface roughness, the flow rate is found either higher or lower than the baseline case, depending on the shape and size of cross-sections. High wettability of the solid surface to water is essential for achieving high recovery of heptane, regardless of surface roughness. When the solid surface is less wetting or non-wetting to water, the existence of an optimal waterflooding driven power is found to allow for the highest possible recovery. A detailed analysis on the evolution of the transient water-heptane interface in those nanochannels is presented to elucidate the underlying mechanisms that impact or dictate the multiphase flow behaviors.

58 GEOSCIENCES↗
DOE OSTI2024

Molecular Mechanisms of Sorbed Ion Effects during Boehmite Particle Aggregation

Classical theories of particle aggregation, such as Derjaguin–Landau–Verwey–Overbeek (DLVO), do not explain recent observations of ion-specific effects or the complex concentration dependence for aggregation. Thus, here, we probe the molecular mechanisms by which selected alkali nitrate ions (Na + , K + , and NO 3 –) influence aggregation of the mineral boehmite (γ-AlOOH) nanoparticles. Nanoparticle aggregation was analyzed using classical molecular dynamics (CMD) simulations coupled with the metadynamics rare event approach for stoichiometric surface terminations of two boehmite crystal faces. Calculated free energy landscapes reveal how electrolyte ions alter aggregation on different crystal faces relative to pure water. Consistent with experimental observations, we find that adding an electrolyte significantly reduces the energy barrier for particle aggregation (~3–4×). However, in this work, we show this is due to the ions disrupting interstitial water networks, and that aggregation between stoichiometric (010) basal–basal surfaces is more favorable than between (001) edge–edge surfaces (~5–6×) due to the higher interfacial water densities on edge surfaces. The interfacial distances in the interlayer between aggregated particles with electrolytes (~5–10 Å) are larger than those in pure water (a few Ångströms). Together, aggregation/disaggregation in salt solutions is predicted to be more reversible due to these lower energy barriers, but there is uncertainty on the magnitudes of the energies that lead to aggregation at the molecular scale. By analyzing the peak water densities of the first monolayer of interstitial water as a proxy for solvent ordering, we find that the extent of solvent ordering likely determines the structures of aggregated states as well as the energy barriers to move between them. Further, the results suggest a path for developing a molecular-level basis to predict the synergies between ions and crystal faces that facilitate aggregation under given solution conditions. Such fundamental understanding could be applied extensively to the aggregation and precipitation utilization in the biological, pharmaceutical, materials design, environmental remediation, and geological regimes.

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

Molecular Origin of Wettability Alteration of Subsurface Porous Media upon Gas Pressure Variations

Upon extraction/injection of a large quantity of gas from/into a subsurface system in shale gas production or carbon sequestration, the gas pressure varies remarkably, which may significantly change the wettability of porous media involved. Mechanistic understanding of such changes is critical for designing and optimizing a related subsurface engineering process. Using molecular dynamics simulations, in this work we have calculated the contact angle of a water droplet on various solid surfaces (kerogen, pyrophyllite, calcite, gibbsite, and montmorillonite) as a function of CO 2 or CH 4 gas pressure up to 200 atm at a temperature of 300 K. The calculation reveals a complex behavior of surface wettability alteration by gas pressure variation depending on surface chemistry and structure, and molecular interactions of fluid molecules with surfaces. As the CO 2 gas pressure increases, a partially hydrophilic kerogen surface becomes highly hydrophobic, while a calcite surface becomes more hydrophilic. Considering kerogen and calcite being the major components of a shale formation, we postulate that the wettability alteration of a solid surface induced by a gas pressure change may play an important role in fluid flows in shale gas production and geological carbon sequestration.

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

Iron Impurity Impairs the CO 2 Capture Performance of MgO: Insights from Microscopy and Machine Learning Molecular Dynamics

Magnesium oxide (MgO) is a promising sorbent for direct air capture (DAC) of carbon dioxide. Iron (Fe) is a common impurity in naturally occurring MgO and minerals used to produce MgO, yet a molecular-scale understanding of Fe-doping effects on carbonation is lacking. Here, in this study, we observed reduced carbonation performance in Fe-doped MgO experimentally. The energetics of adsorbing a (bi)carbonate ion on pristine and Fe-doped MgO(001) surfaces were further investigated using ab initio and machine learning potential molecular dynamics coupled with metadynamics simulations. Both pristine and Fe-doped surfaces exhibited a basic (OH – ) hydration layer, where the (bi)carbonate ion adsorption is thermodynamically favorable. However, the dissolution of surface Fe had smaller energy barriers and was more favorable than Mg. Leached Fe likely neutralized the near-surface basicity, yielding reduced reactivity on Fe-doped MgO. Our observations offer critical insights for material selection and emphasize the importance of evaluating the geologic origin of earth materials used for DAC.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Role of Polymer Architecture in CO 2 Capture from Air Using Supported Poly(alkylenimine)s: Linear vs Branched Polymers

Direct air capture (DAC) of CO 2 coupled with geologic storage is a promising climate change mitigation strategy, with some applications employing amines supported on porous solids as CO 2 sorbents. While branched poly(ethylenimine) (PEI) is the standard benchmark amine material, it suffers from limited oxidative stability. Poly(propylenimine) (PPI), as an alternative, has previously demonstrated improved resistance to degradation under harsh oxidative conditions. Linear and branched PEI are commercially available, though at different molecular weights, while PPI is not commercially available. For this reason, a comparative study of all four polymers (linear PEI, branched PEI, linear PPI, branched PPI) has not been reported for DAC. In this study, we synthesize and compare low-molecular-weight (∼800 g/mol) linear (L) and branched (B) PEI and PPI supported on a model support, SBA-15 silica. These materials are evaluated for CO 2 adsorption under dry, DAC-relevant conditions (400 ppm of CO 2 , 30 °C). LPPI exhibited the highest amine efficiency at all loadings, reaching a maximum of 0.14 mmol CO 2 /mmol N, outperforming BPEI, while LPEI consistently showed the lowest uptake capacity. Temperature-programmed desorption reveals that the structure of the amine polymer impacts the CO 2 binding strength, with branched polymers displaying higher desorption energies of 102−111 kJ/mol. In situ infrared spectroscopy experiments show that all sorbents preferentially capture CO 2 as ammonium carbamate. Isobaric CO 2 uptake studies further underscore the influence of polymer mobility and support pore crowding on performance, while demonstrating the sorbents’ performance at elevated temperatures and CO 2 concentrations. All materials demonstrated good stability over 25 adsorption−desorption cycles using thermal regeneration in an inert gas purge, with only BPPI displaying a 10−11% decrease in capacity/amine efficiency during cycling, possibly due to the loss of low molecular weight, oligomeric amines. This is the first side-by-side comparison of the CO 2 sorption properties of linear and branched PEI and PPI with similar molecular weights. These findings highlight the significant role of polymer architecture in CO 2 capture efficiency and inform future designs of durable, high-performance DAC sorbents.

adsorbents↗
DOE OSTI2020

Simulations and Experiments Reveal Effect of Nanopores on Helium Diffusion in Quartz

The diffusion properties of noble gases in minerals are widely used to reconstruct the thermal histories of rocks. In this work, we combine density functional theory (DFT) calculations with laboratory experiments to investigate controls on helium diffusion in quartz. DFT calculations for perfect α-quartz predict substantially lower activation energies and frequency factors for helium diffusion than observed in laboratory experiments, especially in the [001] direction. These results imply that no helium could be retained in quartz at Earth surface temperatures, which conflicts with observations of partial cosmogenic 3 He retention over geologic time scales. Here, we implement a model of helium diffusion in α-quartz modulated by nanopore defects that disrupt energetically favorable diffusion pathways. In this model, we find that laboratory-determined diffusivities can be most closely reproduced when a helium atom encounters ~70 nanopore sites per million interstitial sites. The results of our model indicate that diffusion of helium in natural quartz, like other noble gases in other minerals, can be significantly modulated by extended defects.

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

High-Temperature Thermodynamics of Cerium Silicates, A-Ce 2 Si 2 O 7 , and Ce 4.67 (SiO 4 ) 3 O

Lanthanide disilicates and oxyapatites have potential roles in high temperature applications as thermal (TBC) and environmental barrier coatings (EBC), or possible alteration phases in geological nuclear waste repositories. However, those Ce 3+ -bearing silicates have only been limitedly studied. In this work, we performed detailed structural and thermodynamic investigations on A-Ce 2 Si 2 O 7 (tetragonal, P4 1 ) and Ce 4.67 (SiO 4 ) 3 O (hexagonal, P6 3 /m). The high temperature structural behaviors and coefficients of thermal expansion were determined by in situ high temperature synchrotron X-ray diffraction (HT-XRD) implemented with Rietveld analysis and thermogravimetric analysis coupled with differential scanning calorimetry (TGA-DSC). A-Ce 2 Si 2 O 7 was found to be stable in N 2 and air up to ~1483 K with anisotropic thermal expansion along the a and c axes (α a = 12.3 × 10 -6 K -1 and α c = 12.4 × 10 -6 K -1 ). Ce 4.67 (SiO 4 ) 3 O had a slow partial oxidation between 533 K and 873 K to a new nonstoichiometric phase Ce 3+ 1.67- xCe4 + x Ce 3+ 3 (SiO 4 ) 3 O 1+0.5x , followed by a thermal decomposition to CeO 2 and SiO 2 at ~1000 K in air. By using high temperature oxide melt solution calorimetry at 973 K with lead borate as the solvent, the standard enthalpy of formation was determined for A-Ce 2 Si 2 O 7 (-3825.1 ± 6.0 kJ/mol) and Ce 4.67 (SiO 4 ) 3 O (-7391.3 ± 9.5 kJ/mol). Finally, these thermodynamic parameters were compared with those of CeO 2 , CeSiO 4 , and other silicate oxyapatites for examining their chemical stability in high temperature environments relevant for aeronautical applications, mineral formation, and nuclear fuel cycle.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Photocatalytic Oxidation of Dissolved Mn 2+ by TiO 2 and the Formation of Tunnel Structured Manganese Oxides

The redox reaction of manganese (Mn) is of great environmental, geological, and public health significance, as Mn oxides control the distribution and electron flow of numerous nutrients and contaminants in natural and engineered environments. Current understanding on the oxidation pathways of Mn(II) to Mn(III/IV) mainly focuses on biotic processes due to their much higher oxidation rates than those of abiotic processes. This study demonstrates rapid photocatalytic oxidation of Mn 2+ (aq) under circumneutral conditions catalyzed by naturally abundant semiconducting TiO 2 minerals. Notably, the photocatalytic oxidation rates are comparable to or even higher than those of reported biotic/abiotic processes. In addition, the rapid photocatalytic oxidation leads to the formation of large tunnel structured Mn oxides (todorokite and romanechite) on the surface of TiO 2 . These findings suggest that photocatalytic oxidation of Mn 2+ (aq) by natural semiconducting minerals is likely an important yet previously overlooked pathway for understanding the occurrence of natural Mn oxide coatings on rock surfaces. In addition, considering the increasing input of photoreactive engineered nanoparticles into environmental systems, this work shows the potential impacts of nanoparticles on influencing natural redox cycles.

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

Plutonium Co-precipitation with Calcite

The mobility of plutonium (Pu) in the environment is affected by Pu–mineral interactions, such as adsorption–desorption and structural incorporation. Calcite (CaCO 3 ) is a common secondary phase in near surface environments and a major component of many rocks and soils and is expected to form as an alteration product of cement-based materials planned for use in geological repositories. The reactivity of the calcite surface and its ability to tolerate significant variations in its chemical composition through substitution of Ca for other cations make calcite a potentially important sink for environmental contaminants. Here, single crystals of calcite were synthesized from aqueous solutions in equilibrium with air containing Pu as either Pu(VI) or Pu(IV) and characterized using a combination of laser ablation inductively coupled plasma mass spectrometry (LA–ICP–MS) and X-ray absorption spectroscopy (XAS). These data are used to assess the amount, structure, and oxidation state of Pu co-precipitated into calcite, providing insight into the potential for Pu sequestration in calcite precipitates. Overall, the XAS and LA–ICP–MS data support the co-precipitation of plutonyl [Pu(VI/V)] in the bulk calcite, although the exact nature of the co-precipitated Pu complex is difficult to elucidate in the synthesized material. Co-precipitated plutonyl could be incorporated in either distorted Ca lattice sites or defect sites, and we provide evidence to suggest that Pu(VI) is reduced mainly to Pu(V) in the precipitated solid. LA–ICP–MS additionally shows that the co-precipitation of Pu(VI/V) is favored over the co-precipitation of Pu(IV). Altogether, our results suggest that Pu sequestration in calcite under environmental conditions could immobilize Pu and isolate it from groundwater interactions in contaminated environments.

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

Structure–Composition Relationships for Mg–Ni and Mg–Fe Olivine

Olivine is a dynamic and important mineral in the crust and mantle with relevance to processes important to climate change technology, such as geologic carbon storage and critical mineral recovery. In this work, we critically evaluated and compiled a new database of olivine diffraction data, lattice parameters, and composition to enable rapid Ni-Mg-Fe olivine composition determination. A compilation of olivine X-ray diffraction data and chemical compositions from both the literature and the International Centre for Diffraction Data (ICDD) powder database was assembled to plot both the forsterite-fayalite and forsterite-liebenbergite solid solution lines. Here we present an expanded dataset to delineate equations and relationships used for quantifying the correlations between olivine lattice parameters and chemical compositions in Mg 2 SiO 4 -Fe 2 SiO 4 (forsterite-fayalite) and Mg 2 SiO 4 -Ni 2 SiO 4 (forsterite-liebenbergite) olivine solid solution series.

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

Ab Initio Calculation of Equilibrium Isotopic Fractionations of Potassium and Rubidium in Minerals and Water

Here, we used first-principle approaches to calculate the equilibrium isotopic fractionation factors of potassium (K) and rubidium (Rb) in a variety of minerals of geological relevance (orthoclase, albite, muscovite, illite, sylvite, and phlogopite). We also used molecular dynamics simulations to calculate the equilibrium isotopic fractionation factors of K in water. Our results indicate that K and Rb form bonds of similar strengths and that the ratio between the equilibrium fractionations of K and Rb is approximately 3–4. Under low-temperature conditions relevant to weathering of continents or alteration of seafloor basalts (~25 °C), the K isotopic fractionation between solvated K+ and illite (a proxy for K-bearing clays) is +0.24‰, exceeding the current analytical precision, so equilibrium isotopic fractionation can induce measurable isotopic fractionations for this system at low temperature. These findings, however, cannot easily explain why the δ 41 K value of seawater is shifted by +0.6‰ relative to igneous rocks. Our results indicate that part of the observed fractionation is most likely due to kinetic effects. The narrow range of mean force constants for K and Rb in silicate minerals suggests that phase equilibrium is unlikely to create large K and Rb isotopic fractionations at magmatic temperatures (at least in silicate systems). Kinetic effects associated with diffusion can, however, produce large K and Rb isotopic fractionations in igneous rocks.

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

Enhanced Diffusive Transport in Fluctuating Porous Media

Mass transport within porous structures is a ubiquitous process in biological, geological, and technological systems. Despite the importance of these phenomena, there is no comprehensive theory that describes the complex and diverse transport behavior within porous environments. While the porous matrix itself is generally considered a static and passive participant, many porous environments are in fact dynamic, with fluctuating walls, pores that open and close, and dynamically changing cross-links. While diffusion has been measured in fluctuating structures, notably in model biological systems, it is rarely possible to isolate the effect of fluctuations because of the absence of control experiments involving an identical static counterpart, and it is generally impossible to observe the dynamics of the structure. In this paper, we present a direct comparison of the diffusion of nanoparticles of various sizes within a trackable, fluctuating porous matrix and a geometrically equivalent static matrix, in conditions spanning a range of regimes from obstructed to highly confined. The experimental system comprised a close-packed layer of colloidal spheres that were either immobilized to a planar surface or allowed to fluctuate locally, within the space defined by their nearest neighbors. Interestingly, the effective long-time diffusion coefficient was approximately 35–65% greater in the fluctuating porous matrix than in the static one (depending on the size of the nanoparticle probes), regardless of the geometric regime. This was explained by considering the enhancing effects of matrix fluctuations on the short-time diffusion coefficient and cooperative “gate-opening” motions of matrix particles and nanoparticle probes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗