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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

A spatially resolved spectral analysis of giant radio galaxies with MeerKAT

ABSTRACT In this study we report the spatially resolved, wideband spectral properties of three giant radio galaxies (GRGs) in the COSMOS field: MGTC J095959.63+024608.6, MGTC J100016.84+015133.0, and MGTC J100022.85+031520.4. One of these galaxies, MGTC J100022.85+031520.4, is reported here for the first time, with a projected linear size of 1.29 Mpc at a redshift of 0.1034. Unlike the other two, it is associated with a brightest cluster galaxy (BCG), making it one of the few GRGs known to inhabit cluster environments. We examine the spectral age distributions of the three GRGs using new MeerKAT UHF-band (544–1088 MHz) observations, and L-band (900–1670 MHz) data from the MeerKAT International GHz Tiered Extragalactic Exploration (MIGHTEE) survey. We test two models of spectral ageing, the Jaffe–Perola and Tribble models, using the Broadband Radio Astronomy Tools (brats) software, and find that they agree well with each other. We estimate the Tribble spectral age for MGTC J095959.63+024608.6 as 68 Myr, for MGTC J100016.84+015133.0 as 47 Myr, and for MGTC J100022.85+031520.4 as 67 Myr. We find significant disagreements between these spectral age estimates and the estimates of the dynamical ages of these GRGs, modelled in cluster and group environments. Our results highlight the need for additional processes that are not accounted for in either the dynamic age or the spectral age estimations.

Charlton, K. K. L. (ORCID:0000000229252047)↗
NASA NTRS1989

Dust continuum spectra from model HII regions

The infrared spectrum emitted by nebular dust, heated by the ionizing stars in H II blisters and spherical H II regions, is calculated for various model parameters. Absorption of the non-ionizing radiation in a neutral layer is included. Heating by the Lyman alpha photon field is taken into account. The dust is composed of silicate and graphite grains, and evaporation of the grains in the inner region is considered. The models are presented with a view to interpretation of infrared observations of dusty H II regions and can be applied directly to the infrared astronomy satellite survey data. The continuum emission is compared with calculated fine structure line emission.

Aannestad, P. A.↗
DOE OSTI2024

Neural Network Prediction of Strong Lensing Systems with Domain Adaptation and Uncertainty Quantification

Modeling strong gravitational lenses is computationally expensive for the complex data from modern and next-generation cosmic surveys. Deep learning has emerged as a promising approach for finding lenses and predicting lensing parameters, such as the Einstein radius. Mean-variance Estimators (MVEs) are a common approach for obtaining aleatoric (data) uncertainties from a neural network prediction. However, neural networks have not been demonstrated to perform well on out-of-domain target data successfully - e.g., when trained on simulated data and applied to real, observational data. In this work, we perform the first study of the efficacy of MVEs in combination with unsupervised domain adaptation (UDA) on strong lensing data. The source domain data is noiseless, and the target domain data has noise mimicking modern cosmology surveys. We find that adding UDA to MVE increases the accuracy on the target data by a factor of about two over an MVE model without UDA. Including UDA also permits much more well-calibrated aleatoric uncertainty predictions. Advancements in this approach may enable future applications of MVE models to real observational data.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

A Compound Poisson Generator Approach to Point-source Inference in Astrophysics

Abstract The identification and description of point sources is one of the oldest problems in astronomy, yet even today the correct statistical treatment for point sources remains one of the field’s hardest problems. For dim or crowded sources, likelihood-based inference methods are required to estimate the uncertainty on the characteristics of the source population. In this work, a new parametric likelihood is constructed for this problem using compound Poisson generator (CPG) functionals that incorporate instrumental effects from first principles. We demonstrate that the CPG approach exhibits a number of advantages over non-Poissonian template fitting (NPTF)—an existing method—in a series of test scenarios in the context of X-ray astronomy. These demonstrations show that the effect of the point-spread function, effective area, and choice of point-source spatial distribution cannot, generally, be factorized as they are in NPTF, while the new CPG construction is validated in these scenarios. Separately, an examination of the diffuse-flux emission limit is used to show that most simple choices of priors on the standard parameterization of the population model can result in unexpected biases: when a model comprising both a point-source population and diffuse component is applied to this limit, nearly all observed flux will be assigned to either the population or to the diffuse component. A new parameterization is presented for these priors that properly estimates the uncertainties in this limit. In this choice of priors, CPG correctly identifies that the fraction of flux assigned to the population model cannot be constrained by the data.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1990

Observations of the galactic plane by the zodiacal infrared project

The two rocket flights of the Zodiacal Infrared Project (ZIP), flown 18 August 1980 and 31 July 1981, were intended to provide data on the near-infrared thermal emission of the interplanetary dust cloud over a broad range of ecliptic coordinates (latitudes -60 to +85 degrees, solar elongation angles 22 to 90 degrees and 140 to 180 degrees). In addition, their multiple crossings of the Galactic plane provided low resolution spectral data (delta lambda/lambda ranging from 1. to 0.1, for effective wavelengths from 3 to 30 microns) for most of the first quadrant (longitudes 30 to 100 degrees). Examples are displayed. Having made a thorough reanalysis of the calibration of the ZIP database, researchers present the salient features of the Galactic plane as observed by ZIP. The binned, in-plane data, corrected for zodiacal emission, generally show an exponential decrease with increasing longitude. The fitted exponential scale-length is 0.038/degree, and can be inverted to derive a radial density profile. Channel ratios are converted to temperatures by using model spectra in which thermal emitters with emissivity approx. 1/lambda are convolved with the filter responses. The results for channels 5 (11 microns) and 12 (21 microns) are shown, along with similarly derived temperatures from Infrared Astronomy Satellite (IRAS) 12 microns and 25 microns data. The ZIP data show little variation with longitude, consistent with IRAS results. A narrow spectral feature at 13 microns appears consistently in data for the plane (uncorrected for zodiacal emission). However, this is strongly contaminated by calibration problems for channel 8. Researchers suggest that residual emission at 13 microns arises from the (NeII) line at 12.8 microns.

Rickard, L. J.↗
DOE OSTI2021

Singling out modified gravity parameters and data sets reveals a dichotomy between Planck and lensing

ABSTRACT An important route to testing general relativity (GR) at cosmological scales is usually done by constraining modified gravity (MG) parameters added to the Einstein perturbed equations. Most studies have analysed so far constraints on pairs of MG parameters, but here, we explore constraints on one parameter at a time while fixing the other at its GR value. This allows us to analyse various models while benefiting from a stronger constraining power from the data. We also explore which specific data sets are in tension with GR. We find that models with (μ = 1, η) and (μ, η = 1) exhibit a 3.9σ and 3.8σ departure from GR when using Planck18 + Supernovae type Ia (SNe) + Baryon acoustic oscillation (BAO), while (μ, η) shows a tension of 3.4σ. We find no tension with GR for models with the MG parameter Σ fixed to its GR value. Using a Bayesian model selection analysis, we find that some one-parameter MG models are moderately favoured over ΛCDM when using all data set combinations except Planck cosmic microwave background lensing and dark energy survey data. Namely, Planck18 shows a moderate tension with GR that only increases when adding any combination of redshift space distortion, SNe, or BAO. However, adding lensing diminishes or removes these tensions, which can be attributed to the ability of lensing in constraining the MG parameter Σ. The two overall groups of data sets are found to have a dichotomy when performing consistency tests with GR, which may be due to systematic effects, lack of constraining power, or modelling. These findings warrant further investigation using more precise data from ongoing and future surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Sapsan: Framework for Supernovae Turbulence Modeling with Machine Learning

Sapsan is a framework designed to make Machine Learning (ML) more accessible in the study of turbulence, with a focus on astrophysical applications. Sapsan includes modules to load, filter, subsample, batch, and split the data from hydrodynamic (HD) simulations for training and validation. Next, the framework includes built-in conventional and physically-motivated estimators that have been used for turbulence modeling. This ties into Sapsan’s custom estimator module, aimed at designing a custom ML model layer-by-layer, which is the core benefit of using the framework. To share your custom model, every new project created via Sapsan comes with pre-filled, ready-for-release Docker files. Furthermore, training and evaluation modules come with Sapsan as well. The latter, among other features, includes the construction of power spectra and comparison to established analytical turbulence closure models, such as a gradient model. Thus, Sapsan attempts to minimize the hard work required for data preparation and analysis, leaving one to focus on the ML model design itself.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Analyzing Black Hole X-Rays with Simulations

We generate synthetic X-ray spectra of black holes from large-scale 3D GRMHD simulations performed on HPC machines. Here are a sampling of model spectra, characterized by the parameters indicated in the diagram below. The dashed red lines show power law fits to the parts of the spectrum observed by the X-ray telescope NuSTAR. How well these model spectra compare to real data helps us to characterize the physical realism of the underlying simulation techniques.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1986

Modeling and investigative studies of Jovian low frequency emissions

Jovian decametric (DAM) and hectometric (HOM) emissions were first observed over the entire spectrum by the Voyager 1 and 2 flybys of the planet. They display unusual arc-like structures on frequency-versus-time spectrograms. Software for the modeling of the Jovian plasma and magnetic field environment was performed. In addition, an extensive library of programs was developed for the retrieval of Voyager Planetary Radio Astronomy (PRA) data in both the high and low frequency bands from new noise-free, recalibrated data tapes. This software allows the option of retrieving data sorted with respect to particular sub-Io longitudes. This has proven to be invaluable in the analyses of the data. Graphics routines were also developed to display the data on color spectrograms.

Menietti, J. D.↗
NASA NTRS1992

Fortran for the nineties

Fortran has largely enjoyed prominence for the past few decades as the computer programming language of choice for numerically intensive scientific, engineering, and process control applications. Fortran's well understood static language syntax has allowed resulting parsers and compiler optimizing technologies to often generate among the most efficient and fastest run-time executables, particularly on high-end scalar and vector supercomputers. Computing architectures and paradigms have changed considerably since the last ANSI/ISO Fortran release in 1978, and while FORTRAN 77 has more than survived, it's aged features provide only partial functionality for today's demanding computing environments. The simple block procedural languages have been necessarily evolving, or giving way, to specialized supercomputing, network resource, and object-oriented paradigms. To address these new computing demands, ANSI has worked for the last 12-years with three international public reviews to deliver Fortran 90. Fortran 90 has superseded and replaced ISO FORTRAN 77 internationally as the sole Fortran standard; while in the US, Fortran 90 is expected to be adopted as the ANSI standard this summer, coexisting with ANSI FORTRAN 77 until at least 1996. The development path and current state of Fortran will be briefly described highlighting the many new Fortran 90 syntactic and semantic additions which support (among others): free form source; array syntax; new control structures; modules and interfaces; pointers; derived data types; dynamic memory; enhanced I/O; operator overloading; data abstraction; user optional arguments; new intrinsics for array, bit manipulation, and system inquiry; and enhanced portability through better generic control of underlying system arithmetic models. Examples from dynamical astronomy, signal and image processing will attempt to illustrate Fortran 90's applicability to today's general scalar, vector, and parallel scientific and engineering requirements and object oriented programming paradigms. Time permitting, current work proceeding on the future development of Fortran 2000 and collateral standards will be introduced.

Himer, J. T.↗
DOE OSTI2023

Core States of Neutron Stars from Anatomizing Their Scaled Structure Equations

Given an Equation of State (EOS) for neutron star (NS) matter, there is a unique mass–radius sequence characterized by a maximum mass ${M}_{\mathrm{NS}}^{\max}$ at radius $R$ max . We first show analytically that the ${M}_{\mathrm{NS}}^{\max}$ and $R$ max scale linearly with two different combinations of the NS central pressure $P$ $c$ and energy density $ε$ $c$ , by dissecting perturbatively the dimensionless Tolman–Oppenheimer–Volkoff (TOV) equations governing NS internal variables. The scaling relations are then verified via 87 widely used and rather diverse phenomenological as well as 17 microscopic NS EOSs with/without considering hadron–quark phase transitions and hyperons, by solving numerically the original TOV equations. The EOS of the densest NS matter allowed before it collapses into a black hole is then obtained. Using the universal ${M}_{\mathrm{NS}}^{\max}$ and $R$ max scalings and Neutron Star Interior Composition Explorer and XMM-Newton mass–radius observational data for PSR J0740+6620, a very narrow constraining band on the NS central EOS is extracted directly from the data for the first time, without using any specific input EOS model.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1971

New astronomy space experiments with television scanning.

Application of data from the Celescope Catalog of Ultraviolet Observations to various problems of stellar and interstellar astrophysics. These include refinements in the incorporation of line blanketing in theoretical stellar-atmosphere models, variations in the law of interstellar extinction from one region of the sky to another, and selection and identification of various types of peculiar stars for further investigation. In addition, the data were analyzed to determine the photometric and astrometric accuracy of the Celescope equipment.

Davis, R. J.↗
DOE OSTI2025

The Importance of Being Adaptable: An Exploration of the Power and Limitations of Domain Adaptation for Simulation-Based Inference with Galaxy Clusters

The application of deep machine learning methods in astronomy has exploded in the last decade, with new models showing remarkably improved performance on benchmark tasks. Not nearly enough attention is given to understanding the models' robustness, especially when the test data are systematically different from the training data, or "out of domain." Domain shift poses a significant challenge for simulation-based inference, where models are trained on simulated data but applied to real observational data. In this paper, we explore domain shift and test domain adaptation methods for a specific scientific case: simulation-based inference for estimating galaxy cluster masses from X-ray profiles. We build datasets to mimic simulation-based inference: a training set from the Magneticum simulation, a scatter-augmented training set to capture uncertainties in scaling relations, and a test set derived from the IllustrisTNG simulation. We demonstrate that the Test Set is out of domain in subtle ways that would be difficult to detect without careful analysis. We apply three deep learning methods: a standard neural network (NN), a neural network trained on the scatter-augmented input catalogs, and a Deep Reconstruction-Regression Network (DRRN), a semi-supervised deep model engineered to address domain shift. Although the NN improves results by 17% in the Training Data, it performs 40% worse on the out-of-domain Test Set. Surprisingly, the Scatter-Augmented Neural Network (SANN) performs similarly. While the DRRN is successful in mapping the training and Test Data onto the same latent space, it consistently underperforms compared to a straightforward Yx scaling relation. These results serve as a warning that simulation-based inference must be handled with extreme care, as subtle differences between training simulations and observational data can lead to unforeseen biases creeping into the results.

Ntampaka, Michelle [Baltimore, Space Telescope Sci↗
DOE OSTI2023

When Spectral Modeling Meets Convolutional Networks: A Method for Discovering Reionization-era Lensed Quasars in Multiband Imaging Data

Over the last two decades, around 300 quasars have been discovered at z ≳ 6, yet only one has been identified as being strongly gravitationally lensed. We explore a new approach—enlarging the permitted spectral parameter space, while introducing a new spatial geometry veto criterion—which is implemented via image-based deep learning. We first apply this approach to a systematic search for reionization-era lensed quasars, using data from the Dark Energy Survey, the Visible and Infrared Survey Telescope for Astronomy Hemisphere Survey, and the Wide-field Infrared Survey Explorer. Our search method consists of two main parts: (i) the preselection of the candidates, based on their spectral energy distributions (SEDs), using catalog-level photometry; and (ii) relative probability calculations of the candidates being a lens or some contaminant, utilizing a convolutional neural network (CNN) classification. The training data sets are constructed by painting deflected point-source lights over actual galaxy images, to generate realistic galaxy–quasar lens models, optimized to find systems with small image separations, i.e., Einstein radii of θ E ≤ 1''. Visual inspection is then performed for sources with CNN scores of P lens > 0.1, which leads us to obtain 36 newly selected lens candidates, which are awaiting spectroscopic confirmation. These findings show that automated SED modeling and deep learning pipelines, supported by modest human input, are a promising route for detecting strong lenses from large catalogs, which can overcome the veto limitations of primarily dropout-based SED selection approaches.

High-redshift galaxies↗
DOE OSTI2023

The star-formation history in the last 10 billion years from CIB cross-correlations

ABSTRACT The cosmic infrared background (CIB) traces the emission of star-forming galaxies throughout all cosmic epochs. Breaking down the contribution from galaxies at different redshifts to the observed CIB maps would allow us to probe the history of star formation. In this paper, we cross-correlate maps of the CIB with galaxy samples covering the range z ≲ 2 to measure the bias-weighted star-formation rate (SFR) density 〈bρSFR〉 as a function of time in a model independent way. This quantity is complementary to direct measurements of the SFR density ρSFR, giving a higher weight to more massive haloes, and thus provides additional information to constrain the physical properties of star formation. Using cross-correlations of the CIB with galaxies from the DESI Legacy Survey and the extended Baryon Oscillation Spectroscopic Survey, we obtain high signal-to-noise ratio measurements of 〈bρSFR〉, which we then use to place constraints on halo-based models of the star-formation history. We fit halo-based SFR models to our data and compare the recovered ρSFR with direct measurements of this quantity. We find a qualitatively good agreement between both independent data sets, although the details depend on the specific halo model assumed. This constitutes a useful robustness test for the physical interpretation of the CIB, and reinforces the role of CIB maps as valuable astrophysical probes of the large-scale structure. We report our measurements of 〈bρSFR〉 as well as a thorough account of their statistical uncertainties, which can be used to constrain star-formation models in combination with other data.

Astronomy & Astrophysics↗
DOE OSTI2024

Toward distinguishing Dirac from Majorana neutrino mass with gravitational waves

We propose a new method toward distinguishing the Dirac versus Majorana nature of neutrino masses from the spectrum of gravitational waves (GWs) associated with neutrino mass genesis. Motivated by the principle of generating small neutrino masses without tiny Yukawa couplings, we assume generic seesaw mechanisms for both Majorana and Dirac neutrino masses. For Majorana neutrinos, we further assume a spontaneously broken gauged U ( 1 ) B − L symmetry, independently of the type of Majorana seesaw mechanism, which gives a cosmic string induced GW signal flat over a wide range of frequencies. For Dirac neutrinos, we assume the spontaneous breaking of a Z 2 symmetry, the minimal symmetry choice associated with all Dirac seesaw mechanisms, which is softly broken, generating a peaked GW spectrum from the annihilation of the resulting domain walls. In fact, the GW spectra for all types of Dirac seesaws with such a broken Z 2 symmetry are identical, subject to a mild caveat. As an illustrative example, we study the simplest respective type-I seesaw mechanisms, and show that the striking difference in the shapes of the GW spectra can help differentiate between these Dirac and Majorana seesaws, complementing results of neutrinoless double beta decay experiments. We also discuss detailed implications of the recent NANOGrav data for Majorana and Dirac seesaw models. Published by the American Physical Society 2024

Astronomy & Astrophysics↗