SEARCH · Search NASA

Results for “Structural Model Decomposition”

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.

SourceResearch connectionFrom yearThrough year

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

At least 289 records · Page 16

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Leveraging High-resolution Molecular Composition of Soil Organic Matter to Enhance Carbon Cycling Modeling

Soils store more carbon than the atmosphere and vegetation combined, yet Earth system models still struggle to predict how this vast reservoir will respond to environmental change. A central limitation is that most soil biogeochemical models represent organic matter using bulk conceptual pools or chemically homogeneous fractions, preventing direct use of rapidly expanding molecular-scale datasets. Here we develop and test a new soil decomposition framework that explicitly integrates high-resolution information on organic matter composition. First, we construct a molecularly informed litter decomposition module in which plant inputs are partitioned into five functional compound classes—carbohydrates, proteins, lignin-like aromatics, lipids, and carbonyls—using a molecular mixing model calibrated to solid-state 13 C Nuclear Magnetic Resonance (NMR) spectra. Class-specific kinetics, lignin-dependent physical protection, and substrate-driven microbial carbon use efficiency allow the module to capture metabolic tradeoffs associated with enzyme production and nutrient limitation. We then embed this litter module within a microbially explicit whole-soil model that tracks the transformation of these compound classes through particulate organic matter, dissolved organic matter, mineral-associated organic matter, and microbial biomass. High-resolution Fourier Transform Ion Cyclotron Resonance mass spectrometry (FTICR-MS) data are used to link internal pools to measurable soil organic matter fractions and to constrain key process parameters. Applications at soil-core and ecosystem scales demonstrate that the new model reproduces observed soil respiration dynamics while providing mechanistic attribution of CO 2 fluxes to specific chemical classes and pools. Compared to existing frameworks such as the Community Land Model soil biogeochemistry module and the Millennial model, our approach maintains competitive predictive skill while substantially improving interpretability and opportunities for data–model integration. This work illustrates a viable pathway for leveraging molecular-scale observations to reduce structural uncertainty in soil carbon–climate feedback projections.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Warming and Increased Respiration Have Transformed an Alpine Steppe Ecosystem on the Tibetan Plateau From a Carbon Dioxide Sink Into a Source

Cold region ecosystems store vast amounts of soil organic carbon (C), which upon warming and decomposition can affect the C balance and potentially change these ecosystems from C sinks to carbon dioxide (CO 2 ) sources. We quantified the decadal year-round CO 2 flux from an alpine steppe-ecosystem on the Tibetan Plateau using eddy covariance and automatic chamber approaches during a period of significant warming (0.13°C per 10 years; and 0.18°C in the non-growing season alone: 1st October to next 30th April). The results showed that ongoing climate change, mainly warming within the topsoil layers, is the main reason for the site’s change from a sink for to a source of CO 2 in the atmosphere. Non-growing-season ecosystem respiration accounted for 51% of the annual ecosystem respiration and has increased significantly. The growing seasons (1st May to 30th September) were consistent CO 2 sink periods without significant changes over the study period. Observations revealed high-emission events from the end of the non-growing season to early in the growing season (1st March to fifteenth May), which significantly (p < 0.01) increased at a rate of 22.6 g C m -2 decade -1 , ranging from 14.6 ± 10.7 g C m -2 yr -1 in 2012 to 35.3 ± 12.1 g C m -2 yr -1 in 2017. Further, structural equation modeling suggested that active layer warming was the key factor in explaining changes in ecosystem respiration, leading to significant changes in net ecosystem exchange over the period 2011–2020 and indicated that these changes have already transformed the ecosystem from a CO 2 sink into a source. These results can be used to improve our understanding of the sensitivity of ecosystem respiration to increased warming during the non-growing period.

58 GEOSCIENCES↗

DOE OSTI2021

Diol isomer revealed as a source of methyl ketene from propionic acid unimolecular decomposition

Carboxylic acids are important combustion intermediates, especially regarding the combustion of oxygenates such as ethyl esters. The purpose of this study is to discern the unimolecular decomposition pathways of one such carboxylic acid, propionic acid, using microreactor flow experiments with line-tunable photoionization mass spectrometry and infrared spectroscopy, along with high-level electronic structure calculations (CCSD(T)/cc-pV∞Z//M06-2X/cc-pVTZ level of theory) and master equation theory. Here, microreactor experiments were performed at 300–1500 K, pressures decreasing from roughly 300 Torr to vacuum pressure, and around 100 μs residence times. Primary products are methyl ketene, ketene, ethylene, and methyl radical. Theory suggests the two lowest energy bond fissions are active at these conditions and responsible for the production of ketene, methyl radical, and some ethylene. Importantly, theory revealed the significance of the isomerization of propionic acid to propene-1,1-diol, and the subsequent dehydrogenation reaction of the diol as an alternative explanation for the unimolecular formation of methyl ketene. This is predicted to be the dominant thermal decomposition pathway at lower temperatures (up to ~850 K). Line-tunable VUV photons allowed for isomer resolution of the products and showed no evidence of the diol of propionic acid surviving until detection, suggesting propene-1,1-diol decomposes either to methyl ketene and water rapidly, fragments upon ionization to a distonic ion inseparable from methyl ketene through our detection methods, or never stabilizes as propionic acid well-skips directly to methyl ketene. Infrared spectroscopy showed no evidence of the decarboxylation of propionic acid at these conditions. Updated rate constants and branching ratios for propionic acid decomposition are calculated and provided for future modeling studies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

The evolution of coal porosity during pyrolysis

Gasification of coal, municipal waste, or other organic materials is a potential hydrogen source that entails complex thermal decomposition and transport processes. This study provides a multiscale analysis of these processes for sub-bituminous (Usibelli, Healy, Alaska) and lignite (Center, North Dakota) coals and provides data useful for process design. The chemistry, mineralogy, and pore structures of pyrolyzed coal and their evolution with thermal decomposition are discussed. Samples pyrolyzed at 200–1000 °C were analyzed by small-angle neutron scattering; ultra-small, small-, and wide-angle X-ray scattering; and other complementary techniques. Scanning electron microscopy showed new pores in the high-temperature-pyrolyzed material. Upon heating, the coals became progressively denser, and the concentration of hydrogen decreased. Changes in pore volume fell into three temperature ranges: an initial, low-temperature range that, for the Usibelli coal, involved an increase in overall porosity; a mid-temperature range associated with pore volume loss; and a high-temperature range associated with significant porosity increase and char formation. This transformation was paralleled by changes in fractal dimension and correlation length. The higher the pyrolysis temperature the greater the small-pore-volume fraction and overall surface area became. Pyrolysis increased the lateral size of coal crystallites, decreased the amorphous fraction, and increased the aromatics fraction and overall coal rank. Comparisons of neutron and X-ray scattering data and subsequent water uptake studies showed that pre-dried coals can re-hydrate relatively rapidly upon exposure to air, which can significantly affect the porosity calculated from small-angle-scattering data. Fits to the cumulative porosity curves provide a method for modeling the physical and chemical transformation of hydrogen-containing feedstock during gasification.

Anovitz, Lawrence {Larry} [ORNL] (ORCID:0000000226↗

DOE OSTI2022

Energetics of hydroxylbastnäsite solid solutions, La 1-$\chi$ Nd $\chi$ CO 3 OH

Bastnäsites (LnCO 3 (F,OH)) are a group of common rare earth elements (REE)-bearing minerals and are one of the primary global sources of REE. Due to the chemical similarities among REE, bastnӓsites tend to occur as solid solutions instead of end members in REE containing ores. To better understand the processes and the mechanisms of formation of such deposits, it is essential to determine the thermodynamic properties of bastnӓsites, including hydroxylbastnӓsite (LnCO 3 OH) solid solutions. In this work, we performed detailed structural and calorimetric investigations on synthetic hexagonal La–Nd hydroxylbastnӓsite (La 1-$\chi$ NdCO 3 OH, $\chi$ = 0, 0.25, 0.5, 0.75, 1) solid solutions. X-ray diffraction confirms the crystal structure of the solid solution series in the $\ P$6 space group, and pair distribution function (PDF) analysis reveals local bonding environments characterized by three different types of 9-coordinated metal-oxygen polyhedra. Unit cell parameters of La1–xNdxCO3OH exhibit a nearly linear relation with the Nd content x, suggesting a random distribution of La and Nd in the structure. Their standard enthalpies of formation (Δ$\ H$° f ) were determined by high temperature oxide melt drop solution calorimetry, from which the enthalpies of mixing (Δ$\ H$ mix ) were derived. The Δ$\ H$ mix can be fitted by a regular solution model with an interaction parameter of 12.58 ± 0.16 kJ/mol, suggesting enthalpic metastability of La1–xNdxCO3OH relative to the two endmembers. Combining entropy and enthalpy, we further estimated the Gibbs free energies of mixing (Δ$\ G$ mix ) at relevant temperatures, revealing favorable temperatures under which the intermediate La 1-$\chi$ NdCO 3 OH phases can be stabilized. Such entropy-driven stabilization, as is consistent with our geochemical modeling results, may explain the enhancement of thermal stability of the solid solutions in nature. Additionally, the temperature range constrained from this study may be used to estimate the thermal history of REE bastnӓsite deposit.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Sooting tendencies of 20 bio-derived fuels for advanced spark-ignition engines

The sooting tendencies of 20 bio-derived fuels indicated as potential blendstocks for spark-ignition engines by the Co-Optimization of Fuels & Engines (Co-Optima) Initiative are studied. The Yield Sooting Index (YSI) is used to quantify the sooting tendencies. The Co-Optima Chemical Model is used to predict the numerical YSIs for all test fuels, and these results are compared with measurements. Experimental YSIs are newly measured for 2 furans and are taken from our previous work for the other compounds. Overall, the predicted YSIs agree well with measurements, within the experimental uncertainties, except for 2,5-dimethylfuran (2,5-DMF). It is found that Polycyclic Aromatic Hydrocarbons (PAH) growth reactions have little influence on the relative soot production from each fuel and therefore on YSI predictions. A modified sensitivity coefficient formulation is proposed to evaluate reaction sensitivities specifically for YSI. This formulation is applied to the 2,5-DMF-doped flame, and six 2,5-DMF initial decomposition reactions, which lead to different major soot precursors, are identified as the most sensitive reactions. The impact of the chemical kinetic uncertainties embedded in these reactions is quantified by randomly perturbing their reaction rates within a factor of two. The resulting prediction uncertainty in the 2,5-DMF-doped flame is found to be ± 8 YSI units (±11%), which demonstrates that the test-fuel-specific decomposition reactions indeed have considerable influence on its YSI prediction. Our research suggests that more accurate kinetic parameters for fuel initial decomposition reactions can potentially improve YSI predictions effectively, without altering the YSI predictions of other fuels with significantly different molecular structures.

09 BIOMASS FUELS↗

DOE OSTI2021

Combined modeling and experiments to predict corrosion and embrittlement in dual-phase stainless steels within the MARMOT framework

Extending service of Light Water Reactors (LWRs) to beyond 60 years will demand high integrity of materials and components in the reactors. The accurate evaluation and prediction of materials performance under anticipated operating conditions are of particular importance for ensuring the safety of nuclear power plants over their extended lifetime. Stainless steels are extensively used as structural alloys in light water reactor (LWR) systems because of their excellent combination of mechanical properties and corrosion resistance. Austenitic stainless steel welds and cast austenitic stainless steels (CASS) contain significant amounts of ferrite ~5-30% and are ubiquitous in LWR piping (elbows, pump casings, valves) and internal structural components. During their lifetime, these components are subjected to thermal aging at the temperature range of 288-327°C for pressurized water reactors (PWR). A phase transformation within the ferrite phase of these materials leads to embrittlement and degradation in corrosion resistance, which undermines the materials selection criteria.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2022

Mineral Protection and Resource Limitations Combine to Explain Profile-Scale Soil Carbon Persistence

The fate of soil carbon (C) is largely controlled by microbial oxidation of organic matter (OM), which is constrained by a variety of mechanisms. OM association with soil minerals provides pronounced protection against microbial decomposition. However, factors such as climate, occlusion, and resource limitations also contribute to OM preservation. Here, we explore the factors explaining C distribution and age within an upland rainforest soil in Hawai'i, a site with abundant preferential flow paths (PFPs) and high short-range order (SRO) mineral content. We characterized lateral and vertical changes in Δ 14 C, SRO mineral content, C-functional group chemistry, and microbial community composition to elucidate the contributions of multiple protection mechanisms to OM preservation. Consistent with our expectation, SRO mineral content and Δ 14 C were strongly correlated (R 2 = 0.95), indicating strong mineral protection of OM throughout the profile. Surprisingly, distance from PFP was also a significant predictor of Δ 14 C and improved model fit, particularly in the shallow horizons (R 2 = 0.97). Elevated C/N ratios, decreased microbial abundance, and greater SRO mineral content suggest nitrogen limitations and enhanced mineral protection constrain OM turnover with distance from PFPs in deep, subsurface mineral horizons. Steady microbial abundance, increasing putative anaerobe abundance, and changes in C-functional group chemistry indicate oxygen limitations constrain OM turnover in the matrix of shallow mineral horizons. Given that oxygen and nutrient limitations contribute to OM preservation in this high SRO system—an exemplar of mineral protection—resource limitations may play an even more important role in OM preservation in other well-structured soils.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2022

Soil-related developments of the Biome-BGCMuSo v6.2 terrestrial ecosystem model

Terrestrial biogeochemical models are essential tools to quantify climate–carbon cycle feedback and plant–soil relations from local to global scale. In this study, a theoretical basis is provided for the latest version of the Biome-BGCMuSo biogeochemical model (version 6.2). Biome-BGCMuSo is a branch of the original Biome-BGC model with a large number of developments and structural changes. Earlier model versions performed poorly in terms of soil water content (SWC) dynamics in different environments. Moreover, lack of detailed nitrogen cycle representation was a major limitation of the model. Since problems associated with these internal drivers might influence the final results and parameter estimation, additional structural improvements were necessary. In this paper the improved soil hydrology as well as the soil carbon and nitrogen cycle calculation methods are described in detail. Capabilities of the Biome-BGCMuSo v6.2 model are demonstrated via case studies focusing on soil hydrology, soil nitrogen cycle, and soil organic carbon content estimation. Soil-hydrology-related results are compared to observation data from an experimental lysimeter station. The results indicate improved performance for Biome-BGCMuSo v6.2 compared to v4.0 (explained variance increased from 0.121 to 0.8 for SWC and from 0.084 to 0.46 for soil evaporation; bias changed from -0.047 to -0.007 m 3 m -3 for SWC and from -0.68 to -0.2 mm d -1 for soil evaporation). Simulations related to nitrogen balance and soil CO 2 efflux were evaluated based on observations made in a long-term field experiment under crop rotation. The results indicated that the model is able to provide realistic nitrate content estimation for the topsoil. Soil nitrous oxide (N 2 O) efflux and soil respiration simulations were also realistic, with overall correspondence with the observations (for the N 2 O efflux simulation bias was between -0.13 and -0.1 mgNm -2 d -1 , and normalized root mean squared error (NRMSE) was 32.4 %–37.6 %; for CO 2 efflux simulations bias was 0.04–0.17 gCm -2 d -1 , while NRMSE was 34.1 %–40.1 %). Sensitivity analysis and optimization of the decomposition scheme are presented to support practical application of the model. The improved version of Biome-BGCMuSo has the ability to provide more realistic soil hydrology representation as well as nitrification and denitrification process estimation, which represents a major milestone.

58 GEOSCIENCES↗

DOE OSTI2021

Efficient Space–Time Reduced Order Model for Linear Dynamical Systems in Python Using Less than 120 Lines of Code

A classical reduced order model (ROM) for dynamical problems typically involves only the spatial reduction of a given problem. Recently, a novel space–time ROM for linear dynamical problems has been developed [Choi et al., Space–tume reduced order model for large-scale linear dynamical systems with application to Boltzmann transport problems, Journal of Computational Physics, 2020], which further reduces the problem size by introducing a temporal reduction in addition to a spatial reduction without much loss in accuracy. The authors show an order of a thousand speed-up with a relative error of less than 10−5 for a large-scale Boltzmann transport problem. In this work, we present for the first time the derivation of the space–time least-squares Petrov–Galerkin (LSPG) projection for linear dynamical systems and its corresponding block structures. Utilizing these block structures, we demonstrate the ease of construction of the space–time ROM method with two model problems: 2D diffusion and 2D convection diffusion, with and without a linear source term. For each problem, we demonstrate the entire process of generating the full order model (FOM) data, constructing the space–time ROM, and predicting the reduced-order solutions, all in less than 120 lines of Python code. We compare our LSPG method with the traditional Galerkin method and show that the space–time ROMs can achieve O(10−3) to O(10−4) relative errors for these problems. Depending on parameter–separability, online speed-ups may or may not be achieved. For the FOMs with parameter–separability, the space–time ROMs can achieve O(10) online speed-ups. Finally, we present an error analysis for the space–time LSPG projection and derive an error bound, which shows an improvement compared to traditional spatial Galerkin ROM methods.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

Promise of Graph Sparsification and Decomposition for Noise Reduction in QAOA: Analysis for Trapped-Ion Compilations

We develop new approximate compilation schemes that significantly reduce the expense of compiling the Quantum Approximate Optimization Algorithm (QAOA) for solving the Max-Cut problem. Our main focus is on compilation with trapped-ion simulators using Pauli-X operations and all-to-all Ising Hamiltonian HIsing evolution generated by Molmer-Sorensen or optical dipole force interactions, though some of our results also apply to standard gate-based compilations. Our results are based on principles of graph sparsification and decomposition; the former reduces the number of edges in a graph while maintaining its cut structure, while the latter breaks a weighted graph into a small number of unweighted graphs. Though these techniques have been used as heuristics in various hybrid quantum algorithms, there have been no guarantees on their performance, to the best of our knowledge. This work provides the first provable guarantees using sparsification and decomposition to improve quantum noise resilience and reduce quantum circuit complexity. For quantum hardware that uses edge-by-edge QAOA compilations, sparsification leads to a direct reduction in circuit complexity. For trapped-ion quantum simulators implementing all-to-all HIsing pulses, we show that for a (1−ϵ) factor loss in the Max-Cut approximation (ϵ>0), our compilations improve the (worst-case) number of HIsing pulses from O(n2) to O(nlog(n/ϵ)) and the (worst-case) number of Pauli-X bit flips from O(n2) to O(nlog(n/ϵ)ϵ2) for n-node graphs. This is an asymptotic improvement for any constant ϵ>0. We demonstrate that significant improvements to the approximation ratio are obtained using decomposition in simulated trapped-ion experiments with dephasing noise. We further present a generic argument showing that sparsification results in an exponentially improved circuit fidelity lower bound in digital computing schemes based on one- and two-qubit gates, which are relevant to a wide variety of hardwares such as superconducting qubits and certain neutral atom or trapped ion setups, and more sophisticated noise models. We anticipate these approximate compilation techniques will be useful tools in a variety of future quantum computing experiments.

Moondra, Jai [Georgia Institute of Technology]↗

DOE OSTI2021

A fast multi-fidelity method with uncertainty quantification for complex data correlations: Application to vortex-induced vibrations of marine risers

Here we develop a fast multi-fidelity modeling method for very complex correlations between high- and low-fidelity data by working in modal space to extract the proper correlation function. We apply this method to infer the amplitude of motion of a flexible marine riser in cross-flow, subject to vortex-induced vibrations (VIV). VIV are driven by an absolute instability in the flow, which imposes a frequency (Strouhal) law that requires a matching with the impedance of the structure; this matching is easily achieved because of the rapid parametric variation of the added mass force. As a result, the wavenumber of the riser spatial response is within narrow bands of uncertainty. Hence, an error in wavenumber prediction can cause significant phase-related errors in the shape of the amplitude of response along the riser, rendering correlation between low- and high-fidelity data very complex. Working in modal space as outlined herein, dense data from low-fidelity data, provided by the semi-empirical computer code VIVA, can correlate in modal space with few high-fidelity data, obtained from experiments or fully-resolved CFD simulations, to correct both phase and amplitude and provide predictions that agree very well overall with the correct shape of the amplitude response. We also quantify the uncertainty in the prediction using Bayesian modeling and exploit this uncertainty to formulate an active learning strategy for the best possible location of the sensors providing the high fidelity measurements.

42 ENGINEERING↗

DOE OSTI2024

Galaxy bispectrum in the spherical Fourier-Bessel basis

The bispectrum, the three-point correlation in Fourier space, is a crucial statistic for studying many effects targeted by the next-generation galaxy surveys, such as primordial non-Gaussianity (PNG) and general relativistic (GR) effects on large scales. In this work we develop a formalism for the bispectrum in the spherical Fourier-Bessel (SFB) basis—a natural basis for computing correlation functions on the curved sky, as it diagonalizes the Laplacian operator in spherical coordinates. Working in the SFB basis allows for line-of-sight effects such as redshift space distortions and GR to be accounted for exactly, i.e., without having to resort to perturbative expansions to go beyond the plane-parallel approximation. Only analytic results for the SFB bispectrum exist in the literature given the intensive computations needed. We numerically calculate the SFB bispectrum for the first time, enabled by a few techniques: We implement a template decomposition of the redshift-space kernel Z 2 into Legendre polynomials, and separately treat the PNG and velocity-divergence terms. We derive an identity to integrate a product of three spherical harmonics connected by a Dirac delta function as a simple sum and use it to investigate the limit of a homogeneous and isotropic Universe. Furthermore, we present a formalism for convolving the signal with separable window functions and use a toy spherically symmetric window to demonstrate the computation and give insights into the properties of the observed bispectrum signal. While our implementation remains computationally challenging, it is a step toward a feasible full extraction of information on large scales via a SFB bispectrum analysis.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Multiscale Flow for robust and optimal cosmological analysis

We propose Multiscale Flow, a generative Normalizing Flow that creates samples and models the field-level likelihood of two-dimensional cosmological data such as weak lensing. Multiscale Flow uses hierarchical decomposition of cosmological fields via a wavelet basis and then models different wavelet components separately as Normalizing Flows. The log-likelihood of the original cosmological field can be recovered by summing over the log-likelihood of each wavelet term. This decomposition allows us to separate the information from different scales and identify distribution shifts in the data such as unknown scale-dependent systematics. The resulting likelihood analysis can not only identify these types of systematics, but can also be made optimal, in the sense that the Multiscale Flow can learn the full likelihood at the field without any dimensionality reduction. We apply Multiscale Flow to weak lensing mock datasets for cosmological inference and show that it significantly outperforms traditional summary statistics such as power spectrum and peak counts, as well as machine learning–based summary statistics such as scattering transform and convolutional neural networks. We further show that Multiscale Flow is able to identify distribution shifts not in the training data such as baryonic effects. Finally, we demonstrate that Multiscale Flow can be used to generate realistic samples of weak lensing data.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

PRISMS-PF: A general framework for phase-field modeling with a matrix-free finite element method

Abstract A new phase-field modeling framework with an emphasis on performance, flexibility, and ease of use is presented. Foremost among the strategies employed to fulfill these objectives are the use of a matrix-free finite element method and a modular, application-centric code structure. This approach is implemented in the new open-source PRISMS-PF framework. Its performance is enabled by the combination of a matrix-free variant of the finite element method with adaptive mesh refinement, explicit time integration, and multilevel parallelism. Benchmark testing with a particle growth problem shows PRISMS-PF with adaptive mesh refinement and higher-order elements to be up to 12 times faster than a finite difference code employing a second-order-accurate spatial discretization and first-order-accurate explicit time integration. Furthermore, for a two-dimensional solidification benchmark problem, the performance of PRISMS-PF meets or exceeds that of phase-field frameworks that focus on implicit/semi-implicit time stepping, even though the benchmark problem’s small computational size reduces the scalability advantage of explicit time-integration schemes. PRISMS-PF supports an arbitrary number of coupled governing equations. The code structure simplifies the modification of these governing equations by separating their definition from the implementation of the numerical methods used to solve them. As part of its modular design, the framework includes functionality for nucleation and polycrystalline systems available in any application to further broaden the phenomena that can be used to study. The versatility of this approach is demonstrated with examples from several common types of phase-field simulations, including coarsening subsequent to spinodal decomposition, solidification, precipitation, grain growth, and corrosion.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Microscale spatiotemporal patterns of water, soil organic carbon, and enzymes in plant litter detritusphere

Understanding biophysical and biochemical processes in detritusphere is critical to quantifying and modeling plant residue decomposition dynamics and subsequent soil organic carbon (C) accrual. The objectives of the study were to explore (i) micro-environmental conditions within the detritusphere formed around decomposing corn and soybean leaves, and (ii) the relationships between moisture distribution, enzyme activity, and C dynamics during decomposition. We assessed spatial and temporal dynamics of moisture distribution using X-ray and neutron computed tomography, and activities of β-glucosidase and chitinase, two enzymes involved in soil C and N processing, using zymography. We used 13 C labeled residue to track residue contribution to atmospheric CO 2 and soil organic C. Moisture redistribution pattern varied depending on the residue type. While the water was immediately absorbed by the corn leaves and maintained afterward, switchgrass leaves absorbed water more slowly and created water-deficient zones within ~1mm from the residue. This initial moisture depletion led to lower chitinase activity and residue-derived CO 2 emissions in switchgrass. In contrast to chitinase, β-glucosidase activity was influenced by a combination of vegetation history and residue type, and it was higher when the origin of the residue matched the vegetation history of the soil. Pore size had an opposite impact on the studied enzymes, supporting the notion that contrasting soil pore architecture can stimulate activities of different enzymes through a selection of dominant enzyme producers. We concluded that the decomposition dynamics of plant residues is not only a simple function of residue chemistry, but rather a combined effect of the vegetation history, in part through its effect on microbial community composition, the plant residue chemical and likely physical characteristics, and the soil pore structure in the detritusphere. Together, they create temporally dynamic micro-environmental conditions influencing decomposition. Specifically, our study demonstrated that the initial micro-environment formulated in detritusphere can play an important role in enzyme activities and consequent C dynamics.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2020

Towards a self-consistent analysis of the anisotropic galaxy two- and three-point correlation functions on large scales: application to mock galaxy catalogues

We establish a practical method for the joint analysis of anisotropic galaxy two- and three-point correlation functions (2PCF and 3PCF, respectively) on the basis of the decomposition formalism of the 3PCF using tripolar spherical harmonics. We perform such an analysis with MultiDark-Patchy mock catalogues to demonstrate and understand the benefit of the anisotropic 3PCF. We focus on scales above $80\, h^{-1}\, {\rm Mpc}$, and use information from the shape and the baryon acoustic oscillation (BAO) signals of the 2PCF and 3PCF. We also apply density field reconstruction to increase the signal-to-noise ratio of BAO in the 2PCF measurement, but not in the 3PCF measurement. In particular, we study in detail the constraints on the angular diameter distance and the Hubble parameter. We build a model of the bispectrum or 3PCF that includes the non-linear damping of the BAO signal in redshift space. We carefully account for various uncertainties in our analysis including theoretical models of the 3PCF, window function corrections, biases in estimated parameters from the fiducial values, the number of mock realizations to estimate the covariance matrix, and bin size. The joint analysis of the 2PCF and 3PCF monopole and quadrupole components shows a $30{{\ \rm per\ cent}}$ and $20{{\ \rm per\ cent}}$ improvement in Hubble parameter constraints before and after reconstruction of the 2PCF measurements, respectively, compared to the 2PCF analysis alone. Furthermore, this study clearly shows that the anisotropic 3PCF increases cosmological information from galaxy surveys and encourages further development of the modelling of the 3PCF on smaller scales than we consider.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

Primordial black hole dark matter: A quantitative parameter sensitivity comparison across formation mechanisms and particle candidates

Primordial black holes (PBHs) in the asteroid-mass window ( 10 17 – 10 22 g ) can account for all of the dark matter without violating any observational constraint, yet are routinely dismissed as fine-tuned. I put that dismissal to the test by applying three complementary sensitivity measures uniformly across a broad landscape: three noninflationary PBH production mechanisms, six classes of inflationary PBH models, and seven particle dark matter benchmarks, all evaluated against the same observable target. Three distinct naturalness universality classes emerge, determined entirely by the analytic structure of the abundance map rather than by the nature of the dark matter candidate. Biased-domain-wall PBHs, in their least model-dependent (free- V b ) form, have the same low sensitivity, Δ = 4.5 , as off-resonance weakly interacting massive particles and freeze-in particles ( Δ = 2 ), a sensitivity that, because it is constant over the entire parameter space of the construction, also coincides trivially with its own Wilson-normalized average within that parameter space (Section Definition and conventions), an equivalence that concerns only the space over which Δ is computed and is not a naturalness statement about the construction as a whole; a further reduction to Δ = 2 is possible only under the additional, independently motivated but not required, assumption that the domain-wall bias is generated by Planck-suppressed operators; early matter-domination PBHs occupy an intermediate tier alongside coannihilating weakly interacting massive particles (WIMPs), unified by a structural identity in which the sensitivity measure equals the logarithm of the ratio of the formation scale to the matter–radiation equality scale; first-order phase transition PBHs, once the more accurate super-exponential collapse probability is used in place of the single-exponential approximation, instead belong to the same highly sensitive tier as resonant WIMP annihilation and single-field inflationary collapse, for a structurally distinct reason; single-field ultraslow-roll inflationary collapse is severely tuned for a distinct reason: a double exponential in which the power spectrum amplitude is itself exponentially sensitive to the inflaton potential coefficients, on top of the exponential collapse sensitivity of the abundance map. My main conclusion is that the claim that PBH dark matter is generically fine-tuned conflates the worst case with a landscape spanning every naturalness tier. The Barbieri-Giudice sensitivity computed here and the Wilson-normalized measure of Iovino and Riotto answer distinct and mutually consistent questions about the same construction, a distinction I clarify within the two-layer decomposition.

Profumo, Stefano [University of California, Santa ↗