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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

The γ -process nucleosynthesis in core-collapse supernovae

The γ-process nucleosynthesis in core-collapse supernovae is generally accepted as a feasible process for the synthesis of neutron-deficient isotopes beyond iron. However, crucial discrepancies between theory and observations still exist: the average yields of γ-process nucleosynthesis from massive stars are still insufficient to reproduce the solar distribution in galactic chemical evolution calculations, and the yields of the Mo and Ru isotopes are a factor of ten lower than the yields of the other γ-process nuclei. We investigate the γ-process in five sets of core-collapse supernova models published in the literature with initial masses of 15, 20, and 25 M ⊙ at solar metallicity. We compared the γ-process overproduction factors from the different models. To highlight the possible effect of nuclear physics input, we also considered 23 ratios of two isotopes close to each other in mass relative to their solar values. Further, we investigated the contribution of C–O shell mergers in the supernova progenitors as an additional site of the γ-process. Our analysis shows that a large scatter among the different models exists for both the γ-process integrated yields and the isotopic ratios. We find only ten ratios that agree with their solar values, all the others differ by at least a factor of three from the solar values in all the considered sets of models. The γ-process within C–O shell mergers mostly influences the isotopic ratios that involve intermediate and heavy proton-rich isotopes with A > 100. We conclude that there are large discrepancies both among the different data sets and between the model predictions and the solar abundance distribution. More calculations are needed; particularly updating the nuclear network, because the majority of the models considered in this work do not use the latest reaction rates for the γ-process nucleosynthesis. Moreover, the role of C–O shell mergers requires further investigation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

H- and Dissociation in Ultra-hot Jupiters: A Retrieval Case Study of WASP-18b

Atmospheres of a number of ultra-hot Jupiters (UHJs) with temperatures ≳2000 K have been observed recently. Many of these planets show largely featureless thermal spectra in the near-infrared observed with the HST WFC3 spectrograph (1.1–1.7 μm) even though this spectral range contains strong H{sub 2}O opacity. Recent works have proposed the possibility of H- opacity masking the H{sub 2}O feature and/or thermal dissociation of H{sub 2}O causing its apparent depletion at the high temperatures of UHJs. In this work, we test these hypotheses using observations of the exoplanet WASP-18b as a case study. We report detailed atmospheric retrievals of the planet using the HyDRA retrieval code, extended to include the effects of H- opacity and thermal dissociation. We report constraints on the H{sub 2}O, CO, and H- abundances as well as the pressure–temperature profile of the dayside atmosphere for retrievals with and without H-/dissociation for each data set. We find that the H{sub 2}O and H- abundances are relatively unconstrained given the featureless WFC3 spectra. We do not conclusively detect H- in the planet, contrary to previous studies that used equilibrium models to infer its presence. The constraint on the CO abundance depends on the combination of WFC3 and Spitzer data, ranging from solar to super-solar CO values. We additionally see signs of a thermal inversion from two of the data sets. Our study demonstrates the potential of atmospheric retrievals of UHJs, including the effects of H- and thermal dissociation of molecules.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Domain Adaptation for Measurements of Strong Gravitational Lenses

Upcoming surveys are predicted to discover galaxy-scale strong lenses on the order of $10^5$, making deep learning methods necessary in lensing data analysis. Currently, there is insufficient real lensing data to train deep learning algorithms, but the alternative of training only on simulated data results in poor performance on real data. Domain Adaptation may be able to bridge the gap between simulated and real datasets. We utilize domain adaptation for the estimation of Einstein radius ($\Theta_E$) in simulated galaxy-scale gravitational lensing images with different levels of observational realism. We evaluate two domain adaptation techniques - Domain Adversarial Neural Networks (DANN) and Maximum Mean Discrepancy (MMD). We train on a source domain of simulated lenses and apply it to a target domain of lenses simulated to emulate noise conditions in the Dark Energy Survey (DES). We show that both domain adaptation techniques can significantly improve the model performance on the more complex target domain dataset. This work is the first application of domain adaptation for a regression task in strong lensing imaging analysis. Our results show the potential of using domain adaptation to perform analysis of future survey data with a deep neural network trained on simulated data.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

A Conditional Autoencoder for Galaxy Photometric Parameter Estimation

Astronomical photometric surveys routinely image billions of galaxies, and traditionally infer the parameters of a parametric model for each galaxy. This approach has served us well, but the computational expense of deriving a full posterior probability distribution function is a challenge for increasingly ambitious surveys. In this paper, we use deep learning methods to characterize galaxy images, training a conditional autoencoder on mock data. The autoencoder can reconstruct and denoise galaxy images via a latent space engineered to include semantically meaningful parameters, such as brightness, location, size, and shape. Our model recovers galaxy fluxes and shapes on mock data with a lower variance than the Hyper Suprime-Cam photometry pipeline, and returns reasonable answers even for inputs outside the range of its training data. When applied to data in the training range, the regression errors on all extracted parameters are nearly unbiased with a variance near the Cramr-Rao bound.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Dynamical inflection point inflation

Abstract We provide a mechanism to obtain multiple inflection points for the inflaton potential based on the Coleman-Weinberg potential. The key ingredient is the existence of zeros of the inflaton quartic coupling beta function, which can be simply realized by a sign change of a gauge coupling beta function at a massive threshold scale in gauge-Yukawa models. A universal inflaton potential emerges for a wide range of models and can accommodate the observational data in the small-field inflation scenario. While the ratio of tensor over scalar amplitudes is predicted to be small, the running of the scalar spectral index, α s , has branches with both positive and negative values. The positive branch with α s ≳ 0.002 will be tested by the coming cosmic microwave background observations.

Astronomy & Astrophysics↗
DOE OSTI2021

NuSTAR Observations of Abell 2163: Constraints on Non-thermal Emission

Since the first non-thermal reports of inverse Compton (IC) emission from the intracluster medium (ICM) of galaxy clusters at hard X-ray energies, we have yet to unambiguously confirm IC emission in observations with newer facilities. RXTE detected IC emission in one of the hottest known clusters, Abell 2163 (A2163), a massive merging cluster with a giant radio halo—the presumed source of relativistic electrons IC scattering CMB photons to X-ray energies. The cluster's redshift (z ~ 0.2) allows its thermal and non-thermal radio emission to fit the field of view of the Nuclear Spectroscopic Telescope Array (NuSTAR), permitting a deep observation capable of confirming or ruling out the RXTE report. The IC flux provides constraints on the average magnetic field strength in a cluster. To determine the global diffuse IC emission in A2163, we fit its global NuSTAR spectrum with four models: single- (1T) and two-temperature (2T), 1T+power-law component (T+IC), and multi-temperature+power law (9T+IC). Each represents different characterizations of the thermal ICM emission, with power-law components added to represent IC emission. We find that the 3–30 keV spectrum can be described by purely thermal emission, with a global average temperature of kT = (11.8 ± 0.2) keV. The IC flux is constrained to <4.0 × 10 –12 erg s –1 cm –2 using the 1T+IC model and <1.6 × 10 –12 erg s –1 cm –2 with the more physical 9T+IC model, both to 90% confidence levels. Here, combining these limits with 1.4 GHz diffuse radio data from the VLA, we find the average magnetic field strength to be >0.22 μG and >0.35 μG, respectively, providing the strongest constraints on these values in A2163 to date.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Modeling Dust Production, Growth, and Destruction in Reionization-era Galaxies with the CROC Simulations: Methods and Parameter Exploration

We introduce a model for the explicit evolution of interstellar dust in a cosmological galaxy formation simulation. We post-process a simulation from the Cosmic Reionization on Computers project, integrating an ordinary differential equation for the evolution of the dust-to-gas ratio along pathlines in the simulation sampled with a tracer particle technique. This model incorporates the effects of dust grain production in asymptotic giant branch star winds and supernovae, grain growth due to the accretion of heavy elements from the gas phase of the interstellar medium, and grain destruction due to thermal sputtering in the high-temperature gas of supernova remnants due to thermal sputtering. A main conclusion of our analysis is the importance of a carefully chosen dust destruction model, for which different reasonable parameterizations can predict very different values at the ∼100 pc resolution of the interstellar medium (ISM) in our simulations. We run this dust model on the single most massive galaxy in a 10 h$^{−1}$ comoving megaparsec box, which attains a stellar mass of ∼2 × 10$^{9}$ M $_{⊙}$ by z = 5. We find that the model is capable of reproducing dust masses and dust-sensitive observable quantities broadly consistent with existing data from high-redshift galaxies. The total dust mass in the simulated galaxy is somewhat sensitive to parameter choices for the dust model, especially the timescale for grain growth due to accretion in the ISM. Consequently, observable quantity observations that can constrain galaxy dust masses at these epochs are potentially useful for placing constraints on dust physics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Rapid 3D nanoscale coherent imaging via physics-aware deep learning

Phase retrieval, the problem of recovering lost phase information from measured intensity alone, is an inverse problem that is widely faced in various imaging modalities ranging from astronomy to nanoscale imaging. The current process of phase recovery is iterative in nature. As a result, the image formation is time consuming and computationally expensive, precluding real-time imaging. Here, we use 3D nanoscale X-ray imaging as a representative example to develop a deep learning model to address this phase retrieval problem. We introduce 3D-CDI-NN, a deep convolutional neural network and differential programing framework trained to predict 3D structure and strain, solely from input 3D X-ray coherent scattering data. Our networks are designed to be “physics-aware” in multiple aspects; in that the physics of the X-ray scattering process is explicitly enforced in the training of the network, and the training data are drawn from atomistic simulations that are representative of the physics of the material. We further refine the neural network prediction through a physics-based optimization procedure to enable maximum accuracy at lowest computational cost. 3D-CDI-NN can invert a 3D coherent diffraction pattern to real-space structure and strain hundreds of times faster than traditional iterative phase retrieval methods. Our integrated machine learning and differential programing solution to the phase retrieval problem is broadly applicable across inverse problems in other application areas.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2023

Observational constraints on warm natural inflation

Warm natural inflation is studied for the case of the original cosine potential. The radiation bath during inflation induces a dissipation (friction) rate in the equation of motion for the inflation field, which can potentially reduce the field excursion needed for an observationally viable period of inflation. Here we examine if the dissipation thus provides a mechanism to avoid the large decay constant f ≳ M pl of cold cosine natural inflation. Whereas temperature independent dissipation has previously been shown to alleviate the need for a trans-Planckian decay constant f, we illustrate here the difficulties of accommodating a significantly sub-Planckian decay constant (f < 10 -1 M pl ) for the case of a temperature dependent dissipation rate in the form Γ ∝ T c , with c = {1,3}. Such dissipation rates represent physically well-motivated constructions in the literature. For each model, we map its location in the r-n s plane and compare with cosmic microwave background (CMB) data. For c = 1 (c = 3), we find that agreement with CMB data requires that dissipation be in the weak (moderate) regime and that the minimum allowed value of the decay constant in the potential is f min = 0.3 (0.8) M pl , respectively.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Modeling neutrino-induced scale-dependent galaxy clustering for photometric galaxy surveys

Abstract The increasing statistical precision of photometric redshift surveys requires improved accuracy of theoretical predictions for large-scale structure observables to obtain unbiased cosmological constraints. In ΛCDM cosmologies, massive neutrinos stream freely at small cosmological scales, suppressing the small-scale power spectrum. In massive neutrino cosmologies, galaxy bias modeling needs to accurately relate the scale-dependent growth of the underlying matter field to observed galaxy clustering statistics. In this work, we implement a computationally efficient approximation of the neutrino-induced scale-dependent bias (NISDB). Through simulated likelihood analyses of Dark Energy Survey Year 3 (DESY3) and Legacy Survey of Space and Time Year 1 (LSSTY1) synthetic data that contain an appreciable NISDB, we examine the impact of linear galaxy bias and neutrino mass modeling choices on cosmological parameter inference. We find model misspecification of the NISDB approximation and neutrino mass models to decrease the constraining power of photometric galaxy surveys and cause parameter biases in the cosmological interpretation of future surveys. We quantify these biases and devise mitigation strategies.

Astronomy & Astrophysics↗
DOE OSTI2022

Novel Rosenbluth extraction framework for Compton form factors from deeply virtual exclusive experiments

We use a generalization of the Rosenbluth separation method for a model independent simultaneous extraction of the Compton Form Factors $\mathscr{H}$ and $\mathscr{E}$ from virtual Compton scattering data on an unpolarized target. A precise evaluation of $\mathscr{H}$ and $\mathscr{E}$ enabled by the proposed method, is the first step towards pinning down the distribution of angular momentum inside the proton.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Neutrino-assisted early dark energy: theory and cosmology

Abstract The tension between measurements of the Hubble constant obtained at different redshifts may provide a hint of new physics active in the relatively early universe, around the epoch of matter-radiation equality. A leading paradigm to resolve the tension is a period of early dark energy, in which a scalar field contributes a subdominant part of the energy budget of the universe at this time. This scenario faces significant fine-tuning problems which can be ameliorated by a non-trivial coupling of the scalar to the standard model neutrinos. These become non-relativistic close to the time of matter-radiation equality, resulting in an energy injection into the scalar that kick-starts the early dark energy phase, explaining its coincidence with this seemingly unrelated epoch. We present a minimal version of this neutrino-assisted early dark energy model, and perform a detailed analysis of its predictions and theoretical constraints. We consider both particle physics constraints — that the model constitute a well-behaved effective field theory for which the quantum corrections are under control, so that the relevant predictions are within its regime of validity — and the constraints provided by requiring a consistent cosmological evolution from early through to late times. Our work paves the way for testing this scenario using cosmological data sets.

Astronomy & Astrophysics↗
DOE OSTI2023

DSPS: Differentiable stellar population synthesis

ABSTRACT Models of stellar population synthesis (SPS) are the fundamental tool that relates the physical properties of a galaxy to its spectral energy distribution (SED). In this paper, we present DSPS: a python package for SPS. All of the functionality in DSPS is implemented natively in the JAX library for automatic differentiation, and so our predictions for galaxy photometry are fully differentiable, and directly inherit the performance benefits of JAX, including portability onto GPUs. DSPS also implements several novel features, such as i) a flexible empirical model for stellar metallicity that incorporates correlations with stellar age, ii) support for the Diffstar model that provides a physically-motivated connection between the star formation history of a galaxy (SFH) and the mass assembly of its underlying dark matter halo. We detail a set of theoretical techniques for using autodiff to calculate gradients of predictions for galaxy SEDs with respect to SPS parameters that control a range of physical effects, including SFH, stellar metallicity, nebular emission, and dust attenuation. When forward modelling the colours of a synthetic galaxy population, we find that DSPS can provide a factor of 5 speed-up over standard SPS codes on a CPU, and a factor of 300-400 on a modern GPU. When coupled with gradient-based techniques for optimization and inference, DSPS makes it practical to conduct expansive likelihood analyses of simulation-based models of the galaxy–halo connection that fully forward model galaxy spectra and photometry.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Kepler Data Analysis: Non-Gaussian Noise and Fourier Gaussian Process Analysis of Stellar Variability

We develop a statistical analysis model of Kepler stellar flux data in the presence of planet transits, non-Gaussian noise, and stellar variability. We first develop a model for the Kepler noise probability distribution in the presence of outliers, which make the noise probability distribution non-Gaussian. We develop a signal likelihood analysis based on this probability distribution, in which we model the signal as a sum of the star variability and planetary transits. We argue that these components need to be modeled together if optimal signal is to be extracted from the data. For the stellar variability model we develop an optimal Gaussian process analysis using a Fourier-based Wiener filter approach, where the power spectrum is non-parametric and learned from the data. We develop high dimensional optimization of the objective function, where we jointly optimize all the model parameters, including thousands of star variability modes, and planet transit parameters. We apply the method to Kepler-90 data and show that it gives a better match to the stellar variability than the existing methods, and robustly handles noise outliers. As a consequence, the planet radii have a higher value than what the existing methods give, including splines and celerite.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

HOLISMOKES IX. Neural network inference of strong-lens parameters and uncertainties from ground-based images

Modeling of strong gravitational lenses is a necessity for further applications in astrophysics and cosmology. With the large number of detections in current and upcoming surveys, such as the Rubin Legacy Survey of Space and Time (LSST), it is pertinent to investigate automated and fast analysis techniques beyond the traditional and time-consuming Markov chain Monte Carlo sampling methods. Building upon our (simple) convolutional neural network (CNN), we present here another CNN, specifically a residual neural network (ResNet), that predicts the five mass parameters of a singular isothermal ellipsoid (SIE) profile (lens center x and y, ellipticity e x and e y , Einstein radius θ E ) and the external shear (γ ext, 1 , γ ext, 2 ) from ground-based imaging data. In contrast to our previous CNN, this ResNet further predicts the 1σ uncertainty for each parameter. To train our network, we use our improved pipeline to simulate lens images using real images of galaxies from the Hyper Suprime-Cam Survey (HSC) and from the Hubble Ultra Deep Field as lens galaxies and background sources, respectively. We find very good recoveries overall for the SIE parameters, especially for the lens center in comparison to our previous CNN, while significant differences remain in predicting the external shear. From our multiple tests, it appears that most likely the low ground-based image resolution is the limiting factor in predicting the external shear. Given the run time of milli-seconds per system, our network is perfectly suited to quickly predict the next appearing image and time delays of lensed transients. Therefore, we use the network-predicted mass model to estimate these quantities and compare to those values obtained from our simulations. Unfortunately, the achieved precision allows only a first-order estimate of time delays on real lens systems and requires further refinement through follow-up modeling. Nonetheless, our ResNet is able to predict the SIE and shear parameter values in fractions of a second on a single CPU, meaning that we are able to efficiently process the huge amount of galaxy-scale lenses expected in the near future.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Hierarchical Inference with Bayesian Neural Networks: An Application to Strong Gravitational Lensing

In the past few years, approximate Bayesian Neural Networks (BNNs) have demonstrated the ability to produce statistically consistent posteriors on a wide range of inference problems at unprecedented speed and scale. However, any disconnect between training sets and the distribution of real-world objects can introduce bias when BNNs are applied to data. This is a common challenge in astrophysics and cosmology, where the unknown distribution of objects in our universe is often the science goal. In this work, we incorporate BNNs with flexible posterior parameterizations into a hierarchical inference framework that allows for the reconstruction of population hyperparameters and removes the bias introduced by the training distribution. We focus on the challenge of producing posterior PDFs for strong gravitational lens mass model parameters given Hubble Space Telescope–quality single-filter, lens-subtracted, synthetic imaging data. We show that the posterior PDFs are sufficiently accurate (statistically consistent with the truth) across a wide variety of power-law elliptical lens mass distributions. We then apply our approach to test data sets whose lens parameters are drawn from distributions that are drastically different from the training set. We show that our hierarchical inference framework mitigates the bias introduced by an unrepresentative training set's interim prior. Simultaneously, we can precisely reconstruct the population hyperparameters governing our test distributions. Our full pipeline, from training to hierarchical inference on thousands of lenses, can be run in a day. The framework presented here will allow us to efficiently exploit the full constraining power of future ground- and space-based surveys (https://github.com/swagnercarena/ovejero).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Generation of realistic input parameters for simulating atmospheric point-spread functions at astronomical observatories

High-fidelity simulated astronomical images are an important tool in developing and measuring the performance of image-processing algorithms, particularly for high precision measurements of cosmic shear – correlated distortions of images of distant galaxies due to weak gravitational lensing caused by the large-scale mass distribution in the Universe. For unbiased measurements of cosmic shear, all other sources of correlated image distortions must be modeled or removed. One such source is the correlated blurring of images due to optical turbulence in the atmosphere, which dominates the point-spread function (PSF) for ground-based instruments. In this work, we leverage data from weather forecasting models to produce wind speeds and directions, and turbulence parameters, that are realistically correlated with altitude. To study the resulting correlations in the size and shape of the PSF, we generate simulated images of the PSF across a ~10 square-degree field of view – the size of the camera focal plane for the Vera C. Rubin Observatory in Chile – using weather data and historical seeing for a geographic location near the Observatory. We make quantitative predictions for two-point correlation functions (2PCF) that are used in analyses of cosmic shear. We observe a strong anisotropy in the two-dimensional 2PCF, which is expected based on observations in real images, and study the dependence of the orientation of the anisotropy on dominant wind directions near the ground and at higher altitudes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Standardizing reverberation-measured Mg II time-lag quasars, by using the radius–luminosity relation, and constraining cosmological model parameters

ABSTRACT We use 78 reverberation-measured Mg II time-lag quasars (QSOs) in the redshift range 0.0033 ≤ z ≤ 1.89 to constrain cosmological parameters in six different cosmological models. The basis of our method is the use of the radius–luminosity or R−L relation to standardize these 78 Mg II QSOs. In each cosmological model, we simultaneously determine R−L relation and cosmological model parameters, thus avoiding the circularity problem. We find that the R−L relation parameter values are independent of the cosmological model used in the analysis, thus, establishing that current Mg II QSOs are standardizable candles. Cosmological constraints obtained using these QSOs are significantly weaker than, but consistent with, those obtained from a joint analysis of baryon acoustic oscillation (BAO) observations and Hubble parameter [H(z)] measurements. So, we also analyse these QSOs in conjunction with the BAO + H(z) data and find cosmological constraints consistent with the standard spatially flat ΛCDM model as well as with mild dark energy dynamics and a little spatial curvature. A larger sample of higher quality reverberation-measured QSOs should have a smaller intrinsic dispersion and so should provide tighter constraints on cosmological parameters.

79 ASTRONOMY AND ASTROPHYSICS↗