Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 307 records · Page 17

Uncertainty quantification of the virial black hole mass with conformal prediction

Precise measurements of the black hole mass are essential to gain insight on the black hole and host galaxy co-evolution. A direct measure of the black hole mass is often restricted to nearest galaxies and instead, an indirect method using the single-epoch virial black hole mass estimation is used for objects at high redshifts. However, this method is subjected to biases and uncertainties as it is reliant on the scaling relation from a small sample of local active galactic nuclei. In this study, we propose the application of conformalized quantile regression (CQR) to quantify the uncertainties of the black hole predictions in a machine learning setting. We compare CQR with various prediction interval techniques and demonstrated that CQR can provide a more useful prediction interval indicator. In contrast to baseline approaches for prediction interval estimation, we show that the CQR method provides prediction intervals that adjust to the black hole mass and its related properties. That is it yields a tighter constraint on the prediction interval (hence more certain) for a larger black hole mass, and accordingly, bright and broad spectral line width source. Using a combination of neural network model and CQR framework, the recovered virial black hole mass predictions and uncertainties are comparable to those measured from the Sloan Digital Sky Survey. The code is publicly available.

79 ASTRONOMY AND ASTROPHYSICS↗

Nuclear data for space exploration

Understanding the harmful effects of galactic cosmic rays (GCRs) on space exploration requires a substantial amount of nuclear data. Specifically, the interaction of energetic GCR charged particles with spacecraft materials generates secondary radiations that, through energy deposition, can harm astronauts and electronic systems. By identifying the gaps in our knowledge of the relevant nuclear data—such as interaction cross sections—and identifying ways to fill those gaps—with measurements, compilations, evaluations, dissemination, reaction modeling, sensitivity studies, and uncertainty quantification—the safety and viability of space exploration can be improved. This work surveys the state of the art in this interdisciplinary field and identifies promising collaborative research topics that have significant potential to advance our understanding of the effects of the space radiation environment on space exploration.

42 ENGINEERING↗

Sweeping Horndeski canvas: new growth-rate parameterization for modified-gravity theories

We propose and numerically validate a new fitting formula that is sufficiently accurate to model the growth of structure in Horndeski theories of modified gravity for upcoming Stage IV and V large-scale structure surveys. Based on an analysis of more than 18,000 Horndeski models and adopting the popular parameterization of the growth rate f(z) = Ω M (z) γ , we generalize the constant growth index γ to a two-parameter redshift-dependent quantity, γ(z), that more accurately fits these models. We demonstrate that the functional form γ(z) = γ 0 + γ 0 z 2 /(1 + z) improves the median χ2 of the fit to viable Horndeski models by a factor of ~40 relative to that of a constant γ, and is sufficient to obtain unbiased results even for precise measurements expected in Stage IV and V surveys. Finally, we constrain the parameters of the new fitting formula using current cosmological data.

79 ASTRONOMY AND ASTROPHYSICS↗

Modeling Electron Acceleration and Transport in the Early Impulsive Phase of the 2017 September 10th Solar Flare

The X8.2-class limb flare on 2017 September 10 is among the best studied solar flare events owing to its great similarity to the standard flare model and the broad coverage by multiple spacecraft and ground-based observations. These multiwavelength observations indicate that electron acceleration and transport are efficient in the reconnection and flare looptop regions. However, there lacks a comprehensive model for explaining and interpreting the multi-faceted observations. In this work, we model the electron acceleration and transport in the early impulsive phase of this flare. We solve the Parker transport equation that includes the primary acceleration mechanism during magnetic reconnection in the large-scale flare region modeled by MHD simulations. We find that electrons are accelerated up to several MeV and fill a large volume of the reconnection region, similar to the observations shown in microwaves. The electron spatial distribution and spectral shape in the looptop region agree well with those derived from the microwave and hard X-ray emissions before magnetic islands grow large and dominate the acceleration. Future emission modelings using the electron maps will enable direct comparison with microwave and hard X-ray observations. These results shed new light on the electron acceleration and transport in a broad region of solar flares within a data-constrained realistic flare geometry.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA Computational Case Study: Spectral Energy Distribution Fitting

The need for faster, more efficient algorithms is an important aspect of scientific computing. Generally, scientists are only exposed to computational issues that arise in their field. Thus, collaboration between a numerical analyst and a scientist is becoming necessary for scientific computing. The purpose of this case study is to expose computer scientists to processes that an astronomer would use to obtain useful results from raw data. For example, astronomers are interested in determining the properties of galaxies and measuring changes in those properties as a function of time throughout cosmic history. To do so, they use certain models that are designed and refined over time via observations at different wavelengths of the light spectrum. The process of matching these models with observed data from studying celestial bodies is referred to as Spectral Energy Distribution (SED) fitting. In this case study, we learn how to perform the SED fit. This process requires knowledge of both the astronomical and computational issues involved when fitting flux, the total energy from a source as seen from Earth, to a set of physical templates. Once a fit is complete, one can classify the source and estimate a number of physical parameters. The goal is to demonstrate how the computer science skill set can be used in the scientific community and to possibly improve one or more of the computational aspects of this problem. The following provides some background on the astronomy issues, including information on the spectral energy distribution and related physical parameters, and the computational issues, including the fitting procedure..

modeling↗

Radio interferometric detection of a traveling ionospheric disturbance excited by the explosion of Mount St. Helens

A large-amplitude traveling ionospheric disturbance (TID) was detected over Owens Valley, California, on May 18, 1980, by a highly sensitive very long baseline interferometry (VLBI) radio astronomy experiment. This TID is interpreted as the response of the ionosphere to a gravity wave excited in the neutral atmosphere by the explosion of Mount St. Helens that took place at 1532 UT on that day. A model, invoking the point-excitation of internal gravity waves in an isothermal atmosphere, which fits observations of the TID at several other stations, leads to identification of the features observed in the VLBI data. Small-amplitude higher-frequency changes in the ionosphere were detected for several hours after the passage of the large-amplitude Mount St. Helens TID, but it is not clear whether these were excited by the passage of the gravity wave or were background fluctuations.

Roberts, D. H.↗

CMB power spectra and cosmological parameters from Planck PR4 with CamSpec

ABSTRACT We present angular power spectra and cosmological parameter constraints derived from the Planck PR4 (NPIPE) maps of the cosmic microwave background. NPIPE, released by the Planck Collaboration in 2020, is a new processing pipeline for producing calibrated frequency maps from Planck data. We have created new versions of the CamSpec likelihood using these maps and applied them to constrain the lambda cold dark matter (ΛCDM) model and single-parameter extensions. We find excellent consistency between NPIPE and the Planck 2018 maps at the parameter level, showing that the Planck cosmology is robust to substantial changes in the mapmaking. The lower noise of NPIPE leads to ~10 per cent tighter constraints, and we see both smaller error bars and a shift towards the ΛCDM values for beyond-ΛCDM parameters including ΩK and AL.

79 ASTRONOMY AND ASTROPHYSICS↗

Advances in constraining intrinsic alignment models with hydrodynamic simulations

We use galaxies from the illustristng, massiveblack-ii, and illustris-1 hydrodynamic simulations to investigate the behaviour of large scale galaxy intrinsic alignments. Our analysis spans four redshift slices over the approximate range of contemporary lensing surveys z = 0–1. We construct comparable weighted samples from the three simulations, which we then analyse using an alignment model that includes both linear and quadratic alignment contributions. Our data vector includes galaxy–galaxy, galaxy–shape, and shape–shape projected correlations, with the joint covariance matrix estimated analytically. In all of the simulations, we report non-zero IAs at the level of several σ. For a fixed lower mass threshold, we find a relatively strong redshift dependence in all three simulations, with the linear IA amplitude increasing by a factor of ~2 between redshifts z = 0 and z = 1. We report no significant evidence for non-zero values of the tidal torquing amplitude, A 2 , in TNG, above statistical uncertainties, although MBII favours a moderately negative A 2 ~ –2. Examining the properties of the TATT model as a function of colour, luminosity and galaxy type (satellite or central), our findings are consistent with the most recent measurements on real data. We also outline a novel method for constraining the TATT model parameters directly from the pixelized tidal field, alongside a proof-of-concept exercise using TNG. This technique is shown to be promising, although comparison with previous results obtained via other methods is non-trivial.

79 ASTRONOMY AND ASTROPHYSICS↗

Constraints on f ( R ) and normal-branch Dvali-Gabadadze-Porrati modified gravity model parameters with cluster abundances and galaxy clustering

In this report we present forecasted cosmological constraints from combined measurements of galaxy cluster abundances from the Simons Observatory and galaxy clustering from a DESI-like experiment on two well-studied modified gravity models, the chameleon-screened Hu-Sawicki f(R) model and the nDGP braneworld Vainshtein model. A Fisher analysis is conducted using σ 8 constraints derived from thermal Sunyaev-Zel’dovich (tSZ) selected galaxy clusters as well as linear and quasilinear redshift-space 2-point galaxy correlation functions. We find that the cluster abundances drive the constraints on the nDGP model while f(R) constraints are led by galaxy clustering. The two tracers of the cosmological gravitational field are found to be complementary, and their combination significantly improves constraints on the f(R) in particular in comparison to each individual tracer alone. For a model of f(R) with a general relativity (GR) fiducial case (f R0 = 0), we find a 2-σ upper limit of f R0 ≤ 5.68 x10 -7 . For the well-studied log-based fiducial parameter value in f(R), log 10 (f R0 ) = -5, paired with the parameter value n = 1, we find combined 1-σ constraints of σ(log 10 (f R0 )) = 0.12 and σ(n) = 0.36. For the nDGP model with fiducial n nDGP = 1 we find σ( n nDGP ) = 0.087. Our results present the exciting potential to utilize upcoming galaxy and CMB survey data available in the near future to discern and/or constrain cosmic deviations from GR.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Cloud Atlas: Unraveling the Vertical Cloud Structure with the Time-series Spectrophotometry of an Unusually Red Brown Dwarf

Rotational modulations of emission spectra in brown dwarf and exoplanet atmospheres show that clouds are often distributed non-uniformly in these ultracool atmospheres. The spatial heterogeneity in cloud distribution demonstrates the impact of atmospheric dynamics on cloud formation and evolution. In this study, we update the Hubble Space Telescope (HST) time-series data analysis of the previously reported rotational modulations of WISEP J004701+680352—an unusually red late-L brown dwarf with a spectrum similar to that of the directly imaged planet HR 8799e. We construct a self-consistent spatially heterogeneous cloud model to explain the HST and the Spitzer time-series observations, as well as the time-averaged spectra of WISE 0047. Additionally, in the heterogeneous cloud model, a cloud thickness variation of around one pressure scale height explains the wavelength dependence in the HST near-IR spectral variability. By including disequilibrium CO/CH 4 chemistry, our models also reproduce the redder $J-{K}_{{\rm{s}}}$ color of WISE 0047 compared to that of field brown dwarfs. We discuss the impact of vertical cloud structure on atmospheric profile and estimate the minimum eddy diffusivity coefficient for other objects with redder colors. Our data analysis and forward modeling results demonstrate that time-series spectrophotometry with a broad wavelength coverage is a powerful tool for constraining heterogeneous atmospheric structure.

79 ASTRONOMY AND ASTROPHYSICS↗

Definitive orbit determination for the HEAO-2 spacecraft

Precise ephemerides for the High Energy Astronomy Observatory-2 (HEAO-2) were computed to assist in the Charles Stark Draper Laboratory development of an onboard orbit determination technique. Weighted least-squares, batch orbit solutions were calculated using a high-precision earth gravity model and an approximate model for intermittent spacecraft thrusting. With these improvements, orbit solution consistencies at the 50- to 100-meter level were attained using the Goddard Trajectory Determination System and NASA S-band tracking data.

Smith, R. L.↗

Solar maximum mission: Ground support programs at the Harvard Radio Astronomy Station

Observations of the spectral characteristics of solar radio bursts were made with new dynamic spectrum analyzers of high sensitivity and high reliability, over the frequency range 25-580 MHz. The observations also covered the maximum period of the current solar cycle and the period of international cooperative programs designated as the Solar Maximum Year. Radio data on shock waves generated by solar flares were combined with optical data on coronal transients, taken with equipment on the SMM and other satellites, and then incorporated into computer models for the outward passage of fast-mode MHD shocks through the solar corona. The MHD models are non-linear, time-dependent and for the most recent models, quasi-three-dimensional. They examine the global response of the corona for different types of input pulses (thermal, magnetic, etc.) and for different magnetic topologies (for example, open and closed fields). Data on coronal shocks and high-velocity material ejected from solar flares have been interpreted in terms of a model consisting of three main velocity regimes.

Maxwell, A.↗

Weak Lensing And Cosmology In The Dark Energy Survey

In this PhD Thesis, we present studies which aim at understanding and testing the fundamental assumptions of the $\Lambda$CDM cosmological model. We initially show that weak gravitational lensing can be derived as a consequence of General Relativity, and that it can be applied to astronomical data from the Dark Energy Survey (DES) in order to constrain the amplitude of matter fluctuations in the universe. This constraint can be used as an end-to-end cosmological test when it is compared to inferences from the Cosmic Microwave Background. We further extend the cosmological model to account for non-standard parameters which can point to novel physics, and use DES data to assess their statistical significance. We finally show that the Cold Dark Matter (CDM) paradigm can also be fundamentally tested with observational data using morphological signatures of disk galaxies in cluster environments if CDM is self-interacting.

79 ASTRONOMY AND ASTROPHYSICS↗

Study of solid state photomultiplier

Available solid state photomultiplier (SSPM) detectors were tested under low-background, low temperature conditions to determine the conditions producing optimal sensitivity in a space-based astronomy system such as a liquid cooled helium telescope in orbit. Detector temperatures varied between 6 and 9 K, with background flux ranging from 10 to the 13th power to less than 10 to the 6th power photons/square cm-s. Measured parameters included quantum efficiency, noise, dark current, and spectral response. Experimental data were reduced, analyzed, and combined with existing data to build the SSPM data base included herein. The results were compared to analytical models of SSPM performance where appropriate models existed. Analytical models presented here were developed to be as consistent with the data base as practicable. Significant differences between the theory and data are described. Some models were developed or updated as a result of this study.

Hays, K. M.↗

Probing the Extragalactic Mid-infrared Background with HAWC

The extragalactic background light (EBL) contains all the radiation emitted by nuclear and accretion processes in stars and compact objects since the epoch of recombination. Measuring the EBL density directly is challenging, especially in the near-to-far-infrared wave band, mainly due to the zodiacal light foreground. Instead, gamma-ray astronomy offers the possibility to indirectly set limits on the EBL by studying the effects of gamma-ray absorption in the very high energy (VHE: >100 GeV) spectra of distant blazars. The High Altitude Water Cherenkov Gamma Ray Observatory (HAWC) is one of the few instruments sensitive to gamma rays with energies above 10 TeV. This offers the opportunity to probe the EBL in the near/mid-IR region: λ = 1–100 μm. In this study, we fit physically motivated emission models to Fermi-LAT gigaelectronvolt data to extrapolate the intrinsic teraelectronvolt spectra of blazars. We then simulate a large number of absorbed spectra for different randomly generated EBL model shapes and calculate Bayesian credible bands in the EBL intensity space by comparing and testing the agreement between the absorbed spectra and HAWC extragalactic observations of two blazars. The resulting bands are in agreement with current EBL lower and upper limits, showing a downward trend toward higher wavelength values λ > 10 μm also observed in previous measurements.

79 ASTRONOMY AND ASTROPHYSICS↗

Gamma-ray burst data strongly favour the three-parameter fundamental plane (Dainotti) correlation over the two-parameter one

ABSTRACT Gamma-ray bursts (GRBs), observed to redshift z = 9.4, are potential probes of the largely unexplored z ∼ 2.7–9.4 part of the early Universe. Thus, finding relevant relations among GRB physical properties is crucial. We find that the Platinum GRB data compilation, with 50 long GRBs (with relatively flat plateaus and no flares) in the redshift range 0.553 ≤ z ≤ 5.0, and the LGRB95 data compilation, with 95 long GRBs in 0.297 ≤ z ≤ 9.4, as well as the 145 GRB combination of the two, strongly favour the 3D Fundamental Plane (Dainotti) correlation (between the peak prompt luminosity, the luminosity at the end of the plateau emission, and its rest-frame duration) over the 2D one (between the luminosity at the end of the plateau emission and its duration). The 3D Dainotti correlations in the three data sets are standardizable. We find that while LGRB95 data have ∼50 per cent larger intrinsic scatter parameter values than the better-quality Platinum data, they provide somewhat tighter constraints on cosmological-model and GRB-correlation parameters, perhaps solely due to the larger number of data points, 95 versus 50. This suggests that when compiling GRB data for the purpose of constraining cosmological parameters, given the quality of current GRB data, intrinsic scatter parameter reduction must be balanced against reduced sample size.

79 ASTRONOMY AND ASTROPHYSICS↗

Neural Networks for Nuclear Reactions in MAESTROeX

We demonstrate the use of neural networks to accelerate the reaction steps in the MAESTROeX stellar hydrodynamics code. A traditional MAESTROeX simulation uses a stiff ODE integrator for the reactions; here, we employ a ResNet architecture and describe details relating to the architecture, training, and validation of our networks. Our customized approach includes options for the form of the loss functions, a demonstration that the use of parallel neural networks leads to increased accuracy, and a description of a perturbational approach in the training step that robustifies the model. We test our approach on millimeter-scale flames using a single-step, 3-isotope network describing the first stages of carbon fusion occurring in Type Ia supernovae. We train the neural networks using simulation data from a standard MAESTROeX simulation, and show that the resulting model can be effectively applied to different flame configurations. This work lays the groundwork for more complex networks, and iterative time-integration strategies that can leverage the efficiency of the neural networks.

79 ASTRONOMY AND ASTROPHYSICS↗