Search NASA⌕ Search

SEARCH · Search NASA

Results for “STEP FUNCTION”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

A spline-based method to obtain spatially dependent viscosity in confined flows

Coupling chemical physics to continuum theories is a critical step to understanding multi-scale phenomena. This paper will connect non-equilibrium molecular dynamics simulations to a continuum-based Navier-Stokes equation that has relaxed the assumption of spatial uniformity in viscosity. Using a form for viscosity based on spline interpolation, viscosity as a function of position is obtained from the least squares fit of the velocity profile measured from molecular simulations of flow in a nanochannel. Viscosity can vary widely, particularly near the channel boundaries, indicating that uniform viscosity is no longer appropriate. Variations of the viscosity near the channel surfaces imply that considering solution and surface chemistry could be necessary to rigorously understand molecular-scale flows in nanochannels.

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

The safety of magnetic resonance imaging contrast agents

Gadolinium-based contrast agents are increasingly used in clinical practice. While these pharmaceuticals are verified causal agents in nephrogenic systemic fibrosis, there is a growing body of literature supporting their role as causal agents in symptoms associated with gadolinium exposure after intravenous use and encephalopathy following intrathecal administration. Gadolinium-based contrast agents are multidentate organic ligands that strongly bind the metal ion to reduce the toxicity of the metal. The notion that cationic gadolinium dissociates from these chelates and causes the disease is prevalent among patients and providers. We hypothesize that non-ligand-bound (soluble) gadolinium will be exceedingly low in patients. Soluble, ionic gadolinium is not likely to be the initial step in mediating any disease. The Kidney Institute of New Mexico was the first to identify gadolinium-rich nanoparticles in skin and kidney tissues from magnetic resonance imaging contrast agents in rodents. In 2023, they found similar nanoparticles in the kidney cells of humans with normal renal function, likely from contrast agents. We suspect these nanoparticles are the mediators of chronic toxicity from magnetic resonance imaging contrast agents. This article explores associations between gadolinium contrast and adverse health outcomes supported by clinical reports and rodent models.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

HfZr_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys hafnium-zirconium (Hf-Zr). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Hf and Zr. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Physics-informed latent neural operator for real-time predictions of time-dependent parametric PDEs

Deep operator network (DeepONet) has shown significant promise as surrogate models for systems governed by partial differential equations (PDEs), enabling accurate mappings between infinite-dimensional function spaces. However, when applied to systems with high-dimensional input-output mappings arising from large numbers of spatial and temporal collocation points, these models often require heavily overparameterized networks, leading to long training times. Latent DeepONet addresses some of these challenges by introducing a two-step approach: first learning a reduced latent space using a separate model, followed by operator learning within this latent space. While efficient, this method is inherently data-driven and lacks mechanisms for incorporating physical laws, limiting its robustness and generalizability in data-scarce settings. Here, in this work, we propose PI-Latent-NO, a physics-informed latent neural operator framework that integrates governing physics directly into the learning process. Our architecture features two coupled DeepONets trained end-to-end: a Latent-DeepONet that learns a low-dimensional representation of the solution, and a Reconstruction-DeepONet that maps this latent representation back to the physical space. By embedding PDE constraints into the training via automatic differentiation, our method eliminates the need for labeled training data and ensures physics-consistent predictions. The proposed framework is both memory and compute-efficient, exhibiting near-constant scaling with problem size and demonstrating significant speedups over traditional physics-informed operator models. We validate our approach on a range of parametric PDEs, showcasing its accuracy, scalability, and suitability for real-time prediction in complex physical systems.

Latent representations↗
DOE OSTI2025

Ammonium-coordinated exchanger (ACE) functionalized silica sorbents for recovering/removing aqueous anionic contaminants

There are limited studies of functionalized silica anion exchange sorbents used for critical/heavy metal recovery/removal relative to polymeric and other inorganic materials. This work features ammonium-coordinated exchanger (ACE) anion exchange particle sorbents prepared by either acid-washing epoxy-crosslinked polyethylenimine (PEI) hydrogen bonded within/to a silica particle sorbent (two-step method) or reacting a di-chlorinated crosslinker, α,α-dichloro-p-xylene (DPX), with PEI within silica (single-step method). Energy dispersive X-ray spectroscopy (EDS) and infrared spectroscopy confirmed the presence of -NH 2 + ···Cl - and -NH 3 + ···Cl - groups, which removed oxyanionic species –arsenate, selenate, chromate, sulfate, phosphate, and nitrate– plus bromide from ideal solutions, authentic acid mine drainage (AMD), and authentic flue gas desulfurization (FGD) wastewater. Affinity of the anions for ACE varied across single- and mixed-element solutions. However, affinity for CrO 4 2- was among the highest in both cases. Total anion uptake by the optimized ACE, PEI-E3-HCl_1.1, reached 1.2 mmol anion/g-sorb. (0.56 mmol CrO 4 2- /g), or ∼2.3 mmol negative charge/g-sorb. This was close to the 0.52 mmol CrO 4 2- /g of a commercial anion exchange resin. Near-consistent removal of 20–80 % of each anion from FGD during an eight-cycle adsorption-desorption (1 M NaCl) test predicted good ACE viability for testing under practical conditions at larger scale.

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

Fundamental benchmarking of the discharge properties of negative electrodes in lead acid batteries

The unprecedented scale of energy storage deployment needed to allow high penetration of intermittent renewable energy sources into the electrical grid places significant economic cost and performance targets on battery technologies. Lead batteries, owing to their highly abundant and inexpensive raw materials, can be considered an important solution to grid storage as long as they achieve high material utilization, fast recharge rates, and long cycle life. To address the limitations in material utilization, we need to revisit the electrochemical processes during the discharge step and answer the following question: what are the fundamental limits of the discharge reaction? In this work we discuss how well-defined lead metal surfaces with nanometer scale roughness allow us to resolve the accessible discharge capacity for both faradaic and non-faradaic processes as a function of electrolyte concentration and in the presence of traditional additives. The direct connection between PbSO 4 layer morphology and discharge rates gives us important insights on the mechanism of discharge as related to the dissolution and passivation events. In conclusion, this work sets the stage for reimagining the design of lead batteries with implications to grid storage applications.

25 ENERGY STORAGE↗
DOE OSTI2024

Heterogenous catalysis for oxygen tolerant photoredox atom transfer radical polymerization and small-molecule dehalogenation

Heterogeneous photocatalysts (PCs) have garnered attention for their sustainability and cost-effectiveness. Despite the existence of various types of these PCs, their synthesis often involves complex, multi-step procedures and laborious purification. Herein, we propose a simple method for attaching small-molecule photocatalytic species onto crosslinked 3-D polymer networks as insoluble scaffolds to create robust heterogeneous PCs. The highly swellable poly(ethylene glycol)-based ChemMatrix (CM) resin, known for its amphiphilic properties and high functional group loading, facilitated the covalent immobilization of the photoredox dye Eosin Y (EY), but also streamlined functionalization with Ir( III ) complexes. The resulting heterogeneous CM-EY demonstrated efficient photocatalytic performance in open-to-air dual photoredox catalysis of atom transfer radical polymerization (photo-ATRP) under green light. This was confirmed by the well-controlled synthesis of polymers with molecular masses ranging from 20 kDa to 300 kDa and low dispersities. Furthermore, CM-EY exhibited excellent photostability and recyclability over multiple cycles of ATRP. The heterogeneous catalysis of photo-ATRP provided high temporal control and enabled benign conditions for synthesizing protein-polymer hybrids (PPH). When combined with the initiator-modified CM (CM-BIB), CM-EY facilitated the solid-phase synthesis of homopolymers and block copolymers with recyclable performance. However, the coordinatively bound Ir@CM showed decreased catalytic activity and efficiency toward photoredox dehalogenation due to the leaching of active species during recycling. This study highlights the advantages of the covalent linking of catalysts to solid supports over non-covalent interactions, underscoring the potential of functionalized polymer resin as a promising scaffold. Such an approach offers customization and tunability, presenting opportunities for innovation in green chemistry.

Kapil, Kriti↗
DOE OSTI2026

Ultrafast solvent-to-solute proton transfer mediated by intermolecular coherent vibrations

Ultrafast photoinduced excited-state proton transfer (ESPT) plays a crucial role in protecting biomolecules and functional materials from photodamage. However, the influence of solute-solvent interactions on ESPT dynamics remains under active investigation. Here, we present an ultrafast spectroscopic study of ESPT in the photobase 2-(2´-pyridyl)benzimidazole (PBI) in methanol. Ultrafast absorption spectroscopy, supported by quantum chemical calculations, reveals three distinct kinetic steps: (1) a 2.2 ps solvent-to-solute proton transfer, (2) subsequent nonradiative relaxation to the ground state within 31 ps, producing a vibrationally hot ensemble with substantial excess kinetic energy, and (3) equilibration as this energy dissipates into the surrounding solvent bath over 186 ps. Femtosecond-resolved dynamics exhibit oscillatory signals indicative of coherent wavepacket motion on the S 1 potential energy surface. A phase flip in the excited-state absorption maximum confirms this assignment. Fourier analysis resolves two dominant periods (∼117 fs and ∼340 fs), corresponding to in-plane and out-of-plane vibrational modes coupled between PBI and the hydrogen-bonded methanol molecule. The rapid dephasing ( < 300 fs) suggests that the nuclear wavefunction evolves on an anharmonic potential energy surface while traversing the ESPT reaction coordinate.

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

HfTa_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys hafnium-tantalum (Hf-Ta). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Hf and Ta. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect BCC lattice sites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,040 randomized atomic structures over 31 chemical compositions. Further methodological and structural information is contained in the dataset README.txt file.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Clustering and Cliques in Preferential Attachment Random Graphs with Edge Insertion

In this paper, we investigate the global clustering coefficient (a.k.a transitivity) and clique number of graphs generated by a preferential attachment random graph model with an additional feature of allowing edge connections between existing vertices. Specifically, at each time step t, either a new vertex is added with probability f(t), or an edge is added between two existing vertices with probability 1 – f(t). We establish concentration inequalities for the global clustering and clique number of the resulting graphs under the assumption that f(t) is a regularly varying function at infinity with index of regular variation –$\gamma$, where $\gamma$ $\in$ [0, 1). Finally, we also demonstrate an inverse relation between these two statistics: the clique number is essentially the reciprocal of the global clustering coefficient.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

PROACTIVE Focus Area 4: Structured Decision Metrics Analysis (Final Report)

A structured decision metrics analysis was developed as one of the tasks under PROACTIVE’s Focus Area 4 (FA4) and used structured decision metrics for processes, items, and facilities (PIF) to determine the efficacy of an M&V system in meeting treaty goals and technical objectives. The method starts with the functional decomposition of a “treaty” with the M&V system overlaid on it. A Functional Decomposition Rubric (FDR) for each PIF is used to determine a score ranging from 0 (weak) to 4 (strong) for assurance, security, and burden. The individual scores are combined to provide an overall Verification System Score (VSS) which assesses the overall verification system’s suitability given goals; scores for each step of the process are also given. The outcome of the evaluation is to identify needs or modifications to the enterprise model, verification system, or proposed testbed capabilities. VeriScore was used to evaluate FA4’s Spiral 0 exercise; those results are included in this report. The VeriScore and FDR framework can also be used in a non-treaty environment, as any goal/objective/method structure will work, thus expanding its applicability beyond traditional arms control structures.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2025

LDRD conclusion poster - Synthesizing Heterometallic Uranium Single Crystals to Understand the Influence of the Secondary Metals on Uranyl Axial Bond Strength

Understanding how transition metals influenced the chemistry of lanthanide and actinide (f-element) materials is critical for advancing separation technologies, materials design, and coordination chemistry. This project examined how incorporating first-row transition metals affected the structural and spectroscopic properties of f-element coordination polymers. In uranium(VI)-based systems synthesized with 2,6-pyridinedicarboxylic acid (PDC) ligands, single-crystal X-ray diffraction and Raman spectroscopy revealed that the presence of transition metals shortened the uranyl axial bond and induced a blue shift in its symmetric stretching vibration—evidence of increased bond strength. Electronic structure analysis, including Density of States (DOS) calculations using density functional theory (DFT), revealed altered orbital overlaps and highlighted the role of transition metal d-orbitals in modulating bonding. Raman modes were modeled using truncated structural fragments in collaboration with the University of Notre Dame, and although the predicted frequencies were lower than experimental values, they remained within expected ranges. In parallel, similar experiments with cerium (Ce) in the presence of cobalt (Co) and PDC demonstrated multi-step single-crystal-to-single-crystal transformations—behavior not observed in the uranium systems. Initial products included light yellow, orange, and polycrystalline materials. Single-crystal X-ray diffraction studies, conducted in collaboration with the Colorado School of Mines, identified the yellow phase as monometallic Ce(PDC)2(H2O)2·4H2O and the orange phase as heterometallic Ce2Co(PDC)4(H2O)6. After standing in solution for one week, both phases fully transformed into a dark yellow crystalline phase, [Ce3(PDC)5(H2O)8].6(H2O). Remarkably, this transformation was reversible—disturbing the equilibrium by removing some crystals caused reversion to the initial Ce(PDC)2(H2O)2·4H2O phase, highlighting dynamic behavior. All three structures were previously unreported. Solid-state UV-visible and Raman spectroscopy further distinguished these phases, revealing ligand-to-metal charge transfer involving Ce and characteristic d–d transitions from Co(II). The precise mechanism driving these transformations remained unclear; however, pH-dependent experiments confirmed that the transformation did not occur when the pH decreased. Overall, the project demonstrated that transition metals could be employed to tune bonding interactions, structural dimensionality, and optical properties in f-element materials, establishing new pathways for designing functional heterometallic systems. The work resulted in several novel structural discoveries and fostered productive collaborations with the University of Notre Dame and the Colorado School of Mines.

37 - INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL C↗
DOE OSTI2026

Astrocyte Reactivity by Alcohol Dependence in the Central Amygdala

Astrocytes play essential roles in the brain, but prolonged insult can result in the loss or diminishment of homeostatic functions and increased neuroinflammatory response. The central amygdala (CeA) is an important limbic region in reward and stress responsivity as well as in the development of alcohol dependence. Using a multi-omics approach with Aldh1l1-EGFP/Rpl10a mice and the chronic intermittent ethanol – two-bottle choice exposure model, we have characterized the translational response of CeA astrocytes, as well as the proteomic and phosphoproteomic changes in dependent, non-dependent, and naïve mice. We identified astrocyte-specific alterations in neuroimmune functions and antioxidant/oxidative stress pathways in dependent animals while non-dependent EtOH drinking altered cytoskeletal plasticity related pathways. Proteomic analysis showed down-regulation of astrocyte physiological functions in dependent animals while phosphoproteomic analysis identified cytoskeleton remodeling in both dependent and non-dependent animals suggesting the development of a reactive astrocyte phenotype. Astrocyte morphological reconstruction showed increased CeA astrocyte complexity in dependent and non-dependent groups compared to naïve mice. The astrocyte-specific increase in neuroimmune functions, down-regulation of astrocyte homeostatic functions, alteration in protein phosphorylation-mediated cytoskeleton remodeling, and increased astrocyte complexity demonstrate EtOH induced astrocyte activation in the CeA and suggest the presence of both adaptive and maladaptive reactive astrocytes. These findings highlight the roles CeA astrocytes play in the progression of alcohol intake to dependence and are the first step in the identification of novel astrocyte-specific therapeutic approaches for alcohol use disorder through the potentiation of adaptive changes and inhibition of maladaptive changes in astrocytes.

Hashimoto, Joel↗
DOE OSTI2024

Active learning for SNAP interatomic potentials via Bayesian predictive uncertainty

Bayesian inference with a simple Gaussian error model is used to efficiently compute prediction variances for energies, forces, and stresses in the linear SNAP interatomic potential. Here, the prediction variance is shown to have a strong correlation with the absolute error over approximately 24 orders of magnitude. Using this prediction variance, an active learning algorithm is constructed to iteratively train a potential by selecting the structures with the most uncertain properties from a pool of candidate structures. The relative importance of the energy, force, and stress errors in the objective function is shown to have a strong impact upon the trajectory of their respective net error metrics when running the active learning algorithm. Batched training of different batch sizes is also tested against singular structure updates, and it is found that batches can be used to significantly reduce the number of retraining steps required with only minor impact on the active learning trajectory.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

On the numerical sensitivity of cellular automata grain structure predictions to large thermal gradients and cooling rates

Cellular automata (CA) models of as-solidified grain structure, originally developed and applied to casting, have become a common means of predicting grain structure resulting from Additive Manufacturing (AM) processes. The majority of these models are based on the decentered octahedron approach, which attempts to correct for the effect of grid anisotropy on the prediction of competitive solidification of dendritic grains. However, AM solidification occurs under cooling rates ($\dot{T}$) and thermal gradients (G) that are orders of magnitude larger than those encountered in casting, and no systematic investigation on the effect of the CA model cell size (Δx) and time step (Δt) on AM microstructure predictions has been performed. Here, in this study, such an investigation is first performed via simulation of individual grains of various crystallographic orientations with a fixed, unidirectional G, showing that CA prediction of the steady-state undercooling matched the expected values based on the interfacial response function at small G and deviated from the expected values at large G. Simulation of competitive growth of multiple grains showed a weakening of the predicted texture as G and Δx became large. Simulation of solidification under AM conditions, where G and $\dot{T}$ vary spatially across the melt pools, showed that not only does grain selection weaken and deviate from expectations at large Δx, but grains with crystallographic $\langle$100$\rangle$ aligned with the grid directions are more adversely affected by the temperature field discontinuities than grains with other crystallographic orientations. Despite the fact that the exact grain competition results depended on Δt, the overall texture development was notably less sensitive to Δt than Δx, provided that a reasonable value of Δt is selected based on the ratio of Δx to the maximum local solidification velocity in the simulation domain. Finally, from the directional solidification and AM simulation results, an analysis of computational cost compared to simulation resolution is performed based on an equation derived to quantify the relatively inaccuracy in grain selection based on the model and temperature field inputs. From this analysis, it is concluded that there is a need for algorithmic improvements to improve CA grain competition accuracy for large G processing conditions as sufficiently small Δx to resolve the necessary competition is intractable for many AM processing conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Dual atom catalysts for rapid electrochemical reduction of CO to ethylene

Strong CO adsorption and facile CO dimerization are the key challenges in electrochemical CO2 reduction towards multi-carbon (C2+) products. We recently showed that CoPc immobilized on a single-walled carbon nanotube can selectively reduce CO2 to methanol. This is enabled through molecular strain, which dramatically improves the CO adsorption energy to CoPc, which in turn facilitates methanol formation. We now examine the extended Phthalocyanine (PcEx) dual atom catalyst (DAC), which is intrinsically strained and contains two catalyst centers, making it a candidate for reducing CO to C2+ products. Using Quantum Mechanics (QM), we screened 20 elements embedded in the PcEx, seeking catalysts with weak hydrogen binding, strong CO binding, and facile CO dimerization. We identi>ied Fe, Ru, Co, and Ir as the best performers and subsequently evaluated the entire CO to C2H4 mechanism (9 steps) using each of these elements as catalysts. In terms of limiting potential and overall exergonicity, we identi>ied CoPcEx as the best catalyst, followed by IrPcEx. We then examined the full CO to C2H4 mechanism on the bimetallic IrCoPcEx catalyst using grand canonical QM to obtain the reaction energetics as a function of applied potential. We conclude that the bimetallic IrCoPcEx is most promising for ef>iciently converting CO to ethylene.

Musgrave, Charles B.↗
DOE OSTI2025

Feasibility of carbon foam-based sorbents for the abatement of gaseous mercury and iodine

The U.S. Department of Energy Hanford Site in Washington State is in the process of commissioning the Waste Treatment and Immobilization Plant to process a portion of the 54 million gallons of radioactive and chemical waste from cold war weapon production. Technologies for the capture of volatile species of concern are still being assessed, and new methods and materials are developed as operational and flowsheet mission risks are identified. One such area still being assessed is the abatement efficacy of the Carbon Adsorber units to retain gaseous mercury and 129 I released during processing. It is challenging to predict the mercury chemistry due to the variability of the feed, and different methods/materials are required for the capture of gaseous Hg 0 and Hg II compounds. In this study, the feasibility of using developmental carbon foam (CF) sorbents for the capture of iodine and mercury was assessed using static and dynamic flow testing and compared against a commercially available sorbent, BATII-37. Both CF and CF functionalized with bismuth particles (CF-Bi) chemisorbed iodine, and CF-Bi had similar mercury capture performance to BATII-37 in dynamic flow tests. While species loading concentrations were measurable, limitations in achieving a mass balance prevented a full evaluation of capture efficacy. Nonetheless, the results serve as an important first step in demonstrating the potential for simultaneous iodine and mercury capture.

Hanford site↗
DOE OSTI2026

NbTi_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys niobium-titanium (Nb-Ti). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP), https://github.com/ORNL/ScalableWorkflow_VASP_Calculations. Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Nb and Ti. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions. The calculations have been collected on Air Force HPC11 cluster using the VASP 6.5.1. Additional methodology and file structure information is available in the dataset README.txt file.

36 MATERIALS SCIENCE↗