Search NASA⌕ Search

SEARCH · Search NASA

Results for “astronomy data modeling”

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 19 records
All sourcesNASA NTRSDOE OSTI
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 NTRS1992

C++, objected-oriented programming, and astronomical data models

Contemporary astronomy is characterized by increasingly complex instruments and observational techniques, higher data collection rates, and large data archives, placing severe stress on software analysis systems. The object-oriented paradigm represents a significant new approach to software design and implementation that holds great promise for dealing with this increased complexity. The basic concepts of this approach will be characterized in contrast to more traditional procedure-oriented approaches. The fundamental features of objected-oriented programming will be discussed from a C++ programming language perspective, using examples familiar to astronomers. This discussion will focus on objects, classes and their relevance to the data type system; the principle of information hiding; and the use of inheritance to implement generalization/specialization relationships. Drawing on the object-oriented approach, features of a new database model to support astronomical data analysis will be presented.

Farris, A.↗
DOE OSTI2020

SN 2018fif: The Explosion of a Large Red Supergiant Discovered in Its Infancy by the Zwicky Transient Facility

High-cadence transient surveys are able to capture supernovae closer to their first light than ever before. Applying analytical models to such early emission, we can constrain the progenitor stars' properties. Here, we present observations of SN 2018fif (ZTF 18abokyfk). The supernova was discovered close to first light and monitored by the Zwicky Transient Facility (ZTF) and the Neil Gehrels Swift Observatory. Early spectroscopic observations suggest that the progenitor of SN 2018fif was surrounded by relatively small amounts of circumstellar material compared to all previous cases. This particularity, coupled with the high-cadence multiple-band coverage, makes it a good candidate to investigate using shock-cooling models. We employ the SOPRANOS code, an implementation of the model by Sapir & Waxman and its extension to early times by Morag et al. Compared with previous implementations, SOPRANOS has the advantage of including a careful account of the limited temporal validity domain of the shock-cooling model as well as allowing usage of the entirety of the early UV data. We find that the progenitor of SN 2018fif was a large red supergiant with a radius of $R={744.0}_{-128.0}^{+183.0}\,{R}_{\odot }$ and an ejected mass of ${M}_{\mathrm{ej}}={9.3}_{-5.8}^{+0.4}\,{M}_{\odot }$. Our model also gives information on the explosion epoch, the progenitor's inner structure, the shock velocity, and the extinction. The distribution of radii is double-peaked, with smaller radii corresponding to lower values of the extinction, earlier recombination times, and a better match to the early UV data. If these correlations persist in future objects, denser spectroscopic monitoring constraining the time of recombination, as well as accurate UV observations (e.g., with ULTRASAT), will help break the extinction/radius degeneracy and independently determine both.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2020

SN 2018fif: The Explosion of a Large Red Supergiant Discovered in its Infancy by the Zwicky Transient Facility

High-cadence transient surveys are able to capture supernovae closer to their first light than ever before. Applying analytical models to such early emission, we can constrain the progenitor stars' properties. In this paper, we present observations of SN 2018fif (ZTF 18abokyfk). The supernova was discovered close to first light and monitored by the Zwicky Transient Facility (ZTF) and the Neil Gehrels Swift Observatory. Early spectroscopic observations suggest that the progenitor of SN 2018fif was surrounded by relatively small amounts of circumstellar material compared to all previous cases. This particularity, coupled with the high-cadence multiple-band coverage, makes it a good candidate to investigate using shock-cooling models. We employ the SOPRANOS code, an implementation of the model by Sapir & Waxman and its extension to early times by Morag et al. Compared with previous implementations, SOPRANOS has the advantage of including a careful account of the limited temporal validity domain of the shock-cooling model as well as allowing usage of the entirety of the early UV data. We find that the progenitor of SN 2018fif was a large red supergiant with a radius of R = 744.0 (sup +183.0) (sub -128.0) R⨀ and an ejected mass of M(sub ej) = 9.3 (sup +0.4)(sub -5.8) M⨀. Our model also gives information on the explosion epoch, the progenitor's inner structure, the shock velocity, and the extinction. The distribution of radii is double-peaked, with smaller radii corresponding to lower values of the extinction, earlier recombination times, and a better match to the early UV data. If these correlations persist in future objects, denser spectroscopic monitoring constraining the time of recombination, as well as accurate UV observations (e.g., with ULTRASAT), will help break the extinction/radius degeneracy and independently determine both.

Supernovae↗
DOE OSTI2022

Far-ultraviolet to FIR Spectral-energy Distribution Modeling of the Stellar Formation History of the M31 Bulge

M31 is being surveyed at far- and near-ultraviolet with the UVIT telescope on AstroSat. The central bulge of M31 was observed in the N279N (275–280 nm), N219M (200–240 nm), F172M (160–185 nm), F169M (145–175 nm), and F148W (120–180 nm) filters. These images are made publicly available here. The UVIT data are supplemented with Sloan Digital Sky Survey data in optical, Spitzer data in near-infrared and Herschel data in mid- and far-infrared. The resulting far-ultraviolet to IR spectral-energy distributions for the bulge and for 10 subregions, are modeled using combinations of simple stellar populations and with CIGALE models. We find a dominant old (10–12 Gyr) metal-rich ([Z/H] ~ 0.3) population and a younger (600 Myr) solar abundance ([Z/H] ~ 0) population throughout the bulge. For the innermost 120" we find an additional very young (25 Myr) metal-poor ([Z/H] ~ –0.7) population. The results are consistent with the most recent stellar population studies of the bulge, which find the two populations for the whole bulge and a third young population in the innermost bulge.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

A Reanalysis of Public Galactic Bulge Gravitational Microlensing Events from OGLE-III and -IV

Abstract Modern surveys of gravitational microlensing events have progressed to detecting thousands per year, and surveys are capable of probing Galactic structure, stellar evolution, lens populations, black hole physics, and the nature of dark matter. One of the key avenues for doing this is the microlensing Einstein radius crossing time (t E) distribution. However, systematics in individual light curves as well as oversimplistic modeling can lead to biased results. To address this, we developed a model to simultaneously handle the microlensing parallax due to Earth's motion, systematic instrumental effects, and unlensed stellar variability with a Gaussian process model. We used light curves for nearly 10,000 OGLE-III and -IV Milky Way bulge microlensing events and fit each with our model. We also developed a forward model approach to infer the t E distribution by forward modeling from the data rather than using point estimates from individual events. We find that modeling the variability in the baseline removes a source of significant bias in individual events, and the previous analyses overestimated the number of t E > 100 day events due to their oversimplistic model ignoring parallax effects. We use our fits to identify the hundreds filling a regime in the microlensing parameter space that are 50% pure of black holes. Finally, we have released the largest-ever catalog of Markov Chain Monte Carlo parameter estimates for microlensing events.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Gravitational Microlensing Event Statistics for the Zwicky Transient Facility

Microlensing surveys have discovered thousands of events, with almost all events discovered within the Galactic bulge or toward the Magellanic Clouds. The Zwicky Transient Facility (ZTF), while not designed to be a microlensing campaign, is an optical time-domain survey that observes the entire northern sky, including the Galactic plane, every few nights. The ZTF observes ~10 9 stars in the g and r bands and can significantly contribute to the observed microlensing population. In this paper, we predict that the ZTF will observe ~1100 microlensing events in 3 yr of observing within 10° latitude of the Galactic plane, with ~500 events in the outer Galaxy (ℓ ≥ 10°). This yield increases to ~1400 (~800) events by combining every three ZTF exposures, ~1800 (~900) events if the ZTF observes for a total of 5 yr, and ~2400 (~1300) events for a 5 yr survey with postprocessing image stacking. Using the microlensing modeling software PopSyCLE, we compare the microlensing populations in the Galactic bulge and the outer Galaxy. We also present an analysis of the microlensing event ZTF18abhxjmj to demonstrate how to leverage these population statistics in event modeling. The ZTF will constrain Galactic structure, stellar populations, and primordial black holes through photometric microlensing.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Methods for Incorporating Model Uncertainty into Exoplanet Atmospheric Analysis

A key goal of exoplanet spectroscopy is to measure atmospheric properties, such as abundances of chemical species, in order to connect them to our understanding of atmospheric physics and planet formation. In this new era of high-quality JWST data, it is paramount that these measurement methods are robust. When comparing atmospheric models to observations, multiple candidate models may produce reasonable fits to the data. Typically, conclusions are reached by selecting the best-performing model according to some metric. This ignores model uncertainty in favor of specific model assumptions, potentially leading to measured atmospheric properties that are overconfident and/or incorrect. In this paper, we compare three ensemble methods for addressing model uncertainty by combining posterior distributions from multiple analyses: Bayesian model averaging, a variant of Bayesian model averaging using leave-one-out predictive densities, and stacking of predictive distributions. We demonstrate these methods by fitting the Hubble Space Telescope (HST) + Spitzer transmission spectrum of the hot Jupiter HD 209458b using models with different cloud and haze prescriptions. All of our ensemble methods lead to uncertainties on retrieved parameters that are larger but more realistic and consistent with physical and chemical expectations. Since they have not typically accounted for model uncertainty, uncertainties of retrieved parameters from HST spectra have likely been underreported. We recommend stacking as the most robust model combination method. Our methods can be used to combine results from independent retrieval codes and from different models within one code. They are also widely applicable to other exoplanet analysis processes, such as combining results from different data reductions.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2023

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↗
DOE OSTI2022

From the Fire: A Deeper Look at the Phoenix Stream

We use 6 yr of data from the Dark Energy Survey to perform a detailed photometric characterization of the Phoenix stellar stream, a 15° long, thin, dynamically cold, low-metallicity stellar system in the Southern Hemisphere. We use natural splines, a nonparametric modeling technique, to simultaneously fit the stream track, width, and linear density. This updated stream model allows us to improve measurements of the heliocentric distance (17.4 ± 0.1 (stat.) ±0.8 (sys.) kpc) and distance gradient (−0.009 ± 0.006 kpc deg$^{−1}$) of Phoenix, which corresponds to a small change of 0.13 ± 0.09 kpc in heliocentric distance along the length of the stream. We measure linear intensity variations on degree scales, as well as deviations in the stream track on ∼2° scales, suggesting that the stream may have been disturbed during its formation and/or evolution. We recover three peaks and one gap in linear intensity along with fluctuations in the stream track. Compared to other thin streams, the Phoenix stream shows more fluctuations and, consequently, the study of Phoenix offers a unique perspective on gravitational perturbations of stellar streams. We discuss possible sources of perturbations to Phoenix, including baryonic structures in the Galaxy and dark matter subhalos.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2022

Measurement of Galactic 26Al with the Compton Spectrometer and Imager

The Compton Spectrometer and Imager (COSI) is a balloon-borne compact Compton telescope de-17signed to survey the 0.2–5 MeV sky. COSI’s energy resolution of ∼0.2% at 1.8 MeV, single-photon reconstruction, and wide field of view make it capable of studying astrophysical nuclear lines, particularly the 1809 keVγ-ray line from decaying Galactic 26Al. Most 26Al originates in massive stars and core-collapse supernova nucleosynthesis, but the path from stellar evolution models to Galaxy-wide emission remains unconstrained. In 2016, COSI had a successful 46-day flight on a NASA super pressure balloon. Here, we detail the first search for the 1809 keV26Al line in the COSI 2016 balloon flight using a maximum likelihood analysis. We find a Galactic 26Al flux of (8.6±2.5) ×10−4ph cm−2s−124 within the Inner Galaxy (|ℓ|≤30◦,|b|≤10◦) with 3.7σ significance above background. Within uncertainties, this flux is consistent with expectations from previous measurements by SPI and COMPTEL. This analysis demonstrates COSI’s powerful capabilities for studies of γ-ray lines and underscores the scientific potential of future compact Compton telescopes. In particular, the next iteration of COSI as a NASA Small Explorer satellite has recently been approved for launch in 2025.

Gamma-ray lines↗
DOE OSTI2021

First M87 Event Horizon Telescope Results. VII. Polarization of the Ring

In 2017 April, the Event Horizon Telescope (EHT) observed the near-horizon region around the supermassive black hole at the core of the M87 galaxy. These 1.3 mm wavelength observations revealed a compact asymmetric ring-like source morphology. This structure originates from synchrotron emission produced by relativistic plasma located in the immediate vicinity of the black hole. Here we present the corresponding linear-polarimetric EHT images of the center of M87. We find that only a part of the ring is significantly polarized. The resolved fractional linear polarization has a maximum located in the southwest part of the ring, where it rises to the level of ~15%. The polarization position angles are arranged in a nearly azimuthal pattern. We perform quantitative measurements of relevant polarimetric properties of the compact emission and find evidence for the temporal evolution of the polarized source structure over one week of EHT observations. The details of the polarimetric data reduction and calibration methodology are provided. We carry out the data analysis using multiple independent imaging and modeling techniques, each of which is validated against a suite of synthetic data sets. The gross polarimetric structure and its apparent evolution with time are insensitive to the method used to reconstruct the image. These polarimetric images carry information about the structure of the magnetic fields responsible for the synchrotron emission. Their physical interpretation is discussed in an accompanying publication.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Building Trustworthy Machine Learning Models for Astronomy

Astronomy is entering an era of data-driven discovery, due in part to modern machine learning (ML) techniques enabling powerful new ways to interpret observations. This shift in our scientific approach requires us to consider whether we can trust the black box. Here, we overview methods for an often-overlooked step in the development of ML models: building community trust in the algorithms. Trust is an essential ingredient not just for creating more robust data analysis techniques, but also for building confidence within the astronomy community to embrace machine learning methods and results.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS2017

Update on Bayesian Blocks: Segmented Models for Sequential Data

The Bayesian Block algorithm, in wide use in astronomy and other areas, has been improved in several ways. The model for block shape has been generalized to include other than constant signal rate - e.g., linear, exponential, or other parametric models. In addition the computational efficiency has been improved, so that instead of O(N**2) the basic algorithm is O(N) in most cases. Other improvements in the theory and application of segmented representations will be described.

Scargle, Jeff↗
DOE OSTI2023

Towards an astronomical foundation model for stars with a transformer-based model

ABSTRACT Rapid strides are currently being made in the field of artificial intelligence using transformer-based models like Large Language Models (LLMs). The potential of these methods for creating a single, large, versatile model in astronomy has not yet been explored. In this work, we propose a framework for data-driven astronomy that uses the same core techniques and architecture as used by LLMs. Using a variety of observations and labels of stars as an example, we build a transformer-based model and train it in a self-supervised manner with cross-survey data sets to perform a variety of inference tasks. In particular, we demonstrate that a single model can perform both discriminative and generative tasks even if the model was not trained or fine-tuned to do any specific task. For example, on the discriminative task of deriving stellar parameters from Gaia XP spectra, we achieve an accuracy of 47 K in Teff, 0.11 dex in log g, and 0.07 dex in [M/H], outperforming an expert XGBoost model in the same setting. But the same model can also generate XP spectra from stellar parameters, inpaint unobserved spectral regions, extract empirical stellar loci, and even determine the interstellar extinction curve. Our framework demonstrates that building and training a single foundation model without fine-tuning using data and parameters from multiple surveys to predict unmeasured observations and parameters is well within reach. Such ‘Large Astronomy Models’ trained on large quantities of observational data will play a large role in the analysis of current and future large surveys.

Leung, Henry W. (ORCID:0000000200362752)↗
NASA NTRS1988

Polar cap emission model of Uranian kilometric radiation

Theoretical radiation patterns of Uranian kilometric radiation, consisting of a planet-wide hollow emission cone whose axis of symmetry is along the dipole axis, have been constructed for several frequencies. A comparison of model results with planetary radio astronomy data shows that the pulse profile is influenced by the geometry of the spacecraft relative to the magnetic field. The introduction of a small longitudinal variation of the emitting L-shell is found to eliminate the 10-deg phase discrepancy noted between theory and observations.

Barbosa, D. D.↗
NASA NTRS2007

Status on Iterative Transform Phase Retrieval Applied to the GBT Data

This slide presentation reviews the use of iterative transform phase retrieval in the analysis of the Green Bank Radio Telescope (GBT) Data. It reviews the NASA projects that have used phase retrieval, and the testbed for the algorithm to be used for the James Webb Space Telescope. It shows the comparison of phase retrieval with an interferometer, and reviews the two approaches used for phase retrieval, iterative transform (ITA) or parametric (non-linear least squares model fitting). The concept of ITA Phase Retrieval is reviewed, and the application to Radio Antennas is reviewed. The presentation also examines the National Radio Astronomy Observatory (NRAO) data from the GBT, and the Fourier model that NRAO uses to analyze the data. The challenge for ITA phase retrieval is reviewed, and the coherent approximation for incoherent data is shown. The validity of the approximation is good for a large tilt. There is a review of the proof of concept of the Phase Review simulation using the input wavefront, and the initial sampling parameters estimate from the focused GBT data.

Dean, Bruce↗
NASA NTRS2000

Chandra and XMM-Newton: Atomic Data Needs For X-ray Astronomy

With the launches of the Chandra X-ray Observatory and XM M-Newton, high resolution X-ray spectra of cosmic sources are broadening our understanding of the physical conditions, such as temperature, density, Ionization state, and elemental abundances. X-ray emitting astrophysical plasmas can be generally classified by their dominant ionization mechanism, either collisional ionization or X-ray photoionization. The atomic data needs are significantly different for these two cases; however, for both cases it is important that we identify robust and accurate diagnostics and that we verify completeness of the broadband models. We discuss the status of tile atomic data currently used in atomic databases for X-ray astronomy, in view of theoretical and experimental atomic physics considerations. We will also discuss the application of these models to new astrophysical data.

Brickhouse, Nancy↗