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 433 records · Page 24

Data-constrained Solar Modeling with GX Simulator

To facilitate the study of solarflares and active regions, we have created a modeling framework, the freelydistributed GX Simulator IDL package, that combines 3D magnetic and plasma structures with thermal andnonthermal models of the chromosphere, transition region, and corona. Its object-based modular architecture,which runs on Windows, Mac, and Unix/Linux platforms, offers the ability to either import 3D density andtemperature distribution models, or to assign numerically defined coronal or chromospheric temperatures anddensities, or their distributions, to each individual voxel. GX Simulator can apply parametric heating modelsinvolving average properties of the magneticfield lines crossing a given voxel, as well as compute and investigatethe spatial and spectral properties of radio,(sub)millimeter, EUV, and X-ray emissions calculated from the model,and quantitatively compare them with observations. The package includes a fully automatic model productionpipeline that, based on minimal users input, downloads the required SDO/HMI vector magneticfield data,performs potential or nonlinear force-freefield extrapolations, populates the magneticfield skeleton withparameterized heated plasma coronal models that assume either steady-state or impulsive plasma heating, andgenerates non-LTE density and temperature distribution models of the chromosphere that are constrained byphotospheric measurements. The standardized models produced by this pipeline may be further customizedthrough specialized IDL scripts, or a set of interactive tools provided by the graphical user interface. Here, wedescribe the GX Simulator framework and its applications.

Solar active regions↗

Physics-inspired spatiotemporal-graph AI ensemble for the detection of higher order wave mode signals of spinning binary black hole mergers

We present a new class of AI models for the detection of quasi-circular, spinning, non-precessing binary black hole mergers whose waveforms include the higher order gravitational wave modes ($\ell$, |m|) = {(2,2), (2,1), (3,3), (3,2), (4,4)}, and mode mixing effects in the $\ell$ = 3, |m| = 2 harmonics. These AI models combine hybrid dilated convolution neural networks to accurately model both short- and long-range temporal sequential information of gravitational waves; and graph neural networks to capture spatial correlations among gravitational wave observatories to consistently describe and identify the presence of a signal in a three detector network encompassing the Advanced LIGO and Virgo detectors. We first trained these spatiotemporal-graph AI models using synthetic noise, using 1.2 million modeled waveforms to densely sample this signal manifold, within 1.7 h using 256 NVIDIA A100 GPUs in the Polaris supercomputer at the Argonne Leadership Computing Facility. This distributed training approach exhibited optimal classification performance, and strong scaling up to 512 NVIDIA A100 GPUs. With these AI ensembles we processed data from a three detector network, and found that an ensemble of 4 AI models achieves state-of-the-art performance for signal detection, and reports two misclassifications for every decade of searched data. We distributed AI inference over 128 GPUs in the Polaris supercomputer and 128 nodes in the Theta supercomputer, and completed the processing of a decade of gravitational wave data from a three detector network within 3.5 h. Finally, we fine-tuned these AI ensembles to process the entire month of February 2020, which is part of the O3b LIGO/Virgo observation run, and found 6 gravitational waves, concurrently identified in Advanced LIGO and Advanced Virgo data, and zero false positives. This analysis was completed in one hour using one NVIDIA A100 GPU.

79 ASTRONOMY AND ASTROPHYSICS↗

Derivations of the Total Radiation Belt Electron Content

Abstract We present multiple derivations of the Total Radiation Belt Electron Content (TRBEC), an indicator of the global number of electrons that instantaneously occupy the radiation belts. Derived from electron flux measurements, the TRBEC reduces the spatial information into a scalar quantity that concisely describes global aspects of the system. This index provides a simple, global, and long‐term assessment of the radiation belts that enables systematic analysis. In this work, we examine the TRBEC using the adiabatic invariants of which has been used in previous articles as this coordinate system removes reversible adiabatic effects. We then introduce a new expression to compute the TRBEC using the non‐adiabatic coordinates of , relevant in the contexts of energetic electron precipitation, chorus, and hiss scattering where adiabatic invariant quantities are no longer conserved. From both expressions of the TRBEC we demonstrate that an erroneous factor of that appeared in previous works using the adiabatic derivation led to an overestimate of the reported electron populations. In addition, we quantify electron loss in the outer radiation belt via a case study using the Van Allen Probes data over a 20‐day period from March 2013 specifying particle populations both in terms of the aforementioned adiabatic and non‐adiabatic variables. The total number of electrons in the outer radiation belt reached upwards of electrons at the peak of the storm, a rest mass of roughly 10 g.

Astronomy & Astrophysics↗

Damping Wing-Like Features in the Spectra of High Redshift Quasars: a Challenge for Fully-Coupled Simulations

Recently, several observational detections of damping-wing-like features at the edges of “dark gaps” in the spectra of distant quasars (the “Malloy-Lidz effect”) have been reported, rendering strong support for the existence of “neutral islands” in the universe at redshifts as low as z<5.5. We apply the procedure from one of these works, Zhu et al 2024, to the outputs of fully coupled cosmological simulations from two recent large projects, “Cosmic Reionization On Computers” (CROC) and “Thesan”. Synthetic spectra in both simulations have statistics of dark gaps similar to observations, but do not exhibit the damping wing features. Moreover, a toy model with neutral islands added “by hand” only reproduces the observational results when the fraction of neutral islands among all dark gaps approaches 90%. I.e., simulations and observations appear to produce two distinct “populations” of dark gaps. In addition, in the simulations, the neutral islands at z=5.9 should be short-lived and should not extend to z<5.5. A plausible explanation for this discrepancy is that both simulations underestimate the fluctuations in the photoionization rate and, hence, are missing a population of long-lived neutral islands, located in the large downward fluctuations of the photoionization rate.

79 ASTRONOMY AND ASTROPHYSICS↗

The R -Process Alliance: Exploring the cosmic scatter among ten r -process sites with stellar abundances

Context. The astrophysical origin of the rapid neutron-capture process (r-process), responsible for producing roughly half of the elements heavier than iron, remains uncertain. Detailed chemical signatures from the oldest, most metal-poor stars, which act as fossil records of the earliest nucleosynthesis events, can be used to identify the dominant r-process sites. Aims. We present a homogeneous chemical abundance analysis of ten r-process element-enhanced stars. These old and metal-poor stars are strongly enriched in r-process elements with minimal contamination from other nucleosynthetic sources. By focusing on this chemically pure sample, we aim to investigate intrinsic variations in the r-process abundance patterns and explore their implications for the nature and potential diversity of r-process sites. Methods. We performed a detailed chemical abundance analysis of high-resolution, high-signal-to-noise spectra. For each star, we inspected over 1400 individual absorption lines using a combination of equivalent width measurements and spectral synthesis. The analysis was conducted under the assumption of 1D local thermodynamic equilibrium and employing the MOOG radiative transfer code. Results. We derived abundances for 54 chemical species, including 29 neutron-capture (n-capture) elements, covering the full mass range of the r-process abundance pattern. A kinematic analysis reveals that stars likely originated from ten kinematically distinct systems. Based on this assumption, we used the sample to probe the maximum variation expected from ten independent r-process nucleosynthesis events and computed the intrinsic dispersion of each element relative to Zr and Eu for the light and heavy r-process elements, respectively. This exercise resulted in a remarkably low cosmic scatter across the ten r-process sites enriching these stars; for the rare earth and third peak elements, for example, we find σ [La/Eu] = 0.08 and σ [Os/Eu] = 0.11 dex, while the scatter between light and heavy elements, σ [Zr/Eu] , is slightly higher at 0.18 dex. Conclusions. The elemental abundance patterns across the ten independent r-process sites show remarkably small cosmic dispersions. This minimal dispersion suggests a high degree of uniformity in r-process yields across diverse astrophysical environments.

Astronomy and AstroPhysics↗

Distinguishing Orbiting and Infalling Dark Matter Particles with Machine Learning

Dark matter halos are typically defined as spheres that enclose some overdensity, but these sharp, somewhat arbitrary boundaries introduce nonphysical artifacts such as backsplash halos, pseudo-volution, and an incomplete accounting of halo mass. A more physically motivated alternative is to define halos as the collection of particles that are physically orbiting within their potential well. However, existing methods to classify particles as orbiting or infalling suffer from trade-offs between accuracy, computational cost, and generalizability across cosmologies. We present an efficient, yet accurate, supervised machine learning approach using decision trees. The classification is based on only the particle radii and velocities at two epochs. Compared to detailed analysis of particle trajectories, we find that our model matches the classification of 97% of particles. Consequently, we are able to quickly and accurately reproduce the density profiles of the orbiting and infalling components out to many virial radii. We demonstrate that our model generalizes to a significantly different cosmology that lies outside the training data set. We make publicly available both our final model and the code to train similar models.

79 ASTRONOMY AND ASTROPHYSICS↗

On the Physical Nature of Lyα Transmission Spikes in High-redshift Quasar Spectra

We investigate Lyman-alpha (Lyα) transmission spikes at 5.2 < z < 6.8 using synthetic quasar spectra from the “Cosmic Reionization on Computers” simulations. We focus on understanding the relationship between these spikes and the properties of the intergalactic medium (IGM). Disentangling the complex interplay between IGM physics and the influence of galaxies on the generation of these spikes presents a significant challenge. To address this, we employ Explainable Boosting machines, an interpretable machine learning algorithm, to quantify the relative impact of various IGM properties on the Lyα flux. Our findings reveal that gas density is the primary factor influencing absorption strength, followed by the intensity of background radiation and the temperature of the IGM. Ionizing radiation from local sources (i.e., galaxies) appears to have a minimal effect on Lyα flux. The simulations show that transmission spikes predominantly occur in regions of low gas density. Our results challenge recent observational studies suggesting the origin of these spikes in regions with enhanced radiation. We demonstrate that Lyα transmission spikes are largely a product of the large-scale structure, of which galaxies are biased tracers.

79 ASTRONOMY AND ASTROPHYSICS↗

Cosmic structure strikes back: The elimination of vector-mediated nonstandard interaction models as a mechanism for sterile neutrino dark matter production

We revisit sterile neutrino production enabled by nonstandard interactions (NSIs) among active neutrinos mediated by new bosons. We focus on vector mediators, including neutrinophilic, gauged 𝐿 𝜇 −𝐿 𝜏 , and 𝐵−𝐿 realizations, that modify in-medium dispersion and scattering, thereby altering the active-sterile conversion history. Building on a novel production framework with NSI thermal potentials and collision integrals, we compute nonthermal phase-space distributions across sterile neutrino mixing and NSI parameters and map each point to an equivalent thermal warm dark matter particle mass 𝑚 th via linear theory transfer function fitting with the cosmological structure formation Boltzmann solver. This enables a direct reinterpretation of state-of-the-art structure formation limits from Milky Way satellites, strong lensing, and the Lyman-𝛼 forest. These limits, in conjunction with x-ray decay searches, as well as results from a wide variety of particle physics experiments allow for a more complete examination of these models. We find that these vector-mediated models are ruled out when the full combination of current constraints, listed above, are taken into account. NSI scalar-mediated models and models with low reheating temperatures remain viable.

cosmology↗

Cryogenic optical beam steering for superconducting device calibration

We have developed a calibration system based on a micro-electromechanical systems (MEMS) mirror that is capable of delivering an optical beam over a wavelength range of 180 -- 2000 nm (0.62 -- 6.89 eV) in a sub-Kelvin environment. This portable, integrated system can steer the beam over a $\sim$3 cm $\times$ 3 cm area on the surface of any sensor with a precision of $\sim$100 $μ$m, enabling characterization of device response as a function of position. This fills a critical need in the landscape of calibration tools for sub-Kelvin devices, including those used for dark matter detection and quantum computing. These communities have a shared goal of understanding the impact of ionizing radiation on device performance, which can be pursued with our system. This paper describes the design of the first-generation calibration system and the results from successfully testing its performance at room temperature and 20 mK.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

SIM Planetquest Science and Technology: A Status Report

Optical interferometry will open new vistas for astronomy over the next decade. The Space Interferometry Mission (SIM-PlanetQuest), operating unfettered by the Earth's atmosphere, will offer unprecedented astrometric precision that promises the discovery of Earth-analog extra-solar planets as well as a wealth of important astrophysics. Results from SIM will permit the determination of stellar masses to accuracies of 2% or better for objects ranging from brown dwarfs through main sequence stars to evolved white dwarfs, neutron stars, and black holes. Studies of star clusters will yield age determinations and internal dynamics. Microlensing measurements will present the mass spectrum of the Milky Way internal to the Sun while proper motion surveys will show the Sun's orbital radius and speed. Studies of the Galaxy's halo component and companion dwarf galaxies permit the determination of the Milky Way's mass distribution, including its Dark Matter component and the mass distribution and Dark Matter component of the Local Group. Cosmology benefits from precision (1-2%) determination of distances to Cepheid and RR Lyrae standard candles. The emission mechanism of supermassive black holes will be investigated. Finally, radio and optical celestial reference frames will be tied together by an improvement of two orders of magnitude. Optical interferometers present severe technological challenges. The Jet Propulsion Laboratory, with the support of Lockheed Martin Advanced Technology Center (LM ATC) and Northrop Grumman Space Technology (NGST), has addressed these challenges with a technology development program that is now complete. The requirements for SIM have been satisfied, based on outside peer review, using a series of laboratory tests and appropriate computer simulations: laser metrology systems perform with 10 picometer precision; mechanical vibrations have been controlled to nanometers, demonstrating orders of magnitude disturbance rejection; and knowledge of component positions throughout the whole test assembly has been demonstrated to the required picometer level. Technology transfer to the SIM flight team is now well along.

astrometry↗

Cosmic Reionization on Computers: Statistical Properties of the Distributions of Mean Opacities

Quasar absorption lines provide a unique window to the relationship between galaxies and the intergalactic medium during the Epoch of Reionization. In particular, high redshift quasars enable measurements of the neutral hydrogen content of the universe. However, the limited sample size of observed quasar spectra, particularly at the highest redshifts, hampers our ability to fully characterize the intergalactic medium during this epoch from observations alone. In this work, we characterize the distributions of mean opacities of the intergalactic medium in simulations from the Cosmic Reionization on Computers (CROC) project. We find that the distribution of mean opacities along sightlines follows a non-trivial distribution that cannot be easily approximated by a known distribution. When comparing the cumulative distribution function of mean opacities measurements in subsamples of sample sizes similar to observational measurements from the literature, we find consistency between CROC and observations at redshifts $z\lesssim 5.7$. However, at higher redshifts ($z\gtrsim5.7$), the cumulative distribution function of mean opacities from CROC is notably narrower than those from observed quasar sightlines implying that observations probe a systematically more opaque intergalactic medium at higher redshifts than the intergalactic medium in CROC boxes at these same redshifts. This is consistent with previous analyses that indicate that the universe is reionized too early in CROC simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

Oxygen Opacity Measurements at High-Energy-Density Conditions

We present the first benchmark measurements of oxygen opacity at high-energy-density, crucial for understanding stellar interiors. A 3–5−μm thick SiO 2 sample was heated using a Z facility x-ray source, and opacity was measured in the 5–19.5 Å range with multiple spectrometers. The electron temperature and density inferred from Si 𝐾-shell spectra were 148 ± 4 eV and 8.6 ±1.4 × 10 21 cm −3 , respectively. Six experiments showed opacity reproducibility within ±15% and 5%–25% opacity uncertainty. Measured oxygen opacity and seven different opacity models exhibit extraordinary overall agreement. Discrepancies near the H-like photon ionization threshold and in certain spectral lines suggest the need for model refinements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

An analytically tractable marked power spectrum

The increasing precision of cosmology data in the modern era is calling for methods to allow the extraction of non-Gaussian information using tools beyond two-point statistics. The marked power spectrum has the potential to extract beyond two-point information in a computationally efficient way while using much of the infrastructure already available for the power spectrum. In this work we explore the marked power spectrum from an analytical perspective. In particular, we explore a low-order polynomial for the mark that allows us to better control the theoretical uncertainties and we show that with minimal new degrees of freedom the analytical results match measurements from N-body simulations for both the matter field and biased tracers in redshift space. Finally, we show that even within the limited forms of mark that we consider, there are degeneracies that can be broken by inclusion of the marked auto-spectrum or the cross-spectrum with the unmarked field. I n conclusion, we discuss future theoretical developments that would enable us to apply this approach to survey data.

79 ASTRONOMY AND ASTROPHYSICS↗

LUVOIR Primary Mirror Segment Alignment Control With Joint Laser Metrology and Segment Edge Sensing

An approach is developed for the alignment and stability maintenance of the LUVOIR segmented primary mirror using a segment state estimation and wavefront control method based on a hybrid segment motion sensing architecture of laser truss metrology and segment edge sensors. Our current computer model was generated for LUVOIR Architecture Option A with a 15-meter aperture, 120segment primary mirror. The methodology and simulation results will be presented and analyzed. JPL has a long history of technology development in laser metrology and edge sensors, including work in SIM [7], Keck and TMT [8], CCAT [3] and LUVOIR [1]. We will discuss our current efforts of LUVOIR laser metrology and edge-sensor models development, showing sensitivities of sensor measurements to various mirror eigenmodes, removing global modes and strengthening weak modes by performing joint (hybrid) laser-metrology and edge sensing. We will define and derive an important performance metric called wavefront error multiplier (WEM), and show that WEM provides a simple link between sensor errors and the closed-loop (controlled) system wavefront error. We will show WEM values for several hybrid sensor configuration options studied. We will discuss an algorithm for mirror shape control and maintenance through segment state and wavefront estimations using joint edge-metrology sensing. We will compare simulated performance of mirror state estimation, wavefront estimation and wavefront control based on joint edge-metrology sensing among several sensor configurations, and show the impact of sensor error distributions on the segmented mirror alignment performance. Mirror shape control performance will be also be evaluated in terms of imaging contrast between inner working angles (IWA) and outer working angles (OWA) of a LUVOIR coronagraph. The research was carried out at the Jet Propulsion Laboratory, California Institute of Technology, under a contract with the National Aeronautics and Space Administration. LUVOIR is a mission concept under study for the next astronomy Decadal Survey—the information presented in this paper is pre-decisional and is provided for planning and discussion purposes only.

Shelton, Chris↗

Toward Accelerated Nuclear-physics Parameter Estimation from Binary Neutron Star Mergers: Emulators for the Tolman–Oppenheimer–Volkoff Equations

Abstract Gravitational-wave observations of binary neutron-star (BNS) mergers have the potential to revolutionize our understanding of the nuclear equation of state (EOS) and the fundamental interactions that determine its properties. However, Bayesian parameter estimation frameworks do not typically sample over microscopic nuclear-physics parameters that determine the EOS. One of the major hurdles in doing so is the computational cost involved in solving the neutron-star structure equations, known as the Tolman–Oppenheimer–Volkoff (TOV) equations. In this paper, we explore approaches to emulating solutions for the TOV equations: multilayer perceptrons (MLPs), Gaussian processes, and a data-driven variant of the reduced basis method (RBM). We implement these emulators for three different parameterizations of the nuclear EOS, each with a different degree of complexity represented by the number of model parameters. We find that our MLP-based emulators are generally more accurate than the other two algorithms, whereas the RBM results in the largest speedup with respect to the full high-fidelity TOV solver. We employ these emulators for a simple parameter inference using a potentially loud BNS observation and show that the posteriors predicted by our emulators are in excellent agreement with those obtained from the full TOV solver.

79 ASTRONOMY AND ASTROPHYSICS↗

Force-free Wave Interaction in Magnetar Magnetospheres: Computational Modeling in Axisymmetry

Abstract Crustal quakes of highly magnetized neutron stars can disrupt their magnetospheres, triggering energetic phenomena like X-ray and fast radio bursts. Understanding plasma wave dynamics in these extreme environments is vital for predicting energy transport across scales to the radiation length. This study models relativistic plasma wave interaction in magnetar magnetospheres with force-free electrodynamics simulations. For propagation along curved magnetic field lines, we observe the continuous conversion of Alfvén waves to fast magnetosonic (FMS) waves. The conversion efficiency can be up to three times higher when counter-propagating Alfvén waves interact in the equatorial region. Alfvén waves generate FMS waves of twice their frequency during their first crossing of the magnetosphere. After the initial transient burst of FMS waves, Alfvén waves convert to FMS waves periodically, generating variations on timescales of the magnetospheric Alfvén wave crossing time. This decaying FMS wave tail carries a significant portion (half) of the total energy emitted. Plastic damping of “bouncing” Alfvén waves by the magnetar crust has minimal impact on the FMS efficiency. We discuss the implications of the identified wave phenomena for magnetar observations. Outgoing FMS waves can develop electric zones, potential sources of coherent radiation. Long wavelength FMS waves could generate FRBs through reconnection beyond the light cylinder.

Astronomy & Astrophysics↗

In-depth analysis of solar models with high-metallicity abundances and updated opacity tables

Context. As a result of the high-quality constraints available for the Sun, we are able to carry out detailed combined analyses using neutrino, spectroscopic, and helioseismic observations. These studies lay the ground for future improvements of the key physical components of solar and stellar models because ingredients such as the equation of state, the radiative opacities, or the prescriptions for macroscopic transport processes of chemicals are then used to study other stars in the Universe. Aims. We study the existing degeneracies in solar models using the recent high-metallicity spectroscopic abundances by comparing them to helioseismic and neutrino data and discuss the effect on their properties of changes in the micro and macro physical ingredients. Methods. We carried out a detailed study of solar models computed with a high-metallicity composition from the literature based on averaged 3D models that were claimed to resolve the solar modelling problem. We compared these models to helioseismic and neutrino constraints. Results. The properties of the solar models are significantly affected by the use of the recent OPLIB opacity tables and the inclusion of macroscopic transport. The properties of the standard solar models computed using the OPAL opacities are similar to those for which the OP opacities were used. We show that a modification of the temperature gradient just below the base of the convective zone is required to remove the discrepancies in solar models, particularly in the presence of macroscopic mixing. This can be simulated by a localised increase in the opacity of a few percent. Conclusions. We conclude that the existing degeneracies and issues in solar modelling are not removed by using an increase in the solar metallicity, in contradiction to what has been suggested in the recent literature. Therefore, standard solar models cannot be used as an argument for a high-metallicity composition. While further work is required to improve solar models, we note that direct helioseismic inversions indicate a low metallicity in the convective envelope, in agreement with spectroscopic analyses based on full 3D models.

74 ATOMIC AND MOLECULAR PHYSICS↗

Interacting dark sector within ETHOS: Cosmological constraints from SPT cluster abundance with DES and HST weak lensing data

We use galaxy cluster abundance measurements from the South Pole Telescope enhanced by multicomponent matched filter confirmation and complemented with mass information obtained using weak-lensing data from Dark Energy Survey Year 3 (DES Y3) and targeted Hubble Space Telescope observations for probing deviations from the cold dark matter paradigm. Concretely, we consider a class of dark sector models featuring interactions between dark matter (DM) and a dark radiation (DR) component within the framework of the effective theory of structure formation (ETHOS). We focus on scenarios that lead to power suppression over a wide range of scales, and thus can be tested with data sensitive to large scales, as realized, for example, for DM–DR interactions following from an unbroken non-Abelian 𝑆⁢𝑈⁡(𝑁) gauge theory (interaction rate with power-law index 𝑛 = 0 within the ETHOS parametrization). Cluster abundance measurements are mostly sensitive to the amount of DR interacting with DM, parametrized by the ratio of DR temperature to the cosmic microwave background (CMB) temperature, 𝜉 DR = 𝑇 DR /𝑇 CMB . We find an upper limit 𝜉 DR < 17% at 95% credibility. When the cluster data are combined with Planck 2018 CMB data along with baryon acoustic oscillation (BAO) measurements we find 𝜉 DR < 10%, corresponding to a limit on the abundance of interacting DR that is around 3 times tighter than that from CMB + BAO data alone. We also discuss the complementarity of weak lensing informed cluster abundance studies with probes sensitive to smaller scales, explore the impact on our analysis of massive neutrinos, and comment on a slight preference for the presence of a nonzero interacting DR abundance, which enables a physical solution to the 𝑆 8 tension.

79 ASTRONOMY AND ASTROPHYSICS↗