Search NASA⌕ Search

SEARCH · Search NASA

Results for “reference values”

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

Quantitative infrared spectroscopy of aerosols: Mie theory modeling with experimental validation

Infrared spectroscopy is a well-established method for identifying solid, liquid, and gas-phase chemicals. Accurate infrared spectroscopic analysis requires reference libraries where library endmembers reflect all optical phenomena contributing to the observed spectra. Traditional spectral libraries most often contain molecular-based absorption spectra, but these do not account for the complex scattering effects that become significant when measuring aerosols. In this work, we combine the laboratory-derived, wavelength-dependent complex optical vectors ( n / k ) of liquid dioctyl sebacate (DOS) with Mie scattering theory and the Beer–Lambert law to generate synthetic infrared transmission spectra of aerosolized DOS. Additionally, we record experimental infrared transmission spectra using an FTIR spectrometer coupled to a simple aerosol chamber filled with a quantified number size distribution of aerosolized DOS. The modeled and measured spectra show strong agreement, with Mie scattering effects clearly altering the overall spectral shape as well as the positions and profiles of absorption features. The results demonstrate that synthetic spectra generated from n / k values can reliably capture aerosol-specific spectral behavior and thus serve as a foundation for building scalable, physics-based aerosol reference libraries to enable infrared spectroscopic detection of aerosols.

Salcido, Jessica M. O. [Pacific Northwest National↗
DOE OSTI2024

Composition Quantification of SiGeSn Alloys Through Time-of-Flight Secondary Ion Mass Spectrometry: Calibration Methodologies and Validation With Atom Probe Tomography

Time-of-Flight Secondary Ion Mass Spectrometry (ToF-SIMS) is a powerful technique for elemental compositional analysis and depth profiling of materials. However, it encounters the problem of matrix effects that hinder its application. In this work, we introduce a pioneering ToF-SIMS calibration method tailored for SixGeySnz ternary alloys. SixGe1-x and Ge1-zSnz binary alloys with known compositions are used as calibration reference samples. Through a systematic SIMS quantification study of SiGe and GeSn binary alloys, we unveil a linear correlation between secondary ion intensity ratio and composition ratio for both SiGe and GeSn binary alloys, effectively mitigating the matrix effects. Extracted relative sensitivity factor (RSF) value from SixGe1-x (0.07 < x < 0.83) and Ge1-zSnz (0.066 < z < 0.183) binary alloys are subsequently applied to those of SixGeySnz (0.011 < x < 0.113, 0.863 < y < 0.935 and 0.023 < z < 0.103) ternary alloys for elemental compositions quantification. These values are cross-checked by Atom Probe Tomography (APT) analysis, an indication of the great accuracy and reliability of as-developed ToF-SIMS calibration process. Furthermore, the proposed method and its reference sample selection strategy in this work provide a low-cost as well as simple-to-follow calibration route for SiGeSn composition analysis, thus driving the development of next-generation multifunctional SiGeSn-related semiconductor devices.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Technical and Economic Assessment and Gap Analysis of Advanced Nuclear Reactor Integration with a Reference Methanol Synthesis Plant

Efforts continue to identify the most-economic methods to decarbonize several sectors of the United States (U.S.) economy. Industrial processes such as synfuel synthesis and high value commodity chemicals rely heavily on energy-dense and easily stored and transported fossil fuels, which power and feed their operations. Steam methane reforming (SMR) is a widely used process for producing methanol. In this process, methane (CH 4 ) from natural gas (NG) reacts with steam (H 2 O) over a catalyst at high temperatures (700-1,000°C) to produce syngas, a mixture of hydrogen (H 2 ) and carbon monoxide (CO). The syngas is then converted into methanol (CH 3 OH) through a second catalytic reaction. This method is known for being an efficient and commonly employed pathway for industrial methanol production. The high-temperature heat needed for SMR, which is currently used in the natural-gas-to-methanol process, cannot be supplied by small modular nuclear reactor (SMNR) direct heating; the temperatures required for the SMR process exceed those of the main steam produced by near-market high-temperature gas reactors (HTGRs). For the conventional methanol process, this leaves possible nuclear-integration opportunities that include: (1) blending nuclear hydrogen into the SMR NG fuel, or (2) assessing alternative synthesis routes leveraging nuclear capabilities and steam electrolysis outputs. In the reference methanol plant, SMR provides the methanol-synthesis reactor with H 2 and co. In Case (2), the state-of-the-art reverse water gas shift (RWGS) pathway achieves the same, sourcing carbon from an industrial CO 2 source.

08 HYDROGEN↗
DOE OSTI2025

Technical and Economic Assessment and Gap Analysis of Advanced Nuclear Reactor Integration with a Reference Methanol Synthesis Plant

Efforts continue to identify the most-economic methods to decarbonize several sectors of the United States (U.S.) economy. Industrial processes such as synfuel synthesis and high value commodity chemicals rely heavily on energy-dense and easily stored and transported fossil fuels, which power and feed their operations. Steam methane reforming (SMR) is a widely used process for producing methanol. In this process, methane (CH 4 ) from natural gas (NG) reacts with steam (H 2 O) over a catalyst at high temperatures (700°1,000°C) to produce syngas, a mixture of hydrogen (H 2 ) and carbon monoxide (CO). The syngas is then converted into methanol (CH 3 OH) through a second catalytic reaction. This method is known for being an efficient and commonly employed pathway for industrial methanol production. The high-temperature heat needed for SMR, which is currently used in the natural-gas-to-methanol process, cannot be supplied by small modular nuclear reactor (SMNR) direct heating; the temperatures required for the SMR process exceed those of the main steam produced by near-market high-temperature gas reactors (HTGRs). For the conventional methanol process, this leaves possible nuclear-integration opportunities that include: (1) blending nuclear hydrogen into the SMR NG fuel, or (2) assessing alternative synthesis routes leveraging nuclear capabilities and steam electrolysis outputs. In the reference methanol plant, SMR provides the methanol-synthesis reactor with H 2 and co. In Case (2), the state-of-the-art reverse water gas shift (RWGS) pathway achieves the same, sourcing carbon from an industrial CO 2 source.

08 - HYDROGEN↗
DOE OSTI2024

Transgenic Sugarcane–Oilcane: An Alternative Feedstock for the Production of Drop-in Fuel and Value-Added Bioproducts

The utilization of plants and other agricultural produce can partly offset petroleum dependency for energy requirements and can potentially provide sustainable solutions to global environmental problems. To this end, synthetic biology has shown great potential in developing transgenic bioenergy grasses such as sugarcane, sorghum, miscanthus, and energy cane that hyperaccumulate energy-rich lipid molecules in their vegetative tissues, such as leaves, stems, and roots. These perennial high-biomass transgenic C4 grasses are not targeted to grow on prime agricultural land and can be dedicatedly used to produce biofuels and other value-added bioproducts. Recently, sugarcane has been metabolically engineered to sequestrate carbon from juice towards biosynthesis of lipid molecules in the vegetative tissues. Transgenic sugarcane is referred to as “oilcane.” Transgenic sugarcane–oilcane has improved energy density due to an elevated lipid content in the vegetative tissues. Transgenic bioenergy crops cater to both cellulosic sugars and vegetative lipids, hence resulting in higher biofuel yield (biodiesel/renewable diesel and bioethanol) per unit area of cultivable land. In this book chapter, the genetic engineering of sugarcane and its bioprocessing are discussed to illustrate its development and use as an alternative feedstock for the production of biofuels (bio-jet fuel, biodiesel, and bioethanol) and value-added bioproducts.

Maitra, Shraddha↗
DOE OSTI2026

Hanford Site Post-Fire Vegetation Monitoring Report for Calendar Year 2025

The Hanford Site and the surrounding greater Columbia Basin is comprised of an expanse of shrub-steppe habitats that provides exceptional ecological value to plants and animals. The U.S. Department of Energy, Hanford Field Office (HFO)-managed portion of the Hanford Site, referred to from herein as Central Hanford, has experienced repeated wildfires resulting in profound, often adverse, effects on the landscape. This report summarizes these fires and fire monitoring efforts that occurred in calendar year (CY) 2025 and analyzes the post-fire vegetation recovery in restored areas to inform future post-fire response actions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Central Hanford Ecological Integrity Assessments: Final Report

The Hanford Site is comprised of an expanse of shrub-steppe habitats that provides exceptional ecological value to plants and animals located on the site and in the surrounding greater Columbia Basin. The U.S. Department of Energy, Hanford Field Office (HFO)-managed portion of the Hanford Site, referred to herein as Central Hanford, has been the focus of various ecological monitoring efforts, such as vegetation monitoring. The scope and goals of vegetation surveys have varied greatly since the Hanford Site was established, but studies have documented a rapidly changing landscape, making it clear that routine vegetation monitoring is integral to understanding ecological changes and preserving the ecological value of the Hanford Site. A new vegetation monitoring effort was initiated in calendar year (CY) 2023 and continued through CY 2025 using methods based on ecological integrity assessments (EIA) (NHR-2024-04, Field Manual for Applying Rapid Ecological Integrity Assessments in Upland Plant Communities of Washington State ) developed by the Washington Natural Heritage Program (WNHP), a division of the Department of Natural Resources (DNR), and NatureServe.® The methods were modified and supplemented to meet monitoring goals at the Hanford Site. The EIA monitoring effort consisted of field surveys to evaluate vegetation and soil conditions for upland habitats throughout Central Hanford. Vegetation cover estimates were used to score a variety of metrics for vegetation condition. This report summarizes monitoring methods, monitoring results, and provides management recommendations.

54 ENVIRONMENTAL SCIENCES↗
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

High‐Speed Slot‐Die Coating with Donor‐Priority Rapid Aggregation Kinetics for Improved Morphology and Efficiency in Ecofriendly Organic Solar Cells

Abstract Solution‐processable organic solar cells (OSCs) represent a promising renewable photovoltaic technology with significant potential for eco‐compatible production. While high power conversion efficiencies (PCEs) have been achieved in OSCs, scaling this technology for high‐throughput manufacturing remains challenging. Key reason lies in the lack of efficient control strategies for the complex and long‐duration morphology evolution during high‐speed coating process with ecofriendly solvents. Here, a donor‐priority rapid aggregation process (DP‐RAP) scheme is proposed to solve this issue by adjusting the aggregation kinetics of donor and acceptor components. DP‐RAP enables blends with a nanoscale fiber network structure and favorable crystallinity, which contributes to balanced carrier transport and reduced recombination losses. As a result, the PCE is improved from 14.3% (reference) to 17.4% (DP‐RAP) for ultra‐high speed coated PM6:BTP‐eC9 devices in atmosphere, which is one of the highest values for non‐halogenated solvent‐processed solar cells at coating speeds of 500 mm s −1 . Moreover, the DP‐RAP based devices remain a stable PCE of approximately 17.4% across a broad range of coating speeds (20–500 mm s −1 ), illustrating its tolerance to the varied manufacturing conditions. This work highlights a promising avenue for the high‐speed, ecofriendly production of efficient OSCs, pushing the boundaries of practical manufacturing in renewable energy technologies.

Chemistry↗
DOE OSTI2024

Utah FORGE: 2024 Discrete Fracture Network Model Data

The Utah FORGE 2024 Discrete Fracture Network (DFN) Model dataset provides a set of files representing discrete fracture network modeling for the FORGE site near Milford, Utah. The dataset includes four distinct DFN model file sets, each corresponding to different time frames and modeling approaches in 2024. These models characterize both natural and induced fractures in the geothermal reservoir, which consists of crystalline granitic and metamorphic rock approximately 8,000 feet below the ground surface. The dataset includes a reference DFN model from February 2024 that incorporates planar fractures and well trajectories, as well as upscaled permeability, porosity, compressibility, and storage values on specified grids. Additionally, there are models based on new microseismic (MEQ) data from May and July 2024, including fracture planes fitted to the latest MEQ catalog datasets, tensile fractures from hydraulic stimulation, and an alternative connected DFN for modeling purposes. Coordinate data is provided in both global and local frames, with detailed instructions on the transformations used to align with principal stress orientations. The dataset also includes notes and calculation files for estimating fracture sizes and differences between various fracture sets. There are subfolders for Global Coordinates and Local Coordinates. To move from the global to the local coordinate frame, fractures and wells were a) rotated 20 degrees counterclockwise looking down about the global point (335376.400482041, 4263189.99998761, 250.093546450195) to better align with the principal stresses; and b) translated by (-335408.68, -4263010.9, 1150). Upscaled permeability values using the _XYZ suffix show directions with respect to the global XYZ coordinate frame, while those using the _IJK suffix are aligned with local coordinate frame.

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

Core Model Proposal 401: Ukraine as an independent region in GCAM

The goal of this core model proposal (CMP) is to break out Ukraine from the Europe_Eastern region. This work aims to establish Ukraine as an independent region in the GCAM core (region 14) while moving Belarus and Moldova to region 15 (Europe_Non_EU). We have: 1) Updated several mappings to recode region 14 (formerly Europe_Eastern) as Ukraine and moved Belarus and Moldova to region 15 (Europe_Non_EU); 2) Updated several assumptions in the raw data files which provide information by region to reframe Ukraine as the 14th region, including coefficients, base year values, share weight interpolation values and rules, pipeline networks for gas trade, elasticities, shares, etc. 3) Changed documentation and in-code comments at several places referring to fixed 32 regions in GCAM to indicate that GCAM can have any number of regions; 4) Updated code base in gcamdata to dynamically process data for Ukraine given special cases.

Global Change Analysis Model (GCAM)↗
DOE OSTI2026

A Material Frame: Hard Recoils from Slow Force Carriers

We modify the Standard Model by making the longitudinal component of the photon physical, propagating with small speed $c_L \ll 1$. This modification breaks gauge invariance and Lorentz symmetry, defining a preferred reference frame. No new fields or scales are introduced, and the photon mass is protected by a Galilean higher-form symmetry. The theory has a smooth $c_L\to0$ limit: the new modes decouple from standard matter and ordinary gauge theory predictions are recovered. Experimental constraints are severe. Charged particles Cherenkov-emit these slow modes, experiencing hard recoils against the preferred frame. The recoil rate decreases with $c_L$, but the momentum transfer remains of the order of the particle's momentum. Dark matter detectors, sensitive to keV-scale nuclear recoils, imply an order-of-magnitude bound $c_L \lesssim 10^{-60}$, tens of orders of magnitude stronger than speed-difference bounds. At such values, $c_L$ is better understood as a recoil-rate factor rather than a physically relevant speed. The slow modes, practically fixed in space, make the reference frame a material medium capable of absorbing momentum from charged particles.

Serra, Francesco [Johns Hopkins U.]↗
DOE OSTI2026

PIPS-2025: an updated practical pressure scale for the large-volume presses

In 1997, AIRAPT recommended a set of room-temperature pressure reference points (PRPs) for the large-volume press (LVP) community, known as the Practical International Pressure Scale (PIPS-97). With recent AIRAPT-recommended ruby scale Ruby2020 for diamond-anvil cell (DAC) experiments, it is now necessary to re-examine PIPS-97 for consistency with the DAC pressure scale. Here, we reconstruct equations of state (EOSs) for NaCl and Au that are explicitly tied to Ruby2020 and recalculate the PRPs using these EOSs. The revised PRP values are given. GaAs is removed from the PRP list because of the complexity of its phase transitions. The ZnS value, although retained, requires further investigation. Overall, the revised pressures are systematically higher than PIPS-97 by up to 4% at 35 GPa. We propose that this updated set of PRPs, referred to as PIPS-2025, be adopted for LVP experiments to improve pressure consistency between the LVP and DAC communities.

Large volume press↗
DOE OSTI2024

A semi–automatic analytical methodology for characterizing the energy consumption of MRI systems using load duration curves

Background and purpose: Magnetic resonance imaging (MRI) scanners are a major contributor to greenhouse gas emissions from the healthcare sector, and efforts to improve energy efficiency and reduce energy consumption rely on quantification of the characteristics of energy consumption. The purpose of this work was to develop a semi-automatic analytical methodology for the characterization of the energy consumption of MRI systems using only the load duration curve (LDC). LDCs are a fundamental tool used across various fields to analyze and understand the behavior of loads over time. Methods: An electric current transformer sensor and data logger were installed on two 3T MRI scanners from two vendors, termed M1 (outpatient scanner) and M2 (inpatient/emergency scanner). Data was collected for 1 month (7/11/2023 to 8/11/2023). Active power was calculated, assuming a balanced three-phase system, using the average current measured across all three phases, a 480 V reference voltage for both machines, and vendor-provided power factors. An LDC was constructed for each system by sorting the active power values in descending order and computing the cumulative time (in units of percentage) for each data point. The first derivative of the LDC was then computed (LDC’), smoothed by convolution with a window function (sLDC’), and used to detect transitions between different system modes including (in descending power levels): scan, prepared-to-scan, idle, low-power, and off. The final, segmented LDC was used to measure time (% total time), total energy (kWh), and mean power (kW) for each system mode on both scanners. The method was validated by comparing mean power values, computed using the segmented 1-month LDC, for each nonproductive system mode (i.e., prepared-to-scan, idle, lower-power, and off) against power levels measured after a deliberate system shutdown was performed for each scanner (1 day worth of data). Results: The validation revealed differences in mean power values <1.4% for all nonproductive modes and both scanners. In the scan system mode, the mean power values ranged from 29.8 to 37.2 kW and the total energy consumed for 1 month ranged from 11 106 to 14 466 kWh depending on the scanner. Over the course of 1 month, the portion of time the scanners were in nonproductive modes ranged from 76% to 80% across scanners and the nonproductive energy consumption ranged from 8010 to 6722 kWh depending on the scanner. The M1 (outpatient) scanner consumed 99.9 and 183.9 kWh/day in idle mode for weekdays and weekends, respectively, because the scanner spent 23% more time proportionally in idle mode on the weekends. Conclusions: A semi-automatic method for quantifying energy consumption characteristics of MRI scanners was introduced and validated. This method is relatively simple to implement as it requires only power data from the scanners and avoids the technical challenges associated with extracting and processing scanner log files. Finally, the methodology enables quantitative evaluation of the power, time, and energy characteristics of MRI scanners in scan and nonproductive system modes, providing baseline data and the capability of identifying potential opportunities for enhancing the energy efficiency of MRI scanners.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

Findings of the African Combustion Aerosol Collaborative Intercomparison Analysis (ACACIA) Pilot Project to Understand the Optical Properties of Biomass Burning Smoke

Africa is a critical source of biomass burning (BB) aerosols, and its importance is increasing. The African Combustion Aerosol Collaborative Intercomparison Analysis (ACACIA) Pilot Project set to optically characterize BB aerosol generated from sub-Saharan African fuels. We used a photoacoustic spectrometer as a reference instrument to determine the multiple-scattering correction factor C λ for an AE33 aethalometer at three wavelengths, which produced weighted mean values of C 370 =3.69, C 470 =5.65, and C 520 =6.39. C λ increased with wavelength and C 370 was statistically independent of the others, suggesting a single C λ is insufficient, especially in BB scenarios. While a dependence of C λ on burning state was not found, C λ was shown to strongly relate to particle single scattering albedo (SSA, ω). When Cλ was plotted against SSA, values slowly rose at low SSA values, followed by a sharp rise around an SSA of ∼ 0.9; indicating a larger correction needed for less absorbing aerosol. A number of functions operating on either SSA or C λ were explored and the best function was -C λ /(1-C λ )=Aω+B. This is an important parametrization of C λ specifically geared towards BB aerosol from African fuels under different aging states, and is of particular importance for future field work in that continent. An Ångström matrix plot shows that African BB aerosol can have values more akin to dust, which demonstrates that these fuels are distinct in their wavelength dependence from more typical BB aerosol. Lastly, we examined the mass extinction and absorbance cross sections for BB aerosol generated for the same fuels with two different tube furnace setups. Not only is this combustion method flexible, it was found to be reproducible between labs.

47 OTHER INSTRUMENTATION↗
DOE OSTI2026

Neural network potentials with effective charge separation for non-equilibrium dynamics of ionic solids: a ZnO case study

Developing neural network potentials (NNPs) accurate under non-equilibrium dynamics is challenging, as such systems require extensive sampling beyond equilibrium phases. Here we construct high-fidelity NNPs for zinc oxide (ZnO), a polymorphic ionic solid, using density functional theory (DFT) reference data. To efficiently capture transitional configurations, we combine enhanced-sampling molecular dynamics with empirical potentials, data distillation, and pretraining on short-range atomic energies (A-Train), followed by transfer learning with DFT-relabeled datasets. This hierarchical approach improves transferability across polymorphs and stress states. We further introduce effective charge separation, treating long-range Coulombic terms analytically while short-range residual interactions are learned by the NNP. The optimal effective charges fall in the range 0.5–1.0 q e , consistent with dielectric-screened values derived from formal charges but distinct from Bader estimates. Motivated by this observation, we propose a simple data-driven protocol in which effective charges are optimized by comparing DFT reference energies with explicit Coulomb calculations, without additional NNP training. This strategy improves accuracy and transferability in DFT-level predictions of energies, forces, and stress. Together, these results provide a practical charge-selection framework for robust NNP development in ionic solids, enabling reliable simulation of polymorphic phase transformations and non-equilibrium dynamics.

Chemistry↗
DOE OSTI2024

Pressure Transducer Measurement Variability in Deep Wells Screened Across the Water Table

Abstract Automated water level measurements collected using vented pressure transducers in deep wells screened across the water table may exhibit a greater response to barometric pressure changes than the true water level. The cause was hypothesized to be disequilibrium in barometric pressure between the wellbores and land surface due to air exchange with the deep vadose zone. In this study, vented and nonvented pressure transducers were installed and operated simultaneously in two deep wells screened across the water table. A vent tube open to the atmosphere at land surface allowed for barometric compensation of the vented transducers. Two nonvented transducers were installed in each well, one submerged in the water and one above the water surface. The difference in readings allowed for barometric compensation. Manual measurements were also collected. It was confirmed that measurements from the vented transducers exhibited greater variability in response to barometric pressure changes than the nonvented transducers and manual measurements. Comparison of the downhole barometric pressure measurements to values from a nearby meteorology station showed the response in the wells to changes in barometric pressure was time‐lagged and attenuated. Thus, the reference pressure from land surface supplied to the vented transducers was not representative of the air pressure within the wells. This caused fluctuations of the transducer readings in response to barometric pressure changes to be greater than the true water level change. This issue can be resolved by the use of nonvented pressure transducers.

McDonald, John P.↗