Search NASA⌕ Search

SEARCH · Search NASA

Results for “astronomy computing”

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

Not-quite-transcendental Functions for Logarithmic Interpolation of Tabulated Data

From tabulated nuclear and degenerate equations of state to photon and neutrino opacities and nuclear reaction rates, tabulated data is ubiquitous in computational astrophysics. The dynamic range that must be covered by these tables typically spans many orders of magnitude. Here we present a novel strategy for accurately and performantly interpolating tabulated data that spans these large dynamic ranges. We demonstrate the efficacy of this strategy in tabulated lookups for nuclear and terrestrial equations of state. We show that this strategy is a faster drop-in replacement for linear interpolation of logarithmic grids.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

AthenaK: A Performance-portable Version of the Athena++ Adaptive Mesh Refinement Framework

We describe AthenaK: a new implementation of the Athena++ block-based adaptive mesh refinement framework using the Kokkos programming model. Finite volume methods for Newtonian, special relativistic, and general relativistic (GR) hydrodynamics and magnetohydrodynamics (MHD), and GR-radiation hydrodynamics and MHD, as well as a module for evolving Lagrangian tracer or charged test particles (e.g., cosmic rays) are implemented using the framework. In two companion papers, we describe (1) a new solver for the Einstein equations based on the Z4c formalism, and (2) a GRMHD solver in dynamical spacetimes also implemented using the framework, enabling new applications in numerical relativity. By adopting Kokkos, the code can be run on virtually any hardware, including CPUs, GPUs from multiple vendors, and emerging Advanced RISC Machine processors. AthenaK shows excellent performance and weak scaling, achieving over 1 billion cell updates per second for hydrodynamics in three dimensions on a single NVIDIA Grace Hopper processor. It does this with a typical parallel efficiency of 80% on 65,536 AMD GPUs on the OLCF Frontier system. Such performance portability enables AthenaK to leverage modern exascale computing systems for challenging applications in astrophysical fluid dynamics, numerical relativity, and multimessenger astrophysics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Constituent model of light hybrid meson decays

A model of light hybrid mesons and their strong decays is developed. The model employs a gluonic quasiparticle to describe low energy gluodynamics and uses the QCD Hamiltonian in Coulomb gauge to guide the construction of states and decay amplitudes. We compute the partial widths of the twelve low lying isovector and vector hybrids. Implications of these results on hybrid searches are also made, with the chief conclusions being that direct observation of the vector states will be difficult, that a hybrid π ( 1800 ) has distinctive decay characteristics, a narrow η ( 1900 ) hybrid should exist, an η 1 ( 1750 ) should be sought, and that the exotic nature of J P C = 2 − + hybrid mesons should be discernible with sufficient data. We argue that the isovector π 2 hybrid has been discovered, giving a total of four possible hybrid mesons, π 1 ( 1600 ) , η 1 ( 1855 ) , π ( 1800 ) , and π 2 ( 2360 ) , which appear to be filling out the low lying hybrid supermultiplet in the expected fashion. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
NASA NTRS2018

The Apollo Sample Suite: 50 Years of Solar System Insight

The Apollo program was undoubtable a crowning achievement in human history. In addition to the obvious cultural significance, scientific results from the Apollo program had a lasting impression on a range of scientific fields, none more so that the effect the samples had on the fields of geology and cosmochemistry. The six Apollo missions collected 382 kg of rock, regolith, and core samples from geologically diverse locations on the Moon. In the nearly 50 years since the first samples were returned, there have been over 3000 different requests for samples, each yielding insights into fields as disparate as biology, medicine, astronomy, engineering, material science, and of course geology. Early studies of the Apollo samples revealed primary insights into the origin and evolution of the Moon, and of the Earth-Moon system, but the results also had implications for bodies throughout the solar system, e.g., defining crater counting rates. Over the decades, continued study of the Apollo samples by new generations of scientists using new instruments have continued to yield significant new discoveries, including the presence of endogenous water in the Moon and the possible presence of a lunar cataclysm, that in turn has contributed to new models of solar system formation and evolution. The Apollo samples have often been used as a proxy for studying other bodies like Mercury or asteroids. The Apollo samples have also directly contributed to the interpretation of remotely sensed data sets, including their use as ground truth for both Clementine and Lunar Prospector global geochemical maps. Despite the Apollo samples being a static collection, recent efforts will ensure that investigators continue to have access to new samples. For example, there was a recent solicitation for study of previously unopened Apollo samples in vacuum-sealed containers, as well as new access to samples stored frozen or in a He atmosphere. Similarly, the use of X-ray computed tomography as part of the curation process is identifying new clasts within polymict breccias that are available for study. Finally, the MoonDB project is putting all previously published lunar geochemical analyses into a searchable database, which should facilitate new investigations.

Zeigler, Ryan↗
DOE OSTI2024

Anharmonic effects in nuclear recoils from sub-GeV dark matter

Direct detection experiments are looking for nuclear recoils from scattering of sub-GeV dark matter (DM) in crystals, and have thresholds as low as ∼ 10 eV or DM masses of ∼ 100 MeV . Future experiments are aiming for even lower thresholds. At such low energies, the free nuclear recoil prescription breaks down, and the relevant final states are phonons in the crystal. Scattering rates into single as well as multiple phonons have already been computed for a harmonic crystal. However, crystals typically exhibit some anharmonicity, which can significantly impact scattering rates in certain kinematic regimes. In this work, we estimate the impact of anharmonic effects on scattering rates for DM in the mass range ∼ 1 – 10 MeV , where the details of multiphonon production are most important. Using a simple model of a nucleus in a bound potential, we find that anharmonicity can modify the scattering rates by up to two orders of magnitude for DM masses of O ( MeV ) . However, such effects are primarily present at high energies where the rates are suppressed, and thus only relevant for very large DM cross sections. We show that anharmonic effects are negligible for masses larger than ∼ 10 MeV . Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2023

RadioGalaxyNET: Dataset and novel computer vision algorithms for the detection of extended radio galaxies and infrared hosts

Abstract Creating radio galaxy catalogues from next-generation deep surveys requires automated identification of associated components of extended sources and their corresponding infrared hosts. In this paper, we introduce RadioGalaxyNET, a multimodal dataset, and a suite of novel computer vision algorithms designed to automate the detection and localization of multi-component extended radio galaxies and their corresponding infrared hosts. The dataset comprises 4 155 instances of galaxies in 2 800 images with both radio and infrared channels. Each instance provides information about the extended radio galaxy class, its corresponding bounding box encompassing all components, the pixel-level segmentation mask, and the keypoint position of its corresponding infrared host galaxy. RadioGalaxyNET is the first dataset to include images from the highly sensitive Australian Square Kilometre Array Pathfinder (ASKAP) radio telescope, corresponding infrared images, and instance-level annotations for galaxy detection. We benchmark several object detection algorithms on the dataset and propose a novel multimodal approach to simultaneously detect radio galaxies and the positions of infrared hosts.

Astronomy & Astrophysics↗
DOE OSTI2025

$\mathrm{SageNet}$: Fast Neural Network Emulation of the Stiff-amplified Gravitational Waves from Inflation

Accurate modeling of the inflationary gravitational waves (GWs) requires time-consuming, iterative numerical integrations of differential equations to take into account their backreaction on the expansion history. To improve computational efficiency while preserving accuracy, we present the Stiff-amplified Gravitational-wave Emulator Network (SageNet), a deep learning framework designed to replace conventional numerical solvers (code available at https://github.com/YifangLuo/SageNet). SageNet employs a long short-term memory architecture to emulate the present-day energy density spectrum of the inflationary GWs with possible stiff amplification, Ω GW (f). Trained on a data set of 25,689 numerically generated solutions, SageNet allows accurate reconstructions of Ω GW (f) and generalizes well to a wide range of cosmological parameters; 90.9% of the test emulations with randomly distributed parameters exhibit errors of under 4%. In addition, SageNet demonstrates its ability to learn and reproduce the artificial, adaptive sampling patterns in numerical calculations, which implement denser sampling of frequencies around changes in spectral indices in Ω GW (f). The dual capability of learning both physical and artificial features of the numerical GW spectra establishes SageNet as a robust alternative to exact numerical methods. Finally, our benchmark tests show that SageNet reduces the computation time from tens of seconds to milliseconds, achieving a speedup of ∼10 4 times over standard CPU-based numerical solvers with the potential for further acceleration on GPU hardware. These capabilities make SageNet a powerful tool for accelerating Bayesian inference procedures for extended cosmological models. In a broad sense, the SageNet framework offers a fast, accurate, and generalizable solution to modeling cosmological observables whose theoretical predictions demand costly differential equation solvers.

Astronomy data modeling↗
NASA NTRS2002

Optical Testing of Diamond Machined, Aspheric Mirrors for Groundbased, Near-IR Astronomy

The Infrared Multi-Object Spectrometer (IRMOS) is a facility-class instrument for the Kitt Peak National Observatory 4 and 2.1 meter telescopes. IRMOS is a near-IR (0.8-2.5 micron) spectrometer and operates at approximately 80 K. The 6061-T651 aluminum bench and mirrors constitute an athermal design. The instrument produces simultaneous spectra at low- to mid-resolving power (R=lambda/delta lambda= 300-3000) of approximately 100 objects in its 2.8 x 2.0 arcmin field. We describe ambient and cryogenic optical testing of the IRMOS mirrors across a broad range in spatial frequency (figure error, mid-frequency error, and microroughness). The mirrors include three rotationally symmetric, off-axis conic sections, one off-axis biconic, and several flat fold mirrors. The symmetric mirrors include convex and concave prolate and oblate ellipsoids. They range in aperture from 94x86 mm to 286x269 mm and in f-number from 0.9 to 2.4. The biconic mirror is concave and has a 94x76 mm aperture, R(sub x)=377 mm, k(sub x)=0.0778, R(sub y)=407 mm, and k(sub y)=0.1265 and is decentered by -2 mm in X and 227 mm in Y. All of the mirrors have an aspect ratio of approximately 6:1. The surface error fabrication tolerances are less than 10 nm RMS microroughness, 'best effort' for mid-frequency error, and less than 63.3 nm RMS figure error. Ambient temperature (approximately 293 K) testing is performed for each of the three surface error regimes, and figure testing is also performed at approximately 80 K. Operation of the ADE Phaseshift MicroXAM white light interferometer (micro-roughness) and the Bauer Model 200 profilometer (mid-frequency error) is described. Both the sag and conic values of the aspheric mirrors make these tests challenging. Figure testing is performed using a Zygo GPI interferometer, custom computer generated holograms (CGH), and optomechanical alignment fiducials. Cryogenic CGH null testing is discussed in detail. We discuss complications such as the change in prescription with temperature and thermal gradients. Correction for the effect of the dewar window is also covered. We discuss the error budget for the optical test and alignment procedure. Data reduction is accomplished using commercial optical design and data analysis software packages. Results from CGH testing at cryogenic temperatures are encouraging thus far.

Chambers, V. John↗
DOE OSTI2024

CIRCLEZ : Reliable photometric redshifts for active galactic nuclei computed solely using photometry from Legacy Survey Imaging for DESI

Photometric redshifts for galaxies hosting an accreting supermassive black hole in their center, known as active galactic nuclei (AGNs), are notoriously challenging. At present, they are most optimally computed via spectral energy distribution (SED) fittings, assuming that deep photometry for many wavelengths is available. However, for AGNs detected from all-sky surveys, the photometry is limited and provided by a range of instruments and studies. This makes the task of homogenizing the data challenging, presenting a dramatic drawback for the millions of AGNs that wide surveys such as SRG/eROSITA are poised to detect. This work aims to compute reliable photometric redshifts for X-ray-detected AGNs using only one dataset that covers a large area: the tenth data release of the Imaging Legacy Survey (LS10) for DESI. LS10 provides deep grizW1-W4 forced photometry within various apertures over the footprint of the eROSITA-DE survey, which avoids issues related to the cross-calibration of surveys. We present the results from CIRCLEZ, a machine-learning algorithm based on a fully connected neural network. CIRCLEZ is built on a training sample of 14 000 X-ray-detected AGNs and utilizes multi-aperture photometry, mapping the light distribution of the sources. The accuracy (σNMAD) and the fraction of outliers (η) reached in a test sample of 2913 AGNs are equal to 0.067 and 11.6%, respectively. The results are comparable to (or even better than) what was previously obtained for the same field, but with much less effort in this instance. We further tested the stability of the results by computing the photometric redshifts for the sources detected in CSC2 and Chandra-COSMOS Legacy, reaching a comparable accuracy as in eFEDS when limiting the magnitude of the counterparts to the depth of LS10. The method can be applied to fainter samples of AGNs using deeper optical data from future surveys (for example, LSST, Euclid), granting LS10-like information on the light distribution beyond the morphological type. Along with this paper, we have released an updated version of the photometric redshifts (including errors and probability distribution functions) for eROSITA/eFEDS.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Robust Dark Energy Constraints with the Dark Energy Spectroscopic Survey (Final Technical Report)

This project developed and applied advanced theoretical, computational, and data-analysis methodologies to extract robust and precise cosmological constraints from the Dark Energy Spectroscopic Instrument (DESI). The work focused on maximizing the scientific return of DESI through optimized survey strategy, novel higher-order clustering statistics, improved modeling of small-scale structure, and rigorous mitigation of observational systematics. Over the award period, the project made substantial contributions to DESI science planning, produced new methods for bispectrum and three-point correlation function analyses, advanced constraints on primordial non-Gaussianity, and delivered widely used software tools. The project also played a major role in training graduate students and a postdoctoral researcher who contributed directly to DESI key projects. The results have significantly enhanced the cosmological reach of DESI and provide a strong foundation for future surveys such as DESI-II and Stage-V experiments.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2016

Systematic Comparison of Photoionized Plasma Codes with Application to Spectroscopic Studies of AGN in X-Rays

Atomic data and plasma models play a crucial role in the diagnosis and interpretation of astrophysical spectra, thus influencing our understanding of the Universe. In this investigation we present a systematic comparison of the leading photoionization codes to determine how much their intrinsic differences impact X-ray spectroscopic studies of hot plasmas in photoionization equilibrium. We carry out our computations using the Cloudy, SPEX, and XSTAR photoionization codes, and compare their derived thermal and ionization states for various ionizing spectral energy distributions. We examine the resulting absorption-line spectra from these codes for the case of ionized outflows in active galactic nuclei. By comparing the ionic abundances as a function of ionization parameter, we find that on average there is about 30 deviation between the codes in where ionic abundances peak. For H-like to B-like sequence ions alone, this deviation in is smaller at about 10 on average. The comparison of the absorption-line spectra in the X-ray band shows that there is on average about 30 deviation between the codes in the optical depth of the lines produced at log 1 to 2, reducing to about 20 deviation at log 3. We also simulate spectra of the ionized outflows with the current and upcoming high-resolution X-ray spectrometers, on board XMM-Newton, Chandra, Hitomi, and Athena. From these simulations we obtain the deviation on the best-fit model parameters, arising from the use of different photoionization codes, which is about 10 to40. We compare the modeling uncertainties with the observational uncertainties from the simulations. The results highlight the importance of continuous development and enhancement of photoionization codes for the upcoming era of X-ray astronomy with Athena.

plasmas – atomic processes – atomic data – t↗
DOE OSTI2024

A Full Accounting of the Visible Mass in SDSS MaNGA Disk Galaxies

We present a study of the ratio of visible mass to total mass in spiral galaxies to better understand the relative amount of dark matter present in galaxies of different masses and evolutionary stages. Using the velocities of the Hα emission line measured in spectroscopic observations from the Sloan Digital Sky Survey (SDSS) MaNGA Data Release 17 (DR 17), we evaluate the rotational velocity of over 5500 disk galaxies at their 90% elliptical Petrosian radii, R 90 . We compare this to the velocity expected from the total visible mass, which we compute from the stellar, H i, H 2 , and heavy metals and dust masses. H 2 mass measurements are available for only a small subset of galaxies observed in SDSS MaNGA DR17, so we derive a parameterization of the H 2 mass as a function of absolute magnitude in the r band using galaxies observed as part of SDSS DR7. With these parameterizations, we calculate the fraction of visible mass within R 90 that corresponds to the observed velocity. Based on statistically analyzing the likelihood of this fraction, we conclude that the null hypothesis (no dark matter) cannot be excluded at a confidence level better than 95% within the visible extent of the disk galaxies. We also find that when all mass components are included, the ratio of visible to total mass within the visible extent of star-forming disk galaxies increases with galaxy luminosity.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Performance of high-order Godunov-type methods in simulations of astrophysical low Mach number flows

High-order Godunov methods for gas dynamics have become a standard tool for simulating different classes of astrophysical flows. Their accuracy is mostly determined by the spatial interpolant used to reconstruct the pair of Riemann states at cell interfaces and by the Riemann solver that computes the interface fluxes. In most Godunov-type methods, these two steps can be treated independently, so that many different schemes can in principle be built from the same numerical framework. Because astrophysical simulations often test out the limits of what is feasible with the computational resources available, it is essential to find the scheme that produces the numerical solution with the desired accuracy at the lowest computational cost. However, establishing the best combination of numerical options in a Godunov-type method to be used for simulating a complex hydrodynamic problem is a nontrivial task. In fact, formally more accurate schemes do not always outperform simpler and more diffusive methods, especially if sharp gradients are present in the flow. For this work, we used our fully compressible Seven-League Hydro (SLH) code to test the accuracy of six reconstruction methods and three approximate Riemann solvers on two- and three-dimensional (2D and 3D) problems involving subsonic flows only. We considered Mach numbers in the range from 10 −3 to 10 −1 , which are characteristic of many stellar and geophysical flows. In particular, we considered a well-posed, 2D, Kelvin–Helmholtz instability problem and a 3D turbulent convection zone that excites internal gravity waves in an overlying stable layer. Although the different combinations of numerical methods converge to the same solution with increasing grid resolution for most of the quantities analyzed here, we find that (i) there is a spread of almost four orders of magnitude in computational cost per fixed accuracy between the methods tested in this study, with the most performant method being a combination of a low-dissipation Riemann solver and a sextic reconstruction scheme; (ii) the low-dissipation solver always outperforms conventional Riemann solvers on a fixed grid when the reconstruction scheme is kept the same; (iii) in simulations of turbulent flows, increasing the order of spatial reconstruction reduces the characteristic dissipation length scale achieved on a given grid even if the overall scheme is only second order accurate; (iv) reconstruction methods based on slope-limiting techniques tend to generate artificial, high-frequency acoustic waves during the evolution of the flow; and (v) unlimited reconstruction methods introduce oscillations in the thermal stratification near the convective boundary, where the entropy gradient is steep.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

JAXtronomy: A JAX port of lenstronomy

Gravitational lensing is a phenomenon where light bends around massive objects, resulting in distorted images seen by an observer. Studying gravitationally lensed systems provides insights into cosmology and astrophysics, including constraints of the expansion rate of the Universe and the distribution of dark matter. Thus, we introduce JAXtronomy, a re-implementation of the gravitational lensing software package lenstronomy (Birrer, 2021; Birrer & Amara, 2018) using JAX (Bradbury et al., 2018). JAX is a Python library that uses an accelerated linear algebra (XLA) compiler to improve the performance of computing software. Our core design principle of JAXtronomy is to maintain an identical API to that of lenstronomy. The main JAX features utilized in JAXtronomy are just-in-time compilation, which can lead to significant reductions in execution time, and automatic differentiation, which allows for the implementation of gradient-based algorithms that were previously impossible. Additionally, JAX allows code to be run on GPUs or parallelized across CPU cores, further boosting the performance of JAXtronomy.

astronomy↗
NASA NTRS1992

The SERENDIP piggyback SETI project

The SERENDIP project is an ongoing program of monitoring and processing broadband radio signals acquired by existing radio astronomy observatories. SERENDIP operates in a piggyback mode: it makes use of whatever observing plan (sequence of frequencies, sky coordinates, and polarizations) is under way at its host observatory. Moreover, the SERENDIP data acquisition system, once installed, operates autonomously. This approach makes it possible to obtain large amounts of high quality observing time in a manner that is economical and that does not adversely affect ongoing radio astronomy survey work. The SERENDIP II system has been installed at the NRAO 300-foot telescope at Green Bank, West Virginia, and has operated there for several thousand hours. In this report, we summarize our findings from these observations and describe the present status of the project. Two key elements of SERENDIP are the automated data acquisition system that uses adaptive thresholds and logs only statistically significant peaks in the real-time power spectra, and the subsequent off-line analysis programs that identify and reject a variety of interference signals. Several specific correlations have been identified that offer promise. At present, the development and testing of these interference rejection algorithms is the main thrust of our work.

NASA Discipline Number 52-60↗
DOE OSTI2025

Quantum‐Classical Simulation of Quantum Field Theory by Quantum Circuit Learning

Abstract Quantum circuit learning is employed to simulate quantum field theories (QFTs). Typically, when simulating QFTs with quantum computers, significant challenges are encountered due to the technical limitations of quantum devices when implementing the Hamiltonian using Pauli spin matrices. To address this challenge, quantum circuit learning is leveraged, employing a compact configuration of qubits and low‐depth quantum circuits to predict real‐time dynamics in quantum field theories. The key advantage of this approach is that a single‐qubit measurement can accurately forecast various physical parameters, including fully‐connected operators. To demonstrate the effectiveness of this method, it is used to predict quench dynamics, chiral dynamics and jet production in a 1+1‐dimensional model of quantum electrodynamics. It is found that our predictions closely align with the results of rigorous classical calculations, exhibiting a high degree of accuracy. This hybrid quantum‐classical approach illustrates the feasibility of efficiently simulating large‐scale QFTs on cutting‐edge quantum devices.

Ikeda, Kazuki [Department of Physics University of↗
DOE OSTI2025

Shaping dark photon spectral distortions

The cosmic microwave background (CMB) spectrum is an extraordinary tool for exploring physics beyond the Standard Model. The exquisite precision of its measurement makes it particularly sensitive to small effects caused by hidden sector interactions. In particular, CMB spectral distortions can unveil the existence of dark photons which are kinetically coupled to the standard photon. In this work, we use the COBE-FIRAS dataset to derive accurate and robust limits on photon-to-dark-photon oscillations for a large range of dark photon masses, from 10$^{−10}$ to 10$^{−4}$ eV. We consider in detail the redshift dependence of the bounds, computing CMB distortions due to photon injection/removal using a Green’s function method. Our treatment improves on previous results, which had set limits studying energy injection/removal into baryons rather than photon injection/removal, or ignoring the redshift evolution of distortions. The difference between our treatment and previous ones is particularly noticeable in the predicted spectral shape of the distortions, a smoking gun signature for photon-to-dark-photon oscillations. The characterization of the spectral shape is crucial for future CMB missions, which could improve the present sensitivity by orders of magnitude, exploring regions of the dark photon parameter space that are otherwise difficult to access [inline-graphic not available: see fulltext].

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

A comparison of the dielectric function for p-point water models

We compute the longitudinal dielectric function ɛl(k) for p-point water models, p = 3 (SPC/E, TIP3P), p = 4 (TIP4P-EW, OPC), and p = 5 (TIP5P-E) from the charge–charge fluctuation function SZZ(k), obtained from two methods: average of the charge in k-space (method:sum) and from the pair distribution function in Fourier space (method:integral). The latter requires a continuation to small k-values, which we thoroughly discuss. We conclude with a detailed comparison of the longitudinal dielectric function for the different p-point models and its main characteristics: small k and large k limits, position of the poles, and its (complex) zeros.

Shevy, Noah [Ames National Laboratory, Department ↗