Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Shape function”

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,567 records · Page 87
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Response of soil nutrient pools and microbiomes to recurrent wildfire disturbance and varying burn severities in a mixed conifer forest

Wildfire is a pervasive disturbance in mixed-conifer forests, yet the relative influence of fire recurrence versus burn severity on soil biogeochemistry and microbial communities remains poorly quantified. We examined a natural gradient of fire history (0–3 prior fires) and burn severity (low–high) spanning 50 yr in a mixed-conifer ecosystem to assess how repeated fire shapes soil carbon (C) and nitrogen (N) pools, their isotopic signatures, mineral and particulate fractions, microbial community composition, carbon-use, CO₂ fluxes, and vegetation cover. Successive fires produced progressively higher bare-ground percentages and lower tree cover, which were tightly linked to declines in microbial diversity and reductions bulk %C, and %N. δ 13 C increased with fire frequency, indicating preferential loss of labile C through combustion or enhanced microbial oxidation, thereby explaining the observed net soil-C decline. Conversely, δ 15 N decreased and pH increased as tree density declined, reflecting altered N cycling and reduced acidification in post-fire soils. Fire recurrence, more than severity, corresponded with a marked shift in the bacterial community: for example, Xanthobacteraceae—key N-fixers and C-cyclers—diminished, while N-fixing Bacillaceae increased, underscoring the tightly coupled nature of soil nutrient dynamics and microbiome composition after repeated burns. Our results demonstrate that fire recurrence appears to be a stronger driver of post-fire soil ecosystem responses in this mixed-conifer forest, influencing both abiotic nutrient pools and the functional potential of the soil microbiome. These findings provide a more enhanced assessment and understanding to date of the biogeochemical consequences of repeated wildfire disturbance that can be used to inform management strategies aimed at preserving soil health in fire-prone landscapes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Advanced Functional Materials ZnO@ZIF-8 Heterostructures for Neuromorphic Synaptic Devices

Metal-Organic Frameworks have shown great promise for memristive synaptic devices because of their adaptable electrical properties and inherent porosity, which help in efficient ion transport and storage of functional molecules. However, the limited conductivity and sensitivity to moisture associated with most MOF materials hinder their use in electronic applications. In this study, MOFs derived from metabolites have been synthesized, using ZnO as a self-sacrificial metallic source to grow Zeolitic Imidazolate Framework-8 (ZIF-8), forming ZnO@ZIF-8 heteronanostructures for neuromorphic applications. These hybrids ZnO@ZIF-8 act as promising candidates for advanced electronic applications. Through systematic investigations, including a comparison of ZnO@ZIF-8 bilayer devices with ZnO single-layer counterparts, a significant discovery emerged. The bilayer synaptic devices outperformed their single layer counterparts in electrical characteristics, as evidenced by Current-Voltage curve analysis and Long-Term Potentiation/Long Term Depression (LTP/LTD) measurements. Additionally, the ZnO layer thickness (10, 30, 50 nm) has been emphasized and the effects of annealing at 300°C in a vacuum were carefully studied for Neuromorphic Synaptic Devices. Remarkably, the 10 nm ZnO layer exhibited the lowest set voltage and a consistent On/Off ratio. This work marks a significant advancement in understanding bilayer synaptic devices, highlighting the crucial role of ZIF-8 in shaping the future of neuromorphic systems

Kolleboyina, Jayaramulu↗
DOE OSTI2020

Low–severity fire as a mechanism of organic matter protection in global peatlands: Thermal alteration slows decomposition

Worldwide, regularly-recurring wildfires shape many peatland ecosystems to the extent that fire-adapted species often dominate plant communities, suggesting that wildfire is an integral part of peatland ecology rather than an anomaly. The most destructive blazes are smoldering fires that are usually initiated in periods of drought and can combust entire peatland carbon stores. However, peatland wildfires more typically occur as low-severity surface-burns that arise in the dormant season when vegetation is desiccated, and soil moisture is high. In such low-severity fires, surface layers experience flash heating, but there is little loss of underlying peat to combustion. This study examines the potential importance of such processes in several peatlands that span a gradient from hemiboreal to tropical ecozones and experience a wide range of fire return intervals. We show that low-severity fires can increase the pool of stable soil carbon by thermally altering the chemistry of soil organic matter (SOM), thereby reducing rates of microbial respiration. Using XPS and FTIR, we demonstrate that low-severity fires significantly increase the degree of carbon-condensation and aromatization of SOM functional groups, particularly on the surface of peat aggregates. Laboratory incubations show lower CO 2 emissions from peat subjected to low-severity fire, and predict lower cumulative CO 2 emissions from burned peat after one to three years. Also, low-severity fires reduce the temperature sensitivity (Q 10 ) of peat, indicating that these fires can inhibit microbial access to SOM. The increased stability of thermally-altered SOM may allow a greater proportion of organic matter to survive vertical migration into saturated and anaerobic zones of peatlands where environmental conditions physiochemically protect carbon stores from decomposition for thousands of years. Furthermore, across latitudes, low-severity fire is an overlooked factor influencing carbon cycling in peatlands, which is relevant to global carbon budgets as climate-change alters fire regimes worldwide.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

De novo design of small beta barrel proteins

Small beta barrel proteins are attractive targets for computational design because of their considerable functional diversity despite their very small size (<70 amino acids). However, there are considerable challenges to designing such structures, and there has been little success thus far. Because of the small size, the hydrophobic core stabilizing the fold is necessarily very small, and the conformational strain of barrel closure can oppose folding; also intermolecular aggregation through free beta strand edges can compete with proper monomer folding. Here, we explore the de novo design of small beta barrel topologies using both Rosetta energy–based methods and deep learning approaches to design four small beta barrel folds: Src homology 3 (SH3) and oligonucleotide/oligosaccharide-binding (OB) topologies found in nature and five and six up-and-down-stranded barrels rarely if ever seen in nature. Both approaches yielded successful designs with high thermal stability and experimentally determined structures with less than 2.4 Å rmsd from the designed models. Using deep learning for backbone generation and Rosetta for sequence design yielded higher design success rates and increased structural diversity than Rosetta alone. The ability to design a large and structurally diverse set of small beta barrel proteins greatly increases the protein shape space available for designing binders to protein targets of interest.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Forest Carbon Storage in the Western United States: Distribution, Drivers, and Trends

Abstract Forests are a large carbon sink and could serve as natural climate solutions that help moderate future warming. Thus, establishing forest carbon baselines is essential for tracking climate‐mitigation targets. Western US forests are natural climate solution hotspots but are profoundly threatened by drought and altered disturbance regimes. How these factors shape spatial patterns of carbon storage and carbon change over time is poorly resolved. Here, we estimate live and dead forest carbon density in 19 forested western US ecoregions with national inventory data (2005–2019) to determine: (a) current carbon distributions, (b) underpinning drivers, and (c) recent trends. Potential drivers of current carbon included harvest, wildfire, insect and disease, topography, and climate. Using random forests, we evaluated driver importance and relationships with current live and dead carbon within ecoregions. We assessed trends using linear models. Pacific Northwest (PNW) and Southwest (SW) ecoregions were most and least carbon dense, respectively. Climate was an important carbon driver in the SW and Lower Rockies. Fire reduced live and increased dead carbon, and was most important in the Upper Rockies and California. No ecoregion was unaffected by fire. Harvest and private ownership reduced carbon, particularly in the PNW. Since 2005, live carbon declined across much of the western US, likely from drought and fire. Carbon has increased in PNW ecoregions, likely recovering from past harvest, but recent record fire years may alter trajectories. Our results provide insight into western US forest carbon function and future vulnerabilities, which is vital for effective climate change mitigation strategies.

Environmental Sciences & Ecology↗
DOE OSTI2022

A method to predict texture effect on ion beam channeling analysis of polycrystals and the application to study the mosaic spreading effect in highly oriented pyrolytic graphite

We propose a method to convert the channeling Rutherford backscattering spectrum yield map of a single crystal to a polycrystal through a matrix rotation technique. The rotation matrix is determined by the deviation of the crystal axial direction from the original z axis. The final yield map is created after averaging the rotated yields using the texture function as the weight factor. For highly oriented pyrolytic graphite (HOPG) exhibiting mosaic spread, the method leads to a Gaussian kernel averaging of the map obtained from a single crystal. The yield map of a single crystal is obtained by a simulation of ion trajectories in a potential field described by Moliere screened Coulomb potentials. Yield maps are calculated under various σ values (standard deviations of mosaic spread). The simulated results are compared with experimental results obtained using 1.2 MeV alpha particle. σ is extracted through the best fitting, demonstrating that the method can be used to obtain texture details. The effects of mosaic spread on minimum yield χ min and the half-width at half maximum of angular scans ψ ½ are systematically modeled and compared with previous theoretical equations. The study also shows that previous theoretical equations are valid only at small σ values. The proposed method can be applied to any type of polycrystal and is not limited to HOPG. It provides near-surface mosaic spread and crystallography information with a longitudinal depth resolution of tens of nanometers and is not influenced by grain shapes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Skyrme–Hartree–Fock–Bogoliubov mass models on a 3D mesh: IIb. Fission properties of BSkG2

Large-scale models of nuclear structure are currently the only way to provide consistent datasets for the many properties of thousands of exotic nuclei that are required by nucleosynthesis simulations. In [W. Ryssens et al., Eur. Phys. J. A 58, 246 (2022)], we recently presented the new BSkG2 model based on an energy density functional of the Skyrme type. Relying on a flexible three-dimensional coordinate representation of the nucleus, the model takes into account both triaxial deformation and time-reversal symmetry breaking. BSkG2 achieves a state-of-the-art global description of nuclear ground state (g.s.) properties and reproduces in particular the known masses with a root-mean-square (rms) deviation of 678 keV. Moving beyond g.s. properties, the model also reproduces all empirical values for the primary and secondary barriers as well as isomer excitation energies of actinide nuclei with rms deviations below 500 keV, i.e. with unprecedented accuracy. Here we discuss in detail the extension of our framework to the calculation of the fission barriers of 45 actinide nuclei, including odd-mass and odd-odd systems. We focus in particular on the impact of symmetry breaking which is key to the accuracy of the model: we allow systematically for axial, reflection and time-reversal symmetry breaking. The effect of the latter on the fission properties of odd-mass and odd-odd nuclei is small, but we find that allowing for shapes with triaxial or octupole deformation, as well as shapes with both, is crucial to achieving this accuracy. Here, the numerical accuracy of our coordinate space approach, the variety of nuclear configurations explored and the simultaneous successful description of fission properties and known masses makes BSkG2 the tool of choice for the large-scale study of nuclear structure.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Rational Design of Novel Biomimetic Sequence-Defined Polymers for Mineralization Applications

Silica biomineralization is a naturally occurring process, wherein organisms use proteins and other biological structures to direct the formation of complex, hierarchical nanostructures. Discovery and characterization of such proteins and their underlying mechanisms spurred significant efforts to identify routes for biomimetic mineralization that reproduce the exquisite shapes and size selectivities found in nature. A common strategy has been the use of short peptide sequences with chemistry mimicking those found in natural systems, such as the use of the silaffin-derived R5 peptide. While progress has been made using this approach, there are many limitations that have prevented breakthroughs in biomimicry. To advance our ability to use charged macromolecules for silica formation, we propose to use sequence-defined synthetic polymers known as peptoids, or N-substituted polyglycines, which present significant capability for the precise tuning of sequence and structure beyond what can often be achieved with peptides alone. This study presents a computationally predicted design of these polymers that leads to the controlled formation of silica nanomaterials. We investigate surface adsorption and the mineralization process through analysis of binding mechanisms and energetics of the R5 system. Next, we synthesized two R5-inspired peptoids and validated our prediction in the design of mineralization polymers through characterization using surface plasmon resonance and electron microscopy. Here, this computationally guided study holds great promise for designing new sequences with unprecedented control of the placement of chemical functional groups, thus allowing for further unraveling of silicification mechanisms and the eventual design of sequence-defined synthetic polymers leading to the predictive synthesis of nanostructured functional materials.

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

First Constraints on Growth Rate from Redshift-space Ellipticity Correlations of SDSS Galaxies at 0.16 < z < 0.70

We report the first constraints on the growth rate of the universe, f(z)σ 8 (z), with intrinsic alignments (IAs) of galaxies. We measure the galaxy density-intrinsic ellipticity cross-correlation and intrinsic ellipticity autocorrelation functions over 0.16 < z < 0.7 from luminous red galaxies (LRGs) and LOWZ and CMASS galaxy samples in the Sloan Digital Sky Survey (SDSS) and SDSS-III BOSS survey. We detect clear anisotropic signals of IA due to redshift-space distortions. By combining measured IA statistics with the conventional galaxy clustering statistics, we obtain tighter constraints on the growth rate. The improvement is particularly prominent for the LRG, which is the brightest galaxy sample and known to be strongly aligned with underlying dark matter distribution; using the measurements on scales above 10 h -1 Mpc, we obtain $f{\sigma }_{8}={0.5196}_{-0.0354}^{+0.0352}$ (68% confidence level) from the clustering-only analysis and $f{\sigma }_{8}={0.5322}_{-0.0291}^{+0.0293}$ with clustering and IA, meaning 19% improvement. The constraint is in good agreement with the prediction of general relativity, f σ 8 = 0.4937 at z = 0.34. For LOWZ and CMASS samples, the improvement of constraints on f σ 8 is found to be 10% and 3.5%, respectively. Our results indicate that the contribution from IA statistics for cosmological constraints can be further enhanced by carefully selecting galaxies for a shape sample.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

The Lyman- α forest from LBGs: First 3D correlation measurement with DESI and prospects for cosmology

The Lyman-α (Lyα) forest is a key tracer of large-scale structure at redshifts z > 2, traditionally studied using the spectra of luminous but relatively rare quasars. In this work, we explore the viability of using the fainter yet significantly more abundant Lyman Break Galaxies (LBGs) as alternative background sources for Lyα forest studies. We analyze 4,151 Lyα forest skewers extracted from LBG spectra obtained in the DESI pilot surveys conducted in the COSMOS and XMM-LSS fields. From this dataset, we present the first measurement of the Lyα forest auto-correlation function derived exclusively from LBG spectra, probing comoving separations up to 48 h -1 Mpc at an effective redshift of z eff = 2.70. The measured LBG Lyα forest auto-correlation is consistent with that derived from DESI DR2 quasar Lyα forest spectra at a comparable redshift, validating the use of LBGs as reliable background sources for Lyα forest analyses. In addition, we measure the cross-correlation between the LBG Lyα forest and the positions of 13,362 galaxies, demonstrating that this observable serves as a sensitive diagnostic for assessing the precision and accuracy of galaxy redshift estimates, and for identifying and correcting systematic offsets. Finally, using both synthetic LBG spectra and Fisher matrix forecasts, we show that a future wide-area survey covering ∼5,000 deg 2 , targeting 1,000 LBGs per square degree at signal-to-noise levels comparable to our sample, could enable LBG-based Lyα forest baryon acoustic oscillation (BAO) measurements with expected uncertainties of σ α ISO = 0.4% (isotropic) and σ α AP = 1.3% (Alcock-Paczynski). This performance is further enhanced when combining the BAO analysis with a Lyα forest Full Shape (FS) approach, yielding a predicted uncertainty of σ α ISO FS = 0.6%. These results open a new avenue for precision cosmology at high redshift using the Lyα forest in dense LBG samples.

Lyman alpha forest↗
DOE OSTI2025

Measuring and unbiasing the BAO shift in the Ly α forest with AbacusSummit

ABSTRACT The Dark Energy Spectroscopic Instrument (DESI) places sub- per cent constraints on measurements of the Baryon Acoustic Oscillation (BAO) scaling parameters from the Ly $\alpha$ forest. However, no systematic error budget stemming from non-linearities in the three-dimensional clustering of the Ly $\alpha$ forest is included in the DESI-Ly $\alpha$ analysis. In this work, we measure the size of the shift of the BAO peak using large Ly $\alpha$ forest mocks produced on the N-body simulation suite AbacusSummit, which adopt the Fluctuating–Gunn–Peterson Approximation (FGPA). Specifically, we measure the Ly $\alpha$ autocorrelation and the Ly $\alpha$-quasar cross-correlation functions. To mitigate the noise, we adopt a linear control variates technique, reducing the error bars by a factor of up to $\sim \sqrt{50}$ on large scales. From the autocorrelation, we detect a small positive shift in radial direction of $\Delta \alpha _{\parallel }= 0.35~{{\ \rm per\ cent}}$ at the 3$\sigma$ level and virtually no shift in the transverse direction, $\alpha _\perp$. From the cross-correlation, we see a similar shift to $\Delta \alpha _\parallel$, albeit with larger error bars, and a small negative shift, $\Delta \alpha _{\perp }=\sim$0.25 per cent, at the 2$\sigma$ level. We also make a connection with the Ly $\alpha$ forest effective field theory (EFT) framework and find that the one-loop EFT power spectrum yields unbiased measurements of the BAO shift parameters in radial and transverse direction for Ly $\alpha$ auto- and the Ly $\alpha$-quasar cross-correlation measurements. When using the one-loop EFT framework, we find that we can recover the BAO parameters without a shift, which has important implications for future Ly $\alpha$ forest analyses based on EFT. This work paves the way for novel full-shape analyses of the currently observing DESI and future surveys such as the PFS, WEAVE-QSO, and 4MOST.

Hadzhiyska, Boryana↗
DOE OSTI2022

Additive manufacturing of LiNi 1/3 Mn 1/3 Co 1/3 O 2 battery electrode material via vat photopolymerization precursor approach

Additive manufacturing, also called 3D printing, has the potential to enable the development of flexible, wearable and customizable batteries of any shape, maximizing energy storage while also reducing dead-weight and volume. In this work, for the first time, three-dimensional complex electrode structures of high-energy density LiNi 1/3 Mn 1/3 Co 1/3 O 2 (NMC 111) material are developed by means of a vat photopolymerization (VPP) process combined with an innovative precursor approach. This innovative approach involves the solubilization of metal precursor salts into a UV-photopolymerizable resin, so that detrimental light scattering and increased viscosity are minimized, followed by the in-situ synthesis of NMC 111 during thermal post-processing of the printed item. The absence of solid particles within the initial resin allows the production of smaller printed features that are crucial for 3D battery design. The formulation of the UV-photopolymerizable composite resin and 3D printing of complex structures, followed by an optimization of the thermal post-processing yielding NMC 111 is thoroughly described in this study. Based on these results, this work addresses one of the key aspects for 3D printed batteries via a precursor approach: the need for a compromise between electrochemical and mechanical performance in order to obtain fully functional 3D printed electrodes. In addition, it discusses the gaps that limit the multi-material 3D printing of batteries via the VPP process.

25 ENERGY STORAGE↗
DOE OSTI2025

Next generation Arctic vegetation maps: Aboveground plant biomass and woody dominance mapped at 30 m resolution across the tundra biome

The Arctic is warming faster than anywhere else on Earth, placing tundra ecosystems at the forefront of global climate change. Plant biomass is a fundamental ecosystem attribute that is sensitive to changes in climate, closely tied to ecological function, and crucial for constraining ecosystem carbon dynamics. However, the amount, functional composition, and distribution of plant biomass are only coarsely quantified across the Arctic. Therefore, we developed the first moderate resolution (30 m) maps of live aboveground plant biomass (g m −2 ) and woody plant dominance (%) for the Arctic tundra biome, including the mountainous Oro Arctic. We modeled biomass for the year 2020 using a new synthesis dataset of field biomass harvest measurements, Landsat satellite seasonal synthetic composites, ancillary geospatial data, and machine learning models. Additionally, we quantified pixel-wise uncertainty in biomass predictions using Monte Carlo simulations and validated the models using a robust, spatially blocked and nested cross-validation procedure. Observed plant and woody plant biomass values ranged from 0 to ∼6000 g m −2 (mean ≈ 350 g m −2 ), while predicted values ranged from 0 to ∼4000 g m −2 (mean ≈ 275 g m −2 ), resulting in model validation root-mean-squared-error (RMSE) ≈ 400 g m −2 and R 2 ≈ 0.6. Our maps not only capture large-scale patterns of plant biomass and woody plant dominance across the Arctic that are linked to climatic variation (e.g., thawing degree days), but also illustrate how fine-scale patterns are shaped by local surface hydrology, topography, and past disturbance. By providing data on plant biomass across Arctic tundra ecosystems at the highest resolution to date, our maps can significantly advance research and inform decision-making on topics ranging from Arctic vegetation monitoring and wildlife conservation to carbon accounting and land surface modeling.

Climate change↗
DOE OSTI2022

Tunable Solid Acid Catalyst Thin Films Prepared by Atomic Layer Deposition

Solid acid catalysts, including zeolites and amorphous silica-aluminas (ASAs), are industrially important materials widely used in the fuel and petrochemical industries. The versatility of zeolites is due to the Brønsted acidity of the bridging hydroxyl and shape selectivity that can be tailored during and after synthesis. This is in contrast to amorphous silica-alumina, where tailoring acidity is a major challenge as the Brønsted acid structure in ASA is still debated. In both cases, however, the pore size and acidity cannot be tuned independently, and this is particularly limiting in the application of biomass conversion, where zeolite pores are too small for the molecules of interest. Herein, we present a method using atomic layer deposition (ALD) to prepare thin films of solid acid materials where the ratio of Brønsted to Lewis acid sites can be tuned precisely. This capability, combined with the sub-nm pore size control afforded by ALD yields a powerful and flexible method for synthesizing solid acid catalysts inside virtually any mesoporous host. We demonstrate the utility of these materials in two acid-catalyzed reactions relevant to biomass conversion: (1) Meerwein–Ponndorf–Verley–Oppenauer (MPVO) reaction and dehydration of fructose and (2) cascade reaction of glucose to 5-hydroxymethylfurfural. Finally, we propose a plausible structure for the Brønsted acid sites in our materials based on infrared spectroscopy and solid-state nuclear magnetic resonance measurements and density functional theory calculations and argue that this same structure might apply to conventional ASAs as well.

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

Towards the spatial resolution of metalloprotein charge states by detailed modeling of XFEL crystallographic diffraction

Oxidation states of individual metal atoms within a metalloprotein can be assigned by examining X-ray absorption edges, which shift to higher energy for progressively more positive valence numbers. Indeed, X-ray crystallography is well suited for such a measurement, owing to its ability to spatially resolve the scattering contributions of individual metal atoms that have distinct electronic environments contributing to protein function. However, as the magnitude of the shift is quite small, about +2 eV per valence state for iron, it has only been possible to measure the effect when performed with monochromated X-ray sources at synchrotron facilities with energy resolutions in the range 2–3 × 10 −4 (Δ E / E ). This paper tests whether X-ray free-electron laser (XFEL) pulses, which have a broader bandpass (Δ E / E = 3 × 10 −3 ) when used without a monochromator, might also be useful for such studies. The program nanoBragg is used to simulate serial femtosecond crystallography (SFX) diffraction images with sufficient granularity to model the XFEL spectrum, the crystal mosaicity and the wavelength-dependent anomalous scattering factors contributed by two differently charged iron centers in the 110-amino-acid protein, ferredoxin. Bayesian methods are then used to deduce, from the simulated data, the most likely X-ray absorption curves for each metal atom in the protein, which agree well with the curves chosen for the simulation. The data analysis relies critically on the ability to measure the incident spectrum for each pulse, and also on the nanoBragg simulator to predict the size, shape and intensity profile of Bragg spots based on an underlying physical model that includes the absorption curves, which are then modified to produce the best agreement with the simulated data. This inference methodology potentially enables the use of SFX diffraction for the study of metalloenzyme mechanisms and, in general, offers a more detailed approach to Bragg spot data reduction.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

An Eulerian multimaterial framework for simulating high-explosive aquarium tests

Aquarium tests of cylindrical high-explosive charges provide optical data of the detonation front velocity and shape, propagation of the shock wave in the surrounding water, and expansion rates of the detonation products behind the front. Data from aquarium experiments is often used for calibration of reactive burn models based on phenomenological equations of state (EOS) and reaction rate laws. This paper presents a multimaterial numerical modeling framework to solve the 2D axisymmetric reactive Euler equations for high-explosive aquarium tests, in particular for ammonium nitrate - fuel oil (ANFO) explosives. An extension of the Ghost Fluid Method (GFM) is used to handle the dynamic material interfaces for the ANFO explosion products, the charge-confining material (polymethyl methacrylate PMMA), and the surrounding water. This study analyzes the sensitivity of calculations (both computational efficiency and numerical accuracy) to different algorithms for the material interface models including the original GFM versus Riemann solver-based strategies. A novel method for defining the left and right states in the interfacial Riemann problem eliminates the need for sorting or nodal interpolation during the projection along the material interface. Numerical tests indicate that populating the interface node values using the Riemann solution mitigate the overheating error observed in steady-state calculations. Solution convergence and computational efficiency are explored as a function of the spatial and temporal order of the schemes. Results from the computational model with analytical equations of state and fitted reaction rate parameters show very good quantitative agreement with experimentally observed detonation front velocity, reaction products expansion, and shock wave propagation in the surrounding water for a cylindrical ANFO charge encased in PMMA. Finally, the proposed modeling framework, in conjunction with experimental tests, provides a reliable tool to assess equations of state and reaction rate expressions for reactive burn models of confined high explosives.

42 ENGINEERING↗
DOE OSTI2022

A critical comparison of CH∙∙∙π versus π∙∙∙π interactions in the benzene dimer: obtaining benchmarks at the CCSD(T) level and assessing the accuracy of lower scaling methods

In this study, we have established CCSD(T)/CBS (Complete Basis Set) limits for 3 stationary points on the benzene dimer potential energy surface, corresponding to the π∙∙∙π (parallel displaced or PD(C 2h ), minimum) and CH∙∙∙π (T-shaped or T(C 2v ), transition state) and tilted T-shaped (or TT(Cs), minimum) bonding scenarios considering both the structure and binding energy. The CCSD(T)/CBS binding energies are -2.65 ± 0.02 (PD), -2.74 ± 0.03 (T), and -2.83 ± 0.01 kcal mol -1 (TT). To this end, the CH∙∙∙π is ~0.2 kcal mol -1 stronger than the π∙∙∙π interaction, whereas the tilting of the CH donating benzene molecule with respect to the other benzene is worth 0.1 kcal mol -1 . As previously discussed in the literature, the MP2 level of theory does not provide a close match for either the energy or structure, yet the SCS-MP2 yields structures in excellent agreement with respect to the CCSD(T) result. It is found that the SCS-MI-MP2 also gives optimized structures very close to SCS-MP2 (within ~0.01 Å of the benchmark). Despite the closer match in structure, the spin-biased MP2 methods (SCS-, SCS-MI-, and SOS-MP2) incorrectly predict the relative stabilities of the isomers. That said, none of the spin biased MP2 methods offers a good compromise between energy and structure for the systems examined. Finally, the CCSD(T)/CBS benchmarks were used to assess the performance of 13 DFT functionals selected from different rungs of Jacob's ladder. Several functionals such as TPSS-D3, B3LYP-D3, B97-D, B97-D3, and B2PLYP-D3 provided a good description of the binding energies for both CH∙∙∙π and π∙∙∙π interactions, yielding values within 6% of the CCSD(T)/CBS benchmark values. Unlike the MP2 methods, these functionals correctly predict the relative stability of the PD(C 2h ) and T(C 2v ) dimers. Further, we find that there is no systematic improvement as Jacob's ladder is ascended (increased complexity of functional). The best functionals that result in a good compromise between structure and energy accuracy are B97-D3 and B2PLYP-D3 for both the CH∙∙∙π and π∙∙∙π interaction. Despite the impressive performance of these functionals, a challenge that remains is ensuring the transferability of these density functionals in accurately describing the interaction between dimers of larger aromatic molecules, the latter requiring high-level benchmarks for these systems.

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

wavess 1.2: presenting an HLA-aware within-host virus sequence simulation framework

Motivation Understanding how virus sequences are shaped by selection can inform vaccine design and transmission inference. Modeling within-host evolution to interrogate these questions requires a detailed mechanistic framework that accurately captures sequence diversification. The CD8 + cytotoxic T-lymphocyte (CTL) response plays an important role in immune-mediated selection and can leave strong signatures in virus sequences; however, existing sequence-based within-host virus modeling frameworks do not explicitly include a human leukocyte antigen (HLA)-aware CTL response. Results We extended our previously published within-host sequence evolution simulator, wavess, to include an explicit CTL response, and share a method for identifying HLA-specific CTL epitopes given a founder virus sequence. We also updated the model to permit a variable recombination rate, which allows for modeling non-adjacent genes, segmented genomes, and recombination hotspots. These extensions to wavess allow for more accurate simulation of viruses and virus genes, particularly in regions of the genome where the immune response is dominated by CTLs (rather than antibodies). It also provides the foundation for investigations of how these newly-added biological mechanisms influence within-host evolution. Availability and implementation The core of wavess is written in Python 3, with helper functions written in R. It is available at https://github.com/MolEvolEpid/wavess.

60 APPLIED LIFE SCIENCES↗