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 649 records · Page 36

Cosmology from clustering, cosmic shear, CMB lensing, and cross correlations: combining Rubin observatory and Simons Observatory

ABSTRACT In the near future, the overlap of the Rubin Observatory Legacy Survey of Space and Time (LSST) and the Simons Observatory (SO) will present an ideal opportunity for joint cosmological data set analyses. In this paper, we simulate the joint likelihood analysis of these two experiments using six two-point functions derived from galaxy position, galaxy shear, and CMB lensing convergence fields. Our analysis focuses on realistic noise and systematics models and we find that the dark energy Figure of merit (FoM) increases by 53 per cent (92 per cent) from LSST-only to LSST+SO in Year 1 (Year 6). We also investigate the benefits of using the same galaxy sample for both clustering and lensing analyses, and find the choice improves the overall signal to noise by ${\sim}30{-}40{{\ \rm per\ cent}}$, which significantly improves the photo-z calibration and mildly improves the cosmological constraints. Finally, we explore the effects of catastrophic photo-z outliers finding that they cause significant parameter biases when ignored. We develop a new mitigation approach termed ‘island model’, which corrects a large fraction of the biases with only a few parameters while preserving the constraining power.

79 ASTRONOMY AND ASTROPHYSICS↗

Understanding Asteroid 16 Psyche’s composition through 3D impact crater modeling

Asteroid 16 Psyche is the largest M-type (metallic) Main Belt Asteroid (MBA). Radar albedo data indicate Psyche’s surface is rich in metallic content, but estimates for Psyche’s bulk structure vary widely, with bulk density estimates ranging from 1.4 ± 0.3 g/cm 3 to 4.5 ± 1.4 g/cm 3 and porosity estimates ranging from 30% to 70%. A forthcoming NASA discovery mission will visit Psyche. Narrowing the range of estimated compositions of Psyche prior to mission launch can help to determine the necessary tools for asteroid analysis. Psyche has two large impact structures in its Southern hemisphere, with estimated diameters from 50 km to 70 km and estimated depths up to 6.4 km. In this work, we present results from 2D and 3D hydrocode simulations of the largest of these impact structures. Through quantitative comparison of the simulated crater dimensions with measured values, our models suggest that Psyche is largely composed of porous, metallic material. In addition, to match the crater dimensions, our simulations indicate that the angle of impact was significantly altered from vertical by 45–60 degrees.

79 ASTRONOMY AND ASTROPHYSICS↗

MaNGA galaxy properties – II. A detailed comparison of observed and simulated spiral galaxy scaling relations

ABSTRACT We present a catalogue of dynamical properties for 2368 late-type galaxies from the MaNGA survey. The latter complements the catalogue of photometric properties for the same sample based on deep optical dark energy sky instrument legacy imaging survey (DESI) photometry processed with autoprof . Rotation curves (RCs), extracted by model-fitting H α velocity maps from the MaNGA Data Analysis Pipeline, extend out to 1.4 (1.9) Re for the primary (secondary) MaNGA samples, respectively. The RCs and ancillary MaNGA Pipe3D data products were used to construct various fundamental galaxy scaling relations that are also compared uniformly with similar relations from numerical investigation of a hundred astrophysical object (NIHAO) zoom-in simulations. Simulated NIHAO galaxies were found to broadly reproduce the observed MaNGA galaxy population for $\log (M_*/{\rm {\rm M}_{\odot })\, \gt\, 8.5}$. Some discrepancies remain, such as those pertaining to central stellar densities and the diversity of RCs due to strong feedback schemes. Also presented are spatially resolved scatters for the velocity-size-stellar mass (VRM*) structural relations using MaNGA and NIHAO samples. The scatter for these relations in the galaxian interiors is a consequence of the diversity of inner RC shapes, while scatter in the outskirts is dictated by the large range of stellar surface densities, which is itself driven by sporadic star formation. The detailed, spatially resolved scatter analysis highlights the complex interplay between local and global astrophysical processes and provides a strong constraint to numerical simulations.

Astronomy & Astrophysics↗

Dark matter substructure or source model systematics? A case study of cluster lens Abell S1063

Mapping the small-scale structure of the universe through gravitational lensing is a promising tool for probing the particle nature of dark matter. Curved Arc Basis (CAB) has been proposed as a local lensing formalism in galaxy clusters, with the potential to detect low-mass dark matter substructure. In this work, we analyse the cluster lens Abell S1063 in search of dark matter substructure with the CAB formalism, using multiband imaging data from James Webb Space Telescope ( JWST ). We use two different source modelling methods: shapelets and pixel-based source reconstruction based on Delaunay triangulation. We find that source modelling systematics from shapelets result in a disagreement between CAB parameters measured from different filters. Source modelling with Delaunay significantly alleviates this systematic, as seen in the improvement in agreement across filters. We also find that inadequate complexity in source modelling can result in convincing spurious detections of dark matter substructure from strong gravitational lenses, as seen by our $\Delta \text{BIC} > 20$ measurement of a $M \sim 10^{10}$ ${\rm M}_{\odot }$ subhalo with shapelets, a spurious detection that is not reproduced with Delaunay source modelling. We demonstrate that multiband analysis with different JWST filters is key for disentangling source and lens model systematics from dark matter substructure detections.

79 ASTRONOMY AND ASTROPHYSICS↗

On the minimum number of radiation field parameters to specify gas cooling and heating functions

Fast and accurate approximations of gas cooling and heating functions are needed for hydrodynamic galaxy simulations. We use machine learning to analyze atomic gas cooling and heating functions in the presence of a generalized incident local radiation field computed by Cloudy. We characterize the radiation field through binned radiation field intensities instead of the photoionization rates used in our previous work. We find a set of 6 energy bins whose intensities exhibit relatively low correlation. We use these bins as features to train machine learning models to predict Cloudy cooling and heating functions at fixed metallicity. We compare the relative SHapley Additive exPlanation (SHAP) value importance of the features. From the SHAP analysis, we identify a feature subset of 3 energy bins (0.5-1, 1-4, and 13-16Ry) with the largest importance and train additional models on this subset. We compare the mean squared errors and distribution of errors on both the entire training data table and a randomly selected 20% test set withheld from model training. The machine learning models trained with 3 and 6 bins, as well as 3 and 4 photoionization rates, have comparable accuracy everywhere, with errors ≳10 times smaller than for the interpolation table of Gnedin and Hollon (2012). We conclude that 3 energy bins (or 3 analogous photoionization rates: molecular hydrogen photodissociation, neutral hydrogen HI, and fully ionized carbon CVI) are sufficient to characterize the dependence of the gas cooling and heating functions on our assumed incident radiation field model.

79 ASTRONOMY AND ASTROPHYSICS↗

Building an Efficient Cluster Cosmology Software Package for Modeling Cluster Counts and Lensing

We introduce a software suite developed for galaxy cluster cosmological analysis with the Dark Energy Survey Data. Cosmological analyses based on galaxy cluster number counts and weak-lensing measurements need efficient software infrastructure to explore an increasingly large parameter space, and account for various cosmological and astrophysical effects. Our software package is designed to model the cluster observables in a wide-field optical survey, including galaxy cluster counts, their averaged weak-lensing masses, or the cluster's averaged weak-lensing radial signals. To ensure maximum efficiency, this software package is developed in C++ in the CosmoSIS software framework, making use of the CUBA integration library. We also implement a testing and validation scheme to ensure the quality of the package. We demonstrate the effectiveness of this development by applying the software to the Dark Energy Survey Year 1 galaxy cluster cosmological data sets, and acquired cosmological constraints that are consistent with the fiducial Dark Energy Survey analysis.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

3D correlations in the Lyman-α forest from early DESI data

We present the first measurements of Lyman-α (Lyα) forest correlations using early data from the Dark Energy Spectroscopic Instrument (DESI). We measure the auto-correlation of Lyα absorption using 88 509 quasars at z > 2, and its cross-correlation with quasars using a further 147 899 tracer quasars at z ≳ 1.77. Then, we fit these correlations using a 13-parameter model based on linear perturbation theory and find that it provides a good description of the data across a broad range of scales. We detect the BAO peak with a signal-to-noise ratio of 3.8σ, and show that our measurements of the auto- and cross-correlations are fully-consistent with previous measurements by the Extended Baryon Oscillation Spectroscopic Survey (eBOSS). Even though we only use here a small fraction of the final DESI dataset, our uncertainties are only a factor of 1.7 larger than those from the final eBOSS measurement. As a result, we validate the existing analysis methods of Lyα correlations in preparation for making a robust measurement of the BAO scale with the first year of DESI data.

79 ASTRONOMY AND ASTROPHYSICS↗

DMTN-221: Periodicity Analysis in Alert Production

The baselined timeseries features to be computed in Alert Production include a Lomb-Scargle periodogram ran on two classes of variable systems: RR Lyrae and Eclipsing Binaries. Based on a simulated LSST-like cadence light curves taken from the Extended LSST Astronomical Time-series Classification Challenge (ELAsTiCC) we perform an end-to-end test to characterize the periodicity recovery on the Alert Production multi-band light curves. In both variable classes, we found that a single-band Lomb-Scargle implementation yields to a low fraction of recovered periods, with a significant preference on the simple periodic phenomena such as RR Lyrae. We also investigated the results from a multi-band Lomb-Scargle and found an increased fraction of recovered periodicities above 15% for the eclipsing binaries, and over 80$\%$ for the RR Lyrae stars. Our findings suggest that a multi-band Lomb-Scargle should be implemented for searching periodic phenomena through AP. We also asses the computational and scientific performance of several configurations on simulated alert data and find that our current configuration scales linearly with the number of detections while assuming an heuristic frequency grid.

79 ASTRONOMY AND ASTROPHYSICS↗

(Un)Explained EMIC Waves: Understanding Quiet Time EMIC Wave Drivers

Electromagnetic ion cyclotron (EMIC) waves are known to be generated through cyclotron resonance with the local ion particle population and grow when there is a large enough temperature anisotropy. In general, these temperature anisotropies necessary for wave growth are found to be associated with either solar wind pressure pulses or particle injections during geomagnetic storms or substorms. However, some EMIC events do not show any clear association with these known drivers and appear unexplained. In our analysis of high-amplitude (>1 nT) non-storm time EMIC waves, we find that 24 out of the 223 (~11%) EMIC events with peak amplitude greater than 1 nT were not found to be associated with any clear EMIC wave driver. This raises a compelling question: What magnetospheric or solar wind driver provides the free energy to grow these quiet-time EMIC waves? Here, we examine two EMIC events on 13 April 2017 and 13 February 2014, which were excited during extremely quiet solar wind and geomagnetic conditions. Furthermore, an in-depth analysis of field and particle measurements from multiple data sets, including ground and in situ data for these two events, indicates that extremely weak and otherwise insignificant pressure values and/or very weak substorm injections occurring multiple hours before the event play a significant role in quiet time wave generation.

79 ASTRONOMY AND ASTROPHYSICS↗

Two years of optical and NIR observations of the superluminous supernova UID 30901 discovered by the UltraVISTA SN survey

ABSTRACT We present deep optical and near-infrared photometry of UID 30901, a superluminous supernova (SLSN) discovered during the UltraVISTA survey. The observations were obtained with VIRCAM (YJHKs) mounted on the VISTA telescope, DECam (griz) on the Blanco telescope, and SUBARU Hyper Suprime-Cam (HSC; grizy). These multiband observations comprise +700 d making UID 30901 one of the best photometrically followed SLSNe to date. The host galaxy of UID 30901 is detected in a deep HST F814W image with an AB magnitude of 27.3 ± 0.2. While no spectra exist for the SN or its host galaxy, we perform our analysis assuming z = 0.37, based on the photometric redshift of a possible host galaxy found at a projected distance of 7 kpc. Fitting a blackbody to the observations, the radius, temperature, and bolometric light curve are computed. We find a maximum bolometric luminosity of 5.4 ± 0.34 × 1043 erg s−1. A flattening in the light curve beyond 600 d is observed and several possible causes are discussed. We find the observations to clearly favour an SLSN type I, and plausible power sources such as the radioactive decay of 56Ni and the spin-down of a magnetar are compared to the data. We find that the magnetar model yields a good fit to the observations with the following parameters: a magnetic field B = 1.4 ± 0.3 × 1014 G, spin period of P = 6.0 ± 0.1 ms, and ejecta mass $M_{\mathrm{ ej}} = 11.9^{+4.8}_{-6.4} \,\mathrm{ M}_{\odot }$.

79 ASTRONOMY AND ASTROPHYSICS↗

Search for hadronic decays of feebly-interacting particles at NA62

The NA62 experiment at CERN has the capability to collect data in a beam-dump mode, where 400 GeV protons are dumped on an absorber. In this configuration, New Physics particles, including dark photons, dark scalars, and axion-like particles, may be produced in the absorber and decay in the instrumented volume beginning approximately 80 m downstream of the dump. A search for these particles decaying in flight to hadronic final states is reported, based on an analysis of a sample of 1.4 x 10 17 protons on dump collected in 2021. No evidence of a New Physics signal is observed, excluding new regions of parameter spaces of multiple models.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

HATS-47b, HATS-48Ab, HATS-49b, and HATS-72b: Four Warm Giant Planets Transiting K Dwarfs

We report the discovery of four transiting giant planets around K dwarfs. The planets HATS-47b, HATS-48Ab, HATS-49b, and HATS-72b have masses of 0.369{sub −0.021}{sup +0.031} M{sub J}, 0.243{sub −0.030}{sup +0.022} M{sub J}, 0.353{sub −0.027}{sup +0.038} M{sub J}, and 0.1254±0.0039 M{sub J}, respectively, and radii of 1.117±0.014 R{sub J}, 0.800±0.015 R{sub J}, 0.765±0.013 R{sub J}, and 0.7224±0.0032 R{sub J}, respectively. The planets orbit close to their host stars with orbital periods of 3.9228 days, 3.1317 days, 4.1480 days, and 7.3279 days, respectively. The hosts are main-sequence K dwarfs with masses of 0.674{sub −0.012}{sup +0.016} M{sub ⊙}, 0.7279 ±0.0066 M{sub ⊙}, 0.7133±0.0075 M{sub ⊙}, and 0.7311±0.0028, and with V-band magnitudes of V=14.829±0.010, 14.35±0.11, 14.998±0.040 and 12.469±0.010. The super-Neptune HATS-72b (a.k.a. WASP-191b and TOI 294.01) was independently identified as a transiting planet candidate by the HATSouth, WASP, and TESS surveys, and we present a combined analysis of all of the data gathered by each of these projects (and their follow-up programs). An exceptionally precise mass is measured for HATS-72b thanks to high-precision radial velocity (RV) measurements obtained with VLT/ESPRESSO, FEROS, HARPS, and Magellan/PFS. We also incorporate TESS observations of the warm Saturn–hosting systems HATS-47 (a.k.a. TOI 1073.01), HATS-48A, and HATS-49. HATS-47 was independently identified as a candidate by the TESS team, while the other two systems were not previously identified from the TESS data. The RV orbital variations are measured for these systems using Magellan/PFS. HATS-48A has a resolved [[Formula]] neighbor in Gaia DR2, which is a common-proper-motion binary star companion to HATS-48A with a mass of 0.22 M{sub ⊙} and a current projected physical separation of ∼1400 au.

79 ASTRONOMY AND ASTROPHYSICS↗

The Chandra view of the relation between X-ray and UV emission in quasars

Here, we present a study of the relation between X-rays and ultraviolet emission in quasars for a sample of broad-line, radio-quiet objects obtained from the cross-match of the Sloan Digital Sky Survey DR14 with the latest Chandra Source Catalog 2.0 (2332 quasars) and the Chandra COSMOS Legacy survey (273 quasars). The non-linear relation between the ultraviolet (at 2500 Å, L UV ) and the X-ray (at 2 keV, L X ) emission in quasars has been proved to be characterised by a smaller intrinsic dispersion than the observed one, as long as a homogeneous selection, aimed at preventing the inclusion of contaminants in the sample, is fulfilled. By leveraging on the low background of Chandra, we performed a complete spectral analysis of all the data available for the SDSS-CSC2.0 quasar sample (i.e. 3430 X-ray observations), with the main goal of reducing the uncertainties on the source properties (e.g. flux, spectral slope). We analysed whether any evolution of the L X – L UV relation exists by dividing the sample in narrow redshift intervals across the redshift range spanned by our sample, z ≃ 0.5–4. We find that the slope of the relation does not evolve with redshift and it is consistent with the literature value of 0.6 over the explored redshift range, implying that the mechanism underlying the coupling of the accretion disc and hot corona is the same at the different cosmic epochs. We also find that the dispersion decreases when examining the highest redshifts, where only pointed observations are available. These results further confirm that quasars are ‘standardisable candles’, that is we can reliably measure cosmological distances at high redshifts where very few cosmological probes are available.

79 ASTRONOMY AND ASTROPHYSICS↗

The cosmological analysis of the SDSS/BOSS data from the Effective Field Theory of Large-Scale Structure

The Effective Field Theory of Large-Scale Structure is a formalism that allows us to predict the clustering of Cosmological Large-Scale Structure in the mildly non-linear regime in an accurate and reliable way. In this paper, after validating our technique against several sets of numerical simulations, we perform the analysis for the cosmological parameters of the DR12 BOSS data. We assume Λ CDM, a fixed value of the baryon/dark-matter ratio, Ω b /Ω c , and of the tilt of the primordial power spectrum, n s , and no significant input from numerical simulations. By using the one-loop power spectrum multipoles, we measure the primordial amplitude of the power spectrum, A s , the abundance of matter, Ω m , and the Hubble parameter, H 0 , {to about 13% , 3.2% and 3.2% respectively, obtaining ln (10 10 A s )=2.72± 0.13 , 0Ω m =0.309± 0.01 , H 0 =68.5± 2.2 km/(s Mpc) at 68% confidence level. If we then add a CMB prior on the sound horizon, the error bar on H 0 is reduced to 1.6% .} These results are a substantial qualitative and quantitative improvement with respect to former analyses, and suggest that the EFTofLSS is a powerful instrument to extract cosmological information from Large-Scale Structure.

79 ASTRONOMY AND ASTROPHYSICS↗

Measuring H 0 using X-ray and SZ effect observations of dynamically relaxed galaxy clusters

ABSTRACT We use a sample of 14 massive, dynamically relaxed galaxy clusters to constrain the Hubble constant, H0, by combining X-ray and Sunyaev–Zel’dovich (SZ) effect signals measured with Chandra, Planck, and Bolocam. This is the first such analysis to marginalize over an empirical, data-driven prior on the overall accuracy of X-ray temperature measurements, while our restriction to the most relaxed, massive clusters also minimizes astrophysical systematics. For a cosmological-constant model with Ωm = 0.3 and ΩΛ = 0.7, we find $H_0 = 67.3^{+21.3}_{-13.3}\, \mathrm{km}\, \mathrm{s}^{-1}\, \mathrm{Mpc}^{-1}$, limited by the temperature calibration uncertainty (compared to the statistically limited constraint of $H_0 = 72.3^{+7.6}_{-7.6}\, \mathrm{km}\, \mathrm{s}^{-1}\, \mathrm{Mpc}^{-1}$). The intrinsic scatter in the X-ray/SZ pressure ratio is found to be 13 ± 4 per cent (10 ± 3 per cent when two clusters with significant galactic dust emission are removed from the sample), consistent with being primarily due to triaxiality and projection. We discuss the prospects for reducing the dominant systematic limitation to this analysis, with improved X-ray calibration and/or precise measurements of the relativistic SZ effect providing a plausible route to per cent level constraints on H0.

79 ASTRONOMY AND ASTROPHYSICS↗

Constructing angular distributions of neutrinos in core-collapse supernovae from zeroth and first moments calibrated by full Boltzmann neutrino transport

Two-moment neutrino transport methods have been widely used for developing theoretical models of core-collapse supernovae (CCSN), since they substantially reduce the computational burden inherent in the multidimensional neutrino-radiation hydrodynamical simulations. The approximation, however, comes at a price; the detailed structure of angular distribution of neutrinos is sacrificed, that is the main drawback of this approach. In this paper, we develop a novel method by which to construct angular distributions of neutrinos from the zeroth and first angular moments. In our method, the angular distribution is expressed with two quadratic functions of the neutrino angle in a piecewise fashion. We determine the best parameters in the fitting function by comparing to the neutrino data in a spherically symmetric CCSN model with full Boltzmann neutrino transport. We demonstrate the capability of our method by using our recent 2D CCSN model. In this work, we find that the essential features of the angular distributions can be well reconstructed, whereas the angular distributions of incoming neutrinos tend to have large errors that increase with flux factor ($\textit{κ}$). This issue originates from the insensitiveness of incoming neutrinos to $\textit{κ}$, that is an intrinsic limitation in moment methods. Based on the results of the demonstration, we assess the reliability of electron- neutrinos lepton number (ELN)-crossing searches with the two-moment neutrino transport. This analysis is complementary to our [L. Johns and H. Nagakura arXiv:2104.04106] that scrutinizes the limitation of crossing searches with a few moments. We find that the systematic errors of angular distributions for incoming neutrinos lead to misjudgements of the crossing at $\textit{κ}$ ≳ 0.5. This casts doubt on the results of ELN-crossing searches based on two-moment methods in some previous studies.

79 ASTRONOMY AND ASTROPHYSICS↗

Gas and Cosmic-Ray Properties in the MBM 53, 54, and 55 Molecular Clouds and the Pegasus Loops Revealed by H i Line Profiles, Dust, and Gamma-Ray Data

In studying the interstellar medium (ISM) and Galactic cosmic rays (CRs), uncertainty of the interstellar gas density has always been an issue. To overcome this difficulty, we used a component decomposition of the 21 cm H i emission line and used the resulting gas maps in an analysis of γ-ray data obtained by the Fermi Large Area Telescope (LAT) for the MBM 53, 54, and 55 molecular clouds and the Pegasus loop. We decomposed the ISM gas into intermediate-velocity clouds, narrow-line and optically thick H i, broad-line and optically thin H i, CO-bright H 2 , and CO-dark H 2 using detailed correlations with the H i line profiles from the HI4PI survey, the Planck dust-emission model, and the Fermi-LAT γ-ray data. We found the fractions of the optical depth correction to the H i column density and CO-dark H 2 to be nearly equal. We fitted the CR spectra directly measured at/near the Earth and the measured γ-ray emissivity spectrum simultaneously. We obtained a spectral break in the interstellar proton spectrum at ~7 GeV, and found that the γ-ray emissivity normalization agrees with the AMS-02 proton spectrum within 10%, relaxing the tension with the CR spectra previously claimed.

79 ASTRONOMY AND ASTROPHYSICS↗

DESI Emission-line Galaxies: Unveiling the Diversity of [O II ] Profiles and Its Links to Star Formation and Morphology

We study the [O II ] profiles of emission-line galaxies (ELGs) from the Early Data Release of the Dark Energy Spectroscopic Instrument (DESI). To this end, we decompose and classify the shape of [O II ] profiles with the first two eigenspectra derived from principal component analysis. Our results show that DESI ELGs have diverse line profiles, which can be categorized into three main types: (1) narrow lines with a median width of ∼50 km s −1 , (2) broad lines with a median width of ∼80 km s −1 , and (3) two redshift systems with a median velocity separation of ∼150 km s −1 , i.e., double-peak galaxies. To investigate the connections between the line profiles and galaxy properties, we utilize the information from the COSMOS data set and compare the properties of ELGs, including star formation rate (SFR) and galaxy morphology, with the average properties of reference star-forming galaxies with similar stellar mass, sizes, and redshifts. Our findings show that, on average, DESI ELGs have a higher SFR and more asymmetrical/disturbed morphology than the reference galaxies. Moreover, we uncover a relationship between the line profiles, the excess SFR, and the excess asymmetry parameter, showing that DESI ELGs with broader [O II ] line profiles have more disturbed morphology and higher SFR than the reference star-forming galaxies. Finally, we discuss possible physical mechanisms giving rise to the observed relationship and the implications of our findings on the galaxy clustering measurements, including the halo occupation distribution modeling of DESI ELGs and the observed excess velocity dispersion of the satellite ELGs.

79 ASTRONOMY AND ASTROPHYSICS↗