Search NASA⌕ Search

SEARCH · Search NASA

Results for “material point method”

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.

268 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

The corrosion behavior of borosilicate glass in the presence of cementitious waste forms

Borosilicate glasses are widely used for radioactive waste disposal due to their ability to incorporate a variety of contaminants and radionuclides while exhibiting high durability in various disposal scenarios. This research evaluated the dissolution of borosilicate glass using both single-pass-flow-through (ASTM C1662-18) and product consistency test (ASTM C1285-21) methods with different solutions, including a cementitious-contacted water (called grout-contacted, GC, from this point) and solutions with varying levels of dissolved cementitious species such as Si, Ca, Al. Here, the results indicated that the presence of Ca plays a crucial role in suppressing glass corrosion, as evidenced by the slower normalized dissolution rates, which were one order of magnitude lower for boron and two orders of magnitude lower for rhenium, observed in both Ca-amended and GC solutions compared to the pH 12 buffer solution. This effect is attributed to the formation of a dense, low-porosity, and strongly bonded calcium silicate hydrate (CSH) layer on the glass surface, which implies that a glass corrosion process is influenced by ion exchange involving alkali ions Na + , K + , Ca 2+ , and hydrogen-containing species. A small number of glass particles treated in the GC solution showed minor corrosion pits in the form of shallow craters with an average diameter of approximately 500 μm. This observation is correlated with a significant reduction, 2000 to 3000 times lower, in the cumulative volume of glass pores, indicating that smaller pore voids were “sealed” in the presence of Ca 2+ ions, likely attributed to the formation of CSH precipitation or other corrosion products such as calcium carbonate saturated from the grout solution. These findings suggest that the presence of dissolved Ca in the GC solution can slow down the dissolution of borosilicate glass, contrary to the expected trend of higher dissolution rates resulting from exposure to high alkaline and thus higher pH solutions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Thermal conductivity suppression in uranium-doped thorium dioxide due to phonon-spin interactions

In this work, impact of low level of uranium (U) atom substitution on thermal conductivity of thorium dioxide (ThO 2 ) is investigated. ThO 2 is an electronic insulator with a wide optical band-gap and no unpaired electrons whose thermal transport is governed by phonons. U-substitution introduces unpaired f-electrons resulting in paramagnetic behavior of U-ThO 2 at room temperature, which significantly suppresses its thermal conductivity. A single crystal of U-ThO 2 with graded composition of U is grown using a hydrothermal synthesis method, and thermal conductivity measurements are performed in regions with uniform composition of U at levels of 0%, 6%, 9% and 16%. Measured thermal conductivity profiles over 77–300 K temperature range are analyzed using an analytical expression for phonon-mediated thermal transport based on Klemens-Callaway model. Temperature dependent thermal conductivity is found to deviate significantly from the Rayleigh scattering trend expected for a simple substitutional point defect with a small perturbation to mass and interatomic forces. With the resonant scattering term, observed large suppression of thermal conductivity at low temperatures can be closely reproduced. Additionally, the extracted phonon-spin coupling constants imply a nonlinear relation of phonon-spin interaction intensity with respect to U doping percentage. Our study reveals how phonon-spin scattering contributed by unpaired f-electrons in U atoms influences thermal transport in the U-ThO 2 system.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Advanced Computational Modeling of High-Level Waste Vitrification at the Hanford Site

The U.S. Department of Energy (DOE) has selected vitrification for stabilizing legacy tank waste at the Hanford site, where radioactive waste from plutonium production was historically stored in underground tanks. This waste will be separated into low-activity waste (LAW) and high-level waste (HLW) fractions and processed at the Waste Treatment and Immobilization Plant (WTP). At WTP, glass melters are used for the vitrification of radioactive tank waste, transforming it into a stable borosilicate glass form for safe long-term storage. The melter vessel is constructed from highly durable and heat-resistant materials, where the vitrification process occurs. The main regions that are modeled are the melt pool, plenum, cold cap, riser/discharge chamber, and surrounding structure with insulation layers. Forced convection induced by air bubblers at the base of the melter ensure uniform temperature distribution and provide heat to the cold cap layer. The cold cap is a region of reacting batch feed that floats on top of the molten glass and is where the batch-to-glass reactions occur. Joule heating provided by electrodes mounted along the vertical walls of the melter and immersed directly in the glass, generates the necessary heat for the net endothermic conversion processes that occur in the cold cap. The high temperatures, radioactivity, and opaque nature of the glass prevent direct observation inside the melters. Therefore, computational models are essential for providing insight into factors that affect melter throughput. Thermocouples in the plenum provide operators with plenum temperature measurements. Operational adjustments include bubbling rate, voltage supplied to the electrodes, feed adjustments, and glass removal rate. Different computational fluid dynamics (CFD) models have been developed, each serving a specific purpose. There are CFD models of different scale melters, as well as models that capture the two-phase flow interfaces of rising bubbles in the molten glass or models with a simplified molten glass region so that the surrounding structure and plenum can be feasibly incorporated. Pilot-scale melter models have been developed to serve as validation of the methods employed in the simulation of the full-scale WTP melters. Models incorporating resolved bubbling are used to develop momentum source terms to implement into a single phase, multi-region, steady-state flow model that is being validated by measured process parameters such as glass production rate, voltage, input power, plenum temperatures, etc. The resolved bubbling model uses the multiphase volume of fluid approach to model the system with a high-resolution interface capturing scheme to maintain sharp interfaces between the molten glass and the air phase. The suite of CFD models is continually being improved to incorporate more realistic physics and achieve faster turnaround time. For example, an incremental controller is implemented to automatically adjust electrode voltage within the simulation to a molten glass set point temperature of 1150°C. Newer models feature improved meshes to ensure conformal meshes between regions and eliminate unnecessary mesh refinement in areas that are not of interest (such as boundary layers in offgas ports). Instead of explicitly modeling the structural, refractory, and insulation layers of the melter, a thermal resistance approach is used with published correlations used for boundary conditions. The development of robust and efficient CFD models will be instrumental in enabling the WTP to successfully fulfill its mission of safely stabilizing legacy nuclear waste.

12 - MGMT OF RADIOACTIVE AND NON-RADIOACTIVE WASTE↗
DOE OSTI2025

Innovative use of CaO in combination with amino acid salt to convert CO 2 as CaCO 3 nanoparticles under mild pH and low temperature

The increasing demand for sustainable CO 2 management has driven the development of innovative methods that can convert point source CO 2 into value-added products. In this study, CaO in combination with amino acid salt was used to convert CO 2 into CaCO 3 nanoparticles. Different from the conventional method where CO 2 diffusion presents a major challenge in reacting with CaO to form CaCO 3 , amino acid salt solvent was applied to absorb CO 2 first and then rapidly reacted with CaO to form CaCO 3 nanoparticles (∼50 nm) at a low temperature (e.g., 60 °C). Our experiments showed that at a glycine (Gly)/NaOH ratio of 2:1 or 3:1, the solution pH values during the CO 2 absorption and conversion were about 8–9 at 60 °C, while at a ratio of 1:1, the solution pH values were about 9–11; without Gly, the solution pH values were about 12. Gly-optimized solvent substantially reduced corrosion risk to reactors. In addition, the use of amino acid (i.e., Gly) led to much smaller CaCO 3 particles, distinctly different chemical phases, and fundamentally different chemical reactions. Moreover, in the presence of Gly, the solution pH was completely reversed and the solution was regenerated for cyclic use when CaO was added. The solvent was recyclable and reusable, highlighting the cost-effectiveness and sustainability of this approach. The Gly-modulated CaCO 3 nanoparticles may have significant potential for industrial applications in the biomedicine, construction, plastics, and rubber industries.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Aboveground Rather Than Belowground Productivity Drives Variability in Miscanthus × giganteus Net Primary Productivity

Quantifying the carbon (C) uptake of Miscanthus × giganteus ( M × g ) in both aboveground and belowground structures (e.g., net primary productivity (NPP)) and differences among methodological approaches is crucial. Our objectives were to directly measure Mxg NPP and evaluate the effects of nitrogen application, location, and belowground biomass sampling methods. We hypothesize that increased nitrogen application increases the overall NPP of M × g and that quantifying rhizome biomass using excavations will produce the lowest variability between replicates. We collected biomass from mature M × g stands from three locations in Iowa with three nitrogen application rates and one site in Illinois. We destructively sampled at two time points, when rhizome mass is anticipated to be at a minimum (initial) and anticipated to be at its maximum (peak). Biomass was collected from 1 × 1 m quadrats in which one in-clump and one beside-clump cores were collected and then excavated to 30 cm depth to extract all rhizomes. We found that aboveground M × g NPP ranged from 15.4 Mg DM ha –1 year –1 to 36.4 Mg DM ha –1 year–1 and belowground M × g NPP ranged from 4.4 Mg DM ha –1 year –1 to 19.6 Mg DM ha –1 year –1 . M × g NPP varied across sites, fertilization, and calculation assumptions. Aboveground NPP (yield) was on average 68.7% of the total NPP. Root-to-shoot ratios at peak biomass decreased with nitrogen application rate, from an average of 1.9 for 0 N plots to 0.89 for 224 N fertilized plots. There was more variation in core data than from excavations; however, when in-clump and beside-clump cores were averaged together, core and excavation averages were not different. Overall, these results show that the range of mature M × g NPP is driven by aboveground productivity, influenced by nitrogen application and site. Our results provide useful data to constrain agro-ecosystem models and provide crucial insights for future perennial belowground sampling.

Hartman, Theodore [Univ. of Illinois at Urbana-Cha↗
DOE OSTI2026

A systematic analytical framework for multi-source municipal solid waste characterization for energy recovery

Advancing municipal solid waste (MSW) management from disposal-oriented practices toward circular, value-driven systems requires standardized methodologies capable of identifying material composition and resource recoverable potential at the point of generation. Despite extensive research, MSW characterization remains fragmented due to inconsistences in sampling methodologies, waste sorting categories, and temporal coverage across previous studies which limit cross-site comparability, reproducibility, and constrain the reliable evaluation of potential resource recovery pathways. This lack of consistency has hindered the development of a unified framework for MSW characterization and resource assessment. This study introduces a standardized, field-validated protocol for MSW sampling and composition analysis that ensures consistent, traceable data across diverse waste sources. The protocol integrates randomized spatial sampling, systematic material sorting, and controlled subsampling for multi-site and multi-season field campaigns. Validation included MSW collection from residential, grocery, restaurant, and school MSW streams across five U.S. states, including Maryland, Idaho, Virginia, Ohio, and Mississippi, to demonstrate the protocol’s ability to identify source-based composition patterns relevant to resource recovery applications. Grocery and restaurant streams were dominated by food waste and high-moisture organics, while school waste contained higher paper content and residential waste showed greater heterogeneity. Aggregation into energy-relevant fractions highlighted practical recovery pathways via anaerobic digestion or gasification, supporting data-driven planning, policy, and circular economy strategies for sustainable waste management across waste sources.

09 BIOMASS FUELS↗
DOE OSTI2026

Extreme confinement of hydrogen gas within fullerenelike nanoporous carbon

Nanoporous carbons and carbon nanostructures can store hydrogen at cryogenic temperatures but lack the volumetric and gravimetric capacity to be industrially significant. Recent inelastic neutron scattering experiments suggest a highly dense phase of hydrogen at temperatures well above the melting point of solid hydrogen. However, it remains unclear how pore geometry and intermolecular interactions enable these dense phases to exist, with dispersion (van der Waals) or electrostatic/induction suggested to be the key effects in slit and curved pores but their relative contributions have yet to be quantified. In this paper, we perform benchmark electronic structure calculations allowing the interactions between planar and curved aromatic molecules with hydrogen to be accurately determined. Dispersion was found to dominate over electrostatic and inductive effects with some many-body charge transfer (Dobson type-A) effects needed to capture the most highly curved structures. Density functional methods that include type-A many-body effects were found to accurately describe the intermolecular interactions at a fraction of the cost of coupled-cluster simulations and these approaches were used to calculate the energies inside large carbon bowl and slit pores. The interaction energies inside the bowl pores were found to depend on the orientation of the hydrogen molecule. This rotational barrier, modeled as a quantum hindered rotor, could reproduce the peak splitting observed in inelastic neutron scattering experiments, with weak splitting arising from bowl-like fullerene pores and strong splitting from highly confining nanotubelike pores. Increasing the fraction of such curved pores in nanoporous carbons may therefore offer a pathway to enhance their hydrogen-storage capacity. Moreover, the preferential adsorption of ortho hydrogen on nanotubelike pores could enable the storage of high-density hydrogen without the need to remove heat produced during the ortho-para hydrogen conversion.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Woody Feedstock 2022 State of Technology Report

The U.S. Department of Energy promotes production of advanced liquid transportation fuels from lignocellulosic biomass by funding fundamental and applied research that advances the state of technology (SOT). As part of its involvement in this mission, Idaho National Laboratory completes an annual SOT report for n th -plant and 1 st -plant woody biomass feedstock logistics. The purpose of the SOT is to provide the status of feedstock supply system technology development for woody biomass to biofuels relative to technical targets and cost goals from specific design cases, based on data and experimental results. Conventional feedstock supply systems need to be modified to meet the demands of conversion pathways, specifically to have the ability to adjust the quality of the raw biomass materials. Advanced systems incorporate innovative methods of material handling, preprocessing and supply chain configuration. In advanced designs, variability of the raw biomass can be reduced to produce feedstocks of a uniform format, moving toward biomass commoditization. Against this backdrop, the 2022 Woody SOT for low-ash woody feedstocks utilizes feedstock fractionation by incorporating technologies that can separate the biomass into its anatomical fractions (wood, bark, needle, and extrinsic ash) to reduce impurities and attempt to maximize the retention of usable fractions that satisfy downstream quality considerations. By using a series of air classification steps, this strategy can reduce the extrinsic ash in forest residues, separate out a majority of the incoming needles (which can be supplied to alternate markets), and maximize the retention of whitewood in the usable fraction. The fractionated forest residues are then mixed with clean-pine chips in a 50-50 blend to prepare the feedstock for the desired conversion pathway. The n th -plant analysis estimated the delivered cost for the feedstock at $\$$69.23/dry ton (2016$\$$) which represents a $\$$6.64/dry ton decrease compared to the cost estimate of the 2021 Woody SOT supply system for low-ash woody feedstocks. The quality requirements in the 2022 Woody SOT were identical to those of the 2021 Woody SOT at = 1.00 wt % ash and = 50.51 wt% carbon. The cost savings derive primarily from reductions in dry matter losses during air classification. The GHG emissions for the n th -plant analysis were estimated at 178.39 kg CO2e/dry ton compared to 178.71 kg CO 2 e/dry ton in the 2021 Woody SOT, a decrease of 0.32 kg CO2e/dry ton. The small change stems from an increase in emissions attributed to preprocessing and slightly larger savings in emissions from transportation. In the 1 st -plant analysis of the 2022 Woody SOT system, the average throughput was estimated to be approximately 2,128 dry tons/day or 96.51% of the name plate capacity. During the simulation the daily throughput ranged from 1,090 dry tons/day to 2,200 dry tons/day, or 49.43% to 99.75% of the daily nameplate capacity. After the year of operation 722,403 tons of processed feedstock were produced in total without regard to quality considerations (99.64% of the annual nameplate capacity). The variability in throughput was primarily caused by equipment failures in the system. Regular failures, downtime caused by routine maintenance per manufacturer guidelines, contributed to a majority 62.50% of failures and 62.60% of downtime. Failures due to wear were the other cause of disruption within the system, impacting the rotary shear and orbital screen and accounting for 37.50% of the failures and 37.40% of the total downtime. Ultimately the system was on stream for 87.84% during the simulation period, which is only 2.16 percentage points below the nth-plant assumption for on-stream time. The production cost of the system averaged $\$$71.66/dry ton. The costs ranged from a minimum of $\$$71.23/dry ton to a maximum of $\$$2,115.30/dry ton. When dry matter losses (disposed low-quality fractions as well as other losses such as in grinders) were considered the costs increased to an average of $\$$75.11/dry ton with a minimum of $\$$74.69/dry ton and a maximum of $\$$2,136.86/dry ton...

09 BIOMASS FUELS↗
DOE OSTI2026

Pixelated plastic scintillator array manufacturing using fast-, photo-curable resin

Pixelated plastic scintillator arrays can serve as high efficiency and high resolution neutron imaging detectors. Manufacturing these arrays is intensive in both time and labor. This article presents a fabrication method based on additive manufacturing for two-dimensional plastic organic scintillator arrays using a custom-built automated assembly machine and a custom photocurable resin that has significant non-aromatic acrylate oligomer content. The process involves two main stages: fully autonomous production of one-dimensional layered arrays, followed by semi-autonomous cutting and stacking to form two-dimensional pixel arrays. One-dimensional arrays were manufactured at a rate of around 4 layers per hour with minimal defects and tight dimensional tolerances, while two-dimensional arrays up to 7 x 7 pixels and 70 mm in length were completed in approximately 3.5 hours. Final arrays exhibited dimensional deviations of less than 0.5 mm. Two-dimensional arrays read out by a multi-anode photomultiplier tube demonstrated per-pixel position resolution and pulse-shape discrimination, enabling gamma–neutron interaction separation in mixed radiation environments.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Solid State Quantum Refrigeration Superconducting, Absorption and Measurement Based (Final Technical Report)

During this DOE grant, DE-SC0017890, in place for the past six years, all proposed research was carried out and published in peer-reviewed papers, as well as other projects that emerged during the research. In that effort the research team accomplished all proposed research, as well as many closely related research projects discovered and conceived of during the grant. These works included “Efficient Quantum Measurement Engines”, a work published in Physical Review Letters, giving a theory of quantum measurement-based engines, which uses quantum measurement as a resource. These engines are designed to efficiently convert energy from the stochastic quantum measurement process into useful work. Further publications include “Experimental Realization of a Quantum Dot Energy Harvester”, a joint theory and experimental work in collaboration with the group of Charles Smith in Cambridge, UK, as well as long time theoretical collaborators, Rafael Sánchez and Björn Sothmann. This work, featured as an Editor’s Suggestion in Physical Review Letters, realized an earlier theoretical proposal of ours, whereby two resonant tunneling quantum dots are connected to a central electronic cavity that is heated by a hot energy source. We also published “Superconducting Quantum Refrigerator: Breaking and Rejoining Cooper Pairs with Magnetic Field Cycles” a work done in collaboration with experimentalist Francesco Giazotto from ENS Pisa, Italy, which also resulted in a patent. This paper, published in Phys. Rev. Applied, advanced the concept of a cyclic fridge based on the normal/superconducting phase transition together with layered materials separated by tunnel junctions. We also completed the proposed research on a heat transistor, publishing “Thermal transistor and thermometer based on Coulomb-coupled conductors”, carried out as a collaboration between my group and theorists Splettstoesser (Lund U., Sweden), Sothmann (U. Duisburg-Essen, Germany), and Sánchez (U. Autónoma de Madrid, Spain). We carried out an analysis of a quantum coupled to a quantum point contact as a sensitive thermometer and heat transistor. We found the optimal statistical estimator for the temperature and compared it with experiments on the same type of devices. We also investigated autonomous quantum absorption refrigerators using quantum dots to cool by using a very hot thermal reservoir to drive heat between two other reservoirs. In the article “Quantifying the quantum heat contribution from a driven superconducting circuit”, we demonstrated that for a driven superconducting circuit, we showed heat flow provided by a hot source to the qubit can be switched on and off by varying external parameters, the frequency and the intensity of the driving. In the work “Stochastic thermodynamic cycles of a mesoscopic thermoelectric engine”, we reconsidered the autonomous thermoelectric heat engine in terms of underlying cycles. Rather than periodic behavior, the cycles were stochastic in nature. Nevertheless, by undertaking a graph theoretical analysis of the elementary transport processed, great quantitative and qualitative insight could be found. We also considered the quantum measurement process and showed that a quantum version of Maxwell’s demon could be related to the work extraction of a quantum system, closely related to arrow-of-time measures for quantum measurement, as described in our article “Thermodynamics of quantum measurement and Maxwell's demon's arrow of time”. This work was selected in Phys. Rev. A as an Editor’s Suggestion. A recent preprint titled “Cyclic Superconducting Quantum Refrigerators Using Guided Fluxon Propagation” accomplished an important piece of this grant: to propose a new kind of quantum refrigerator using the dynamics of fluxons in a type II superconductor. This invention envisioned a race-track type geometry where fluxons are confined. By applying a gradient of magnetic field together with electric current in a Corbino geometry, the circulating fluxons can actively cool a cold reservoir, realizing a new type of cyclic superconducting refrigerator. We also investigated the possibility of thermal control from different points of view. The application of quantum measurement to the system gives a new kind of control on the system of interest – we have pioneered this approach and shown that measurement can boost the thermal power of quantum engines as described in “Continuous measurement boosted adiabatic quantum thermal machines”. The ability to have heat flows on demand is an outstanding challenge, and we have provided new solutions to this problem in Thermal control across a chain of electronic nanocavities” for a chain of electron cavities using gating voltage control. The control methods using qubit/qubit coupling to create absorption fridges at their most fundamental level have also been developed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Assessment of the Impact of Realistic Sensor Physics and the Integration of Ex-Core Sensors on Reactor Power Synthesis

In the work documented in this report, a weighting function–based core power synthesis method was applied to multiple Monte Carlo N-Particle (MCNP) reactor models, which are informed based on simulated self-powered neutron detector (SPND) responses. The weighting function method used has been coined the point-based iterative (PBI) method. The goal of this application is to assess the impact of considering realistic sensor physics in the generation of the simulated SPND outputs as well as to consider how the synthesis is impacted based on the inclusion of ex-core detectors in the model. The NuScale small modular reactor (SMR) and Westinghouse AP1000 pressurized water reactor (PWR) are the models that served as the testbeds for the assessment of realistic sensor physics; this was achieved by using Geant4 SPND models in comparison with analytical models, such that the effect of electron transport in realistic SPND geometries in the Geant4 model can be understood in terms of synthesis error and convergence time. The comparison was considered for fuel burnup–induced perturbations, for a range of sensor string densities and synthesized power distribution axial fidelities. The Texas A&M Testing, Research, Isotopes, General Atomics Reactor (TAMU TRIGA) reactor MCNP model was used to assess the impact of ex-core sensors; this was done by performing synthesis with and without the ex-core detectors and by quantifying the synthesis error and number of iterations associated with Gaussian-type perturbations in many locations in the core. The TAMU TRIGA model was particularly pertinent for this study because of the interest in future experimental tests with SPNDs in this reactor, as well as the ease of modifying the MCNP model to include ex-core detectors with heterogeneously described response functions. Results from the comparison between the Geant4 and analytical SPND models indicate that similar average and maximum synthesis errors were obtained for burnup-induced perturbations in both the NuScale SMR and the AP1000. This was true for a range of sensor string densities and axial fidelities. However, there were marked differences between both the Geant4 and analytically informed models in terms of the iterations required to converge on the synthesized power distribution. Namely, the Geant4-informed models tended to lead to fewer iterations, except for a few sensor–core configurations that had particularly numerous iterations. Results from the ex-core sensor assessment with the TAMU TRIGA model indicate that the inclusion of ex-core sensors drastically reduces the synthesis error of Gaussian-type perturbations close to the edge of the core, and it slightly reduces synthesis errors for perturbations closer to the center of the core. This was achieved with a minimal increase in computational cost—that is, the number of iterations required for convergence. The errors were identified to be in the same location as the perturbation in the core, indicating that the methodology remains robust for unperturbed regions of the core. A secondary result from this study with the TAMU TRIGA was yielded by analysis of the neutron flux levels in the in-core and ex-core sensor locations of the core; these flux levels indicate that SPNDs could be used as both in-core and ex-core sensors, so long as the emitter material is sensitive to thermal neutrons. The results from these studies provide a quantitative understanding of the importance of considering realistic sensor physics and including ex-core sensors to perform accurate and timely power distribution synthesis of a reactor core.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

NbZr_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-zirconium (Nb-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 \cite{PAW}. 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 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. The calculations have been collected on NERSC-Perlmutter and OLCF-Summit using the VASP 6.3.2. The VASP calculations for every atomic structure have been performed in 2 main steps: 1. Starting from an ideal body-centered-cubic (BCC) structure, geometry optimization with low precision has been executed to perform a preliminary optimization of the atomic structure. The output for this calculations is available in the files 0.CONTCAR, 0.OUTCAR, rlx1.out. 2. Using the atomic structure resulting from the preliminary geometry optimization, a second geometry optimization has been performed using normal precision. The output for this calculations is available in the files CONTCAR, OUTCAR, rlx2.out, vaspout.h5, and vasprun.xml. Cases 1-10 have been run without generating the file 'vaspout.h5'. Every chemical composition sampled across the composition range in the dataset has its own directory. The convention used to name the directories for binary alloys is AXBY, where A and B refer to the constituents, whereas X and Y are positive integers that represent the number of atoms for each constituent and their values still sum up to 128. Each atomic structure associated with a specific chemical composition has its own sub-directory within the directory of the corresponding chemical composition. The sub-directories for each atomic structure for each chemical composition are named 'case-*', where * is a positive integer that spans all the values from 1 through 100, extremes included. The files contained in each sub-directory 'case-*' for each atomic structure are as follows: FILES contained in each subdirectory with name "case-N" where N ranges between 11 and 100, extremes included: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. KPOINTS: input file that specifies the Bloch vectors (k points) used to sample the Brillouin zone 3. 0.POSCAR: input file that defines the atomic structure of a system 4. 0.CONTCAR: output file that provides the atomic positions and cell parameters after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 5. 0.OUTCAR: output file that contains detailed information about the progress of a calculation after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 6. rlx1.out: file with diagnostic information about the execution of the first geometry optimization with precision variable set to PREC=Low in the INCAR file 7. POSCAR: input file that defines the atomic structure of a system after the first geometry optimization has been run at low precision. This represents the input for the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 8. CONTCAR: output file that provides the atomic positions and cell parameters after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 9. OUTCAR: output file that contains detailed information about the progress of a calculation after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 10. rlx2.out: file with diagnostic information about the execution of the second geometry optimization with precision variable set to PREC=Normal in the INCAR file 11. vaspout.h5: hierarchical HDF5 file containing the inputs and outputs of a VASP calculation. To analyze the data in this file we recommend using py4vasp. This file is only produced if the VASP version used is compiled with HDF5 support 12. vasprun.xml: contains similar information to OUTCAR, but in an xml format. 13. CHGCAR: contains the charge density data of a system. This data is crucial for analyzing electronic structures, calculating electrostatic potential, and studying the distribution of charge in a crystal or molecular system FILES contained in each subdirectory with name "case-N" where N ranges between 1 and 10, extremes included: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. KPOINTS: input file that specifies the Bloch vectors (k points) used to sample the Brillouin zone 3. {ID}.POSCAR: input file that defines the atomic structure of a system at the beginning of ID execution of geometry optimization with PREC=LOW 4. {ID}.CONTCAR: output file that provides the atomic positions and cell parameters at the end of ID execution of geometry optimization with PREC=LOW in the INCAR file 5. {ID}.OUTCAR: output file that contains detailed information about the progress of a calculation after the ID execution of geometry optimization that has been run with the precision variable set to PREC=Low in the INCAR file 6. rlx1.{ID}.out: file with diagnostic information about the execution of the ID execution of the geometry optimization with precision variable set to PREC=Low in the INCAR file 7. N{ID}.POSCAR: input file that defines the atomic structure of a system after the geometry optimization run at low precision. This represents the input for the ID execution of the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 8. N{ID}.CONTCAR: output file that provides the atomic positions and cell parameters after the ID execution of the geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 9. N{ID}.OUTCAR: output file that contains detailed information about the progress of a calculation after the ID execution of the geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 10. rlx2.{ID}.out: file with diagnostic information about the ID execution of geometry optimization with precision variable set to PREC=Normal in the INCAR file 11. vaspout.h5: hierarchical HDF5 file containing the inputs and outputs of a VASP calculation. To analyze the data in this file we recommend using py4vasp. This file is only produced if the VASP version used is compiled with HDF5 support 12. vasprun.xml: contains similar information to OUTCAR, but in an xml format. 13. CHGCAR: contains the charge density data of a system. This data is crucial for analyzing electronic structures, calculating electrostatic potential, and studying the distribution of charge in a crystal or molecular system This research is sponsored by the Artificial Intelligence Initiative as part of the Laboratory Directed Research and Development (LDRD) Program of Oak Ridge National Laboratory, managed by UT-Battelle, LLC, for the US Department of Energy under contract DE-AC05-00OR22725. This work used resources of the Oak Ridge Leadership Computing Facility, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC05-00OR22725, under Directorate Discretionary awards MAT025 (Materials Science) and LRN026 (Machine Learning), and INCITE award MAT201. This work also used resources of the National Energy Research Scientific Computing Center, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231, under award ERCAP0025216. REFERENCES (1) Kresse, G. & Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. review B 47, 558 (1993). (2) Kresse, G. & Hafner, J. Ab initio molecular-dynamics simulation of the liquid-metal–amorphous-semiconductor transition in germanium. Phys. Rev. B 49, 14251 (1994) (3) Kresse, G. & Furthmüller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput. materials science 6, 15–50 (1996) (4) Kresse, G. & Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. review B 54, 11169 (1996) (5) Kresse, G. & Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. review b 59, 1758 (1999)

36 MATERIALS SCIENCE↗
DOE OSTI2024

TaZr_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 tantalum-zirconium (Ta-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 \cite{PAW}. 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 Ta 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. The calculations have been collected on NERSC-Perlmutter and OLCF-Summit using the VASP 6.3.2. The VASP calculations for every atomic structure have been performed in 2 main steps: 1. Starting from an ideal body-centered-cubic (BCC) structure, geometry optimization with low precision has been executed to perform a preliminary optimization of the atomic structure. The output for this calculations is available in the files 0.CONTCAR, 0.OUTCAR, rlx1.out. 2. Using the atomic structure resulting from the preliminary geometry optimization, a second geometry optimization has been performed using normal precision. The output for this calculations is available in the files CONTCAR, OUTCAR, rlx2.out, vaspout.h5, and vasprun.xml. Cases 1-10 have been run without generating the file 'vaspout.h5'. Every chemical composition sampled across the composition range in the dataset has its own directory. The convention used to name the directories for binary alloys is AXBY, where A and B refer to the constituents, whereas X and Y are positive integers that represent the number of atoms for each constituent and their values still sum up to 128. Each atomic structure associated with a specific chemical composition has its own sub-directory within the directory of the corresponding chemical composition. The sub-directories for each atomic structure for each chemical composition are named 'case-*', where * is a positive integer that spans all the values from 1 through 100, extremes included. The files contained in each sub-directory 'case-*' for each atomic structure are as follows: FILES contained in each subdirectory with name "case-N" where N ranges between 11 and 80, extremes included: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. KPOINTS: input file that specifies the Bloch vectors (k points) used to sample the Brillouin zone 3. 0.POSCAR: input file that defines the atomic structure of a system 4. 0.CONTCAR: output file that provides the atomic positions and cell parameters after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 5. 0.OUTCAR: output file that contains detailed information about the progress of a calculation after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 6. rlx1.out: file with diagnostic information about the execution of the first geometry optimization with precision variable set to PREC=Low in the INCAR file 7. POSCAR: input file that defines the atomic structure of a system after the first geometry optimization has been run at low precision. This represents the input for the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 8. CONTCAR: output file that provides the atomic positions and cell parameters after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 9. OUTCAR: output file that contains detailed information about the progress of a calculation after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 10. rlx2.out: file with diagnostic information about the execution of the second geometry optimization with precision variable set to PREC=Normal in the INCAR file 11. vaspout.h5: hierarchical HDF5 file containing the inputs and outputs of a VASP calculation. To analyze the data in this file we recommend using py4vasp. This file is only produced if the VASP version used is compiled with HDF5 support 12. vasprun.xml: contains similar information to OUTCAR, but in an xml format. 13. CHGCAR: contains the charge density data of a system. This data is crucial for analyzing electronic structures, calculating electrostatic potential, and studying the distribution of charge in a crystal or molecular system FILES contained in each subdirectory with name "case-N" where N ranges between 1 and 10 and between 81 and 100, extremes included: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. KPOINTS: input file that specifies the Bloch vectors (k points) used to sample the Brillouin zone 3. {ID}.POSCAR: input file that defines the atomic structure of a system at the beginning of ID execution of geometry optimization with PREC=LOW 4. {ID}.CONTCAR: output file that provides the atomic positions and cell parameters at the end of ID execution of geometry optimization with PREC=LOW in the INCAR file 5. {ID}.OUTCAR: output file that contains detailed information about the progress of a calculation after the ID execution of geometry optimization that has been run with the precision variable set to PREC=Low in the INCAR file 6. rlx1.{ID}.out: file with diagnostic information about the execution of the ID execution of the geometry optimization with precision variable set to PREC=Low in the INCAR file 7. N{ID}.POSCAR: input file that defines the atomic structure of a system after the geometry optimization run at low precision. This represents the input for the ID execution of the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 8. N{ID}.CONTCAR: output file that provides the atomic positions and cell parameters after the ID execution of the geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 9. N{ID}.OUTCAR: output file that contains detailed information about the progress of a calculation after the ID execution of the geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 10. rlx2.{ID}.out: file with diagnostic information about the ID execution of geometry optimization with precision variable set to PREC=Normal in the INCAR file 11. vaspout.h5: hierarchical HDF5 file containing the inputs and outputs of a VASP calculation. To analyze the data in this file we recommend using py4vasp. This file is only produced if the VASP version used is compiled with HDF5 support 12. vasprun.xml: contains similar information to OUTCAR, but in an xml format. 13. CHGCAR: contains the charge density data of a system. This data is crucial for analyzing electronic structures, calculating electrostatic potential, and studying the distribution of charge in a crystal or molecular system This research is sponsored by the Artificial Intelligence Initiative as part of the Laboratory Directed Research and Development (LDRD) Program of Oak Ridge National Laboratory, managed by UT-Battelle, LLC, for the US Department of Energy under contract DE-AC05-00OR22725. This work used resources of the Oak Ridge Leadership Computing Facility, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC05-00OR22725, under Directorate Discretionary awards MAT025 (Materials Science) and LRN026 (Machine Learning), and INCITE award MAT201. This work also used resources of the National Energy Research Scientific Computing Center, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231, under award ERCAP0025216. REFERENCES (1) Kresse, G. & Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. review B 47, 558 (1993). (2) Kresse, G. & Hafner, J. Ab initio molecular-dynamics simulation of the liquid-metal–amorphous-semiconductor transition in germanium. Phys. Rev. B 49, 14251 (1994) (3) Kresse, G. & Furthmüller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput. materials science 6, 15–50 (1996) (4) Kresse, G. & Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. review B 54, 11169 (1996) (5) Kresse, G. & Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. review b 59, 1758 (1999)

36 MATERIALS SCIENCE↗
DOE OSTI2025

Observation of Orbital-Selective Dual Modulations in an Anisotropic Antiferromagnetic Kagome Metal TbTi 3 ⁢Bi 4

Orbital selectivity is pivotal in dictating the phase diagrams of multiorbital systems, with prominent examples including the orbital-selective Mott phase and superconductivity. The intercalation of anisotropic layers represents an effective method for enhancing orbital selectivity and thereby shaping the low-energy physics of multiorbital systems. Despite its potential, related experimental studies, especially those elucidating the correlation between orbital selectivity and magnetism, remain limited. In this work, we systematically examine the interplay between orbital selectivity and magnetism in the newly discovered anisotropic kagome TbTi 3 ⁢Bi 4 single crystal, and report the coexistence of orbital-selective dual-band modulations (𝑞 1 ∼ 1/3⁢𝑎*, 𝑞 2 ∼ 0.28⁢𝑏*) within the antiferromagnetic (AFM) state. By combining soft x-ray and vacuum ultraviolet angle-resolved photoemission spectroscopy measurements, neutron powder diffraction, scanning tunneling microscopy, and density-functional-theory calculations, we identify these dual-band reconstructions as manifestations of the AFM order driven by a (approximately 1/3, 0.28, 0) nesting instability of the intercalated Tb 5⁢𝑑 𝑥⁢𝑧 orbitals. These orbital-selective modulations induce unusual momentum-dependent band folding and lead to the emergence of Dirac cones only at the $\bar{M}$ 1 point, signaling a topological phase transition in the AFM state. Importantly, the discovery of orbital-selective (approximately 1/3, 0.28, 0) AFM order offers crucial insights into the mechanism underlying the fractional magnetization plateau in this kagome AFM metal. Our findings not only underscore the essential role of both conducting and localized electrons in determining the magnetic orders of Ln⁢Ti 3 ⁢Bi 4 (Ln = lanthanide) kagome metals but also offer a pathway for manipulating magnetism through selective control of anisotropic electronic structures.

Zhang, Renjie [Shanghai Jiao Tong University (Chin↗
DOE OSTI2025

The UV Photoinduced Ring-Closing Reaction of Cyclopentadiene Probed with Ultrafast Electron Diffraction

Conjugated cyclic organic molecules are common across many fields such as pharmaceuticals, are naturally occurring in biological systems, and are used in synthetic materials. One particular area of interest from a photochemical point of view is the formation of highly strained cyclic organics. We investigate the photoinduced reaction of cyclopentadiene, a five-membered organic ring molecule which can form strained three and four carbon rings after photoexcitation with UV light, with the gas-phase ultrafast electron diffraction instrument at the SLAC MeV-UED facility. Electron diffraction offers a direct probe sensitive to the nuclear geometry during the reaction, allowing for the determination of the distribution of products formed following photoexcitation. We observe the simultaneous formation of the highly strained ring- closed bicyclo[2.1.0]pentene and vibrationally hot cyclopentadiene within the temporal resolution of the experiment, and determine the relative yield of all reaction products. Furthermore, the experimental results are in good agreement with the predictions of trajectory simulations.

Computational Chemistry↗
DOE OSTI2025

Mondo: integrating disease terminology across communities

Precision medicine aims to enhance diagnosis, treatment, and prognosis by integrating multimodal data at the point of care. However, challenges arise due to the vast number of diseases, differing methods of classification, and conflicting terminological coding systems and practices used to represent molecular definitions of disease. This lack of interoperability artificially constrains the potential for diagnosis, clinical decision support, care outcome analysis, as well as data linkage across research domains to support the development or repurposing of therapeutics. There is a clear and pressing need for a unified system for managing disease entities⁠—including identifiers, synonyms, and definitions. To address these issues, we created the Mondo disease ontology—a community-driven, open-source, unified disease classification system that harmonizes diverse terminologies into a consistent, computable framework. Mondo integrates key medical and biomedical terminologies, including Online Mendelian Inheritance in Man (OMIM), Orphanet, Medical Subject Headings (MeSH), National Cancer Institute Thesaurus (NCIt), and more, to provide a comprehensive and accurate representation of disease concepts with fully provenanced and attributed links back to the sources. Mondo can be used as the handle for curation of gene–disease associations utilized in diagnostic applications, research applications such as computational phenotyping, and in clinical coding systems in clinical decision support by pointing the clinician to the numerous knowledge resources linked to the Mondo identifier. Mondo's community-centric approach, stewarded by the Monarch Initiative's expertise in ontologies, ensures that the ontology remains adaptable to the evolving needs of biomedical research and clinical communities, as well as the knowledge providers.

biomedical informatics↗