Search NASA⌕ Search

SEARCH · Search NASA

Results for “Math Model”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Co-design for Particle Applications at Exascale

Co-design across the Exascale Computing Project (ECP) has been critical for both enabling science applications and bringing disparate communities together. Developing and porting applications to the various high-performance computing (HPC) architectures on pre-exascale and exascale computers has been quite challenging due to the diversity of hardware features and software stacks. The Co-design Center for Particle Applications (CoPA) has developed and enhanced the Cabana and PROGRESS/BML libraries to facilitate the creation of new particle applications, make existing particle applications exascale capable, and allow teams to explore new capabilities. Particle methods from atomistic, mesoscale, continuum, through cosmological scales have been built with Cabana, along with new possibilities for application coupling. Similarly, the PROGRESS/BML library has enabled quantum particle applications with linear algebra solvers to use advanced hardware. Across these CoPA-developed libraries, the co-design abstraction layer combines performance portability with math library support to facilitate separation of concerns and directly support science runs.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Scalable DPG multigrid solver for Helmholtz problems: A study on convergence

This paper presents a scalable multigrid preconditioner targeting large-scale systems arising from discontinuous Petrov–Galerkin (DPG) discretizations of high-frequency wave operators. This work is built on previously developed multigrid preconditioning techniques of Petrides and Demkowicz (Comput. Math. Appl. 87 (2021) pp. 12–26) and extends the convergence results from $\mathscr{O}$(10 7 ) degrees of freedom (DOFs) to $\mathscr{O}$(10 9 ) DOFs using a new scalable parallel MPI/OpenMP implementation. Novel contributions of this paper include an alternative definition of coarse-grid systems based on restriction of fine-grid operators, yielding superior convergence results. In the uniform refinement setting, a detailed convergence study is provided, demonstrating h and p robust convergence and linear scaling with respect to the wave frequency. Finally, the paper concludes with numerical results on hp -adaptive simulations including a large-scale seismic modeling benchmark problem with high material contrast.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Lens Model Accuracy in the Expected LSST Lensed AGN Sample

Strong gravitational lensing of active galactic nuclei (AGN) enables measurements of cosmological parameters through time-delay cosmography (TDC). With data from the upcoming LSST survey, we anticipate using a sample of O(1000) lensed AGN for TDC. To prepare for this dataset and enable this measurement, we construct and analyze a realistic mock sample of 1300 systems drawn from the OM10 (Oguri & Marshall 2010) catalog of simulated lenses with AGN sources at $z<3.1$ in order to test a key aspect of the analysis pipeline, that of the lens modeling. We realize the lenses as power law elliptical mass distributions and simulate 5-year LSST i-band coadd images. From every image, we infer the lens mass model parameters using neural posterior estimation (NPE). Focusing on the key model parameters, $θ_E$ (the Einstein Radius) and $γ_{lens}$ (the projected mass density profile slope), with consistent mass-light ellipticity correlations in test and training data, we recover $θ_E$ with less than 1% bias per lens, 6.5% precision per lens and $γ_{lens}$ with less than 3% bias per lens, 8% precision per lens. We find that lens light subtraction prior to modeling is only useful when applied to data sampled from the training prior. If emulated deconvolution is applied to the data prior to modeling, precision improves across all parameters by a factor of 2. Finally, we combine the inferred lens mass models using Bayesian Hierarchical Inference to recover the global properties of the lens sample with less than 1% bias.

Venkatraman, Padmavathi [Illinois U., Urbana; KIPA↗
DOE OSTI2021

xSDK: Building an ecosystem of highly efficient math libraries for exascale

Current efforts to build increasingly powerful computer architectures are opening up new avenues for more complex and higher fidelity simulations coupled with data analytics and learning, leading to new scientific insights and deeper understanding. At one extreme, exascale computers will be much faster than previous computer generations (performing 10 18 operations per second—that is, 1,000 times faster than petascale). To achieve these performance improvements, computer architectures are becoming increasingly complex, with deep memory hierarchies, very high node and core counts, and heterogeneous features such as graphics processing units (GPUs). Such architectural changes impact the full breadth of computing scales, as heterogeneity pervades even current-generation laptops, workstations, and moderate-sized clusters. While emerging advanced architectures provide unprecedented opportunities, they also present significant challenges for developers of scientific applications, such as multiphysics and multiscale codes, who must adapt their software to handle disruptive changes in architectures and new programming models that have not yet stabilized. Developers must consider increasing concurrency while reducing communication and synchronization, and other complexities such as the potential for using mixed precision to leverage the compute power available in low-precision tensor cores. On one hand, developers must implement new scientific capabilities, which in turn increase code complexity. On the other hand, the codes must be ported to new architectures, requiring the inclusion of new programming models and the restructuring of code to achieve good performance. Addressing these issues is beyond the capability of any single person or team—leading to the need for collaboration among many teams, who encapsulate their expertise in reusable software and work together to create sustainable software ecosystems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Full calibration of the tomographic redshift distribution from the HSC PDR3 Shape Catalog with DESI

The calibration of tomographic redshift distributionsis essential for cosmological analysis of weak lensing data.In this work, we calibrate all four tomographic bins of the Hyper Suprime Camera (HSC) weak lensing catalog with the Dark Energy Spectroscopic Instrument (DESI) Data Release 1 and 2 using the clustering redshifts technique. We include z > 1.2 redshift sources such as emission line galaxies (ELG) and quasars (QSO) sources in our calibration, which were not available in the previous HSC calibration (Rau et al. (2022), Mon. Not. Roy. Astron. Soc. 524 (2023) 5109), allowing a complete calibration of all the redshift bins. We find the first tomographic bin exhibits a small shift towards low redshifts. The second bin is in good agreement with the photometric calibration, while third and fourth bin exhibit a shift towards higher redshifts. However, these shifts are considerably smaller than the shifts obtained in the HSC Year 3 cosmic shear analyses. We evaluate the impact of galaxy bias and magnification effects from all the samples on the measurements, finding them to be small, and we propose corrections to reduce them further. Specifically, we relax the assumption of linear bias and only assume no redshift evolution of the cross-correlation coefficient, allowing us to leverage smaller clustering scales. We model the redshift distributions with splines and compare our results to previous analyses as well as to other parameterizations found in literature. For the two high-redshift tomographic bins, we find the shifts to higher redshifts with respect to the measurements performed in Rau+2022 to be Δz$_{3}$ =-0.039$^{+0.020}$$_{-0.021}$ and Δz$_{4}$ = -0.048$^{+0.012}$$_{-0.012}$.

Choppin de Janvry, J. [LBL, Berkeley; UC, Berkeley↗
DOE OSTI2026

A New Window into the Baryon Cycle at Cosmic Noon with Line Intensity Mapping: Forecasts for auto- and cross-correlations in [CII]-158$μ$m, HI 21 cm, CO$_{J+1\rightarrow J}$, and H$α$ galaxies

Across the peak of cosmic star formation at $z\sim1-2$, inflow, processing, and feedback drive rapid changes in the spatial distribution and chemical composition of baryons in galaxies and surrounding reservoirs; this baryon cycle can be tomographically mapped by line intensity mapping (LIM) of atomic hydrogen, ionized carbon, and carbon monoxide. We present a simulation-based forecasting framework for detecting auto- and cross-power spectra between spectroscopic surveys of four such tracers at $z\sim0.5-1.7$ mapping the same deep field - TIM, EoRSpec/FYST, MeerKAT, & Euclid. We forward-model 3-D distributions for these tracers from magnetohydrodynamic simulations, directly capturing the two-halo, one-halo, and shot statistics without relying on analytical decompositions. We further detail a signal-to-noise formalism, tailored to LIM surveys with highly anisotropic geometries and Fourier-space coverage. We demonstrate that galaxy cross-correlations will be the dominant discovery channel for current-generation surveys. These instruments will detect the auto-spectra for CO and HI 21 cm and the CO $\times$ 21 cm cross-spectrum at modest S/N $\sim 1-10$, while placing upper limits on the [CII]-158$μ$m signals. [CII], CO, and HI LIM will be $\sim3-30\times$ ($0.5-1.5$ dex) more sensitive to cross-correlation with the Euclid survey, however, than their respective auto-correlations, constraining all three models of line emission at high significance (S/N $\sim 10-40$) within this decade. Finally, we formulate a staged instrumental trajectory with planned or reasonable improvements, including the as-proposed SKA-Mid. We forecast advancing the per-$k$-mode sensitivities of each auto-, galaxy-line, and line-line spectrum by several orders of magnitude, enabling new percent- and sub-percent level constraints on cosmology and the redshift evolution of star formation and the baryon cycle.

Agrawal, Shubh [Pennsylvania U., Dept. Math.]↗
DOE OSTI2025

Anomalies of global symmetries on the lattice

't Hooft anomalies of global symmetries play a fundamental role in quantum many-body systems and quantum field theory (QFT). In this paper, we make a systematic analysis of lattice anomalies - the analog of 't Hooft anomalies in lattice systems - for which we give a precise definition. Crucially, a lattice anomaly is not a feature of a specific Hamiltonian, but rather is a topological invariant of the symmetry action. The controlled setting of lattice systems allows for a systematic and rigorous treatment of lattice anomalies, shorn of the technical challenges of QFT. We find that lattice anomalies reproduce the expected properties of QFT anomalies in many ways, but also have crucial differences. In particular, lattice anomalies and QFT anomalies are not, contrary to a common expectation, in one-to-one correspondence, and there can be non-trivial anomalies on the lattice that are infrared (IR) trivial: they admit symmetric trivial gapped ground states, and map to trivial QFT anomalies at low energies. Nevertheless, we show that lattice anomalies (including IR-trivial ones) have a number of interesting consequences in their own right, including connections to commuting projector models, phases of many-body localized (MBL) systems, and quantum cellular automata (QCA). We make substantial progress on the classification of lattice anomalies and develop several theoretical tools to characterize their consequences on symmetric Hamiltonians. Our work places symmetries of quantum many-body lattice systems into a unified theoretical framework and may also suggest new perspectives on symmetries in QFT.

Disordered Systems and Neural Networks (cond-mat.d↗
DOE OSTI2025

Dark Energy Survey Year 3 results: $w$CDM cosmology from simulation-based inference with persistent homology on the sphere

We present cosmological constraints from Dark Energy Survey Year 3 (DES Y3) weak lensing data using persistent homology, a topological data analysis technique that tracks how features like clusters and voids evolve across density thresholds. For the first time, we apply spherical persistent homology to galaxy survey data through the algorithm TopoS2, which is optimized for curved-sky analyses and HEALPix compatibility. Employing a simulation-based inference framework with the Gower Street simulation suite, specifically designed to mimic DES Y3 data properties, we extract topological summary statistics from convergence maps across multiple smoothing scales and redshift bins. After neural network compression of these statistics, we estimate the likelihood function and validate our analysis against baryonic feedback effects, finding minimal biases (under $0.3σ$) in the $Ω_\mathrm{m}-S_8$ plane. Assuming the $w$CDM model, our combined Betti numbers and second moments analysis yields $S_8 = 0.821 \pm 0.018$ and $Ω_\mathrm{m} = 0.304\pm0.037$-constraints 70% tighter than those from cosmic shear two-point statistics in the same parameter plane. Our results demonstrate that topological methods provide a powerful and robust framework for extracting cosmological information, with our spherical methodology readily applicable to upcoming Stage IV wide-field galaxy surveys.

Prat, J. [Nordita; Royal Inst. Tech., Sodertalje; ↗
DOE OSTI2023

IRIS-DMEM: Efficient Memory Management for Heterogeneous Computing

This paper proposes an efficient data memory management approach for the Intelligent RuntIme System (IRIS) heterogeneous computing framework along with new data transfer policies. IRIS provides a task-based programming model for extreme heterogeneous computing (e.g., CPU, GPU, DSP, FPGA) with support for today's most important programming languages (e.g., OpenMP, OpenCL, CUDA, HIP, OpenACC). However, the IRIS framework either forces the programmer to introduce data transfer commands for each task or relies on suboptimal memory management for automatic and transparent data transfers. The work described here extends IRIS with novel heterogeneous memory handling and introduces novel data transfer policies by employing the Distributed data MEMory handler (DMEM) for efficient and optimal movement of data among the various computing resources. The proposed approach achieves performance gains of up to 7× for tiled LU factorization and tiled DGEMM (i.e., matrix multiplication) benchmarks. Moreover, this approach also reduces data transfers by up to 71% when compared to previous IRIS heterogeneous memory management handlers. This work compares the performance results of the IRIS framework's novel DMEM with the StarPU runtime and MAGMA math library for GPUs. Experiments show a performance gain of up to 1.95× over StarPU and 2.1× over MAGMA.

Miniskar, Narasinga Rao↗
DOE OSTI2021

Normal stability of slow manifolds in nearly periodic Hamiltonian systems

Kruskal [J. Math. Phys. 3, 806 (1962)] showed that each nearly periodic dynamical system admits a formal U(1) symmetry, generated by the so-called roto-rate. We prove that such systems also admit nearly invariant manifolds of each order, near which rapid oscillations are suppressed. We study the nonlinear normal stability of these slow manifolds for nearly periodic Hamiltonian systems on barely symplectic manifolds—manifolds equipped with closed, non-degenerate 2-forms that may be degenerate to leading order. In particular, we establish a sufficient condition for long-term normal stability based on second derivatives of the well-known adiabatic invariant. We use these results to investigate the problem of embedding guiding center dynamics of a magnetized charged particle as a slow manifold in a nearly periodic system. Here, we prove that one previous embedding and two new embeddings enjoy long-term normal stability and thereby strengthen the theoretical justification for these models.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Summer 2025 SULI: Nucleus ID, TinyTPC, and Scientific Communication

This paper summarizes my work during the Summer 2025 SULI internship, which focused on two main projects and broader scientific development. The first project involved improving the particle identification (PID) of protons, deuterons, and tritons using PIDA distributions and template fitting, with the goal of modeling nuclear final-state interactions (FSI) and testing the robustness of the method against systematic uncertainties. These techniques pave the way for future application to LArTPC data from the ICARUS detector. The second project centered on the optimization and data-taking of the TinyTPC detector, a compact LArTPC used for high-resolution low-energy measurements. I adjusted gain and threshold parameters, performed hardware validation tests, and developed analysis strategies to extract meaningful physics from collected data. Throughout the summer, I also enhanced my scientific communication and mentorship skills through presentations, collaborative analysis, and peer guidance.

McCright, Hannah [Maryland U.]↗
DOE OSTI2026

The 3D clustering of Lyman Alpha Emitters measured with DESI

We present a clustering analysis of Lyman-$α$ emitters (LAEs) using spectroscopic observations from the Dark Energy Spectroscopic Instrument (DESI) of candidates selected from the Blanco/DECam Intermediate-Band Imaging Survey (IBIS). We measure the two-point correlation function and the power spectrum, including cross-correlations with DESI quasars. Using both analytical and halo occupation distribution (HOD) simulation-based modeling, we find a linear bias of $b \sim 2.31$--$2.62$ for LAEs over the redshift range $2.26 < z < 3.41$. The analytical modeling also provides constraints on the strength of radiative transfer effects, while the HOD analysis characterizes the LAE-halo connection across multiple models. Finally, we quantify the magnitude of non-perturbative clustering effects such as Fingers of God in the LAE population, providing essential input for the accurate modeling of LAE-based cosmological analyses in forthcoming high-redshift surveys such as DESI-II.

Ebina, H. [UC, Berkeley; LBL, Berkeley] (ORCID:000↗
DOE OSTI2021

Climate Drives Modeled Forest Carbon Cycling Resistance and Resilience in the Upper Great Lakes Region, USA

Forests dominate the global terrestrial carbon budget, but their ability to continue doing so in the face of a changing climate is uncertain. A key uncertainty is how forests will respond (resistance) and recover from (resilience) to rising levels of disturbance of varying intensities. This knowledge gap can optimally be addressed by integrating manipulative field experiments with ecophysiological modeling. We used the ED-2.2 model to project carbon fluxes for a northern temperate deciduous forest subjected to a real-world disturbance severity manipulation experiment. ED-2.2 was run for 150 years, starting from near bare ground in 1900, and subjected to three disturbance treatments under an ensemble of climate conditions. Both disturbance severity and climate strongly affected carbon fluxes, and interacted with one another. We then calculated metrics of functional ecosystem stability. Modeled gross primary production exhibited a mean resistance of -0.13, -0.23, and -0.42 (unitless) in response to increasing disturbance severity of 45%, 65%, and 85% mortality, respectively; and resilience of 0.002, 0.004, 0.008 (yr-1) respectively. This pattern held true for other carbon fluxes, meaning that in general simulated GPP declined more initially but recovered more quickly from higher mortality. Notably, however, heterotrophic respiration responded more slowly to disturbance, and its highly variable response was affected by different drivers. This work provides insight into how future conditions might affect the functional stability of mature forests in this region under ongoing climate change and changing disturbance regimes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Generalized Wigner theorem for noninvertible symmetries

In this article, we establish the conditions under which a conservation law associated with a non-invertible operator may be realized as a symmetry in quantum physics. As established by Wigner, all quantum symmetries must be represented by either unitary or antiunitary transformations. Relinquishing an implicit assumption of invertibility, we demonstrate that the fundamental invariance of quantum transition probabilities under the application of symmetries mandates that all non-invertible symmetries may only correspond to projective unitary or antiunitary transformations, i.e., partial isometries. This extends the notion of physical states beyond conventional rays in Hilbert space to equivalence classes in an extended, gauged Hilbert space, thereby broadening the traditional understanding of symmetry transformations in quantum theory. Our generalized theorem applies irrespective of the origin of the (non)invertible symmetry, holds in arbitrary spatial dimensions, and is independent of the Hamiltonian or action. We explore its physical consequences and, using simple model systems, illustrate how the distinction between invertible and non-invertible symmetries can sometimes be tied to the choice of boundary conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

The DESI DR1 Peculiar Velocity Survey: Growth Rate Measurements from the Galaxy Power Spectrum

The large-scale structure of the Universe and its evolution encapsulate a wealth of cosmological information. A powerful means of unlocking this knowledge lies in measuring the auto-power spectrum and/or the cross-power spectrum of the galaxy density and momentum fields, followed by the estimation of cosmological parameters based on these spectrum measurements. In this study, we generalize the cross-power spectrum model to accommodate scenarios where the density and momentum fields are derived from distinct galaxy surveys. The growth rate of the large-scale structures of the Universe, commonly represented as $fσ_8$, is extracted by jointly fitting the monopole and quadrupole moments of the auto-density power spectrum, the monopole of the auto-momentum power spectrum, and the dipole of the cross-power spectrum. Our estimators, theoretical models and parameter-fitting framework have been tested using mocks, confirming their robustness and accuracy in retrieving the fiducial growth rate from simulation. These techniques are then applied to analyze the power spectrum of the DESI Bright Galaxy Survey and Peculiar Velocity Survey, and the fit result of the growth rate is $fσ_8=0.440^{+0.080}_{-0.096}$ at effective redshift $z_{\rm eff}=0.07$. By synthesizing the fitting outcomes from correlation functions, maximum likelihood estimation and power spectrum, yields a consensus value of $fσ_8(z_{\rm eff}=0.07) = 0.450 ^{+0.055}_{-0.055}$, and correspondingly we obtain $γ=0.580^{+0.110}_{-0.110}$, $Ω_\mathrm{m}=0.301^{+0.011}_{-0.011}$ and $σ_8=0.834^{+0.032}_{-0.032}$. The measured $fσ_8$ and $γ$ are consistent with the prediction of the $Λ$ Cold Dark Matter Model and General Relativity.

Qin, F. [Marseille, CPPM] (ORCID:0000000179507864)↗
DOE OSTI2025

Dark Energy Survey Year 3 Results: Cosmological Constraints from Cluster Abundances, Weak Lensing, and Galaxy Clustering

Galaxy clusters provide a unique probe of the late-time cosmic structure and serve as a powerful independent test of the $\Lambda$CDM model. This work presents the first set of cosmological constraints derived with ~16,000 optically selected redMaPPer clusters across nearly 5,000 $\rm{deg}^2$ using DES Year 3 data sets. Our analysis leverages a consistent modeling framework for galaxy cluster cosmology and DES-Y3 joint analyses of galaxy clustering and weak lensing (3x2pt), ensuring direct comparability with the DES-Y3 3x2pt analysis. We obtain constraints of $S_8 = 0.864 \pm 0.035$ and $\Omega_{\rm{m}} = 0.265^{+0.019}_{-0.031}$ from the cluster-based data vector. We find that cluster constraints and 3x2pt constraints are consistent under the $\Lambda$CDM model with a Posterior Predictive Distribution (PPD) value of $0.53$. The consistency between clusters and 3x2pt provides a stringent test of $\Lambda$CDM across different mass and spatial scales. Jointly analyzing clusters with 3x2pt further improves cosmological constraints, yielding $S_8 = 0.811^{+0.022}_{-0.020}$ and $\Omega_{\rm{m}} = 0.294^{+0.022}_{-0.033}$, a $24\%$ improvement in the $\Omega_{\rm{m}}-S_8$ figure-of-merit over 3x2pt alone. Moreover, we find no significant deviation from the Planck CMB constraints with a probability to exceed (PTE) value of $0.6$, significantly reducing previous $S_8$ tension claims. Finally, combining DES 3x2pt, DES clusters, and Planck CMB places an upper limit on the sum of neutrino masses of $\sum m_\nu < 0.26$ eV at 95% confidence under the $\Lambda$CDM model. These results establish optically selected clusters as a key cosmological probe and pave the way for cluster-based analyses in upcoming Stage-IV surveys such as LSST, Euclid, and Roman.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Confinement transition to gravitational waves in the one-flavor $SU(4)$ Hyper Stealth Dark Matter theory

The thermodynamics of the $SU(4)$ gauge theory with a single flavor of fundamental quarks is analyzed on the lattice with dynamical fermion simulations, which is the low-energy sector of a realistic, strongly-interacting dark matter model -- the Hyper Stealth Dark Matter. The gravitational wave spectrum from the first-order confinement transition in the early universe is further calculated, where the effect of the dark sea quarks, which decrease the interface tension in the effective potential of the Polyakov loop, is shown numerically to lower the gravitational wave amplitude.

Ayyar, V. [Boston U.]↗
DOE OSTI2025

Announcing the Biomedical Data Translator: Initial Public Release

ABSTRACT The growing availability of biomedical data offers vast potential to improve human health, but the complexity and lack of integration of these datasets often limit their utility. To address this, the Biomedical Data Translator Consortium has developed an open‐source knowledge graph–based system—Translator—designed to integrate, harmonize, and make inferences over diverse biomedical data sources. We announce here Translator's initial public release and provide an overview of its architecture, standards, user interface, and core features. Translator employs a scalable, federated, knowledge graph framework for the integration of clinical, genomic, pharmacological, and other biomedical knowledge sources, enabling query retrieval, inference, and hypothesis generation. Translator's user interface is designed to support the exploration of knowledge relationships and the generation of insights, without requiring deep technical expertise and gradually revealing more detailed evidence, provenance, and confidence information, as needed by a given user. To demonstrate Translator's application and impact, we highlight features of the user interface in the context of three real‐world use cases: suggesting potential therapeutics for patients with rare disease; explaining the mechanism of action of a pipeline drug; and screening and validating drug candidates in a model organism. We discuss strengths and limitations of reasoning within a largely federated system and the need for rich concept modeling and deep provenance tracking. Finally, we outline future directions for enhancing Translator's functionality and expanding its data sources. Translator represents a significant step forward in making complex biomedical knowledge more accessible and actionable, aiming to accelerate translational research and improve patient care.

Research & Experimental Medicine↗