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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2020

Studying the effects of overlapping objects in dark energy

Observing the clustering of galaxies allows us to calculate cosmological parameters necessary for understanding dark energy. However, as the density of observed objects increases, the probability of these objects blending likewise increases, causing multiple galaxies to be observed as one. This affects the inferred values of parameters such as the galaxy bias (b) and the matter energy density (ΩM). To see whether the bias from incorrectly inferring the galaxy count is significant, we compare the correlation functions in simulated data for “true” and “observed” data sets with one-to-one and multiple-to-one correspondences, respectively. For each data set, we create two correlation functions: one “measured” function directly relying on the galaxies’ positions using the TreeCorr python library, and one “model” derived mathematically from the galaxies’ power spectrum using the Cosmological Core Library (CCL). By minimizing the residual between these two functions, we compute the ideal values for b and ΩM across the various possible redshifts that position the galaxies in three dimensional space. This minimization is done with an Markov chain Monte Carlo (MCMC) estimate that finds one value of ΩM and ten values for b corresponding to the ten redshift bins ranging from z = 0.2 to z = 1.2. We find that neither b nor ΩM is particularly affected by inclusion of blended galaxies. Though there is room for improvement, the data suggests that the fluctuations we found are a result of noise or limitations on the modeling rather than blending explicitly.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1994

A reinterpretation of millimeter observations of nearby IRAS excess stars

We analyze new and previously published 1300, 870, and 800 micrometers, single-element bolometer observations of Vega (alpha Lyr), Formalhaut (alpha PsA), epsilon Eri, tau(sup1) Eri, and Beta Leo. We show that these data are consistent with models in which the dust disks around these stars are larger than the radio telescope beams with which they were observed; thus, these disks may be many hundreds of AU in radius or larger. Our interpretation of submillimeter/millimeter measurements of these stars also indicates that for some Infrared Astronomy Satellite (IRAS) IR excess stars the assumption that there is no astrophysical flux in the OFF beams, when using the standard ON-minus-OFF chopping technique and the chop distance is less than 1000AU, may be incorrect. Therefore, for IRAS IR excess stars within approximately 20 pc, virtually all submillimeter/millimeter chopping observations with chop throws less than 100" may have subtracted away some or most of the flux associated with their circumstellar disks. Finally, we present new 1300 micrometers continuum observations of Vega made with chop throws of 500 and 1000 AU. These data are consistent with an interpretation in which Vega has a disk that is at least 1000AU in radius; this disk could have a region with much less material per beam area near 500 AU than at both 100 AU and 1000 AU, corresponding to a gap at the orbital distance of alpha Lyr B. The observations of Vega are also consistent with the assumption that the circumVega dust structure is unresolved. Thus, these new data very effectively illustrate that the 'standard' model of small, unresolved, dust structures around main-sequence IRAS IR excess stars is not unique. Maps of IRAS IR excess stars, which soon will be available from submillimeter/millimeter array detectors, will determine whether the paradigm shift we propose will occur.

Weintraub, David A.↗
NASA NTRS1985

Analysis of Jovian decametric data: Study of radio emission mechanisms

The Voyager 1 and Voyager 2 Planetary Radio Astronomy Experiments (PRA) have produced the finest set of Jovian decametric radio emission data ever obtained. Jovian decametric L-burst and S-burst arcs were characterized and the data reconciled with models for the radio emission geometry and mechanisms. The first major results involve comparisons of the distribution of arc separations with longitudes. The identification and analyses of systematic variations in the PRA data have yielded interesting results, but only the most obvious features of the data were examined. Analyses of the PRA data were extended with the use of new 6-Sec formats that are more sensitive to the S-bursts.

Staelin, D. H.↗
NASA NTRS1985

The peculiar X-ray and radio star AS431

Einstein observations are reported which reveal that the highly reddened emission-line star AS431 has a highly absorbed X-ray spectrum and a relatively strong intrinsic flux of 5 x 10 to the -12th erg/sq cm/s or more. Observations at 20 cm and 6 cm with the National Radio Astronomy Observatory Very Large Array are also reported which show that AS431 is also a moderately strong radio source (about 35 mJy). These data, together with optical and infrared observations, suggest a model in which both the radio and X-ray emissions arise in a chaotic stellar wind emerging from a single luminous Wolf-Rayet star.

Caillault, J.-P.↗
DOE OSTI2021

Model-independent extraction of the proton charge radius from PRad data

The proton radius puzzle has motivated several new experiments that aim to extract the proton charge radius and resolve the puzzle. Recently, PRad, a new electron–proton scattering experiment at Jefferson Lab, reported a proton charge radius of [Formula: see text]. The value was obtained by using a rational function model for the proton electric form factor. We perform a model-independent extraction using [Formula: see text]-expansion of the proton charge radius from PRad data. We find that the model-independent statistical error is more than 50% larger compared to the statistical error reported by PRad.

Astronomy & Astrophysics↗
DOE OSTI2024

Domain-Adaptive Neural Posterior Estimation for Strong Gravitational Lens Analysis

Modeling strong gravitational lenses is prohibitively expensive for modern and next-generation cosmic survey data. Neural posterior estimation (NPE), a simulation-based inference (SBI) approach, has been studied as an avenue for efficient analysis of strong lensing data. However, NPE has not been demonstrated to perform well on out-of-domain target data -- e.g., when trained on simulated data and then applied to real, observational data. In this work, we perform the first study of the efficacy of NPE in combination with unsupervised domain adaptation (UDA). The source domain is noiseless, and the target domain has noise mimicking modern cosmology surveys. We find that combining UDA and NPE improves the accuracy of the inference by 1-2 orders of magnitude and significantly improves the posterior coverage over an NPE model without UDA. We anticipate that this combination of approaches will help enable future applications of NPE models to real observational data.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Efficient learning of accurate surrogates for simulations of complex systems

Machine learning methods are increasingly deployed to construct surrogate models for complex physical systems at a reduced computational cost. However, the predictive capability of these surrogates degrades in the presence of noisy, sparse or dynamic data. Here, we introduce an online learning method empowered by optimizer-driven sampling that has two advantages over current approaches: it ensures that all local extrema (including endpoints) of the model response surface are included in the training data, and it employs a continuous validation and update process in which surrogates undergo retraining when their performance falls below a validity threshold. We find, using benchmark functions, that optimizer-directed sampling generally outperforms traditional sampling methods in terms of accuracy around local extrema even when the scoring metric is biased towards assessing overall accuracy. Finally, the application to dense nuclear matter demonstrates that highly accurate surrogates for a nuclear equation-of-state model can be reliably autogenerated from expensive calculations using few model evaluations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Estimating the color lifetime of energetic quarks

Using a simple geometric framework with a realistic nuclear density distribution, we fit published HERMES data to determine fundamental properties of hadronization using the nuclear medium as a spatial analyzer. Our approach uses a fit to the transverse momentum broadening observable and the hadronic multiplicity ratio; the simultaneous fit to two different observables strongly constrains the outcome. Using the known sizes of the target nuclei, we extract the color lifetime, finding a z h -dependent range of values from 2 to 8 fm/c for these data. We also extract estimates for the qˆ transport coefficient characterizing the strength of the interaction between the quark and the cold nuclear medium, finding an average value of 0.035±0.011 GeV 2 /fm. With a three-parameter model we obtain satisfactory fits to the data with a goodness-of-fit parameter χ 2 /dof of 1.1 or less. In a secondary fit of the results from that model we independently find a value for the Lund String Model string tension of 1.00±0.05 GeV/fm. We evaluated the sensitivity for extracting quark energy loss and effective in-medium hadronic cross sections using four-parameter variants of the model, finding large uncertainties in both cases. Our results suggest that hadronic interaction of forming hadrons in the nuclear medium is the primary dynamical cause of meson attenuation in the HERMES data, with quark energy loss playing a more minor role.

Electron-ion collider↗
NASA NTRS1974

Calculations of the moon's heat history at different concentrations of radioactive elements taking account of the material differentiation with melting

A mathematical procedure for analyzing the heat conductivity of the lunar surface is discussed. The solution is based on homogeneous and laminated moon models and considers the effects of radioactive elements conveyed to the lunar surface by melting. The various parameters which introduce uncertainties into the numerical analysis are identified. The application of data obtained from radio astronomy and from analyses of lunar samples returned by the Apollo flights is explained. Tables of data are included to show the types and amounts of radioactive materials which have been identified.

Arnatskaya, O. I.↗
NASA NTRS1985

Statistical methods for astronomical data with upper limits. I - Univariate distributions

The statistical treatment of univariate censored data is discussed. A heuristic derivation of the Kaplan-Meier maximum-likelihood estimator from first principles is presented which results in an expression amenable to analytic error analysis. Methods for comparing two or more censored samples are given along with simple computational examples, stressing the fact that most astronomical problems involve upper limits while the standard mathematical methods require lower limits. The application of univariate survival analysis to six data sets in the recent astrophysical literature is described, and various aspects of the use of survival analysis in astronomy, such as the limitations of various two-sample tests and the role of parametric modelling, are discussed.

Feigelson, E. D.↗
NASA NTRS1988

Lunar and planetary studies

This grant supports the core program in planetary astronomy at Caltech. The research includes observations in the IR, sub-mm, mm and cm wavelengths at national and Caltech observatories with a strong emphasis on integrating the observations with spacecraft data and with models of atmospheric structure, dynamics and chemistry. Muhleman's group made extensive observations of Saturn, Uranus and Neptune which are being interpreted in terms of deep atmospheric structures which are obvious in the 2 and 6 cm maps of Saturn and Uranus. The microwave measurements are one of the few sources of information below the 2 bar level. Goldreich is investigating the dynamics of narrow rings with postdoctoral fellow, Pierre-Yves Longaretti. Their work has focused on the role of collisional stresses on the precession of the rings, since the Voyager radio science results imply that the previous model based on the ring's self-gravity is not the entire story. In addition Borderies, Goldreich and Tremaine have completed an investigation of the dynamics of the Encke division in Saturn's A ring.

Muhleman, Duane O.↗
DOE OSTI2021

Automatic Search of Cataclysmic Variables Based on LightGBM in LAMOST-DR7

The search for special and rare celestial objects has always played an important role in astronomy. Cataclysmic Variables (CVs) are special and rare binary systems with accretion disks. Most CVs are in the quiescent period, and their spectra have the emission lines of Balmer series, HeI, and HeII. A few CVs in the outburst period have the absorption lines of Balmer series. Owing to the scarcity of numbers, expanding the spectral data of CVs is of positive significance for studying the formation of accretion disks and the evolution of binary star system models. At present, the research for astronomical spectra has entered the era of Big Data. The Large Sky Area Multi-Object Fiber Spectroscopy Telescope (LAMOST) has produced more than tens of millions of spectral data. the latest released LAMOST-DR7 includes 10.6 million low-resolution spectral data in 4926 sky regions, providing ideal data support for searching CV candidates. To process and analyze the massive amounts of spectral data, this study employed the Light Gradient Boosting Machine (LightGBM) algorithm, which is based on the ensemble tree model to automatically conduct the search in LAMOST-DR7. Finally, 225 CV candidates were found and four new CV candidates were verified by SIMBAD and published catalogs. This study also built the Gradient Boosting Decision Tree (GBDT), Adaptive Boosting (AdaBoost), and eXtreme Gradient Boosting (XGBoost) models and used Accuracy, Precision, Recall, the F1-score, and the ROC curve to compare the four models comprehensively. Experimental results showed that LightGBM is more efficient. The search for CVs based on LightGBM not only enriches the existing CV spectral library, but also provides a reference for the data mining of other rare celestial objects in massive spectral data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Beyond Fisher forecasting for cosmology

The planning and design of future experiments rely heavily on forecasting to assess the potential scientific value provided by a hypothetical set of measurements. The Fisher information matrix, due to its convenient properties and low computational cost, provides an especially useful forecasting tool. However, the Fisher matrix only provides a reasonable approximation to the true likelihood when data are nearly Gaussian distributed and observables have nearly linear dependence on the parameters of interest. Also, Fisher forecasting techniques alone cannot be used to assess their own validity. Thorough sampling of the exact or mock likelihood can definitively determine whether a Fisher forecast is valid, though such sampling is often prohibitively expensive. Here we propose a simple test, based on the Derivative Approximation for likelihoods (DALI) technique, to determine whether the Fisher matrix provides a good approximation to the exact likelihood. We show that the Fisher matrix becomes a poor approximation to the true likelihood in regions where two-dimensional slices of level surfaces of the DALI approximation to the likelihood differ from two-dimensional slices of level surfaces of the Fisher approximation to the likelihood. We demonstrate that our method accurately predicts situations in which the Fisher approximation deviates from the true likelihood for various cosmological models and several data combinations, with only a modest increase in computational cost compared to standard Fisher forecasts.

79 ASTRONOMY AND ASTROPHYSICS↗