SEARCH · Search NASA

Results for “pdf”

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.

SourceResearch connectionFrom yearThrough year

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

At least 55 records · Page 3

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Medium-range order and compositional correlation in metallic glasses

The compositional atomic ordering in metallic glasses was studied by simulation focusing on the medium-range order (MRO). Many metallic alloy liquids and glasses show MRO characterized by the oscillations in the atomic pair-distribution function (PDF) beyond the first peak, which decay exponentially with distance. To study the effects of the local chemical order on MRO, we examine the compositionally resolved PDF and its MRO for models of various binary metallic alloy glasses. We show that compositional ordering is limited mostly to the nearest-neighbor atoms and the MRO is largely independent of the compositional order. For some elements that strongly repel each other in the alloy, a second MRO periodicity is observed owing to the distinct correlations among them. These results are discussed in light of the idea that the MRO oscillations in the PDF describe the correlations in the atomic density fluctuations, rather than the detailed local atomic structure.

Atomic structure↗

DOE OSTI2025

Kaon gluon parton distribution and momentum fraction from 2+1+1 lattice QCD with high statistics

We present a high-statistics lattice-QCD determination of the kaon gluon parton distribution function and gluon momentum fraction. We use clover valence fermion action to take 1,296,640 kaon-correlator measurements on a highly improved staggered quark ensemble with 𝑎 ≈ 0.12 fm and 310-MeV pion mass generated by the MILC Collaboration. A detailed investigation into the impact of gauge-link smearing on the gluonic matrix elements indicates that five steps of hypercubic smearing offer an effective balance between signal quality and preservation of long-distance physics. We report a nonperturbatively renormalized kaon gluon momentum fraction of Math output error at 𝜇 = 2 GeV in the Math output error scheme. Using reduced pseudo-Ioffe-time-distribution matrix elements and pseudo-parton-distribution-function (PDF) matching, we extract the kaon gluon PDF and compare with the prediction from the Dyson-Schwinger equation and with the pion PDF obtained from the same ensemble.

First-principles calculations↗

DOE OSTI2025

Assessment of thermally driven local structural phase changes in 1⁢𝑇′−MoTe 2

The role of layer disorder is important in establishing the topological phases of MoTe 2 . A rich tapestry of atomic ordering influences the structural phase transitions (SPTs), but there is little understanding of the mechanistic details of the phase transition. An atomistic level study was conducted to investigate the local structure of the 1⁢𝑇′ and 𝑇 𝑑 phases of MoTe 2 by using the pair distribution function (PDF) technique. While the average structure exhibits an SPT and coexistence of phases as a function of temperature, the local structure showed the suppression of SPT. The sample retained its monoclinic structure at all temperatures in short-range order. A sharp PDF peak observed at short distances indicated a strong atom-atom correlation between the Mo and Te atoms within the Mo octahedra. In addition, a large-box modeling of the PDF data indicated a preferential motion of Te atoms towards 𝑐 axis at all temperatures. Structural defects, such as stacking faults, likely result in the coexistence of phases in the average structure and suppress the local SPT of MoTe 2 . These results are stepping stones to understand the long-debated origins of structural, vibrational, and electronic properties of MoTe 2 and similar transition metal dichalcogenides.

2-dimensional systems↗

DOE OSTI2025

Learning PDFs through interpretable latent representations in Mellin space

Representing the parton distribution functions (PDFs) of the proton and other hadrons through flexible, high-fidelity parametrizations has been a long-standing goal of particle physics phenomenology. This is particularly true since the chosen parametrization methodology can play an influential role in the ultimate PDF uncertainties as extracted in QCD global analyses; these, in turn, are often determinative of the reach of experiments at the LHC and other facilities to nonstandard physics, including at large 𝑥, where parametrization effects can be significant. In this study, we explore a series of encoder-decoder machine-learning (ML) models with various neural-network topologies as efficient means of reconstructing PDFs from meaningful information stored in an interpretable latent space. Given recent effort to pioneer synergies between QCD analyses and lattice-gauge calculations, we formulate a latent representation based on the behavior of PDFs in Mellin space, i.e., their integrated moments, and test the ability of various models to decode PDFs from this information faithfully. We introduce a numerical package, PDFdecoder, which implements several encoder-decoder models to reconstruct PDFs with high fidelity and use this end-to-end tool to explore how such neural-network-based models might connect PDF parametrizations to underlying properties like their Mellin moments. We additionally dissect patterns of learned correlations between encoded Mellin moments and reconstructed PDFs that suggest opportunities for further improvements to ML-based approaches to PDF parametrizations and uncertainty quantification.

Machine learning↗

DOE OSTI2025

Score-Based Physics-Informed Neural Networks for High-Dimensional Fokker–Planck Equations

The Fokker-Planck (FP) equation is a foundational partial differential equation (PDE) in stochastic processes involving Brownian motions. However, the curse of dimensionality (CoD) poses a formidable challenge when dealing with high-dimensional FP equations. Although Monte Carlo simulation and (vanilla) Physics-Informed Neural Networks (PINNs) have shown the potential to tackle CoD, both methods exhibit significant numerical errors in high dimensions when dealing with the probability density function (PDF) associated with Brownian motion. The point-wise PDF values tend to decrease exponentially as dimensionality increases, surpassing the precision of numerical simulations and resulting in substantial errors. In addition, due to its massive sampling, Monte Carlo fails to offer fast sampling. Modeling the logarithm likelihood (LL) via vanilla PINNs transforms the FP equation into a notoriously difficult Hamilton-Jacobi-Bellman (HJB) equation, which is impractical for PINN learning, whose error grows rapidly with dimension. To this end, we propose a novel approach utilizing a score-based solver to fit the score function in stochastic differential equations (SDEs). The score function, defined as the gradient of the LL, plays a fundamental role in inferring LL and PDF and enables fast SDE sampling, offering an effective means to overcome the CoD. Three fitting methods, Score Matching (SM), Sliced Score Matching (SSM), and Score-PINN, are introduced, each contributing unique advantages in computational complexity, accuracy, and generality. The proposed score-based SDE solver operates in two stages: first, employing score matching or Score-PINN to acquire the score function; and second, solving the LL via an ordinary differential equation (ODE) using the obtained score function. Comparative evaluations across these methods showcase varying trade-offs. The proposed methodology is evaluated across diverse SDEs, including anisotropic Ornstein-Uhlenbeck processes, geometric Brownian motion, and Brownian motion with varying eigenspace. We also test various distributions, including Gaussian, Log-normal, Laplace, and Cauchy distributions. The numerical results demonstrate the score-based SDE solver’s stability, speed, and performance across different experimental settings, solidifying its potential as a solution to CoD for high-dimensional FP equations.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

AdaParse

SF-25-127AdaParse (Adaptive Parallel PDF Parsing and Resource Scaling Engine) enables scalable, high-accuracy PDF parsing. AdaParse is a data-driven strategy that assigns an appropriate parser to each document, offering high accuracy for any computational budget. Moreover, it offers a workflow of various PDF parsing software that includes extraction tools: PyMuPDF, pypdf traditional OCR: Tesseract, modern OCR (e.g., Vision Transformers): Nougat and Marker

Siebenschuh, Carlo↗

DOE OSTI2025

Data from: 'Abiotic influences on continuous conifer forest structure across a subalpine watershed'

This package archives the core data used for analysis and inference in 'Abiotic influences on continuous conifer forest structure across a subalpine watershed' (Worsham et al., 2025). All data were collected in the East River, Washington Gulch, Slate River, and Coal Creek watersheds of Colorado. In the paper, we quantified the relative influence of climate, topographic, edaphic, and geologic factors on conifer stand structure and composition, and their functional relationships, at the watershed scale. We used waveform LiDAR data to derive spatially continuous stand structure metrics. We fused these with a species-level classification map to estimate tree species abundance. We applied generalized additive and generalized boosted models to evaluate the covariability of structural and compositional metrics with abiotic variables. The package contains the essential products required for reproducing our analysis and the tables and figures reported in the publication. The products comprise four classes: (1) geospatial data, (2) tabular data used for inferential analysis, (3) tabular data describing analytical results and performance statistics, and (4) a data user guide. (1) includes discretized waveform LiDAR data, locations and attributes of individual tree crowns, sampling locations and domain boundaries, a canopy height model, and raster files of estimated forest structural and compositional metrics at 100 m grid scale. (2) includes all response and explanatory variable values applied in inferential models. Response variables include conifer forest stand density, basal area, 95th percentile height, quadratic mean diameter, and others. Explanatory variables include climatic water deficit, actual evapotranspiration, elevation, heat load, soil available water content, and others. (3) includes results of training and testing several individual tree detection (ITD) algorithms, as well as inferential modeling results. (4) is a PDF user guide for this data package, including detailed descriptions and data dictionaries for all files. The data package root contains 17 assets: 8 compressed tape archive (.tar.gz) files, 5 comma-separated values (.csv) files, 3 Geographic Tagged Image File Format (GeoTIFF) (.tif) files, and 1 Portable Document Format (.pdf) file. The compressed .tar.gz archives contain ESRI shapefiles (.shp) .tif, compressed LASer (.laz), and .csv files. The archives must first be decompressed using the widely distributed command-line software utility TAR. All other files, including constituent files within the .tar.gz archives, can be opened in the open-source R statistical computing environment. Alternatively, .csv files may also be read in any simple text editor software or Microsoft Excel. Geospatial files including .shp and .tif files can also be opened in GIS software, such as QGIS (open-source) or ESRI ArcGIS (proprietary). The .pdf Data User Guide can be read with Adobe Acrobat Reader or other compatible readers.

2018 NEON and 2025 CHESS Campaigns↗

DOE OSTI2026

Data and scripts associated with “Non-random processes impacting organic matter chemistry are maximized in mid-order streams”

NOTE: The manuscript associated with this data package is currently in review. The data may be revised based on reviewer feedback. Upon manuscript acceptance, this data package will be updated with the final dataset and additional metadata. This data package is associated with the publication “Non-random processes impacting organic matter chemistry are maximized in mid-order streams” submitted to Limnology and Oceanography (L&O) by Danczak et al. (in review). This package contains data and scripts used to investigate dissolved organic matter (DOM) molecular chemistry and diversification processes across 47 surface-water sampling sites in the Yakima River Basin, Washington, USA, during an August 2021 sampling campaign. The package contains analyses of ultrahigh-resolution Fourier transform ion cyclotron resonance mass spectrometry (FTICR-MS), geochemical measurements, geospatial attributes, molecular diversity, and meta-metabolome ecological null models needed to reproduce the main manuscript results. The underlying field data were pulled from exising data packages at https://doi.org/10.15485/1892052 (Fulton et al., 2022) and https://doi.org/10.15485/1898914 (Grieger et al., 2022). For details on how to navigate data packages generated by this project, see https://data.ess-dive.lbl.gov/portals/PNNLRiverCorridorSFA/About. We thank the following organizations for providing access to field locations for sample collection: the United States Forest Service, Washington Department of Fish and Wildlife, Washington Department of Natural Resources, the Confederated Tribes and Bands of the Yakama Nation, and the Cowiche Canyon Conservatory. Research was conducted under Washington State Parks and Recreation Commission Scientific Research Permit #210901. We are grateful to the Yakama Nation Tribal Council and Yakama Nation Fisheries for their collaboration in facilitating sample collection and ensuring data usage aligns with their values and worldview. This data package contains an R-Markdown file for analyses and five folders: (1) Data, (2) Geospatial Data, (3) Supplemental_Files, (5) Figures_pdf, (4) and src. The Data folder contains tabular inputs and derived files used in the manuscript analysis. The Geospatial Data folder contains climate and water-balance, hydrologic, land-cover, population/regional water-use, stream, topographic, and stream-order attribute CSV files. The src folder contains scripts used to process data, run analyses, and generate figures. The Figures_pdf folder contains manuscript figure outputs. The Supplemental_Files folder contains supplemental analysis products. All files are .csv, .pdf, .html, .png, .R, .Rmd, .svg, or .tre. This data package is associated with the rcfsa-RC2-SPS_Null_Modeling repository found at https://github.com/river-corridors-sfa/rcfsa-RC2-SPS_Null_Modeling.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Data-Driven Closures and Assimilation for Stiff Multiscale Random Dynamics

Here, we introduce a data-driven and physics-informed framework for propagating uncertainty in stiff, multiscale random ordinary differential equations (RODEs) driven by correlated (colored) noise. Unlike systems subjected to Gaussian white noise, a deterministic equation for the joint probability density function (PDF) of RODE state variables does not exist in closed form. Moreover, such an equation would require as many phase-space variables as there are states in the RODE system. To alleviate this curse of dimensionality, we instead derive exact, albeit unclosed, reduced-order PDF (RoPDF) equations for low-dimensional observables/quantities of interest. The unclosed terms take the form of state-dependent conditional expectations, which are directly estimated from data at sparse observation times. However, for systems exhibiting stiff, multiscale dynamics, data sparsity introduces regression discrepancies that compound during RoPDF evolution. This is overcome by introducing a kinetic-like defect term to the RoPDF equation, which is learned by assimilating in sparse, low-fidelity RoPDF estimates. Two assimilation methods are considered, namely nudging and deep neural networks, which are successfully tested against Monte Carlo simulations.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Size resolved particle volatility for 30, 60, and 90 nm particles collected at the EPCAPE Mount Soledad site from 04/23/2023 to 06/15/2023

VTDMA GF-PDF: Growth factor probability density functions from the Gysel inversion are located here, where probability is out of 200. The instrument measured three different diameters 30 nm, 60nm and 90nm each diameter was selected for an hour. During that hour three different temperatures where scanned (40°C, 80°C, 160°C) in between each temperature a bypass scan was taken at room temperature. Each scan took 10-minutes and scanned growth factors from 0.2-1.2. Calibrations of the instrument with NaCl were performed on 4/26/2023 from 11:20 to 14:20, on 5/16/2023 from 11:40 to 16:00, and on 5/30/2023 from 10:20 to 14:20. Calibrations with (NH4)2SO4 were performed on 5/8/2023 from 11:45 to 14:40, on 5/23/2023 from 11:20 to 3:00, and on 6/14/2023 from 13:00 to 16:00. Time is recorded in seconds since 1/1/1904. GF_Vol_avg: Growth factor volume average calculated by the Gysel inversion. The instrument measured three different diameters 30 nm, 60nm and 90nm each diameter was selected for an hour. During that hour three different temperatures where scanned (40°C, 80°C, 160°C) in between each temperature a bypass scan was taken at room temperature. Each scan took 10-minutes and scanned growth factors from 0.2-1.2. Calibrations of the instrument with NaCl were performed on 4/26/2023 from 11:20 to 14:20, on 5/16/2023 from 11:40 to 16:00, and on 5/30/2023 from 10:20 to 14:20. Calibrations with (NH4)2SO4 were performed on 5/8/2023 from 11:45 to 14:40, on 5/23/2023 from 11:20 to 3:00, and on 6/14/2023 from 13:00 to 16:00. Time is recorded in seconds since 1/1/1904.

30nm growth factor volume average at 160C (GF_vol_↗

DOE OSTI2024

Measurement of the total and differential cross-sections of $t\bar{t}$$W$ production in pp collisions at $\sqrt{s}$ = 13 TeV with the ATLAS detector

Measurements of inclusive and differential production cross-sections of a top-quark-top-antiquark pair in association with a W boson ($t\bar{t}$$W$) are presented. They are performed by targeting final states with two same-sign or three isolated leptons (electrons or muons) and are based on $\sqrt{s}$ = 13 TeV proton-proton collision data with an integrated luminosity of 140 fb −1 , recorded from 2015 to 2018 with the ATLAS detector at the Large Hadron Collider. The inclusive $t\bar{t}$$W$ production cross-section is measured to be 880 ± 80 fb, compared to a reference theoretical prediction of 745 ± 50 (scale) ± 13 (2-loop approx.) ± 19 (PDF, α s ) fb. Differential cross-section measurements characterise this process in detail for the first time. Several particle-level observables are compared with a variety of theoretical predictions, which generally agree well with the normalised differential cross-section results. Additionally, the relative charge asymmetry of $t\bar{t}$$W$ + and $t\bar{t}$$W$ – is measured inclusively to be $A$$^{rel}_{C}$ = 0.33 ± 0.05, in very good agreement with the theoretical prediction of 0.322 ± 0.003 (scale) ± 0.007 (PDF), as well as differentially.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Photometric redshifts probability density estimation from recurrent neural networks in the DECam local volume exploration survey data release 2

Photometric wide-field surveys are imaging the sky in unprecedented detail. These surveys face a significant challenge in efficiently estimating galactic photometric redshifts while accurately quantifying associated uncertainties. In this work, we address this challenge by exploring the estimation of Probability Density Functions (PDFs) for the photometric redshifts of galaxies across a vast area of 17,000 square degrees, encompassing objects with a median 5 σ point-source depth of g = 24.3, r = 23 . 9 , i = 23.5, and z = 22.8 mag. Our approach uses deep learning, specifically integrating a Recurrent Neural Network architecture with a Mixture Density Network, to leverage magnitudes and colors as input features for constructing photometric redshift PDFs across the whole DECam Local Volume Exploration (DELVE) survey sky footprint. Subsequently, we rigorously evaluate the reliability and robustness of our estimation methodology, gauging its performance against other well-established machine learning methods to ensure the quality of our redshift estimations. Our best results constrain photometric redshifts with the bias of − 0 . 0013 , a scatter of 0.0293, and an outlier fraction of 5.1%. These point estimates are accompanied by well-calibrated PDFs evaluated using diagnostic tools such as Probability Integral Transform and Odds distribution. We also address the problem of the accessibility of PDFs in terms of disk space storage and the time demand required to generate their corresponding parameters.We present a novel Autoencoder model that reduces the size of PDF parameter arrays to one-sixth of their original length, significantly decreasing the time required for PDF generation to one-eighth of the time needed when generating PDFs directly from the magnitudes.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Direct numerical simulations of turbulent premixed cool flames: Global and local flame dynamics analysis

The cool flame dynamics, especially in turbulent flows, is of great interest for both practical application and fundamental research. Here, in this study, a series of direct numerical simulations of turbulent premixed n-C 7 H 16 /O 2 /O 3 /N 2 cool flames are performed, with the focus on the influence of turbulence intensity (u'/S L , where S L is the laminar flame speed) on the flame structure as well as the global and local cool flame dynamics. It is found that the cool flame front is considerably wrinkled by turbulence at high u'/S L , leading to significantly thickened turbulent cool flame brush and largely altered local reactivity compared with the reference laminar flame. However, the turbulent flame structure in the temperature space is found to be insensitive to u'/S L . Besides, with increasing u'/S L , the normalized turbulent cool flame speed (S L /S L ) is monotonically increased, attributed to substantial augmentation on the flame surface area (A T /A L ), while the stretching factor (I 0 ) remains almost constant and is smaller than 1. The underlying mechanisms for such variations are revealed through local flame dynamics analysis. Specifically, the local flame displacement speed S d is found to be strongly negatively correlated with flame curvature; meanwhile, such negative correlation and the probability distribution function (PDF) of flame curvature are barely influenced by u'/S L , leading to a weak dependence of I 0 on u'/S L . In contrast, the PDF of the tangential strain rate is found to span a much wider range and shift to the positive side as u'/S L increases, suggesting that the enhanced tangential strain rate is the main cause for the increase in surface area of the turbulent premixed cool flame. Finally, the influence of equivalence ratio on above findings is found to be insignificant, indicating that although the local reactivity of turbulent premixed cool flames is altered due to the differential diffusion, the resultant flame- stretch interaction is insensitive to the equivalence ratio. This study presents some unique cool flame dynamics that are distinct from hot flames, which can help improve the understanding and modeling of turbulent cool flames.

Cool flames↗

DOE OSTI2025

Uncertainty propagation in feed-forward neural network models

We develop new uncertainty propagation methods for feed-forward neural network architectures with leaky ReLU activation functions subject to random perturbations in the input vectors. In particular, we derive analytical expressions for the probability density function (PDF) of the neural network output and its statistical moments as a function of the input uncertainty and the parameters of the network, i.e., weights and biases. A key finding is that an appropriate linearization of the leaky ReLU activation function yields accurate statistical results even for large perturbations in the input vectors. This can be attributed to the way information propagates through the network. We also propose new analytically tractable Gaussian copula surrogate models to approximate the full joint PDF of the neural network output. To validate our theoretical results, we conduct Monte Carlo simulations and a thorough error analysis on a multi-layer neural network representing a nonlinear integro-differential operator between two polynomial function spaces. Our findings demonstrate excellent agreement between the theoretical predictions and Monte Carlo simulations.

MLP networks↗

DOE OSTI2026

Decoding the proton’s gluonic density with lattice QCD-informed machine learning

We present a first machine learning-based decoding of the gluonic structure of the proton from lattice QCD using a variational autoencoder inverse mapper (VAIM). Harnessing the power of generative AI, we predict the parton distribution function (PDF) of the gluon given information on the reduced pseudo-Ioffe-time distributions (RpITDs) as calculated from an ensemble with lattice spacing a ≈ 0.09 fm and a pion mass of M π ≈ 310 MeV. The resulting gluon PDF is consistent with phenomenological global fits within uncertainties, particularly in the intermediate-to-high-x region where lattice data are most constraining. A subsequent correlation analysis confirms that the VAIM learns a meaningful latent representation, highlighting the potential of generative AI to bridge lattice QCD and phenomenological extractions within a unified analysis framework.

Gluon parton distribution function↗

DOE OSTI2026

Simulations of a hypersonic turbulent boundary layer over wavy surfaces

Here, we conduct large-eddy simulations of a Mach 5.84 cold wall turbulent boundary layer over one-dimensional wavy walls with varying amplitudes and wavelengths. Across all wall topologies, a series of alternating shock and expansion waves is shown to influence the entire boundary layer, and generate repeating wave patterns in the turbulent stresses, dispersive stresses, and turbulent kinetic energy budget. The series of alternating shocks and expansions imposes repeating adverse and favourable pressure gradients across the wavy wall, and at sufficient wall amplitude, triggers flow separation in the trough of the wave. Flow separation is demonstrated to influence the behaviour of wall pressure fluctuations over the wavy wall. In attached flows, the prominent frequencies are consistent with integral-scale boundary layer turbulence, whereas in separated flows, a two-decade frequency range is present, akin to two-dimensional shock–boundary layer interactions. Counter-rotating streamwise-oriented structures are observed on the windward side of the wave, which diminish over the wave crest. A conditional analysis demonstrates that these structures are present in the upstream boundary layer, and are amplified with increasing wall amplitude. An examination of the Görtler number and probability density function (PDF) of the fluctuating lateral wall shear stress demonstrates the strong correlation between a large Görtler number and growth of the PDF tail density, suggesting that the amplification of the counter-rotating streamwise-oriented structures are linked to centrifugal instabilities in regions of streamline concavity.

boundary layers↗

DOE OSTI2024

Correlated Anion Disorder in Heteroanionic Cubic TiOF 2

Resolving anion configurations in heteroanionic materials is crucial for understanding and controlling their properties. For anion-disordered oxyfluorides, conventional Bragg diffraction cannot fully resolve the anionic structure, necessitating alternative structure determination methods. We have investigated the anionic structure of anion-disordered cubic (ReO 3 -type) TiOF 2 using X-ray pair distribution function (PDF), 19 F MAS NMR analysis, density functional theory (DFT), cluster expansion modeling, and genetic-algorithm structure prediction. Our computational data predict short-range anion ordering in TiOF 2 , characterized by predominant cis-[O 2 F 4 ] titanium coordination, resulting in correlated anion disorder at longer ranges. To validate our predictions, we generated partially disordered supercells using genetic-algorithm structure prediction and computed simulated X-ray PDF data and 19 F MAS NMR spectra, which we compared directly to experimental data. To construct our simulated 19 F NMR spectra, we derived new transformation functions for mapping calculated magnetic shieldings to predicted magnetic chemical shifts in titanium (oxy)fluorides, obtained by fitting DFT-calculated magnetic shieldings to previously published experimental chemical shift data for TiF 4 . We find good agreement between our simulated and experimental data, which supports our computationally predicted structural model and demonstrates the effectiveness of complementary experimental and computational techniques in resolving anionic structure in anion-disordered oxyfluorides. From additional DFT calculations, we predict that increasing anion disorder makes lithium intercalation more favorable by, on average, up to 2 eV, highlighting the significant effect of variations in short-range order on the intercalation properties of anion-disordered materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Structural incorporation into goethite fractionates rare earth elements

Up to 20% of rare earth elements (REEs) in ion adsorption deposits (IADs) are associated with iron oxide minerals, primarily goethite. Often termed “non-extractable”, goethite-hosted REEs are thought to be structurally incorporated into the mineral lattice. The large mismatch in size and charge density between REEs (ionic radius, r = 0.86–1.03 Å) and Fe 3+ (r = 0.65 Å), however, makes direct substitution energetically unfavorable. To determine REE compatibility with and incorporation into goethite on the atomic level, we used X-ray pair distribution function analysis (PDF) and LIII-edge extended X-ray absorption fine structure (EXAFS) spectroscopy. The compatibility of REEs with goethite and the precursor ferrihydrite (FH) is Lu ≥ Yb ≫ Dy > Nd for both phases. Nd and Dy primarily form secondary amorphous phases, with <30% Nd and Dy incorporated into goethite. In the FH precursor at pH 6.8, Yb and Lu assumed a local REE-OOH like structure with next nearest neighbor Fe. The Yb, Lu, and Nd-FH samples were also matured at ambient conditions for 100 days; despite the presence of only ∼5% goethite, Lu and Yb were 42% and 100% in goethite-like structural environments, respectively, whereas the PDF and EXAFS of Nd showed little evidence of any incorporation. Using ab initio molecular dynamics (AIMD) to model the EXAFS, we determined the presence of protonated Fe vacancies, edge-sharing with structural Lu and Yb, likely helped accommodate these REEs in the goethite structure. Incorporation into Fe oxyhydroxides thus potentially fractionates the REEs during weathering associated with formation of lateritic and ion adsorption deposits.

Mergelsberg, Sebastian T. [Pacific Northwest Natio↗