Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “astronomy data analysis”

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

Integrating the PanDA Workload Management System with the Vera C. Rubin Observatory

The Vera C. Rubin Observatory will produce an unprecedented astronomical data set for studies of the deep and dynamic universe. Its Legacy Survey of Space and Time (LSST) will image the entire southern sky every three to four days and produce tens of petabytes of raw image data and associated calibration data over the course of the experiment’s run. More than 20 terabytes of data must be stored every night, and annual campaigns to reprocess the entire dataset since the beginning of the survey will be conducted over ten years. The Production and Distributed Analysis (PanDA) system was evaluated by the Rubin Observatory Data Management team and selected to serve the Observatory’s needs due to its demonstrated scalability and flexibility over the years, for its Directed Acyclic Graph (DAG) support, its support for multi-site processing, and its highly scalable complex workflows via the intelligent Data Delivery Service (iDDS). PanDA is also being evaluated for prompt processing where data must be processed within 60 seconds after image capture. This paper will briefly describe the Rubin Data Management system and its Data Facilities (DFs). Finally, it will describe in depth the work performed in order to integrate the PanDA system with the Rubin Observatory to be able to run the Rubin Science Pipelines using PanDA.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Investigating the effect of Milky Way dwarf spheroidal galaxies extension on dark matter searches with Fermi-LAT data

Satellite galaxies of the Milky Way with high mass-to-light ratios and little baryon content, i.e., dwarf spheroidal galaxies (dSphs), are among the most promising targets to detect or constrain the nature of dark matter (DM) through its final annihilation products into high-energy photons. Previously, the assumption that DM emission from dSphs is pointlike has been used to set strong constraints on DM candidates using data from the Fermi Large Area Telescope (LAT). However, due to their high DM densities and proximity, dSphs actually have sufficient angular extension to be detected by the Fermi-LAT. Here, we perform a comprehensive analysis about the impact of accounting for angular extension in the search for gamma-ray DM signals toward known dSphs with Fermi-LAT. Further, we show that, depending on the dSph under consideration, limits on the DM cross section can be weakened by up to a factor of 2–2.5, while the impact on the stacked, i.e., combined, limits is at most 1.5–1.8 depending on the annihilation channel. This result is of relevance when comparing dSph limits to other multimessenger DM constraints and for testing the DM interpretation of anomalous “excesses.”

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Subaru Hyper Suprime-Cam excavates colossal over- and underdense structures over 360 deg 2 out to z = 1

Subaru Strategic Program with the Hyper-Suprime Cam (HSC-SSP) has proven to be successful with its extremely wide area coverage in past years. Taking advantages of this feature, we report initial results from exploration and research of expansive over- and underdense structures at z = 0.3–1 based on the second Public Data Release where optical 5-band photometric data for ~ eight million sources with i < 23 mag are available over ~360 deg 2 . We not only confirm known superclusters but also find candidates of titanic over- and underdense regions out to z = 1. The mock data analysis suggests that the density peaks would involve one or more massive dark matter haloes (>1014 M ⊙ ) of the redshift, and the density troughs tend to be empty of massive haloes over >10 comoving Mpc. Besides, the density peaks and troughs at z ≲ 0.6 are in part identified as positive and negative weak lensing signals respectively, in mean tangential shear profiles, showing a good agreement with those inferred from the full-sky weak lensing simulation. The coming extensive spectroscopic surveys will be able to resolve these colossal structures in 3D space. Finally, the number density information over the entire survey field is available as grid-point data on the website of the HSC-SSP data release (https://hsc.mtk.nao.ac.jp/ssp/data-release/).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Harmonic analysis of discrete tracers of large-scale structure

It is commonplace in cosmology to analyze fields projected onto the celestial sphere, and in particular density fields that are defined by a set of points e.g. galaxies. When performing an harmonic-space analysis of such data (e.g. an angular power spectrum) using a pixelized map one has to deal with aliasing of small-scale power and pixel window functions. We compare and contrast the approaches to this problem taken in the cosmic microwave background and large-scale structure communities, and advocate for a direct approach that avoids pixelization. We describe a method for performing a pseudo-spectrum analysis of a galaxy data set and show that it can be implemented efficiently using well-known algorithms for special functions that are suited to acceleration by graphics processing units (GPUs). The method returns the same spectra as the more traditional map-based approach if in the latter the number of pixels is taken to be sufficiently large and the mask is well sampled. The method is readily generalizable to cross-spectra and higher-order functions. It also provides a convenient route for distributing the information in a galaxy catalog directly in harmonic space, as a complement to releasing the configuration-space positions and weights, and a route to spectral apodization. Finally, we make public a code enabling the application of our method to existing and upcoming datasets.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Observing correlations between dark matter accretion and galaxy growth: II. testing the impact of galaxy mass, star formation indicator, and neighbour colours

A crucial question in galaxy formation is what role new accretion has in star formation. Theoretical models have predicted a wide range of correlation strengths between halo accretion and galaxy star formation. Previously, we presented a technique to observationally constrain this correlation strength for isolated Milky Way mass galaxies at z ~ 0.12, based on the correlation between halo accretion and the density profile of neighbouring galaxies. By applying this technique to both observational data from the Sloan Digital Sky Survey and simulation data from the UNIVERSEMACHINE , where we can test different correlation strengths, we ruled out positive correlations between dark matter accretion and recent star formation activity. In this work, we expand our analysis by (1) applying our technique separately to red and blue neighbouring galaxies, which trace different infall populations, (2) correlating dark matter accretion rates with D n 4000 measurements as a longer-term quiescence indicator than instantaneous star-formation rates, and (3) analysing higher-mass isolated central galaxies with 10 11.0 < M * /M ⊙ < 10 11.5 out to z ~ 0.18. Finally, in all cases, our results are consistent with non-positive correlation strengths with ≳ 85 percent confidence, which is most consistent with models where processes such as gas recycling dominate star formation in massive z = 0 galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Laboratory Constraints on the Neutron-Spin Coupling of feV-Scale Axions

Ultralight axionlike particles can contribute to the dark matter near the Sun, leading to a distinct, stochastic signature in terrestrial experiments. We search for such particles through their neutron-spin coupling by reanalyzing approximately 40 days of data from a K- 3 He comagnetometer with a new frequency-domain likelihood-based formalism that properly accounts for stochastic effects over all axion coherence times relative to the experimental time span. Assuming that axions make up all of the dark matter in the Sun’s vicinity, we find a median 95% upper limit on the neutron-spin coupling of 2.4 × 10 −10 GeV −1 for most axion masses from 0.4 to 4 feV, which is about 5 orders of magnitude more stringent than previous laboratory bounds in that mass range. Although several peaks in the experiment’s magnetic power spectrum suggest the rejection of a white-noise null hypothesis, further analysis of their line shapes yields no positive evidence for a dark-matter axion.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Searching for Novel Chemistry in Exoplanetary Atmospheres Using Machine Learning for Anomaly Detection

Abstract The next generation of telescopes will yield a substantial increase in the availability of high-quality spectroscopic data for thousands of exoplanets. The sheer volume of data and number of planets to be analyzed greatly motivate the development of new, fast, and efficient methods for flagging interesting planets for reobservation and detailed analysis. We advocate the application of machine learning (ML) techniques for anomaly (novelty) detection to exoplanet transit spectra, with the goal of identifying planets with unusual chemical composition and even searching for unknown biosignatures. We successfully demonstrate the feasibility of two popular anomaly detection methods (local outlier factor and one-class support vector machine) on a large public database of synthetic spectra. We consider several test cases, each with different levels of instrumental noise. In each case, we use receiver operating characteristic curves to quantify and compare the performance of the two ML techniques.

Astronomy & Astrophysics↗
DOE OSTI2024

A lanthanide-rich kilonova in the aftermath of a long gamma-ray burst

Observationally, kilonovae are astrophysical transients powered by the radioactive decay of nuclei heavier than iron, thought to be synthesized in the merger of two compact objects. Over the first few days, the kilonova evolution is dominated by a large number of radioactive isotopes contributing to the heating rate. On timescales of weeks to months, its behaviour is predicted to differ depending on the ejecta composition and the merger remnant. Previous work has shown that the kilonova associated with gamma-ray burst 230307A is similar to kilonova AT2017gfo, and mid-infrared spectra revealed an emission line at 2.15 micrometres that was attributed to tellurium. Here, in this study, we report a multi-wavelength analysis, including publicly available James Webb Space Telescope data and our own Hubble Space Telescope data, for the same gamma-ray burst. We model its evolution up to two months after the burst and show that, at these late times, the recession of the photospheric radius and the rapidly decaying bolometric luminosity (L bol ∝ t −2.7±0.4 , where t is time) support the recombination of lanthanide-rich ejecta as they cool.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2021

A NICER View of the Massive Pulsar PSR J0740+6620 Informed by Radio Timing and XMM-Newton Spectroscopy

We report on Bayesian estimation of the radius, mass, and hot surface regions of the massive millisecond pulsar PSR J0740+6620, conditional on pulse-profile modeling of Neutron Star Interior Composition Explorer X-ray Timing Instrument event data. We condition on informative pulsar mass, distance, and orbital inclination priors derived from the joint North American Nanohertz Observatory for Gravitational Waves and Canadian Hydrogen Intensity Mapping Experiment/Pulsar wideband radio timing measurements of Fonseca et al. We use XMM-Newton European Photon Imaging Camera spectroscopic event data to inform our X-ray likelihood function. The prior support of the pulsar radius is truncated at 16 km to ensure coverage of current dense matter models. We assume conservative priors on instrument calibration uncertainty. We constrain the equatorial radius and mass of PSR J0740+6620 to be-+12.390.981.30km and-+2.0720.0660.067Me respectively, each reported as the posterior credible interval bounded by the 16% and 84% quantiles, conditional on surface hot regions that are non-overlapping spherical caps of fully ionized hydrogen atmosphere with uniform effective temperature; a posteriori, the temperature is=-+TlogK5.99100.060.05([])for each hot region. All software for the X-ray modeling framework is open-source and all data, model, and sample information is publicly available, including analysis notebooks and model modules in the Python language. Our marginal likelihood function of mass and equatorial radius is proportional to the marginal joint posterior density of those parameters(within the prior support)and can thus be computed from the posterior samples.

Millisecond pulsars↗
DOE OSTI2024

The BOSS bispectrum analysis at one loop from the Effective Field Theory of Large-Scale Structure

Abstract We analyze the BOSS power spectrum monopole and quadrupole, and the bispectrum monopole and quadrupole data, using the predictions from the Effective Field Theory of Large-Scale Structure (EFTofLSS). Specifically, we use the one loop prediction for the power spectrum and the bispectrum monopole, and the tree level for the bispectrum quadrupole. After validating our pipeline against numerical simulations as well as checking for several internal consistencies, we apply it to the observational data. We find that analyzing the bispectrum monopole to higher wavenumbers thanks to the one-loop prediction, as well as the addition of the tree-level quadrupole, significantly reduces the error bars with respect to our original analysis of the power spectrum at one loop and bispectrum monopole at tree level. After fixing the spectral tilt to Planck preferred value and using a Big Bang Nucleosynthesis prior, we measureσ 8 = 0.794 ± 0.037,h = 0.692 ± 0.011, and Ω m = 0.311 ± 0.010 to about 4.7%, 1.6%, and 3.2%, at 68% CL, respectively.This represents an error bar reduction with respect to the power spectrum-only analysis of about 30%, 18%, and 13% respectively. Remarkably, the results are compatible with the ones obtained with a power-spectrum-only analysis, showing the power of the EFTofLSS in simultaneously predicting several observables. We find no tension with Planck.

Astronomy & Astrophysics↗
DOE OSTI2024

Understanding the origin of early-type dwarfs: the spectrophotometric study of CGCG014−074

ABSTRACT Early-type dwarf galaxies constitute a prevalent population in the central regions of rich groups and clusters in the local Universe. These low-luminosity and low-mass stellar systems play a fundamental role in the assembly of the luminous galaxies observed today, according to the Lambda cold dark matter hierarchical theory. The origin of early-type dwarfs has been linked to the transformation of disc galaxies interacting with the intracluster medium, especially in dense environments. However, the existence of low-luminosity early-type galaxies in low-density environments presents a challenge to this scenario. This study presents a comprehensive photometric and spectroscopic analysis of the early-type dwarf galaxy CGCG014−074 using deep Gemini GMOS (Gemini Multi-Object Spectrograph) data, focusing on its peculiarities and evolutionary implications. CGCG014−074 exhibits distinct features, including a rotating inner disc, an extended stellar formation with a quiescent phase since about 2 Gyr ago, and the presence of boxy isophotes. From the kinematic analysis, we confirm CGCG014−074 as a nucleated early-type dwarf galaxy with embedded disc. The study of its stellar population parameters using different methods provides significant insights into the galaxy’s evolutionary history. These results show an old and metal-poor nucleus (${\sim}9.3$ Gyr and $\mathrm{[Z/H]}\sim -0.84$ dex), while the stellar disc is younger (${\sim}4.4$ Gyr) with a higher metallicity ($\mathrm{[Z/H]}\sim -0.40$ dex). These distinctive features collectively position CGCG014−074 as a likely building block galaxy that has evolved passively throughout its history.

Astronomy & Astrophysics↗
NASA NTRS1998

Fifteen Years of Laboratory Astrophysics at Ames

Tremendous strides have been made in our understanding of interstellar material over the past fifteen years thanks to significant, parallel developments in two closely related areas: observational astronomy and laboratory astrophysics. Fifteen years ago the composition of interstellar dust was largely guessed at, the concept of ices in dense molecular clouds ignored, and the notion of large, abundant, gas phase, carbon-rich molecules widespread throughout the interstellar medium (ISM) considered impossible. Today the composition of dust in the diffuse ISM is reasonably well constrained to cold refractory materials comprised of amorphous and crystalline silicates mixed with an amorphous carbonaceous material containing aromatic structural units and short, branched aliphatic chains. In the dense ISM, these cold dust particles are coated with mixed-molecular ices whose compositions are very well known. Lastly, the signature of carbon-rich polycyclic aromatic hydrocarbons (PAHs), shockingly large molecules by early interstellar chemistry standards, is widespread throughout the ISM. This great progress has only been made possible by the close collaboration of laboratory experimentalists with observers and theoreticians, all with the goal of applying their skills to astrophysical problems of direct interest to NASA programs. Such highly interdisciplinary collaborations ensure fundamental, in depth coverage of the wide-ranging challenges posed by astrophysics. These challenges include designing astrophysically focused experiments and data analysis, tightly coupled with astrophysical searches spanning 2 orders of magnitude in wavelength, and detailed theoretical modeling. The impact of our laboratory has been particularly effective as there is constant cross-talk and feedback between quantum theorists; theoretical astrophysicists and chemists; experimental physicists; organic, physical and petroleum chemists; and infrared and UV/Vis astronomers. In this paper, two examples of the Ames Program will be given. We have been involved in identifying 9 out of the 14 interstellar pre-cometary ice species known, determined their abundances and the physical nature of the ice structure. Details on our ice work are given in the paper by Sandford et al. Our group is among the pioneers of the PAH model. We built the theoretical framework, participated in the observations and developed the experimental techniques needed to test the model. We demonstrated that the ubiquitous infrared emission spectrum associated with many interstellar objects can be matched by laboratory spectra of neutral and positively charged PAHs and that PAHs were excellent candidates for the diffuse interstellar band (DIB) carriers. See Salama et al. and Hudgins et al.

Allamandola, L. J.↗
DOE OSTI2024

Assessing and benchmarking the fidelity of posterior inference methods for astrophysics data analysis

In this era of large and complex astronomical survey data, interpreting, validating, and comparing inference techniques becomes increasingly difficult. This is particularly critical for emerging inference methods like Simulation-Based Inference (SBI), which offer significant speedup potential and posterior modeling flexibility, especially when deep learning is incorporated. We present a study to assess and compare the performance and uncertainty prediction capability of Bayesian inference algorithms – from traditional MCMC sampling of analytic functions to deep learning-enabled SBI. We focus on testing the capacity of hierarchical inference modeling in those scenarios. Before we extend this study to cosmology, we first use astrophysical simulation data to ensure interpretability. We demonstrate a probabilistic programming implementation of hierarchical and non-hierarchical Bayesian inference using simulations derived from the DeepBench software library, a benchmarking tool developed by our group that generates simple and controllable astrophysical objects from first principles. This study will enable astronomers and physicists to harness the inference potential of these methods with confidence.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Oscillating scalar fields and the Hubble tension: A resolution with novel signatures

Here, we present a detailed investigation of a subdominant oscillating scalar field [“early dark energy” (EDE)] in the context of resolving the Hubble tension. Consistent with earlier work, but without relying on fluid approximations, we find that a scalar field frozen due to Hubble friction until log 10 ($z_c$) ~ 3.5, reaching $ρ_{\text{EDE}}(z_c)/ρ_{\text{tot}}$ ~ 10 % and diluting faster than matter afterwards, can bring cosmic microwave background (CMB), baryonic acoustic oscillations, supernovae luminosity distances, and the late-time estimate of the Hubble constant from the SH0ES Collaboration into agreement. A scalar field potential that scales as $V (\phi) \propto \phi ^{2n}$ with 2 ≲ $\textit{n}$ ≲ 3.4 around the minimum is preferred at the 68% confidence level, and the Planck polarization places additional constraints on the dynamics of perturbations in the scalar field. In particular, the data prefer a potential that flattens at large field displacements. A Markov-chain Monte Carlo analysis of mock data shows that the next-generation CMB observations (i.e., CMB-S4) can unambiguously detect the presence of the EDE at a very high significance. This projected sensitivity to the EDE dynamics is mainly driven by improved measurements of the E-mode polarization. We also explore new observational signatures of EDE scalar field dynamics: (i) We find that depending on the strength of the tensor-to-scalar ratio, the presence of the EDE might imply the existence of isocurvature perturbations in the CMB. (ii) We show that a strikingly rapid, scale-dependent growth of EDE field perturbations can result from parametric resonance driven by the anharmonic oscillating field for $\textit{n}$ ≈ 2. This instability and ensuing potentially nonlinear, spatially inhomogeneous, dynamics may provide unique signatures of this scenario.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Non-Gaussianity in the weak lensing correlation function likelihood – implications for cosmological parameter biases

ABSTRACT We study the significance of non-Gaussianity in the likelihood of weak lensing shear two-point correlation functions, detecting significantly non-zero skewness and kurtosis in 1D marginal distributions of shear two-point correlation functions in simulated weak lensing data. We examine the implications in the context of future surveys, in particular LSST, with derivations of how the non-Gaussianity scales with survey area. We show that there is no significant bias in 1D posteriors of Ωm and σ8 due to the non-Gaussian likelihood distributions of shear correlations functions using the mock data (100 deg2). We also present a systematic approach to constructing approximate multivariate likelihoods with 1D parametric functions by assuming independence or more flexible non-parametric multivariate methods after decorrelating the data points using principal component analysis (PCA). While the use of PCA does not modify the non-Gaussianity of the multivariate likelihood, we find empirically that the 1D marginal sampling distributions of the PCA components exhibit less skewness and kurtosis than the original shear correlation functions. Modelling the likelihood with marginal parametric functions based on the assumption of independence between PCA components thus gives a lower limit for the biases. We further demonstrate that the difference in cosmological parameter constraints between the multivariate Gaussian likelihood model and more complex non-Gaussian likelihood models would be even smaller for an LSST-like survey. In addition, the PCA approach automatically serves as a data compression method, enabling the retention of the majority of the cosmological information while reducing the dimensionality of the data vector by a factor of ∼5.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2000

Acousto-Optic Imaging Spectrometers for Mars Surface Science

NASA's long term plan for Mars sample collection and return requires a highly streamlined approach for spectrally characterizing a landing site, documenting the mineralogical make-up of the site and guiding the collections of samples which represent the diversity of the site. Ideally, image data should be acquired at hundreds of VIS and IR wavelengths, in order to separately distinguish numerous anticipated species, using principal component analysis and linear unmixing. Cameras with bore-sighted point spectrometers can acquire spectra of isolated scene elements, but it requires 10(exp 2) to 10(exp 2) successive motions and precise relative pointing knowledge in order to create a single data cube which qualifies as a spectral map. These and other competing science objectives have to be accomplished within very short lander/rover operational lifetime (a few sols). True, 2-D imaging spectroscopy greatly speeds up the data acquisition process, since the spectra of all pixels in the scene are collected at once. This task can be accomplished with cameras that use electronically tunable acousto-optic tunable filters (AOTFs) as the optical tuning element. AOTFs made from TeO2 are now a mature technology, and operate at wavelengths from near-UV to about 5 microns. Because of incremental improvements in the last few years, present generation devices are rugged, radiation-hard and operate at temperatures down to at least 150K so they can be safely integrated into the ambient temperature optics of in-situ instruments such as planetary or small-body landers. They have been used for ground-based astronomy, and were also baselined for the ST-4 Champollion IR comet lander experiment (CIRCLE), prior to cancellation of the ST-4 mission last year. AIMS (for Acousto-optic Imaging spectrometer), is a prototype lander instrument which is being built at GSFC with support by the NASA OSS Advanced Technologies and Mission Studies, Mars Instrument Definition and Development Program (MIDP). AIMS is capable of tunable spectroscopic imaging of surface mineralogy, ices and dust between 0.5 and 2.4 microns, at a resolving power (lambda/delta lambda) which is typically several hundred. The design spatial resolution, similar to IMP and SSI, will allow mapping at scales down to about 1 cm.

Glenar, D. A.↗
DOE OSTI2020

MADHAT: Model-Agnostic Dark Halo Analysis Tool

We present the Model-Agnostic Dark Halo Analysis Tool (MADHAT), a numerical tool which implements a Fermi-LAT data-driven, model-independent analysis of gamma-ray emission from dwarf satellite galaxies and dwarf galaxy candidates due to dark matter annihilation, dark matter decay, or other nonstandard or unknown astrophysics. This tool efficiently provides statistical upper bounds on the number of observed photons in excess of the number expected, based on empirical determinations of foregrounds and backgrounds, using a stacked analysis of any selected set of dwarf targets. It also calculates the resulting bounds on the properties of dark matter under any assumptions the user makes regarding dark sector particle physics or astrophysics. As an application, we determine new bounds on Sommerfeld-enhanced dark matter annihilation in a set of eight dwarfs. MADHAT v1.0 includes 58 dwarfs and dwarf candidate targets, and we discuss future planned developments. MADHAT is available and will be maintained at https://github.com/MADHATdm.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

On the Formation of Phantom Electron Phase Space Density Peaks in Single Spacecraft Radiation Belt Data

This study examines the rapid losses and acceleration of trapped relativistic and ultrarelativistic electron populations in the Van Allen radiation belt during the September 7–9, 2017 geomagnetic storm. By analyzing the dynamics of the last closed drift shell (LCDS), and the electron flux and phase space density (PSD), we show that the electron dropouts are consistent with magnetopause shadowing and outward radial diffusion to the compressed LCDS. During the recovery phase, an in-bound pass of Van Allen Probe A shows an apparent local peak in PSD, which does not exist in reality. Here, a careful analysis of the multipoint measurements by the Van Allen Probes reveals instead how the apparent PSD peak arises from aliasing monotonic PSD profiles which are rapidly increasing due to acceleration from very fast inwards radial diffusion. In the absence of such multisatellite conjunctions during fast acceleration events, such peaks might otherwise be associated with local acceleration processes.

79 ASTRONOMY AND ASTROPHYSICS↗