Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Chemical initiation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Utah FORGE 5-2557: Trixial Shear and Wellbore Breakout Tests Report and Data

Included are a report and data from triaxial shear tests on simulated Sierra White fault gouge and a borehole stimulation experiment. This serves as a progress report for the Role of Fluid and Temperature in Fracture Mechanics and Coupled Thermo-Hydro-Mechanical-Chemical (THMC) Processes for Enhanced Geothermal Systems project by Purdue University (Utah FORGE R&D project 5-2557). Tests were done with Sierra White granite (SWG) to investigate slip mechanisms and to evaluate the fracture initiation and propagation from a geothermal borehole. Raw data is provided alongside a report on all experimental testing protocols, experimental design, and results.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

Modeling the Effects of Microwave Heating on Densification in Chemical Vapor Infiltration

Microwave heating has great potential to accelerate the synthesis of ceramic matrix composites (CMCs) by a process called Chemical Vapor Infiltration (CVI). In CVI, reactive gases ingress a porous preform and undergo chemical transformation to deposit solid ceramic phase within the pores at high temperature thus, densifying the preform. However, the competing effects of chemical kinetics and gas transport are known to result in non-uniform densification as the outer surfaces of the preform experience faster depositions compared to the core. Achieving spatial temperature control plays a key role improving the densification quality. Microwave heating can potentially create temperature inversion such that the core of the preform is hotter than the outer surface and subsequently, lead to improved densification while keeping the manufacturing times and costs low. In the present work, a computational modeling strategy has been developed that accounts for the key physical phenomena responsible for densification of porous preforms using microwave heating. Specifically, a chemical kinetics model has been formulated for Silicon Carbide (SiC) deposition from MTS/H2 precursor. The model is implemented in a pore-resolved reactive transport solver, Quilt. The CVI simulations are performed for several conditions. Initially, parameterized temperature control is imposed to identify optimum conditions for good densification quality at a fraction of processing time. Further, simulations of microwave heating of porous SiC preforms are performed using OpenFOAM. A strategy to achieve and enhance temperature inversion is identified. The resulting temperature profiles are used in the pore-resolved densification simulations to analyze the densification behavior. It is observed that the temperature inversion achieved by microwave heating leads to significant improvements in densification quality and at the same time, keeps the manufacturing time low.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Electrodeposition of Tungsten using Hydrotropic Agents

The current work sought to electrodeposit tungsten from water-based solutions. The work’s initial hypothesis was that methoxide reducing agents could be used to generate urea anions that could enable tungsten electrodeposition. Unfortunately, this hypothesis was found to be incorrect. However, the work led to the understanding that chemical reducing agents not only enable the electrodeposition of refractory metals like rhenium from water-based solutions, but that the key to plating tungsten in the future involves chemical reduction followed by stabilization with proper ligands. The work resulted in a manuscript under review at Inorganic Chemistry Communications on the discovered of L-histidine as a suitable reducing agent from rhenium electrodeposition. Rhenium-tungsten alloys with 4% tungsten were deposited. A technical advance was filed for the rhenium chemistry and parts were delivered to an internal Sandia customer that used the chemistry.

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

Ultrathin Boron Growth onto Nanodiamond Surfaces via Electrophilic Boron Precursors

Diamond as a templating substrate is largely unexplored, and the unique properties of diamond, including its large bandgap, thermal conductance, and lack of cytotoxicity, makes it versatile in emergent technologies in medicine and quantum sensing. Surface termination of an inert diamond substrate and its chemical reactivity are key in generating new bonds for nucleation and growth of an overlayer material. Oxidized high-pressure high temperature (HPHT) nanodiamonds (NDs) are largely terminated by alcohols that act as nucleophiles to initiate covalent bond formation when an electrophilic reactant is available. In this work, we demonstrate a templated synthesis of ultrathin boron on ND surfaces using trigonal boron compounds. Boron trichloride (BCl 3 ), boron tribromide (BBr 3 ), and borane (BH 3 ) were found to react with ND substrates at room temperature in inert conditions. BBr 3 and BCl 3 were highly reactive with the diamond surface, and sheet-like structures were produced and verified with electron microscopy. Surface-sensitive spectroscopies were used to probe the molecular and atomic structure of the ND constructs’ surface, and quantification showed the boron shell was less than 1 nm thick after 1–24 h reactions. Observation of the reaction supports a self-terminating mechanism, similar to atomic layer deposition growth, and is likely due to the quenching of alcohols on the diamond surface. X-ray absorption spectroscopy revealed that boron-termination generated midgap electronic states that were originally predicted by density functional theory (DFT) several years ago. DFT also predicted a negative electron surface, which has yet to be confirmed experimentally here. The boron-diamond nanostructures were found to aggregate in dichloromethane and were dispersed in various solvents and characterized with dynamic light scattering for future cell imaging or cancer therapy applications using boron neutron capture therapy (BNCT). The unique templating mechanism based on nucleophilic alcohols and electrophilic trigonal precursors allows for covalent bond formation and will be of interest to researchers using diamond for quantum sensing, additive manufacturing, BNCT, and potentially as an electron emitter.

36 MATERIALS SCIENCE↗
NASA NTRS1971

Orbiting propellant depot safety. Volume 1: Management summary report

Orbital missions are considered which require the use of orbiting vehicles, e.g., a cislunar shuttle, that is either chemically or nuclear propelled, functioning as orbit-to-orbit shuttles combined with space tugs for servicing earth-orbiting payloads. Such vehicles may be stationed in a low earth orbit from which they would initiate and terminate flights. The flight frequency of these vehicles indicates that large quantities of propellants will have to be delivered to them in orbit. Orbiting propellant depots, in both geocentric and selenocentric orbits, are being considered as candidate methods of making the required propellants readily available. An assessment of the potential safety hazards associated with the operation of such a depot is summarized.

Source record↗
NASA NTRS1987

The use of airborne imaging spectrometer data to determine experimentally induced variation in coniferous canopy chemistry

It was experimentally determined whether induced differences in forest canopy chemical composition can be detected using data from the Airborne Imaging Spectrometer (AIS). Treatments were applied to an even-aged forest of Douglas fir trees. Work to date has stressed wet chemical analysis of foilage samples and correction of AIS data. Plot treatments were successful in providing a range of foliar N2 concentrations. Much time was spent investigating and correcting problems with the raw AIS data. Initial problems with groups of drop out lines in the AIS data were traced to the tape recorder and the tape drive. Custom adjustment of the tape drive led to recovery of most missing lines. Remaining individual drop out lines were replaced using average of adjacent lines. Application of a notch filter to the Fourier transform of the image in each band satisfactorily removed vertical striping. The aspect ratio was corrected by resampling the image in the line direction using nearest neighbor interpolation.

Swanberg, Nancy A.↗
NASA NTRS1995

Sensitivity analysis of a model of CO2 exchange in tundra ecosystems by the adjoint method

A model of net primary production (NPP), decomposition, and nitrogen cycling in tundra ecosystems has been developed. The adjoint technique is used to study the sensitivity of the computed annual net CO2 flux to perturbation in initial conditions, climatic inputs, and model's main parameters describing current seasonal CO2 exchange in wet sedge tundra at Barrow, Alaska. The results show that net CO2 flux is most sensitive to parameters characterizing litter chemical composition and more sensitive to decomposition parameters than to NPP parameters. This underlines the fact that in nutrient-limited ecosystems, decomposition drives net CO2 exchange by controlling mineralization of main nutrients. The results also indicate that the short-term (1 year) response of wet sedge tundra to CO2-induced warming is a significant increase in CO2 emission, creating a positive feedback to atmosphreic CO2 accumulation. However, a cloudiness increase during the same year can severely alter this response and lead to either a slight decrease or a strong increase in emitted CO2, depending on its exact timing. These results demonstrate that the adjoint method is well suited to study systems encountering regime changes, as a single run of the adjoint model provides sensitivities of the net CO2 flux to perturbations in all parameters and variables at any time of the year. Moreover, it is shown that large errors due to the presence of thresholds can be avoided by first delimiting the range of applicability of the adjoint results.

Waelbroek, C.↗
NASA NTRS1999

Autonomous Agents: The Origins and Co-Evolution of Reproducing Molecular Systems

The central aim of this award concerned an investigation into, and adequate formulation of, the concept of an "autonomous agent." If we consider a bacterium swimming upstream in a glucose gradient, we are willing to say of the bacterium that it is going to get food. That is, we are willing, and do, describe the bacterium as acting on its own behalf in an environment. All free living cells are, in this sense, autonomous agents. But the bacterium is "just" a set of molecules. We define an autonomous agent as a physical system able to act on its own behalf in an environment, then ask, "What must a physical system be to be an autonomous agent?" The tentative definition for a molecular autonomous agent is that it must be self-reproducing and carry out at least one thermodynamic work cycle. The work carried out in this grant involved, among other features, the development of a detailed model of a molecular autonomous agent, and study of the kinetics of this system. In particular, a molecular autonomous agent must, by the above tentative definition, not only reproduce, but must carry out at least one work cycle. I took, as a simple example of a self-reproducing molecular system, the single-stranded DNA hexamer 3'CCGCGG5' which can line up and ligate its two complementary trimers, 5'CCG3' and 5'CGG3'. But the two ligated trimers constitute the same molecular sequence in the 3' to 5' direction as the initial hexamer, hence this system is autocatalytic. On the other hand the above system is not yet an autonomous agent. At the minimum, autonomous agents, as I have defined them, are a new class of chemical reaction network. At a maximum, they may constitute a proper definition of life itself.

Kauffman, Stuart↗
NASA NTRS2017

Initial SAM Calibration Gas Experiments on Mars: Quadrupole Mass Spectrometer Results and Implications

The Sample Analysis at Mars (SAM) instrument suite of the Mars Science Laboratory (MSL) Curiosity rover is equipped to analyze both martian atmospheric gases and volatiles released by pyrolysis of solid surface materials, with target measurements including chemical and isotopic composition (Mahaffy et al., 2012). To facilitate assessment of instrument performance and validation of results obtained on Mars, SAM houses a calibration cell containing CO2, Ar, N2, Xe, and several fluorinated hydrocarbon compounds (Franz et al., 2014; Mahaffy et al., 2012). This report describes the first two experiments utilizing this calibration cell on Mars and gives results from analysis of data acquired with the SAM Quadrupole Mass Spectrometer (QMS). These data support the accuracy of isotope ratios obtained with the QMS (Conrad et al., 2016; Mahaffy et al., 2013) and provide ground-truth for reassessment of analytical constants required for atmospheric measurements, which were reported in previous contributions (Franz et al., 2015, 2014). The most significant implication of the QMS data involves reinterpretation of pre-launch contamination previously believed to affect only CO abundance measurements (Franz et al., 2015) to affect N2 abundances, as well. The corresponding adjustment to the N2 calibration constant presented here brings the atmospheric volume mixing ratios for Ar and N2 retrieved by SAM into closer agreement with those reported by the Viking mission (Owen et al., 1977; Oyama and Berdahl, 1977).

Sample Analysis at Mars Investigation↗
NASA NTRSDate not supplied

Reversibly Assembled Microparticles for Sustained Applications on the Moon and Mars

A sustained presence in any extra-terrestrial environment will require the capability to generate materials, goods, and ultimately vehicles, construction supplies and habitats to be tenable. As an example, NASA, in collaboration with Made In Space, Inc., recently reported the results of extensive research toward enabling3D printing in space. The extensive characterization and comparison of 3D printed articles suggested that, overall, articles generated in the International Space Station microgravity environment were comparable to those generated on Earth. The benefit of in-space manufacturing can be increased through utilization of recyclable materials. A recent analysis of a hypothetical 1,100 day round trip mission to Mars determined that significant mass savings and increase in mission probability of success could be achieved through in-space manufacturing with recyclable materials. NASA’s project ESPUR (Enabling Sustained Presence Using Recyclables), through support from NASA Langley’s Innovative Research and Development Fund, is investigating novel, polymer-coated, epoxy microparticle systems as an enabling technology to realize in-space manufacturing using recyclable feedstocks. Building from previously reported results, further progress toward synthesis of the polymer coating containing click chemical functionalities will be described here. Finite elemental analysis (FEA) has been utilized, in light of scale-free/small world considerations in an effort to identify compositions most likely to yield robust macroscopic geometries. Initial results will be described here. Finally, potential mass savings and life cycle of this recyclable materials technology will be discussed.

Click Chemistry↗
DOE OSTI2023

Excited-State Dynamics during Primary C–I Homolysis in Acetyl Iodide Revealed by Ultrafast Core-Level Spectroscopy

In typical carbonyl-containing molecules, bond dissociation events follow initial excitation to $nπ_{C=O}$$^*$ states. However, in acetyl iodide, the iodine atom gives rise to electronic states with mixed $nπ_{C=O}$$^*$ and $nπ_{C–I}$$^*$ character, leading to complex excited-state dynamics, ultimately resulting in dissociation. Using ultrafast extreme ultraviolet (XUV) transient absorption spectroscopy and quantum chemical calculations, we present an investigation of the primary photodissociation dynamics of acetyl iodide via time-resolved spectroscopy of core-to-valence transitions of the I atom after 266 nm excitation. The probed I 4d-to-valence transitions show features that evolve on sub-100-fs time scales, reporting on excited-state wavepacket evolution during dissociation. These features subsequently evolve to yield spectral signatures corresponding to free iodine atoms in their spin–orbit ground and excited states with a branching ratio of 1.1:1 following dissociation of the C–I bond. Calculations of the valence excitation spectrum via equation-of-motion coupled cluster with single and double substitutions (EOM-CCSD) show that initial excited states are of spin-mixed character. From the initially pumped spin-mixed state, we use a combination of time-dependent density functional theory (TDDFT)-driven nonadiabatic ab initio molecular dynamics and EOM-CCSD calculations of the N$_{4,5}$ edge to reveal a sharp inflection point in the transient XUV signal that corresponds to rapid C–I homolysis. Here, by examining the molecular orbitals involved in the core-level excitations at and around this inflection point, we are able to piece together a detailed picture of C–I bond photolysis in which d → σ* transitions give way to d → p excitations as the bond dissociates. We also report theoretical predictions of short-lived, weak 4d → 5d transitions in acetyl iodide, validated by weak bleaching in the experimental transient XUV spectra. This joint experimental–theoretical effort has thus unraveled the detailed electronic structure and dynamics of a strongly spin–orbit coupled system.

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

NbTa_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-tantalum (Nb-Ta). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach \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. For Ta, p semi-core states as valence were chosen. 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 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 ternary alloys is AXBYCZ, where A, B, and C refer to the constituents, and X, Y, and Z 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. Subdirectories with name case-N, where N ranges between 1 and 10 (extremes included) contain all the files listed above except 'vaspout.h5'. Subdirectories with name case-N, where N ranges between 41 and 60 (extremes included), contain a duplicate copy of the files listed above except for KPOINTS. The names of the duplicate files end with -bis, and correspond to a second VASP calculation that has converged to a different optimized geometry. 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. and Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. review B 47, 558 (1993). (2) Kresse, G. and 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. and 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. and 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. and Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. review b 59, 1758 (1999)

36 MATERIALS SCIENCE↗
DOE OSTI2023

TaV_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-vanadium (Tanu-V). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the element V. For Ta, p semi-core states as valence were chosen. 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. Subdirectories with name case-N, where N ranges between 1 and 10 (extremes included) contain all the files listed above except 'vaspout.h5'. Subdirectories with name case-N, where N ranges between 41 and 60 (extremes included), contain a duplicate copy of the files listed above except for KPOINTS. The names of the duplicate files end with -bis, and correspond to a second VASP calculation that has converged to a different optimized geometry. 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. and Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. review B 47, 558 (1993). (2) Kresse, G. and 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. and Furthmuller, 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. and Furthmuller, 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. and Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. review b 59, 1758 (1999)

36 MATERIALS SCIENCE↗
DOE OSTI2023

NbV_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-vanadium (Nb-V). 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 V. 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. Subdirectories with name case-N, where N ranges between 1 and 10 (extremes included) contain all the files listed above except 'vaspout.h5'. Subdirectories with name case-N, where N ranges between 41 and 60 (extremes included), contain a duplicate copy of the files listed above except for KPOINTS. The names of the duplicate files end with -bis, and correspond to a second VASP calculation that has converged to a different optimized geometry. 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. and Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. review B 47, 558 (1993). (2) Kresse, G. and 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. and 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. and 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. and Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. review b 59, 1758 (1999)

36 MATERIALS SCIENCE↗
DOE OSTI2023

VIPIR: A High-Throughput Drop-Weight Impact Instrument for Imaging the Initiation and Propagation of Reactions in Energetic Materials

Characterizing the handling safety and sensitivity of explosives has been a challenging area of study for over 60 years. Historically one of the most accessible and widely utilized experiments has been the drop-weight impact test, which involves dropping a weight on a small sample sandwiched between two anvils. Because this experiment generally only utilizes sound thresholds to determine whether or not a sample reacted, the physical and chemical properties governing sensitivity remain convolved. Better understanding of chemical and material characteristics is needed to give the chemistry and engineering communities a predictive tool to determine the handling sensitivity of explosives prior to pursuing expensive and potentially hazardous synthesis and formulation operations. Here, we are developing a high throughput drop tower instrument capable of imaging the deformation and flow of energetic materials during impact and the resulting thermal ignition and propagation events. This instrument is based on previous designs but has been modified for higher throughput and tailorable modifications in the future. Herein, we present key design features that enable high-speed visible and thermal imaging of explosive initiation by sub-shock impacts, as well as preliminary results in which ignition sites were observed in an erythritol tetranitrate sample.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Experimental and computational investigation into the hydrodynamics and chemical dynamics of laser ablation aluminum plasmas

Laser ablation plasma chemistry is governed by a complex interplay between hydrodynamic plasma-gas mixing processes, thermodynamics, and rapid high-temperature chemical reactions. Here, in this work, we investigate the gas-phase oxidation chemistry of ns-laser ablation aluminum plasmas in air using optical spectroscopy combined with advanced multi-physics modeling. Experimental measurements demonstrate the formation of AlO in the plasma plume as early as 1 μs while computational results reveal that several Al x O y species are distributed in the periphery of the plume at even earlier times (<20 ns) in the presence of large temperature gradients and strong shockwaves. Interactions with the ablation crater during rapid plume expansion are shown to initiate vortex formation, followed by mixing dynamics that work to pull AlO into the vortices to react with gas-phase Al to form Al 2 O. Oxygen and several aluminum oxides are simultaneously pulled up through the stem of the fireball, encouraging further intermixing between reacting species and enhanced molecular formation. This work concludes that chemical dynamics in laser ablation plasmas is driven by diffusion processes, concentration gradients, and plume hydrodynamics while strong shockwaves generated during laser ablation do not impede chemical reactions.

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

Research Progress and Perspectives on Pre‐Sodiation Strategies for Sodium‐Ion Batteries

Sodium‐ion batteries (SIBs) with abundant elements have garnered significant attention from researches as a promise compensation to lithium‐ion batteries (LIBs). However, the large‐scale commercial application of SIBs is partially hindered by the limited initial coulombic efficiency (ICE) due to the irreversible formation of solid electrolyte interphase (SEI) and intercalation into the defects in the anode. Similar to pre‐lithiation techniques, pre‐sodiation approaches are considered to be one of the most direct and effective way to compensate for the loss of active sodium at the anode side of SIBs during the initial cycle. In this context, additional sodium ions are pre‐injected to the cathode/anode material by chemical/electrochemical methods, aiming to improve battery span life and energy density. Here, this review delves into the necessity and impact of pre‐sodiation techniques, compiling the latest research progress, for instance, self‐sacrificing cathode additives, over‐sodiated cathode materials, direct contact and solution chemical pre‐sodiation. Notably, the research mechanisms underlying solution chemical pre‐sodiation are highlighted. This comprehensive overview aims to foster a deeper understanding of the pre‐sodiation techniques and expects to provide guidance for realizing the commercial application of high energy density sodium‐ion batteries.

25 ENERGY STORAGE↗
DOE OSTI2024

Weathering Incongruence in Mountainous Mediterranean Climates Recorded by Stream Lithium Isotope Ratios

Lithium isotope ratios (δ 7 Li) of rivers are increasingly serving as a diagnostic of the balance between chemical and physical weathering contributions to overall landscape denudation rates. Here, we show that intermediate weathering intensities and highly enriched stream δ 7 Li values typically associated with lowland floodplains can also describe small upland watersheds subject to cool, wet climates. This behavior is revealed by stream δ 7 Li between +22.4 and +23.5‰ within a Critical Zone observatory located in the Cévennes region of southern France, where dilute stream solute concentrations and significant atmospheric deposition otherwise mask evidence of incongruence. The water-rock reaction pathways underlying this behavior are quantified through a multicomponent, isotope-enabled reactive transport model. Using geochemical characterization of soil profiles, bedrock, and long-term stream samples as constraints, we evolve the simulation from an initially unweathered granite to a steady state weathering profile which reflects the balance between (a) fluid infiltration and drainage and (b) bedrock uplift and soil erosion. Enriched stream δ 7 Li occurs because Li is strongly incorporated into actively precipitating secondary clay phases beyond what prior laboratory experiments have suggested. Chemical weathering incongruence is maintained despite relatively slow reaction rates and moderate clay accumulation due to a combination of two factors. First, reactive primary mineral phases persist across the weathering profile and effectively “shield” the secondary clays from resolubilization due to their greater solubility. Second, the clays accumulating in the near-surface profile are relatively mature weathering byproducts. These factors promote characteristically low total dissolved solute export from the catchment despite significant input of exogenous dust.

54 ENVIRONMENTAL SCIENCES↗