Search NASA⌕ Search

SEARCH · Search NASA

Results for “O(N) model”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

The ANI-1ccx and ANI-1x data sets, coupled-cluster and density functional theory properties for molecules

Abstract Maximum diversification of data is a central theme in building generalized and accurate machine learning (ML) models. In chemistry, ML has been used to develop models for predicting molecular properties, for example quantum mechanics (QM) calculated potential energy surfaces and atomic charge models. The ANI-1x and ANI-1ccx ML-based general-purpose potentials for organic molecules were developed through active learning; an automated data diversification process. Here, we describe the ANI-1x and ANI-1ccx data sets. To demonstrate data diversity, we visualize it with a dimensionality reduction scheme, and contrast against existing data sets. The ANI-1x data set contains multiple QM properties from 5 M density functional theory calculations, while the ANI-1ccx data set contains 500 k data points obtained with an accurate CCSD(T)/CBS extrapolation. Approximately 14 million CPU core-hours were expended to generate this data. Multiple QM calculated properties for the chemical elements C, H, N, and O are provided: energies, atomic forces, multipole moments, atomic charges, etc. We provide this data to the community to aid research and development of ML models for chemistry.

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

The influence of thermochemistry on the reactivity of propane, the pentane isomers and $n$-heptane in the low temperature regime

In this study, the influence of thermochemistry on the reactivity of fuels at low temperatures (600–1000 K) is studied here. Specifically, the effect of different sets of thermochemistry on chemical model predictions is explored, where various sets are calculated at different levels of theory in addition to recently updated group additivity values. Experimentally measured ignition delay times for propane, the pentane isomers and n-heptane are simulated using NUIGMech1.2 and replacing the thermochemistry of the low-temperature species with the calculated values. For propane, three different thermochemistry sets were calculated, namely CCSD(T)-F12/TZ-F12//B2PLYP-D3/TZ//B2PLYP-D3/TZ (QM1), CCSD(T)-F12/TZ-F12//B2PLYP-D3/TZ//ωB97X-D/TZ (QM2) and B2PLYP-D3/TZ/ωB97X-D/6-31G*//ωB97X-D/6-31G* (QM3). The QM2 results provide parameters to optimize new group additivity (NGA) values which are used to calculate the fourth set of thermochemistry. The model predictions using these four sets are compared to those using NUIGMech1.2 for propane. As the QM1 and QM2 calculations are expensive, the thermochemistry calculated from the QM3 and NGA calculations are used in the pentane isomer and n-heptane models. For all of the models, it is found that the thermochemistry of the species involved in the low-temperature reaction sequence ($RH$, $\dot{R}$, $RO_2H$, $R\dot{O}_2$, $\dot{Q}OOH$ and $\dot{O}_2QOOH$ species) significantly affect fuel reactivity. The NGA values were developed based on all of these species except $\dot{O}_2QOOH$ radicals. The thermochemistry of $\dot{O}_2QOOH$ species cannot be accurately calculated with the NGA representations due to the importance of non-next-nearest neighbor interactions of –OOH substitution. Further development of the NGA method to capture such interactions is in progress. Overall, the model developed using the NGA thermochemistry shows better agreement with experimental data than the model using thermochemistry from affordable and prominent QM methods, such as QM3. Based on the results presented for propane, the pentane isomers and n-heptane, the thermochemistry calculated using the NGA method can be used to model the oxidation of higher order hydrocarbons at low temperatures.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

Ion trap long-range XY model for quantum state transfer and optimal spatial search

Abstract Linear ion trap chains are a promising platform for quantum computation and simulation. The XY model with long-range interactions can be implemented with a single side-band Mølmer–Sørensen scheme, giving interactions that decay as 1 / r α , where α parameterises the interaction range. Lower α leads to longer range interactions, allowing faster long-range gate operations for quantum computing. However, decreasing α causes an increased generation of coherent phonons and appears to dephase the effective XY interaction model. We characterise and show how to correct for this effect completely, allowing lower α interactions to be coherently implemented. Ion trap chains are thus shown to be a viable platform for spatial quantum search in optimal O ( N ) time, for N ions. Finally, we introduce a O ( N ) quantum state transfer protocol, with a qubit encoding that maintains a high fidelity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

A combined theoretical and experimental study of the phase coexistence and morphotropic boundaries in ferroelectric-antiferroelectric-antiferrodistortive multiferroics

The physical nature of the ferroelectric (FE), ferrielectric (FEI) and antiferroelectric (AFE) phases, their coexistence and spatial distributions underpins the functionality of antiferrodistortive (AFD) multiferroics in the vicinity of morphotropic phase transitions. In this work, using Landau-Ginzburg-Devonshire (LGD) phenomenology and a semi-microscopic four sublattice model (FSM), we explore the behavior of different AFE, FEI, and FE long-range orderings and their coexistence at the morphotropic phase boundaries in FE-AFE-AFD multiferroics. These theoretical predictions are compared with the experimental observations for dense Bi 1-y R y FeO 3 ceramics, where R is Sm or La atoms with the fraction 0 ≤ y ≤ 0.25, as confirmed by the X-ray diffraction (XRD) and Piezoresponse Force Microscopy (PFM). These complementary measurements were used to study the macroscopic and nanoscopic transformation of the crystal structure with doping. The comparison of the measured and calculated AFE/FE phase fractions demonstrate that the LGD-FSM approach well describes the experimental results obtained by XRD and PFM for Bi 1-y R y FeO 3 . Hence, this combined theoretical and experimental approach provides further insight into the origin of the morphotropic boundaries and coexisting FE and AFE states in model rare-earth doped multiferroics.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Electromagnetic torque measurements in DIII-D using internal/external magnetic field decomposition

Given spatially resolved measurements of normal and tangential components of the magnetic field just outside the surface of a magnetically confined plasma, the field at the measurement location can be uniquely decomposed into contributions from the plasma and from external sources. This principle allows direct measurement of the electromagnetic torque on the plasma without knowledge of the distribution of the internal and external currents, similar to the more well-known formalism using the Maxwell stress tensor. The internal/external field decomposition also enables a mixed approach that incorporates any explicitly known current distributions (e.g., from non-axisymmetric coils). We discuss the requirements and limitations of such an approach to torque measurements. Experimental measurements of the torque evolution as a rotating tearing mode locks to the wall in the DIII-D tokamak are consistent with a simple model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Comprehensive study of radon progeny attachment to surfaces

Low energy, low rate experiments, such as searches for neutrinoless double β decay and dark matter, require unprecedentedly low levels of background in order to deliver their full science potential. 210 Po driven, neutron induced background, caused by nuclear (α, n) reactions on low-Z materials, direct background contributions of the 210 Po α radiation and desorption of the 210 Pb progeny 210 Bi from surfaces into the detector medium are of particular of concern. These backgrounds depend on details of the components’ exposure to radon-loaded laboratory air and, thus, their handling history. The attachment rates of airborne radon progeny to surfaces, needed for the estimation of these background rates, are poorly understood. This article reports the results of a campaign comprising of more than 1200 attachment measurements, performed for nine different materials. Correlations of the attachment with environmental parameters such as air exchange rate, electrical surface potential, temperature, atmospheric pressure, and relative humidity have been studied and found to be significant only in case of the first two. Attachment modeling, using the Jacobi model, is compared to data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

UV-Vis spectrophotometric determination of rare earth elements (REE) speciation at near-neutral to alkaline pH. Part I: m-cresol purple properties from 25-75 °C and Er hydrolysis

The speciation and mobility of rare earth elements (REE) strongly depends on pH which controls the formation of charged aqueous hydroxyl species. The latter potentially play an important role in controlling heavy REE adsorption on clay minerals in near-neutral to alkaline waters such as in regolith-hosted REE mineral deposits. However, accurate REE hydrolysis constants are needed for developing geochemical models that can predict the role of these charged species in natural systems. Here, we develop a robust experimental UV-Vis spectrophotometric method using m-cresol purple to determine in situ pH from 25 to 75 °C. This method is used to derive the average ligand number and hydrolysis constants of erbium (Er) at 25 °C in aqueous solutions with low ionic strength (≤ 0.001 mol/L) at pH from ~7 to 9.5 and in the presence of Er concentrations from 0 to 0.057 mM. The average ligand number ranges between 1 and 3 indicating that Er(OH) 2+ , Er(OH) 2 + and Er(OH) 3 0 control speciation in the experiments. The logarithm of the Er hydrolysis constants (log*β n ° , n= 1 to 3) derived at infinite dilution for the reaction Er 3+ + nH 2 O = Er(OH) n 3-n + nH + are: *β 1 ° = –7.22 ± 0.10, *β 2 ° = –14.52 ± 0.08, *β 3 ° = –23.24 ± 0.04. Implementation of these experimental data into a geochemical model indicates that the Er(OH) 2 + and Er(OH) 3 0 species are both stable in a much wider pH range than previously predicted. Consequently, the positively charged REE hydroxyl complexes can potentially control the fractionation of light vs. heavy REE via adsorption as observed in the formation of certain regolith-hosted REE deposits.

58 GEOSCIENCES↗
DOE OSTI2022

The CGM 2 Survey: Circumgalactic O vι from Dwarf to Massive Star-forming Galaxies

We combine 126 new galaxy-O vι absorber pairs from the CGM 2 survey with 123 pairs drawn from the literature to examine the simultaneous dependence of the column density of O vι absorbers (N O VI ) on galaxy stellar mass, star-formation rate, and impact parameter. The combined sample consists of 249 galaxy-O vι absorber pairs covering z = 0–0.6, with host galaxy stellar masses M * = 10 7.8 –10 11.2 M ⊙ and galaxy-absorber impact parameters R ⊥ = 0–400 proper kiloparsecs. In this work, we focus on the variation of N O VI with galaxy mass and impact parameter among the star-forming galaxies in the sample. We find that the average N O VI within one virial radius of a star-forming galaxy is greatest for star-forming galaxies with M * = 10 9.2 –10 10 M ⊙ . Star-forming galaxies with M * between 10 8 and 10 11.2 M ⊙ can explain most O vι systems with column densities greater than 10 13.5 cm –2 . Sixty percent of the O vι mass associated with a star-forming galaxy is found within one virial radius, and 35% is found between one and two virial radii. In general, we find that some departure from hydrostatic equilibrium in the CGM is necessary to reproduce the observed O vι amount, galaxy mass dependence, and extent. Our measurements serve as a test set for CGM models over a broad range of host galaxy masses.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

An Experimental Study of the Solubility of Rare Earth Chloride Salts (La, Nd, Er) in HCl Bearing Water Vapor from 350 – 425 °C

In this work, the solubilities of the rare earth chlorides REECl 3 , where REE = (La, Nd, Er), were measured in HCl bearing water vapor from 350 – 425°C with water partial pressures ranging from 8 – 170 bar. Solubility data were fit to the Pitzer-Pabalan quasi-chemical model in order to extract thermodynamic parameters for the formation of the gaseous REE-chloride-water clusters REECl 3 (H 2 O) n . The data show that the solubility of the REE chlorides are orders of magnitude higher than salts such as NaCl or CuCl at low water fugacities, despite their sublimation energies being substantially higher. This enhanced solubility is likely due to the high enthalpy associated with binding a single water molecule to form the species REECl 3 (H 2 O), with derived enthalpies ranging from –378 to –465 kJ/mol. Addition of further water molecules to form higher order clusters (n > 1) involves enthalpy changes of ~ -20 kJ/mol, and are in effect thermodynamically suppressed over the temperature range 350 – 425°C. Despite the enhanced solubility of small REECl 3 water clusters, simulations of boiling processes demonstrate that the REE show highly conservative behavior, partitioning strongly into the dense aqueous phase. Not surprisingly, the presence of phosphates in the system makes this effect even more pronounced, completely immobilizing the REE. This would reduce transport in both the vapor and aqueous phase to negligible levels. However, we suggest that vapor phase transport of the REE may play a significant role in systems having a relatively low partial pressure of water (below the saturation point), where the relatively high stability of the first hydrated REE chloride clusters (REECl 3 (H 2 O) and REECl 3 (H 2 O) 2 ) will give a preference for gas transport of the REE relative to other elements. This can likely happen in systems involving a gas/melt exchange in fumarolic exhalations, where water vapor discharges at relatively low (close to atmospheric) pressures.

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

2025 ASMS Investigation of the Collision-Induced Dissociation Mechanism of Protonated TODGA with IRIS

Title (20 words): Investigation of the Collision-Induced Dissociation Mechanism of Protonated TODGA with IRIS Introduction (120 words): One of the challenges facing wide-spread adoption of nuclear power is the development of efficient separation processes for used nuclear fuel. The molecules in separation processes are subjected to an extreme environment due to the high radiation fields from the used fuel and highly acidic media used for fuel dissolution, which results in significant molecular degradation, leading to reduced process efficiency. These degradation products must be identified and studied so mitigation strategies can be developed to maintain process efficiency. However, complex systems can have many degradation products, complicating identification. Untargeted analysis tools could be used to understand radiation chemistry in complex systems. However, this would necessitate improved understanding of the gas-phase fragmentation mechanisms of fuel cycle molecules like tetraoctyldiglycolamide (TODGA). Methods (120 words): The gas-phase fragmentation of protonated TODGA was investigated using collision-induced dissociation (CID), resonance ejection, and infrared ion spectroscopy (IRIS). CID and resonance ejection experiments were conducted using a Bruker Daltonics (Bremen, Gemany) SolariX XR fourier transform ion cyclotron resonance (FT-ICR) mass spectrometer. IRIS spectra of protonated TODGA and its two CID fragmentation products were measured using a modified Bruker amaZon Speed ETD 3D quadrupole ion trap mass spectrometer coupled to the Free Electron Lasers for Infrared eXperiments (FELIX) free electron laser. Measured spectra were compared with density functional theory (DFT) calculations using the Gaussian 16, Revision C.02 software package with the ?B97X-D functional and def2-TZVPP basis sets. Candidate structures were generated using the CREST 3.0 conformational sampling software tool. Preliminary Data (300 words): Collision-induced dissociation of protonated TODGA ([C36H73N2O3]+, m/z=581.562) results two fragment ions, one at m/z=340.285 assigned as [C20H38NO3]+ and the other at m/z=312.290, assigned as [C19H38NO2]+. Based on the assigned formula and the structure of protonated TODGA, the fragment at m/z=340.285 is likely formed from elimination of neutral dioctylamine. Comparison of the IRIS spectrum of m/z=340.285 with DFT predictions suggests it contains a ring structure, and is assigned as N-octyl-N-(6-oxo-1,4-dioxan-2-ylidene)octan-1-aminium. Based on this structure and the structure of protonated TODGA, we hypothesize this fragment formed from elimination of neutral dioctylamine followed by a ring closure mechanism. Comparison of the IRIS spectrum of the fragment at m/z=312.290 with DFT predictions also indicated the presence of a ring structure, assigned as N-(1,3-dioxolan-4-ylidene)-N-octyloctan-1-aminium. This product could be formed from elimination of carbon monoxide from the ring of m/z=340.285 as a sequential fragmentation or formed directly from protonated TODGA via elimination of neutral N,N-dioctylformamide followed by a ring closure. Resonance ejection experiments where m/z=340 was continuously ejected from the IRC cell showed no decrease in intensity of m/z=312.290 across several collision energies, suggesting that the later, direct formation mechanism, dominates. The location of the ionizing proton in protonated TODGA is important for modeling the fragmentation mechanisms. DFT calculations suggested that the position of bands involving the coupled vibrations of the amide C—N and C=O bonds in TODGA are the most sensitive to proton location. Evaluation of the IRIS spectrum of protonated TODGA suggests that the ionizing proton is located between the two amid oxygens. This protonation location was calculated to lie approximately 30 kJ/mol lower in energy than the next lowest energy location, with the proton located solely on one of the amide oxygens. Novel aspect (20 words): Infrared ion spectroscopy combined with resonance ejection experiments and density functional theory to probe the collision-induced dissociation mechanism of tetraoctyldiglycolamide.

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

Machine learning based approach to predict ductile damage model parameters for polycrystalline metals

Damage models for ductile materials typically need to be parameterized, often with the appropriate parameters changing for a given material depending on the loading conditions. This can make parameterizing these models computationally expensive, since an inverse problem must be solved for each loading condition. Using standard inverse modeling techniques typically requires hundreds or thousands of high-fidelity computer simulations to estimate the optimal parameters. Additionally, the time of a human expert is required to set up the inverse model. Machine learning has recently emerged as an alternative approach to inverse modeling in these settings, where the machine learning model is trained in an offline manner and new parameters can be quickly generated on the fly, after training is complete. Here, this work utilizes such a workflow to enable the rapid parameterization of a ductile damage model called TEPLA with a machine learning inverse model. The machine learning model can efficiently estimate the model parameters much faster, as compared to previously employed methods, such as Bayesian calibration. The results demonstrate good accuracy on a synthetic test dataset and is validated against experimental data.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Efficient First-Order Algorithms for Large-Scale, Non-Smooth Maximum Entropy Models with Application to Wildfire Science

Maximum entropy (MaxEnt) models are a class of statistical models that use the maximum entropy principle to estimate probability distributions from data. Due to the size of modern data sets, MaxEnt models need efficient optimization algorithms to scale well for big data applications. State-of-the-art algorithms for MaxEnt models, however, were not originally designed to handle big data sets; these algorithms either rely on technical devices that may yield unreliable numerical results, scale poorly, or require smoothness assumptions that many practical MaxEnt models lack. In this paper, we present novel optimization algorithms that overcome the shortcomings of state-of-the-art algorithms for training large-scale, non-smooth MaxEnt models. Our proposed first-order algorithms leverage the Kullback–Leibler divergence to train large-scale and non-smooth MaxEnt models efficiently. For MaxEnt models with discrete probability distribution of n elements built from samples, each containing m features, the stepsize parameter estimation and iterations in our algorithms scale on the order of O(mn) operations and can be trivially parallelized. Moreover, the strong ℓ1 convexity of the Kullback–Leibler divergence allows for larger stepsize parameters, thereby speeding up the convergence rate of our algorithms. To illustrate the efficiency of our novel algorithms, we consider the problem of estimating probabilities of fire occurrences as a function of ecological features in the Western US MTBS-Interagency wildfire data set. Our numerical results show that our algorithms outperform the state of the art by one order of magnitude and yield results that agree with physical models of wildfire occurrence and previous statistical analyses of wildfire drivers.

Physics↗
DOE OSTI2024

Inverse-Velocity Transformation Wall Model for Reacting Turbulent Hypersonic Boundary Layers

The present study builds on prior work by taking advantage of the novel framework proposed by Griffin et al. (hereafter referred to as the GFM) as a baseline. The model is progressively extended to multicomponent reacting mixtures, accounting for differential diffusion and finite-rate chemistry in a similar fashion to Di Renzo & Urzay and Di Renzo et al. The accuracy of the present approach, as well the prior model of Di Renzo & Urzay, is assessed in an a priori sense for the first time in a turbulent reacting boundary layer, using the boundary-layer data of Williams et al. Five species are included in the present analysis, i.e., N s = 5, namely N 2 , O 2 , NO, N and O, a neutral mixture most representative of dissociation/recombination phenomena for temperatures below 6000 K. A schematic of a flow over a wedge representative of the described configuration is presented in Figure 1. The brief is organized as follows: In Section 2, the wall-model equations and the computational framework are presented. In Section 3, the a priori results of the proposed model are described and compared to the extended EWM. Finally, in Section 4, some conclusions are offered.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Dynamic, symmetry-preserving, and hardware-adaptable circuits for quantum computing many-body states and correlators of the Anderson impurity model

We present a hardware-reconfigurable ansatz on N q -qubits for the variational preparation of many-body states of the Anderson impurity model (AIM) with N imp + N bath = N q /2 sites, which conserves total charge and spin z component within each variational search subspace. The many-body ground state of the AIM is determined as the minimum over all minima of O(N$^2_ q$) distinct charge-spin sectors. Hamiltonian expectation values are shown to require ω(N q ) < N meas. $\leqslant$ O(N imp N bath ) symmetry-preserving, parallelizable measurement circuits, each amenable to postselection. To obtain the one-particle impurity Green’s function we show how initial Krylov vectors can be computed via midcircuit measurement and how Lanczos iterations can be computed using the symmetry-preserving ansatz. For a single-impurity Anderson model with a number of bath sites increasing from one to seven, we show using numerical emulation that the ease of variational ground-state preparation is suggestive of linear scaling in circuit depth and subquartic scaling in optimizer complexity. We therefore expect that, combined with time-dependent methods for Green’s function computation, our ansatz provides a useful tool to account for electronic correlations on early fault-tolerant processors. Finally, with a view towards computing real materials properties of interest like magnetic susceptibilities and electron-hole propagators, we provide a straightforward method to compute many-body, time-dependent correlation functions using a combination of time evolution, midcircuit measurement-conditioned operations, and the Hadamard test.

36 MATERIALS SCIENCE↗
DOE OSTI2023

The kinetics of the PuO 2 to Pu 2 O 3 conversion

Here in an oxidizing environment, the oxide formed on plutonium (Pu) metal is composed of a plutonium dioxide (PuO 2 ) top layer and a thin cubic plutonium sesquioxide (Pu 2 O 3 ) middle layer. In a reducing environment, the PuO 2 layer auto-reduces to cubic Pu 2 O 3 . The speed and extent of this conversion depend on the combination of temperature and time. While PuO 2 provides a strong diffusion barrier against unwanted Pu corrosion by gaseous species (like hydrogen), Pu 2 O 3 does not, since its crystal structure has chains of oxygen vacancies. The kinetics of the PuO 2 reduction are, therefore, of fundamental interest and enable researchers to better protect Pu from corrosion. In this report, the oxygen-diffusion-limited kinetics of the dioxide to sesquioxide conversion were obtained by dynamically heating a PuO 2 -covered Pu sample from 294 to 418 K in a high-vacuum vessel equipped with an in situ spectroscopic ellipsometer. The physical/chemical constraints in the conversion process were combined with the ellipsometry method of multi-sample analysis to track the percentage of PuO 2 and to compute the extent of Pu 2 O 3 formation. The resulting diffusion coefficients were compared against and then combined with complementary literature data to produce a comprehensive set of kinetic parameters for reliably modeling oxide conversion over a larger temperature range than spanned by prior studies. The extracted thermal activation energy barrier (43.7 kJ/mol) and pre-exponential factor (5.0 × 10 -10 cm 2 /s) for the oxygen-diffusion-limited process can be used to accurately model the PuO 2 to Pu 2 O 3 transformation in vacuum and/or inert gas applications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

New thermal decomposition pathway for TATB

Abstract Understanding the thermal decomposition behavior of TATB (1,3,5-triamino-2,4,6-trinitrobenzene) is a major focus in energetic materials research because of safety issues. Previous research and modelling efforts have suggested benzo-monofurazan condensation producing H 2 O is the initiating decomposition step. However, early evolving CO 2 (m/z 44) along with H 2 O (m/z 18) evolution have been observed by mass spectrometric monitoring of head-space gases in both constant heating rate and isothermal decomposition studies. The source of the CO 2 has not been explained, until now. With the recent successful synthesis of 13 C 6 -TATB ( 13 C incorporated into the benzene ring), the same experiments have been used to show the source of the CO 2 is the early breakdown of the TATB ring, not adventitious C from impurities and/or adsorbed CO 2 . A shift in mass m/z 44 (CO 2 ) to m/z 45 is observed throughout the decomposition process indicating the isotopically labeled 13 C ring breakdown occurs at the onset of thermal decomposition along with furazan formation. Partially labeled (N 18 O 2 ) 3 -TATB confirms at least some of the oxygen comes from the nitro-groups. This finding has a significant bearing on decomposition computational models for prediction of energy release and deflagration to detonation transitions, with respect to conditions which currently do not recognize this oxidation step.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Understanding the tool influence function during sub-aperture belt-on-wheel glass polishing

The tool influence function (TIF) during sub-aperture belt-on-wheel polishing has been evaluated as a function of various process conditions (belt use/wear, dwell time, displacement, belt velocity, and wheel modulus and diameter) on fused silica glass workpieces using C e O 2 polishing media. TIF spots are circular or elliptical in shape with a largely flat bottom character. Furthermore, the volumetric removal rate varies significantly with belt use (or wear), stabilizing after ~<!-- ~ --> 15 m i n of use. A modified Preston model, where the pressure dependence is adjusted using a different scaling of the wheel modulus ( E w 0.5 ), largely predicts the volumetric removal rate over the range of process conditions evaluated. The relatively high volumetric removal rate of 30 -<!-- - --> 60 m m 3 / h using a fixed C e O 2 -in-resin-host belt offers a rapid, and hence, more economical, initial polish of aspheric and freeform optics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Solidification cracking of refractory alloys: a computational and machine learning study to investigate composition-dependence for improved weldability and additive manufacturability

Large-batch numerical, CALculation of PHAse Diagrams (CALPHAD)-based solidification cracking calculations are performed and then analyzed with machine learning methods to generate models that relate chemistry of refractory alloys to cracking susceptibility. Kou’s solidification cracking index is used to study the refractory alloys including O, N, C binary mixtures with Mo, Ta, Nb, and W, the molybdenum-based TZM, Niobium-based C103, and Tantalum-based T111 and Ta-10 W, as well as hypothetical refractory ternary alloys. Findings strongly validate Kou’s Crack Susceptibility Index (CSI) against Varestraint test data for Nb- and Ta-based alloys, establishing CSI thresholds where refractory alloys with CSI < 15,000 K are likely weldable, CSI > 15,000 K are prone to cracking, and CSI > 25,000 K are likely unweldable (or unprintable). Furthermore, interstitial elements C, N, and O significantly increase crack susceptibility, with some existing material specifications coinciding with peak cracking susceptibility concentrations. Finally, machine learning-derived elemental potency factors enable rapid prediction of CSI from alloy chemistry for C103, TZM, Ta-10 W, and T-111 alloys. These results provide practical guidance for feedstock selection, powder reuse limits, and alloy specification amendments for welding and additive manufacturing applications.

36 MATERIALS SCIENCE↗