Search NASA⌕ Search

SEARCH · Search NASA

Results for “Divided differences”

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

Near-term heatwave risk in HighResMIP models across different temperature zones of West Africa

This study projects near-future (2031–2050) changes in heatwave (HW) risk across West Africa (WA) using an ensemble of eight high-resolution global climate models from the High-Resolution Model Intercomparison Project under a high-emission scenario. Using K-means clustering, we divided WA into four unique temperature zones and examined projected changes in extreme temperatures, HW occurrence and magnitude. Our results indicate a statistically significant increase in future HW events across most parts of WA, although considerable spread exists over the region and among individual models. The most pronounced increases are evident in the Sahel/Sahara and the Guinea Highlands subregions, with an ensemble mean increase of ∼10 HW events per year. In contrast, the lowest increase in HW events is projected in central WA, with increases ranging between 1 and 5 events per year. Similarly, the magnitude of HW events is projected to increase in most models, with Sahel/Sahara exhibiting the largest increases. Additionally, projections suggest that the strongest HWs will become more frequent, particularly in northern and southwestern WA. These findings highlight significant spatial heterogeneity in future HW risk across WA, emphasizing the need for targeted adaptation strategies.

HighResMIP↗
DOE OSTI2020

Dynamic load balancing with enhanced shared-memory parallelism for particle-in-cell codes

Furthering our understanding of many of today’s interesting problems in plasma physics – including plasma based acceleration and magnetic reconnection with pair production due to quantum electrodynamic effects – requires large-scale kinetic simulations using particle-in-cell (PIC) codes. However, these simulations are extremely demanding, requiring that contemporary PIC codes be designed to efficiently use a new fleet of exascale computing architectures. To this end, the key issue of parallel load balance across computational nodes must be addressed. We discuss the implementation of dynamic load balancing by dividing the simulation space into many small, self-contained regions or ‘‘tiles,’’ along with shared-memory (e.g., OpenMP) parallelism both over many tiles and within single tiles. The load balancing algorithm can be used with three different topologies, including two space-filling curves. Here, we tested this implementation in the code Osiris and show low overhead and improved scalability with OpenMP thread number on simulations with both uniform load and severe load imbalance. Compared to other load-balancing techniques, our algorithm gives order-of-magnitude improvement in parallel scalability for simulations with severe load imbalance issues.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Nuclear speed and cycle length co-vary with local density during syncytial blastoderm formation in a cricket

Abstract The blastoderm is a broadly conserved stage of early animal development, wherein cells form a layer at the embryo’s periphery. The cellular behaviors underlying blastoderm formation are varied and poorly understood. In most insects, the pre-blastoderm embryo is a syncytium: nuclei divide and move throughout the shared cytoplasm, ultimately reaching the cortex. In Drosophila melanogaster , some early nuclear movements result from pulsed cytoplasmic flows that are coupled to synchronous divisions. Here, we show that the cricket Gryllus bimaculatus has a different solution to the problem of creating a blastoderm. We quantified nuclear dynamics during blastoderm formation in G. bimaculatus embryos, finding that: (1) cytoplasmic flows are unimportant for nuclear movement, and (2) division cycles, nuclear speeds, and the directions of nuclear movement are not synchronized, instead being heterogeneous in space and time. Moreover, nuclear divisions and movements co-vary with local nuclear density. We show that several previously proposed models for nuclear movements in D. melanogaster cannot explain the dynamics of G. bimaculatus nuclei. We introduce a geometric model based on asymmetric pulling forces on nuclei, which recapitulates the patterns of nuclear speeds and orientations of both unperturbed G. bimaculatus embryos, and of embryos physically manipulated to have atypical nuclear densities.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

MPACT 4.4 Theory Manual

MPACT is a three-dimensional (3D) full-core neutron transport code capable of calculating subpin power distributions. Calculations are based on the Boltzmann transport equation for neutron fluxes for problems in which the detailed geometrical configuration of fuel components such as the pellet and cladding are explicitly retained. The cross-section data needed for the neutron transport calculation are obtained directly from a multigroup cross section library, which has traditionally been used by lattice physics codes to generate few-group homogenized cross sections for nodal core simulators. Hence, MPACT assumes neither a priori homogenization nor group condensation for the full core spatial solution. The 3D MPACT transport solution can be obtained using the method of characteristics (MOC), which employs discrete ray tracing within each fuel pin. However, for practical reactor applications, the direct application of MOC to 3D core configurations requires an excessive amount of memory and computing time due to the very large number of rays. For practical 3D full-core calculations, MPACT commonly uses an approximate “2D/1D” method that treats the radial (x and y) variables differently from the axial (z) variable. In particular, the radial dependence of the solution is calculated using transport theory, and the axial dependence is calculated using diffusion or P 3 theory. The 2D/1D method requires the core to be divided into a vertical stack of axial slices with a thickness of Δ z ≈ 5–10 cm. Each axial slice is divided radially into coarse spatial cells with boundaries that usually constitute the pin cell boundaries, for which Δ x = Δ y ≈ 1.5 cm. Then, each coarse radial cell (pin cell) is divided into 50–100 fine radial cells, which resolve the angular flux in the fuel, cladding, and moderator regions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Mechanical Loads on Spent Nuclear Fuel in the General 30 cm Package Drop Scenario

The U.S. Department of Energy Office of Nuclear Energy (DOE-NE) Spent Fuel and Waste Science & Technology (SFWST) research program is guided by the high-level goal of closing prioritized knowledge gaps related to spent nuclear fuel (SNF) storage and transportation, which are summarized by Saltzstein et al. (2020). One of the high-priority knowledge gaps is the identification and quantification of mechanical loads that are expected to affect SNF during normal conditions of transportation and storage to inform the range of physical SNF test programs. This report uses modeling and analysis methods to estimate the mechanical loads on spent nuclear fuel (SNF) in the general 30 cm package drop scenario. The drop scenario assumes impact limiters are in place in the transportation configuration and the impact surface is perfectly rigid. The goal of this analysis is to consider the universe of potential mechanical loading conditions that can happen to SNF and present the results in a manner that is useful for materials testing, decision making, and regulatory rule making purposes. This study uses validated finite element models and methods to perform a broad parametric study of key variables that can affect the mechanical loads on SNF during a hypothetical package free drop scenario. Physical drop test data from a cask and fuel assembly drop test campaign is the basis for model validation. Additionally, the results of the parametric study are used to inform a damage model, which uses multiple nonlinear regression to estimate the relationships between input variables and output response. The parametric finite element analyses consider thousands of input variable combinations, while the damage model estimates millions of combinations. The breadth of this study provides confidence that the potential range of mechanical loads that SNF might experience during the general 30 cm package drop scenario are characterized well enough to consider this knowledge gap closed. While this report documents the overall peak values calculated in this study, the 95 th percentile values, the histograms, and the observed trends are equally important. This study covered a large range of SNF temperatures, room temperature to 300°C, and burnups, 10 GWd/MTU to 62 GWd/MTU. Each temperature and burnup combination has a different cladding yield strain, so it is more meaningful to summarize the calculated cladding strain response as its factor of safety, which is defined relative to the yield strain. The factor of safety is calculated as the yield strain divided by the peak cladding strain. A factor of safety greater than unity indicates that the cladding remains below yield, whereas a value less than unity is indicative of plastic deformation. In all cases of this study a safety factor over 1.0 was calculated, although in the most limiting case at 300°C the safety factor was only 1.01, which suggests that yielding could occur when additional loads like rod internal pressure are included. When the temperature is restricted to 200°C the limiting safety factor increases to 1.28, which has significant margin to accommodate internal pressure and potential local cladding defects that could cause a local stress concentration. An important trend in the calculated fuel rod mechanical loads is that the 2 nd highest loaded fuel rod in an assembly tends to be significantly lower than the highest loaded rod. The implication is that even if one rod in an assembly experiences a failure the loads would have to be significantly higher to cause two or more rods to fail.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

An Empirical Quantile Estimation Approach for Chance-Constrained Nonlinear Optimization Problems

We investigate an empirical quantile estimation approach to solve chance-constrained nonlinear optimization problems. Our approach is based on the reformulation of the chance constraint as an equivalent quantile constraint to provide stronger signals on the gradient. In this approach, the value of the quantile function is estimated empirically from samples drawn from the random parameters, and the gradient of the quantile function is estimated via a finite-difference approximation on top of the quantile-function-value estimation. We establish a convergence theory of this approach within the framework of an augmented Lagrangian method for solving general nonlinear constrained optimization problems. The foundation of the convergence analysis is a concentration property of the empirical quantile process, and the analysis is divided based on whether or not the quantile function is differentiable. In contrast to the sampling-and-smoothing approach used in the literature, the method developed in this paper does not involve any smoothing function and hence the quantile-function gradient approximation is easier to implement and there are less accuracy-control parameters to tune. Furthermore, we demonstrate the effectiveness of this approach and compare it with a smoothing method for the quantile-gradient estimation. Numerical investigation shows that the two approaches are competitive for certain problem instances.

Applied Probability↗
DOE OSTI2024

SSTDR and FDR Detection of Un-Energized and Energized Cable Anomalies Including Thermal Degradation Using Machine Learning

Historically, cables are initially qualified for nuclear power plant use for 40 years. As plants extend their operating license to 60 and 80 years, continued use of these cables must shift to a performance-based approach since it is cost prohibitive to completely replace cables that are likely still capable of performing their design function. A variety of cable tests are available and are commonly applied during outages when the cables can be taken out of service. Frequency domain reflectometry (FDR) is one of these test methods that is being more broadly accepted and used because it not only detects anomalies along the cable with a low-voltage signal that does not stress the cable insulation, but the technique also locates the anomalies. This supports follow-up local inspection and local repair or partial replacement of a damaged cable segment. Currently, FDR testing is only applied to cables that are taken out of service since the test instrument would be damaged by operational voltages. A related technology that has found some acceptance in the aircraft and rail industry is spread spectrum time domain reflectometry (SSTDR). This technology has been implemented with a custom commercial instrument by LiveWire Innovation that is designed to operate on live cables up to 1000 volts and with a bandwidth of 48 MHz. Initial evaluation by the Pacific Northwest National Laboratory (PNNL) of the Live Wire system indicated that a broader bandwidth (BW) SSTDR may be better for many kinds of flaws. This led PNNL to develop an SSTDR laboratory instrument suitable for tests up to 500 MHz bandwidth. Testing on energized cables is also desirable for online monitoring systems so an inductive clamshell coupler was developed that allows energized cables to be tested up to at least 5 kV and likely higher voltage levels. Dielectric spectroscopy and tan delta testing plus various laboratory destructive tests were included in this data acquisition campaign directed to feed a machine learning (ML) study. With these kinds of developments, online energized cable tests may be possible with industrial adoption of such hardware advances but it will be completely impractical to have highly skilled data analysts continually examine these complex signals for indications of damage or compromised conditions. If online testing is to be implemented in new test hardware, it must be accompanied by software that can interpret the signals and alert plant operators of changing or degraded conditions. The thermally aged, shielded cable investigated here was separately treated for ML analysis. Visual analysis of electrical data showed generally increasing peaks where the cable entered and exited the oven. These peaks were not exactly aligned with expected locations, but these differences were attributed to velocity of propagation calibration errors. Only supervised ML was applied to the thermally aged data as this data was only available shortly before the committed publication date of this report. The supervised ML was structured to divide the 0 to 70-day responses as ‘normal’ from 0 to 35 days or ‘anomalous’ from 36 to 70 days, based on cable tensile elongation at break (EAB) insulation characterization. Using 80% of the data for training and 20% for testing, the supervised ML predicted normal versus anomalous was 70% accurate. Important conclusions include: • Accuracy to predict the presence of cable damage is improved from the 2023 effort by more training data. Weighted accuracies for comparisons among the instruments ranged from 67 to 89 % for unsupervised ML and 71 to 99% for supervised ML. • Based on the synthetic data tests, the unsupervised models are more generalizable to unseen anomalies. The Multi-Layer Perceptron classifier (MLP) model reported as high as 99.7% accuracy on the test data, but this dropped to 58.3% when tested on the synthetic data. In contrast, the unsupervised Pointwise model only achieved 89.7% accuracy on the experimental data but reported 78.3% accuracy on the synthetic data. • The best anomaly indicators are higher frequency (400 MHz BW) FDR data. Other tests may be interesting but for this study, this was the best predicter.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

A global model–measurement evaluation of particle light scattering coefficients at elevated relative humidity

The uptake of water by atmospheric aerosols has a pronounced effect on particle light scattering properties, which in turn are strongly dependent on the ambient relative humidity (RH). Earth system models need to account for the aerosol water uptake and its influence on light scattering in order to properly capture the overall radiative effects of aerosols. Here we present a comprehensive model–measurement evaluation of the particle light scattering enhancement factor f(RH), defined as the particle light scattering coefficient at elevated RH (here set to 85 %) divided by its dry value. The comparison uses simulations from 10 Earth system models and a global dataset of surface-based in situ measurements. In general, we find a large diversity in the magnitude of predicted f(RH) amongst the different models, which can not be explained by the site types. Based on our evaluation of sea salt scattering enhancement and simulated organic mass fraction, there is a strong indication that differences in the model parameterizations of hygroscopicity and model chemistry are driving at least some of the observed diversity in simulated f(RH). Additionally, a key point is that defining dry conditions is difficult from an observational point of view and, depending on the aerosol, may influence the measured f(RH). The definition of dry also impacts our model evaluation, because several models exhibit significant water uptake between RH = 0 % and 40 %. The multisite average ratio between model outputs and measurements is 1.64 when RH = 0 % is assumed as the model dry RH and 1.16 when RH = 40 % is the model dry RH value. The overestimation by the models is believed to originate from the hygroscopicity parameterizations at the lower RH range which may not implement all phenomena taking place (i.e., not fully dried particles and hysteresis effects). This will be particularly relevant when a location is dominated by a deliquescent aerosol such as sea salt. Our results emphasize the need to consider the measurement conditions in such comparisons and recognize that measurements referred to as dry may not be dry in model terms. Recommendations for future model–measurement evaluation and model improvements are provided.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Integrands of less-supersymmetric Yang-Mills at one loop

We construct a prescriptive, bubble power-counting basis of one-loop integrands suitable for representing amplitude integrands in less-supersymmetric (1 ≤ $\mathcal{N}$ ≤ 4) Yang-Mills theory. With the exception of massless bubbles, all integrands have unambiguous, leading singularities as coefficients defined in field theory; for the massless bubbles on external legs, we find two natural choices which lead to different integrands that highlight distinct aspects of field theory. For concreteness, we give the all-multiplicity integrands for MHV amplitudes, and the split-helicity amplitude integrand for six-particle NMHV. The basis we construct is mostly pure and is divided into to separately UV- and IR-finite sectors of fixed transcendental weight, resulting in UV- and IR-finite ratio functions of n-particle helicity amplitudes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Effect of second Si–O vibrational overtones/combinations on quantifying water in silicate and silica minerals using infrared spectroscopy, and an experimental method for its removal

Infrared spectroscopy (IR) is the most widely used analytical tool to quantify trace water in silicate and silica minerals. A prerequisite for highly accurate IR measurements of trace water is a good understanding of the effect of the second Si–O vibrational overtones/combination bands (2nd Si–O VOCBs) on the water peaks. Silicate and silica minerals can be divided as isolated (Q 0 ), paired (Q 1 ), ring (Q 2 ), chain (Q 1 or Q 2 ), sheet (Q 3 ) and framework (Q 4 ) structures according to the polymerization of their SiO 4 tetrahedral units, and the 2nd Si–O VOCBs of these different structural types attain different vibrational features which are expected to affect the water peaks to different extents. Here, we selected olivine (Q 0 ) and α-quartz (Q 4 ) as two endmember-like structural examples, performed extensive IR measurements on both pristine and heat-treated thin sections prepared for these two minerals, and explored the vibrational features of the 2nd Si–O VOCBs. In this work, we have found that the 2nd Si–O VOCBs are well separated from the water peaks in olivine, but severely overlap with the water peaks in α-quartz, confirming the different roles that the 2nd Si–O VOCBs play in quantifying trace water in silicate and silica minerals with different structural polymerizations. To remove the influence of the 2nd Si–O VOCBs (or any other species rather than water), an experimental protocol has been successfully developed, as approved by some fundamental equations and verified by the data of α-quartz in the literature. This development should lead to significant accuracy improvement in quantifying trace amounts of water in Earth and planetary materials.

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

Explosion behavior investigation and safety assessment of large-format lithium-ion pouch cells

We report large-format lithium-ion (Li-ion) batteries with high energy density for electric vehicles are prone to thermal runaway (or even explosion) under abusive conditions. In this study, overcharge induced explosion behaviors of large-format Li-ion pouch cells with Li[Ni 0.8 Co 0.1 Mn 0.1 ]O 2 cathode at different current rates (C-rates) (0.5C, 1C, 2C) were investigated. The explosion characteristics of the cells were elucidated by discussing the evolution of the cell voltage, the surface temperature and the shock wave pressure. Generally, the whole overcharge process could be divided into four stages according to the evolution of several key parameters and the overcharge behaviors; the overcharge C-rate has a great influence on cells’ thermal behaviors. The experimental results showed that the thermal runaway process of Li-ion cells caused by overcharging consisted of two kinds of explosions, physical explosion and chemical explosion. The existence of observable negative pressure zone in the pressure curves indicated that the Li-ion cells are not a self-supplying oxygen system during the explosion. Further, the explosion dynamics parameters were matched. An explosion TNT-equivalent conversion strategy that depended on the pressure of the shock wave was utilized to evaluate the released energy and its hazards. In addition, with respect to the overcharge of Li-ion pouch cells, a safety assessment method and a safety management method were proposed based on the explosion behaviors. From the perspective of battery safety, this study is of great significance for the safety design of Li-ion cells and can provide guidance for engineers to optimize the safety function of battery packs.

25 ENERGY STORAGE↗
DOE OSTI2024

Quantum Simulation of the First-Quantized Pauli-Fierz Hamiltonian

We provide an explicit recursive divide-and-conquer approach for simulating quantum dynamics and derive a discrete first-quantized nonrelativistic QED Hamiltonian based on the many-particle Pauli-Fierz Hamiltonian. We apply this recursive divide-and-conquer algorithm to this Hamiltonian and compare it to a concrete simulation algorithm that uses qubitization. Our divide-and-conquer algorithm, using lowest-order Trotterization, scales for fixed grid spacing as O ~ ( Λ N 2 η 2 t 2 / ϵ ) for grid size N , η particles, simulation time t , field cutoff Λ , and error ϵ . Our qubitization algorithm scales as O ~ ( N ( η + N ) ( η + Λ 2 ) t log ( 1 / ϵ ) ) . This shows that even a naive partitioning and low-order splitting formula can yield, through our divide-and-conquer formalism, superior scaling to qubitization for large Λ . We compare the relative costs of these two algorithms on systems that are relevant for applications such as the spontaneous emission of photons and the photoionization of electrons. We observe that for different parameter regimes, one method can be favored over the other. Finally, we give new algorithmic and circuit-level techniques for gate optimization, including a new way of implementing a group of multicontrolled- X gates that can be used for better analysis of circuit cost. Published by the American Physical Society 2024

Mukhopadhyay, Priyanka (ORCID:0000000164639100)↗
DOE OSTI2026

BSEC Air Quality: Stationary Black and Brown Carbon Properties and Mass Concentration

Data is 1-min absorption/attenuation of aerosol collected on filter and sampled at seven wavelengths (950, 880, 660, 590, 520, 470, and 370 nm) with a Magee AE-33 dual spot aethalometer. The wavelength value for black carbon, as recommended by the manufacturer, is the 880 nm channel. The derived mass concentrations, using the factory default mass absorption cross-sections, are reported for all seven wavelengths. The absorption at all seven wavelengths are reported by dividing the derived mass concentrations by the factory default mass absorption cross-sections. The estimated mass concentration contribution from biomass burning and fossil fuel black carbon, using factory default values, are reported. The estimated brown carbon mass concentration is derived by finding the difference between the recommended black carbon wavelength (880 nm) and the ultra violet wavelength (370 nm), assuming that brown carbon absorption is strongest at 370 nm and minimal at 880 nm. Finally, the estimated absorption Angstrom exponent for black and brown carbon is derived. For black carbon, the slope between each minute log(AbsXXX), where XXX is the wavelength used versus log(XXX) is used, and the wavelengths are 590, 660, 880, and 950 nm. For brown carbon, the slope between each minute log(AbsXXX), where XXX is the wavelength used versus log(XXX) is used, and the wavelengths are 370, 470, and 520 nm.

aerosol↗
DOE OSTI2024

GPU-enabled extreme-scale turbulence simulations: Fourier pseudo-spectral algorithms at the exascale using OpenMP offloading

Fourier pseudo-spectral methods for nonlinear partial differential equations are of wide interest in many areas of advanced computational science, including direct numerical simulation of three-dimensional (3-D) turbulence governed by the Navier-Stokes equations in fluid dynamics. This paper presents a new capability for simulating turbulence at a new record resolution up to 35 trillion grid points, on the world's first exascale computer, Frontier, comprising AMD MI250x GPUs with HPE's Slingshot interconnect and operated by the US Department of Energy's Oak Ridge Leadership Computing Facility (OLCF). Key programming strategies designed to take maximum advantage of the machine architecture involve performing almost all computations on the GPU which has the same memory capacity as the CPU, performing all-to-all communication among sets of parallel processes directly on the GPU, and targeting GPUs efficiently using OpenMP offloading for intensive number-crunching including 1-D Fast Fourier Transforms (FFT) performed using AMD ROCm library calls. With 99% of computing power on Frontier being on the GPU, leaving the CPU idle leads to a net performance gain via avoiding the overhead of data movement between host and device except when needed for some I/O purposes. Memory footprint including the size of communication buffers for MPI_ALLTOALL is managed carefully to maximize the largest problem size possible for a given node count. Detailed performance data including separate contributions from different categories of operations to the elapsed wall time per step are reported for five grid resolutions, from 2048 3 on a single node to 32768 3 on 4096 or 8192 nodes out of 9408 on the system. Both 1D and 2D domain decompositions which divide a 3D periodic domain into slabs and pencils respectively are implemented. The present code suite (labeled by the acronym GESTS, GPUs for Extreme Scale Turbulence Simulations) achieves a figure of merit (in grid points per second) exceeding goals set in the Center for Accelerated Application Readiness (CAAR) program for Frontier. The performance attained is highly favorable in both weak scaling and strong scaling, with notable departures only for 2048 3 where communication is entirely intra-node, and for 32768 3 , where a challenge due to small message sizes does arise. Communication performance is addressed further using a lightweight test code that performs all-to-all communication in a manner matching the full turbulence simulation code. Performance at large problem sizes is affected by both small message size due to high node counts as well as dragonfly network topology features on the machine, but is consistent with official expectations of sustained performance on Frontier. Overall, although not perfect, the scalability achieved at the extreme problem size of 32768 3 (and up to 8192 nodes — which corresponds to hardware rated at just under 1 exaflop/sec of theoretical peak computational performance) is arguably better than the scalability observed using prior state-of-the-art algorithms on Frontier's predecessor machine (Summit) at OLCF. New science results for the study of intermittency in turbulence enabled by this code and its extensions are to be reported separately in the near future.

3D fast Fourier transform↗
DOE OSTI2022

Formation of S0s in extreme environments Ⅲ: the role of environment in the formation pathways

It is well established that there are at least two main channels to form lenticular (or S0) galaxies. The first, which we name ‘faded spiral’ scenario, includes quenching events that led to consumption or removal of gas from a spiral progenitor. The second, which we call ‘merger’ scenario, includes merger-like events and interactions between galaxies. Each scenario leaves characteristic signatures in the newly formed lenticular galaxy. However, the conditions that trigger one mechanism over another are still unknown. This paper is the third of a series aimed at understanding the role of the environment in the formation of lenticular galaxies. In this study, we combine the kinematics, morphology, and properties of the stellar populations of 329 S0s from the SAMI and MaNGA surveys in order to highlight the role of the environment in the process. We divide the S0s into two classes (A and B) according to their global properties, that we can associate to the products of a faded spiral scenario (class A) or a merger scenario (class B). Here we then study how the various classes are distributed within different environments. Our study reveals that the ‘faded spiral’ pathway is the most efficient channel to produce S0s, and it becomes more efficient as the mass of the group or cluster or local density of galaxies increase. The merger pathway is also a viable channel, and its efficiency becomes higher with decreasing local density or environment mass.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Halo Properties and Mass Functions of Groups/Clusters from the DESI Legacy Imaging Surveys DR9

Abstract Based on a large group/cluster catalog recently constructed from the DESI Legacy Imaging Surveys DR9 using an extended halo-based group finder, we measure and model the group–galaxy weak-lensing signals for groups/clusters in a few redshift bins within redshift range 0.1 ≤ z < 0.6. Here, the background shear signals are obtained based on the DECaLS survey shape catalog, derived with the F ourier _Q uad method. We divide the lens samples into five equispaced redshift bins and seven mass bins, which allow us to probe the redshift and mass dependence of the lensing signals, and hence the resulting halo properties. In addition to these sample selections, we also check the signals around different group centers, e.g., the brightest central galaxy, the luminosity-weighted center, and the number-weighted center. We use a lensing model that includes off-centering to describe the lensing signals that we measure for all mass and redshift bins. The results demonstrate that our model predictions for the halo masses, biases, and concentrations are stable and self-consistent among different samples for different group centers. Taking advantage of the very large and complete sample of groups/clusters, as well as the reliable estimations of their halo masses, we provide measurements of the cumulative halo mass functions up to redshift z = 0.6, with a mass precision at 0.03 ∼ 0.09 dex.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

The maximum extent of the filaments and sheets in the cosmic web: an analysis of the SDSS DR17

Filaments and sheets are striking visual patterns in cosmic web. The maximum extent of these large-scale structures are difficult to determine due to their structural variety and complexity. We construct a volume-limited sample of galaxies in a cubic region from the SDSS, divide it into smaller subcubes and shuffle them around. We quantify the average filamentarity and planarity in the 3D galaxy distribution as a function of the density threshold and compare them with those from the shuffled realizations of the original data. The analysis is repeated for different shuffling lengths by varying the size of the subcubes. The average filamentarity and planarity in the shuffled data show a significant reduction when the shuffling scales are smaller than the maximum size of the genuine filaments and sheets. We observe a statistically significant reduction in these statistical measures even at a shuffling scale of $\sim 130 \, {{\, \rm Mpc}}$, indicating that the filaments and sheets in three dimensions can extend up to this length scale. They may extend to somewhat larger length scales that are missed by our analysis due to the limited size of the SDSS data cube. We expect to determine these length scales by applying this method to deeper and larger surveys in future.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Using respiration quotients to track changing sources of soil respiration seasonally and with experimental warming

Developing a more mechanistic understanding of soil respiration is hampered by the difficulty in determining the contribution of different organic substrates to respiration and in disentangling autotrophic-versus-heterotrophic and aerobic-versus-anaerobic processes. Here, we use a relatively novel tool for better understanding soil respiration: the apparent respiration quotient (ARQ). The ARQ is the amount of CO 2 produced in the soil divided by the amount of O 2 consumed, and it changes according to which organic substrates are being consumed and whether oxygen is being used as an electron acceptor. We investigated how the ARQ of soil gas varied seasonally, by soil depth, and by in situ experimental warming ( +4 ° C) in a coniferous-forest whole-soil-profile warming experiment over 2 years. We then compared the patterns in ARQ to those of soil δ 13 CO 2 . Our measurements showed strong seasonal variations in ARQ, from ≈0.9 during the late spring and summer to ≈0.7 during the winter. This pattern likely reflected a shift from respiration being fueled by oxidized substrates like sugars and organic acids derived from root and root respiration during the growing season to more reduced substrates such as lipids and proteins derived from microbial necromass during the winter. This interpretation was supported by δ 13 CO 2 values, which were lower, like lipids, in the winter and higher, like sugars, in the summer. Furthermore, experimental warming significantly changed how both ARQ and δ 13 CO 2 responded to soil temperature. Wintertime ARQ and δ 13 CO 2 values were higher in heated than in control plots, probably due to the warming-driven increase in microbial activity that may have utilized oxidized carbon substrates, while growing-season values were lower in heated plots. Experimental warming and phenology change the sources of soil respiration throughout the soil profile. The sensitivity of ARQ to these changes demonstrates its potential as a tool for disentangling the biological sources contributing to soil respiration.

59 BASIC BIOLOGICAL SCIENCES↗