Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computational methods”

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

Predicting metal-binding proteins and structures through integration of evolutionary-scale and physics-based modeling

Metals are essential elements in all living organisms, binding to approximately 50% of proteins. They serve to stabilize proteins, catalyze reactions, regulate activities, and fulfill various physiological and pathological functions. While there have been many advancements in determining the structures of protein-metal complexes, numerous metal-binding proteins still need to be identified through computational methods and validated through experiments. Here, to address this need, we have developed the ESMBind workflow, which combines evolutionary scale modeling (ESM) for metal-binding prediction and physics-based protein-metal modeling. Our approach utilizes the ESM-2 and ESM-IF models to predict metal-binding probability at the residue level. In addition, we have designed a metal-placement method and energy minimization technique to generate detailed 3D structures of protein-metal complexes. Our workflow outperforms other models in terms of residue and 3D-level predictions. To demonstrate its effectiveness, we applied the workflow to 142 uncharacterized fungal pathogen proteins and predicted metal-binding proteins involved in fungal infection and virulence.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2019

Type IV failure in weldment of creep resistant ferritic alloys: I. Micromechanical origin of creep strain localization in the heat affected zone

Creep strength enhanced ferritic (CSEF) steels containing 9-12wt% chromium have been extensively used in fossil-fuel-fired power plants. Despite their excellent creep resistance at high temperatures, premature failures (especially Type IV cracking) are often found in the fine-grained heat affected zone (HAZ) or intercritical HAZ of the welded components. This failure mode is preceded by the strain localization in the HAZ, as measured by the Digital Image Correlation (DIC) technique. The present work aims to develop a finite-element based computational method to determine the micromechanical and microstructural origin of the strain localization phenomenon. We construct a two-dimensional digital microstructure based on the actual microstructure of ferritic steel weldments by using the Voronoi-tessellation method, to account for the effects of its large grain-size gradients. A mechanism-based finite element method is developed for modeling the high temperature deformation resulting from a synergy of thermally activated dislocation movements, diffusional flow and grain boundary sliding. The numerical results agree well with the strain measurements by our DIC technique, particularly revealing the effect of pre-welding tempering on the evolution of strain localization in HAZ of creep resistant steel weldments. It is found that the diffusional creep with dependence on grain sizes, dislocation creep with dependence on material strength, and more importantly, grain boundary sliding, contribute synergistically to the creep strain accumulation in the HAZ, and their relative degree of significance is quantified. The creep rupture life will be investigated in the companion paper.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Revealing the complex chemistry of grain boundaries in K-doped BaFe 2 As 2 with atom probe tomography

Iron-based superconductors have attractive properties for high-field applications, but there is a lack of understanding of the effect of grain boundary chemistry on the in-field performance. The near atomic-scale resolution, ppm sensitivity and 3D analysis offered by atom probe tomography make it a powerful tool to investigate the nanoscale structure and chemistry of these defects in fine-grained K-doped BaFe 2 As 2 samples. A computational method to systematically extract and compare the Gibbsian interfacial excess of chemical species across grain boundaries has been explored in this work. The robustness of the method has been tested by evaluating the effects of selected variables on simulated APT datasets. The accuracy and precision of the calculated Gibbsian interfacial excess were found to be stable over a range of analysis conditions: varying grain boundary widths and detection efficiencies, spatial precisions below 1.5 nm, and bin widths between 1.2 and 1.6 nm. For the K-doped BaFe 2 As 2 samples studied, segregation of As, Ba, K and impurities of O, Na, and Sb were found at grain boundaries. The Gibbsian excess values were found to vary widely between different boundaries, showing the complexity of the grain boundary chemistry in this material. Possible links between the observed critical current density (Jc) of these samples and their nano- and micro-structure have also been investigated and discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Integrated simulation of U-10Mo monolithic fuel swelling behavior

Here, a separate computational branch has been implemented within the DART (Dispersion Analysis Research Tool) computational code to simulate the swelling behavior of U-10Mo monolithic fuel under the operating conditions of high-power research and test reactors (RTRs). The monolithic branch of the DART code implements a mechanistic rate-theory-based fission-gas-behavior model for the calculation of fission gas swelling, as well as a suite of thermal, physical, and mechanical models to take into account various processes occurring in RTR fuels during irradiation. In order to accurately simulate and eventually predict U-10Mo monolithic fuel irradiation behavior, the code uses materials properties calculated with lower length-scale computational methods, such as gas atom diffusivity and U-Mo surface energy from atomic simulations and grain-morphology-specific recrystallization kinetics (recrystallized fuel volume fractions vs. fission density) predicted using the phase-field method. The remainder of fission gas behavior parameters used in the model were calibrated with measured intergranular bubble size distributions. With this integrated simulation approach, the swelling behavior of U-10Mo monolithic fuel was simulated for various initial grain sizes at different operating conditions and compared with measured data. Furthermore, because limited experimental data exist for parameter calibration detailed sensitivity studies for the important parameters used in the fission gas behavior model were performed in order to examine their impact on both intergranular gas bubble morphology at low fission density, and on total porosity at high fission density.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Synthesis and characterization of the dinuclear cobalt(III) complex: [(C 2 F 5 ) 3 Co(μ-F)] 2 2–

Here, reaction of the versatile tris(perfluoroethyl) cobalt(III) precursor [fac-(MeCN) 3 Co(C 2 F 5 ) 3 ] with [NMe 4 ]F and [PPh 4 ]Cl in THF formed the unexpected cobalt(III) bridging fluoride dimer ([(C 2 F 5 ) 3 Co(μ-F)] 2 2– . The new fluoro-organometallic cobalt(III) complex was characterized by NMR and UV-vis spectroscopies, X-ray crystallography, cyclic voltammetry, and by computational methods. In the strongly coordinating solvent MeCN, ([(C 2 F 5 ) 3 Co(μ-F)] 2 2– exhibits dynamic processes on the NMR timescale which are consistent with solvent coordination. However, in the weakly coordinating solvent CH 2 Cl 2 , a more static structure is observed suggesting that the dimer retains its structure in solution state. An electrochemical analysis of ([(C 2 F 5 ) 3 Co(μ-F)] 2 2– was performed, and the data are compared to previously reported cobalt(III) perfluoroethyl complexes.

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

Effects of charged interfaces on electrolyte decomposition at the lithium metal anode

Lithium–Sulfur batteries are promising candidates to substitute conventional Li-ion batteries due to their higher energy density and reduced cost. However, several challenges related to the reactivity of the lithium metal anode have prevented this technology from becoming broadly commercialized. Lithium's high reactive nature leads to the continuous decomposition of the electrolyte and the formation of the solid-electrolyte interphase (SEI) layer. Thus, a comprehensive understanding of how the SEI film is formed is crucial to help elucidate improvements for this battery technology. In this work, we use density functional theory (DFT) based computational methods to investigate the effect of charged interfaces on the electrolyte degradation at the lithium anode. Here, several electrolyte mixtures were considered including 1,2-dimethoxyethane (DME) and 1,3-dioxolane (DOL) as solvents, and lithium bis(fluorosulfonyl)imide (LiFSI) and lithium bis(trifluoromethanesulfonyl)imide (LiTFSI) as salts. It is found that the extent of salt decomposition is higher when the interface is charged. In addition, two types of solvent reduction mechanisms are identified: C–O bond cleavage and radical attacks. Charge transfer and evolution analysis are characterized in detail. Finally, solvent decomposition reactions are found to become more thermodynamically favorable and their activation barriers to be diminished under the effect of constant potentials and electric fields.

25 ENERGY STORAGE↗
DOE OSTI2013

Hydrothermal reactions: From the synthesis of ligand to new lanthanide 3D-coordination polymers

The organic ligand 2,5-piperazinedione-1,4-diacetic acid (H{sub 2}PDA) was synthesized under hydrothermal conditions starting from the iminodiacetic acid and catalyzed by oxalic acid. The X-ray powder diffraction data indicates that the compound crystallizes in the P2{sub 1}/c monoclinic system as reported in the literature. The ligand was also characterized by elemental analysis, magnetic nuclear resonance, infrared spectroscopy and thermogravimetric analysis. Two new coordination networks based on lanthanide ions were obtained with this ligand using hydrothermal reaction. In addition to single-crystal X-ray diffraction, the compounds were characterized by infrared spectroscopy, thermogravimetric analysis, scanning electron microscopy and elemental analysis. Single-crystal XRD showed that the compounds are isostructural, crystallizing in P2{sub 1}/n monoclinic system with chemical formula [Ln(PDA){sub 1.5}(H{sub 2}O)](H{sub 2}O){sub 3} (Ln=Gd{sup 3+}(1) and Eu{sup 3+}(2)).The luminescence properties of both compounds were studied. In the compound (1), a broad emission band was observed at 479 nm, redshifted by 70 nm in comparison of the free ligand. In (2), the typical f–f transition was observed with a maximum peak at 618 nm, related with the red emission of the europium ions. Computational methods were performed to simulate the crystal structure of (2). The theoretical calculations of the intensity parameters are in good agreement with the experimental values. - Graphical abstract: Scheme of obtaining the ligand 2,5-piperazinedione-1,4-diacetic acid (H{sub 2}PDA) and two new isostructural 3D-coordination polymers [Ln(PDA){sub 1.5}(H{sub 2}O)](H{sub 2}O){sub 3} (Ln=Gd{sup 3+} and Eu{sup 3+}) by hydrothermal synthesis. Display Omitted - Highlights: • The ligand 2,5-piperazinedione-1,4-diacetic acid was synthetized using the hydrothermic method and characterized. • Two new 3D-coordination polymers with this ligand containing Gd{sup 3+} and Eu{sup 3+} ions were also obtained. • The compounds are isostructural and the typical luminescent properties were observed in both structures.

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

Structures and thermoelectric properties of double-filled (Ca{sub x}Ce{sub 1−x})Fe{sub 4}Sb{sub 12} skutterudites

The structures and thermoelectric properties of the double-filled (Ca{sub x}Ce{sub 1−x})Fe{sub 4}Sb{sub 12} series (x=0, 0.25, 0.5, 0.75, and 1) have been studied using a combined experimental and computational methods. Compounds of (Ca{sub x}Ce{sub 1−x})Fe{sub 4}Sb{sub 12} were obtained only for x=0, 0.5, and 1. Composition with x=0.25 was found to be a mixture of x=0 and 0.5 compounds, and composition with x=0.75 was found to be a mixture of x=1 and 0.5 compounds, respectively. Our conclusions on phase formation are supported by density functional theory (DFT) calculations. In Ca{sub 0.5}Ce{sub 0.5}Fe{sub 4}Sb{sub 12,} Ca substitution in the Ce site of CeFe{sub 4}Sb{sub 12} leads to high hole concentrations, resulting in stronger semimetal transport as compared to CeFe{sub 4}Sb{sub 12}. Ca{sub 0.5}Ce{sub 0.5}Fe{sub 4}Sb{sub 12} yields a slightly higher ZT value than that of CeFe{sub 4}Sb{sub 12}, which is attributed to its lower lattice thermal conductivity. Phonon mode calculations adopting a three-particle bending model suggest that thermal conductivity is reduced upon Ca substitution because of an additional vibration mode which involves both Ca and Ce atoms. - Graphical abstract: Detailed structural information to correlate with thermoelectric properties in a series of double-filled (Ca{sub x}Ce{sub 1−x})Fe{sub 4}Sb{sub 12} skutterudite samples were obtained using synchrotron X-ray diffraction and first principle calculations. - Highlights: • Research focus on phase stability of CaCe double-filled p-type skutterudite. • Stable structure forms for an equiatomic occupancy of the voids by Ca and Ce. • Samples feature two-phase mixtures with variations of filling atoms ratio. • Phase stablility is investigated by DFT total energy calculations. • Sample with Equiatomic occupancy of Ca and Ce show higher ZT.

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

Aqueous-based route toward Fe:ZnSe semiconductor nanocrystals: Synthesis and characterization

Nearly monodisperse and luminescent Fe-doped ZnSe semiconductor nanocrystals were synthesized by an aqueous-based route. The obtained nanocrystals were characterized by X-ray diffraction, transmission electron microscopy, Fourier transform infrared spectroscopy, energy dispersive X-ray spectroscopy, X-ray photoelectron spectroscopy, photoluminescence and ultraviolet-visible absorption spectroscopy. The resulting nanocrystals were well retained in the zinc blende structure, and the Fe dopants were doped into the ZnSe nanocrystals, as confirmed by X-ray photoelectron spectroscopy. From combination of transmission electron microscopy result and computing methods according to ultraviolet-visible absorption spectrum and X-ray diffraction pattern, particle size of the obtained nanocrystals was found to be about 5 nm. Thioglycolic acid as ligand was successfully capped on the surface of the resulting nanocrystals, confirmed by Fourier transform infrared and energy dispersive X-ray spectroscopy. The obtained nanocrystals exhibit an emission peak at 431 nm; and the photoluminescence band is relatively narrow and symmetric, which indicates that the nanocrystals are nearly monodisperse and homogeneous. Absorption shoulder of the resulting nanocrystals is blue-shifted compared with that of corresponding bulk ZnSe, indicating the quantum confinement effect. The possible reaction mechanism of formation of Fe-doped ZnSe nanocrystals in the aqueous solution was also discussed. Our results demonstrate that the aqueous-based route could be successfully adopted for producing high-quality doped semiconductor nanocrystals. - Research Highlights: {yields} An aqueous-based route is demonstrated for synthesizing Fe:ZnSe nanocrystals. {yields} Structure, composition and optical properties of the nanocrystals were investigated. {yields} The synthetic route could be adopted for producing high-quality doped nanocrystals.

36 MATERIALS SCIENCE↗
DOE OSTI2012

CoTiO{sub 3} via cobalt oxalate-TiO{sub 2} precursor. Synthesis and characterization

Purely trigonal ilmenite type powder; CoTiO{sub 3} was successfully synthesized through the thermal decomposition of CoC{sub 2}O{sub 4}{center_dot}2H{sub 2}O-TiO{sub 2} precursor. The decomposition course was followed using differential thermal analysis-thermogravimetry measurements and the intermediate decomposition products as well as the titanate product were characterized using X-ray diffraction, Fourier transform infrared, transmission electron microscopy techniques. Thermal analysis study indicated that, the intermediate decomposition product; Co{sub 3}O{sub 4} decomposed endothermically into CoO at 963 Degree-Sign C with instantaneous formation of CoTiO{sub 3}. X-ray diffraction of the powder calcined at 950 Degree-Sign C showed characteristic lines for Co{sub 2}TiO{sub 4} besides secondary phase of rutile-TiO{sub 2}, which are disappeared by rising the calcination temperature to 1000 Degree-Sign C. Transmission electron microscopy image of CoTiO{sub 3} powder exhibit weak agglomerated particles with heterogeneous morphology in both shape and dimensions. Type II adsorption isotherm estimated using Brunauer-Emmett-Teller procedure indicated a very high specific surface area of 29.7 m{sup 2} g{sup -1}. Electrical properties of CoTiO{sub 3} measured as a function of temperature and frequency showed a semiconducting behavior at high temperatures with conduction activation energy of 0.73 eV. Kinetic analysis of the non-isothermal thermogravimetry curves, using different computational methods, illustrated that the random nucleation and growth mechanisms are the predominant one. - Graphical Abstract: The figure showed the variation of dielectric constant ({epsilon} Prime ) for CoTiO{sub 3} with temperature at different frequencies. {epsilon} Prime appeared to be temperature and frequency independent at low temperatures. By increasing the temperature, the orientation polarization plays a significant role and a gradual increase of {epsilon} Prime occurred. Highlights: Black-Right-Pointing-Pointer CoTiO{sub 3} was synthesized through decomposition of CoC{sub 2}O{sub 4}{center_dot}2H{sub 2}O-TiO{sub 2} precursor. Black-Right-Pointing-Pointer XRD, FT-IR, TEM techniques are used to characterize decomposition course. Black-Right-Pointing-Pointer Specific surface area of 29.7 m{sup 2} g{sup -1} was obtained. Black-Right-Pointing-Pointer The prepared titanate showed semiconducting behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Prediction of perovskite and other ternary oxide multilayers as mirrors for soft X-rays

Highlights: • Theoretical framework to identify optimal multilayers for high x-ray reflectivity. • Perovskite oxide multilayers are proposed as soft x-ray mirrors. • High soft x-ray reflectivity up to 33% are expected in perovskite multilayers. - Abstract: We propose multilayers of perovskite and related ternary oxides with a general formula of ABO{sub 3} as high reflectivity mirrors, especially for use in water window region (2.3–4.4 nm). The high reflectivity combinations of oxide multilayers were deduced using evolutionary search algorithms such as genetic algorithms, and the calculation speed was accelerated via parallel computing methods We also propose a figure of merit for X-ray reflectivity in periodic multilayers systems, which can simplify the future efforts on identifying material combinations, and the search through this multi-dimensional parameter space. The highest reflectivity value was found to be over 33% at 3.1 nm in the water window region. The effect of interface roughness was simulated and the decrease in reflectivity was found to be modest for practically achievable roughness values. This work establishes the foundation for future experimental and theoretical studies towards achieving high reflectivity x-ray mirrors of complex oxide multilayers.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Deep-Learning-Derived Evaluation Metrics Enable Effective Benchmarking of Computational Tools for Phosphopeptide Identification

Tandem mass spectrometry (MS/MS)-based phosphoproteomics is a powerful technology for global phosphorylation analysis. However, applying four computational pipelines to a typical mass spectrometry (MS)-based phosphoproteomic dataset from a human cancer study, we observed a large discrepancy among the reported phosphopeptide identification and phosphosite localization results, underscoring a critical need for benchmarking. While efforts have been made to compare performance of computational pipelines using data from synthetic phosphopeptides, evaluations involving real application data have been largely limited to comparing the numbers of phosphopeptide identifications due to the lack of appropriate evaluation metrics. We investigated three deep learning-derived features as potential evaluation metrics: phosphosite probability, Delta RT and spectral similarity. Predicted phosphosite probability is computed by MusiteDeep, which provides high accuracy as previously reported; Delta RT is defined as the absolute retention time (RT) difference between RTs observed and predicted by AutoRT; and spectral similarity is defined as the Pearson’s correlation coefficient between spectra observed and predicted by pDeep2. Using a synthetic peptide dataset, we found that both Delta RT and spectral similarity provided excellent discrimination between correct and incorrect peptide-spectrum matches (PSMs) both when incorrect PSMs involved wrong peptide sequences and even when incorrect PSMs were caused by only incorrect phosphosite localization. Based on these results, we used all the three deep learning-derived features as evaluation metrics to compare different computational pipelines on diverse set of phosphoproteomic datasets and showed their utility in benchmarking performance of the pipelines. The benchmark metrics demonstrated in this study will enable users to select computational pipelines and parameters for routine analysis of phosphoproteomics data and will offer guidance for developers to improve computational methods.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2019

An investigation of geometrically necessary dislocations and back stress in large grained tantalum via EBSD and CPFEM

This study explores the evolution of GNDs and their effects on back stress through experimental and computational methods. Four large-grained tantalum tensile specimens were strained in uniaxial tension, electron backscatter diffraction (EBSD) data were collected, and geometrically necessary dislocation (GND) maps of the four specimens in the unloaded state were produced. EBSD-based GND maps revealed several types of features with high GND content which caused back stress in the specimens. Correlations between five geometrically-based grain boundary (GB) transmission factors and the GB GND content were evaluated, and statistically significant correlations were found for transmission factors based on Livingston and Chalmer’s N factor, Werner and Prantl’s slip transfer number, and GB misorientation. The sign of individual components of the Nye tensor were used to visually and quantitatively identify clustering of GNDs of the same sign, thus giving additional evidence of increasing back stress due to deformation. Deformation of one of the specimens was simulated using multiple CPFEM based modeling approaches and predicted stress-strain responses are compared. The super dislocation model (SD model) — a crystal plasticity finite element method (CPFEM) which incorporates elastic dislocation interactions — was able to isolate impact of back stress on the overall flow stress. The SD model predicted correct stresses when compared with experimental data; however, when the elastic interactions in the SD model were turned off, stress predictions were 25% too low. Thus, demonstrating the importance of incorporating back stress into the model.

42 ENGINEERING↗
DOE OSTI2016

Neutron-fragment and Neutron-neutron Correlations in Low-energy Fission

A computational method has been developed to simulate neutron emission from thermal-neutron induced fission of {sup 235}U and from spontaneous fission of {sup 252}Cf. Measured pre-emission mass-yield curves, average total kinetic energies and their variances, both as functions of mass split, are used to obtain a representation of the distribution of fragment velocities. Measured average neutron multiplicities as a function of mass split and their dependence on total kinetic energy are used. Simulations can be made to reproduce measured factorial moments of neutron-multiplicity distributions with only minor empirical adjustments to some experimental inputs. The neutron-emission spectra in the rest-frame of the fragments are highly constrained by ENDF/B-VII.1 prompt-fission neutron-spectra evaluations. The n-f correlation measurements of Vorobyev et al. (2010) are consistent with predictions where all neutrons are assumed to be evaporated isotropically from the rest frame of fully accelerated fragments. Measured n-f and n-n correlations of others are a little weaker than the predictions presented here. These weaker correlations could be used to infer a weak scission-neutron source. However, the effect of neutron scattering on the experimental results must be studied in detail before moving away from a null hypothesis that all neutrons are evaporated from the fragments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

The structural and dynamical role of water in natural organic matter: A 2H NMR and XRD study

Natural organic matter (NOM) is an important component in many near-surface geochemical environments, and its properties are greatly affected by the incorporation of water. Because of its importance, the macroscopic behavior and effects of water in NOM and soil organic matter (SOM) have been extensively studied using a wide range of experimental and computational methods. The molecular scale structural and dynamical behavior of water in these materials, however, is less well understood. This paper presents a variable temperature 2 H NMR and XRD study of water in Suwannee River NOM and its fulvic acid (FA) and humic acid (HA) fractions that provides new insight into the dynamical behavior of structurally different types of water and exchangeable hydrogen environments in NOM. The results provide a basis for future studies of more complex natural organic materials and the interaction of organic materials with mineral surfaces. Room temperature 2 H NMR spectra of samples hydrated in 2 H 2 O and then dehydrated, distinguish 2 H 2 O molecules that are in rapid reorientational motion (correlation times, ν c, > 10 5 Hz), 2 H exchanged onto carboxylic sites of the NOM that do not undergo rapid reorientation at frequencies >~10 3 Hz, and 2 H exchanged onto phenolic and possibly other alcohol sites of the NOM that undergoes rapid, but anisotropic, dynamical reorientation. For samples exposed to water and not dried, the XRD results collected at temperatures from 173 to 298 K show the formation of ice-1h in samples exposed to 100% relative humidity (R.H.) but not in samples exposed to 43% R.H. 2 H NMR of those samples collected at temperatures from 313 K to 173 K show the presence of multiple sites. Near room temperature, the spectra contain a narrow resonance for mobile water undergoing rapid isotropic motion, and a broader symmetrical resonance probably due to a combination of more dynamically restricted water molecules and 2 H exchanged onto phenolic and alcohol functional groups undergoing rapid anisotropic motion. The 43% R.H. samples also yield a broader, quadrupole-dominated, resonance for 2 H exchanged onto functional groups of the NOM. With decreasing temperature the resonances for dynamically restricted water molecules and 2 H exchanged onto phenolic and alcohol functional groups become broader, reflecting a decreasing rate of exchange between the water molecules and functional groups and a decreasing rate of reorientation of the 2 H 2 O molecules. The formation of ice-1h is directly reflected in the 2 H spectra of the 100% R.H. samples as a resonance with a quadrupole coupling constant (QCC) of ~180 kHz. For the 43% R.H. samples, there is also a broad, poorly resolved resonance with typical QCCs of ~180 kHz for which the relative signal intensity increases with decreasing temperature. This signal represents 2 H 2 O molecules that are not crystallized in ice-1h but have greatly reduced reorientation frequencies at low temperature and a hydrogen bonding network with hydrogen bond strengths similar to, but somewhat weaker than, ice-1h. Such molecules are also likely to be present in the 100% R.H. samples. At both R.H.s, some of the 2 H 2 O molecules do not freeze and retain their isotropic motion down to 173 K, the lowest temperature investigated.

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

Advances in rapid compression machine studies of low- and intermediate-temperature autoignition phenomena

Rapid compression machines (RCMs) are widely-used to acquire experimental insights into fuel autoignition and pollutant formation chemistry, especially at conditions relevant to current and future combustion technologies. RCM studies emphasize important experimental regimes, characterized by low- to intermediate-temperatures (600–1200 K) and moderate to high pressures (5–80 bar). At these conditions, which are directly relevant to modern combustion schemes including low temperature combustion (LTC) for internal combustion engines and dry low emissions (DLE) for gas turbine engines, combustion chemistry exhibits complex and experimentally challenging behaviors such as the chemistry attributed to cool flame behavior and the negative temperature coefficient regime. Challenges for studying this regime include that experimental observations can be more sensitive to coupled physical-chemical processes leading to phenomena such as mixed deflagrative/autoignitive combustion. Experimental strategies which leverage the strengths of RCMs have been developed in recent years to make RCMs particularly well suited for elucidating LTC and DLE chemistry, as well as convolved physical-chemical processes. Specifically, this work presents a review of experimental and computational efforts applying RCMs to study autoignition phenomena, and the insights gained through these efforts. A brief history of RCM development is presented towards the steady improvement in design, characterization, instrumentation and data analysis. Novel experimental approaches and measurement techniques, coordinated with computational methods are described which have expanded the utility of RCMs beyond empirical studies of explosion limits to increasingly detailed understanding of autoignition chemistry and the role of physical-chemical interactions. Fundamental insight into the autoignition chemistry of specific fuels is described, demonstrating the extent of knowledge of low-temperature chemistry derived from RCM studies, from simple hydrocarbons to multi-component blends and full-boiling range fuels. In conclusion, emerging needs and further opportunities are suggested, including investigations of under-explored fuels and the implementation of increasingly higher fidelity diagnostics.

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

NEML2: An efficient and modular multiphysics constitutive modeling library for hybrid computing environments

This paper presents NEML2, an open-source, high-performance library developed for constitutive material modeling, designed to support the flexible and modular development of models for complex material behavior. Building on the foundational structure of its predecessor, NEML, the NEML2 library introduces significant improvements, including enhanced vectorization, automatic differentiation, and seamless integration with PyTorch, facilitating the application of machine learning techniques in material simulations. NEML2 provides a C++ backend with Python bindings, enabling users to create custom material models that can be executed efficiently on both CPU and GPU platforms. The library also supports coupling with Multiphysics simulation frameworks like MOOSE, making it suitable for realistic simulations involving coupled physical processes. Rigorous quality assurance through unit and regression testing ensures the reliability of results, while the extensible, user-friendly design encourages collaboration and reproducibility across the scientific community. This paper provides an overview of NEML2’s architecture, core features, and applications, highlighting its impact on accelerating material qualification and advancing computational methods in materials science.

GPU↗