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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

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↗
DOE OSTI2020

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↗
NASA NTRS2024

Comparative Analysis of Type III Radio Bursts and Solar Flares: Spatial Localization and Correlation with Solar Flare Intensity

We present a comprehensive study of type III radio bursts and their association with solar flares of magnitude M1.0 and larger, as observed by four widely separated spacecraft (Parker Solar Probe, Solar Orbiter, STEREO-A, and Wind). Our main focus is the introduction and validation of two methods for localizing radio bursts using the available multispacecraft data. The first method utilizes intensity fitting with a circular Gaussian distribution, while the second method is based on the time arrival of radio bursts. We demonstrate the effectiveness of these methods through the analysis of a single type III burst event and compare their results with the traditional radio triangulation technique. Furthermore, we conduct a statistical study of 17 type III bursts associated with M- and X-class solar flares in years 2020–2022. Our findings suggest a possible correlation between solar flare intensities and longitudes, with east limb flares tending to be weaker than west limb flares. We also observe a systematic drift of radio burst longitudes toward the east, potentially explained by a poleward component of the local density gradient. Our results suggest a strong correlation between solar flare intensities and radio burst properties, enhancing our understanding of the relationship between solar flares and type III radio bursts.

Radio astronomy↗
DOE OSTI2021

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↗
DOE OSTI2021

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↗
DOE OSTI2022

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↗
DOE OSTI2024

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↗
DOE OSTI2023

BEYONDPLANCK XIII. Intensity foreground sampling, degeneracies, and priors

We present the intensity foreground algorithms and model employed within the BEYONDPLANCK analysis framework. The BEYONDPLANCK analysis is aimed at integrating component separation and instrumental parameter sampling within a global framework, leading to complete end-to-end error propagation in the Planck Low Frequency Instrument (LFI) data analysis. Given the scope of the BEYONDPLANCK analysis, a limited set of data is included in the component separation process, leading to foreground parameter degeneracies. In order to properly constrain the Galactic foreground parameters, we improve upon the previous Commander component separation implementation by adding a suite of algorithmic techniques. These algorithms are designed to improve the stability and computational efficiency for weakly constrained posterior distributions. These are: (1) joint foreground spectral parameter and amplitude sampling, building on ideas from MIRAMARE; (2) component-based monopole determination; (3) joint spectral parameter and monopole sampling; and (4) application of informative spatial priors for component amplitude maps. We find that the only spectral parameter with a significant signal-to-noise ratio using the current BEYONDPLANCK data set is the peak frequency of the anomalous microwave emission component, for which we find $ν$ p = 25.3 ± 0.5 GHz; all others must be constrained through external priors. Future works will be aimed at integrating many more data sets into this analysis, both map and time-ordered based, thereby gradually eliminating the currently observed degeneracies in a controlled manner with respect to both instrumental systematic effects and astrophysical degeneracies. When this happens, the simple LFI-oriented data model employed in the current work will need to be generalized to account for both a richer astrophysical model and additional instrumental effects. This work will be organized within the Open Science-based COSMOGLOBE community effort.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

S -factor and scattering-parameter extractions from ${}^{3}\mathrm{He}+{}^{4}\mathrm{He}{ \rightarrow }^{7}\mathrm{Be}+\gamma $

Previous studies of the reaction ${}^{3}\mathrm{He}+{}^{4}\mathrm{He}\to {}^{7}\mathrm{Be}+\gamma $ have mainly focused on providing the best central value and error bar for the S factor at solar energies. Experimental measurements of this capture reaction at higher energies, the ${}^{3}\mathrm{He}$–${}^{4}\mathrm{He}$ scattering phase shifts, as well as properties of ${}^{7}\mathrm{Be}$ and its excited state, have been used to constrain the theoretical models employed for this purpose. Here we show that much more information than was previously appreciated can be extracted from angle-integrated capture data alone. We use the next-to-leading-order (NLO) amplitude in an effective field theory (EFT) for ${}^{3}\mathrm{He}+{}^{4}\mathrm{He}\to {}^{7}\mathrm{Be}+\gamma $ to perform the extrapolation. At this order the EFT describes the capture process using an s-wave scattering length and effective range, the asymptotic properties of ${}^{7}\mathrm{Be}$ and its excited state, and short-distance contributions to the E1 capture amplitude. We extract the multi-dimensional posterior of all these parameters via a Bayesian analysis that uses capture data below 2 MeV. We find that properties of the ${}^{7}\mathrm{Be}$ ground and excited states are well constrained. The total S factor $S(0)\,={0.577}_{-0.016}^{+0.015}$ keV b, while the branching ratio for excited- to ground-state capture at zero energy, ${Br}(0)={0.406}_{-0.011}^{+0.013}$, both at 68% degree of belief. This S(0) is broadly consistent with other recent evaluations, and agrees with the previously recommended value $S(0)=0.56\pm 0.03\,\,\mathrm{keV}$ b, but has a smaller error bar. We also find significant constraints on ${}^{3}\mathrm{He}$–${}^{4}\mathrm{He}$ scattering parameters, and we obtain constraints on the angular distribution of capture gamma rays, which is important for interpreting experiments. The path forward for this reaction seems to lie with better measurements of the scattering phase shifts and S(E)'s angular dependence away from zero energy, together with better understanding of the asymptotic normalization coefficients of the ${}^{7}$Be bound states' wave functions. Lastly, data on these could further reduce the uncertainty on S(0).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

A near magnetic-to-kinetic energy equipartition flare from the relativistic jet in AO 0235 + 164 during 2013–2019

ABSTRACT We present the multiwavelength flaring activity of the blazar AO 0235 + 164 during its recent active period from 2013 to 2019. From a discrete correlation function analysis, we find a significant ($\gt 95~{{\ \rm per\ cent}}$) correlation between radio and gamma-ray light curves with flares at longer wavelengths following flares at shorter wavelengths. We identify a new jet component in 43 GHz Very Long Baseline Array data that was ejected from the radio core on MJD $57246^{+26}_{-30}$ (2015 August 12), during the peak of the 2015 radio flare. From the analysis of the jet component, we derived a Doppler factor of δvar = 28.5 ± 8.4, a bulk Lorentz factor of $\Gamma =16.8^{+3.6}_{-3.1}$, and an intrinsic viewing angle of $\theta _{\rm v}=1.42^{+1.07}_{-0.52}\textrm {~degrees}$. Investigation of the quasi-simultaneous radio data revealed a partially absorbed spectrum with the turnover frequency varying in the range of 10−70 GHz and the peak flux density varying in the range of 0.7−4 Jy. We find the synchrotron self-absorption magnetic field strength to be $B_{\rm SSA}=15.3^{+12.6}_{-14.0}\textrm {~mG}$ at the peak of the 2015 radio flare, which is comparable to the equipartition magnetic field strength of $B_{\rm EQ}=43.6^{+10.6}_{-10.4}\textrm {~mG}$ calculated for the same epoch. Additional analysis of the radio emission region in the relativistic jet of AO 0235 + 164 suggests that it did not significantly deviate from equipartition during its recent flaring activity.

Astronomy & Astrophysics↗
DOE OSTI2021

Stability of Cool Cores during Galaxy Cluster Growth: A Joint Chandra/SPT Analysis of 67 Galaxy Clusters along a Common Evolutionary Track Spanning 9 Gyr

We present the results of a joint analysis of Chandra X-ray and South Pole Telescope (SPT) Sunyaev–Zel’dovich (SZ) observations targeting the first sample of galaxy clusters at 0.3 < z < 1.3, selected to be the progenitors of well-studied nearby clusters based on their expected accretion rate. We develop a new procedure in order to tackle the analysis challenge that is estimating the intracluster medium (ICM) properties of low-mass and high-redshift clusters with ~150 X-ray counts. One of the dominant sources of uncertainty on the ICM density profile estimated with a standard X-ray analysis with such shallow X-ray data is due to the systematic uncertainty associated with the ICM temperature obtained through the analysis of the background-dominated X-ray spectrum. We show that we can decrease the uncertainty on the density profile by a factor varying between 2 and 8 with a joint deprojection of the X-ray surface brightness profile measured by Chandra and the SZ-integrated Compton parameter available in the SPT cluster catalog. We apply this technique to the whole sample of 67 clusters in order to track the evolution of the ICM core density during cluster growth. We confirm that the evolution of the gas density profile is well modeled by the combination of a fixed core and a self-similarly evolving non-cool-core profile. We show that the fraction of cool cores in this sample is remarkably stable with redshift although clusters have gained a factor of ~4 in total mass over the past ~9 Gyr.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Comprehensive Analysis of the Neutrino Process in Core-collapsing Supernovae

We investigate the neutrino flavor change effects due to neutrino self-interaction and shock wave propagation, as well as the matter effects on the neutrino process in core-collapsing supernovae (CCSNe). For the hydrodynamics, we use two models: a simple thermal bomb model and a specified hydrodynamics model for SN1987A. For the presupernova model, we take an updated model, adjusted to explain SN1987A, which employs recent developments in the (n, γ) reaction rates for nuclei near the stability line (A ~ 100). As for the neutrino luminosity, we adopt two different models: equivalent neutrino luminosity and nonequivalent luminosity models. The latter is taken from a synthetic analysis of CCSN simulation data, which quantitatively presented the results obtained by various neutrino transport models. Relevant neutrino-induced reaction rates are calculated using a shell model for light nuclei and a quasiparticle random phase approximation model for heavy nuclei. For each model, we present abundances of the light nuclei ( 7 Li, 7 Be, 11 B, and 11 C) and the heavy nuclei ( 92 Nb, 98 Tc, 138 La, and 180 Ta) produced by the neutrino process. The light nuclei abundances turn out to be sensitive to the Mikheyev–Smirnov–Wolfenstein (MSW) region around O-Ne-Mg layer while the heavy nuclei are mainly produced prior to the MSW region. Through detailed analyses, we find that neutrino self-interaction becomes a key ingredient, in addition to the MSW effect, for understanding the neutrino process and the relevant nuclear abundances. The normal mass hierarchy is shown to be more compatible with the meteorite data. The main nuclear reactions for each nucleus are also investigated in detail.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

A More Precise Measurement of the Radius of PSR J0740+6620 Using Updated NICER Data

PSR J0740+6620 is the neutron star with the highest precisely determined mass, inferred from radio observations to be 2.08 ± 0.07 M ⊙ . Measurements of its radius therefore hold promise to constrain the properties of the cold, catalyzed, high-density matter in neutron star cores. Previously, Miller et al. and Riley et al. reported measurements of the radius of PSR J0740+6620 based on Neutron Star Interior Composition Explorer (NICER) observations accumulated through 2020 April 17, and an exploratory analysis utilizing NICER background estimates and a data set accumulated through 2021 December 28 was presented in Salmi et al. Here we report an updated radius measurement, derived by fitting models of X-ray emission from the neutron star surface to NICER data accumulated through 2022 April 21, totaling ~1.1 Ms additional exposure compared to the data set analyzed in Miller et al. and Riley et al., and to data from XMM-Newton observations. We find that the equatorial circumferential radius of PSR J0740+6620 is ${12.92}_{-1.13}^{+2.09}$ km (68% credibility), a fractional uncertainty ~83% the width of that reported in Miller et al., in line with statistical expectations given the additional data. If we were to require the radius to be less than 16 km, as was done in Salmi et al., then our 68% credible region would become $R={12.76}_{-1.02}^{+1.49}$ km, which is close to the headline result of Salmi et al. Our updated measurements, along with other laboratory and astrophysical constraints, imply a slightly softer equation of state than that inferred from our previous measurements.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

A comparison of compact, presumably young with extended, evolved radio active galactic nuclei

Context.The triggering and evolution of active galactic nuclei (AGNs) and the interaction of the AGN with its host galaxy is an important topic in extragalactic astrophysics. Radio sources with peaked spectra (peaked spectrum sources, PSS) and compact symmetric objects (CSO) are powerful, compact, and presumably young AGNs and therefore particularly suitable to study aspects of the AGN-host connection. Aims.We use a statistical approach to investigate properties of a PSS-CSO sample that are related to host galaxies and could potentially shed light on the link between host galaxies and AGNs. The main goal is to compare the PSS-CSO sample with a matching comparison sample of extended sources (ECS) to see if the two have significant differences. Methods.We analysed composite spectra, diagnostic line diagrams, multi-band spectral energy distributions (MBSEDs), star formation (SF) indicators, morphologies, and cluster environments for a sample of 121 PSSs and CSOs for which spectra are available from the Sloan Digital Sky Survey (SDSS). The statistical results were compared with those of the ECS sample, where we generally considered the two subsamples of quasi-stellar objects (QSOs) and radio galaxies separately. The analysis is based on a large set of archival data in the spectral range from the ultraviolet to mid-infrared. Results.We find significant differences between the PSS-CSO and the ECS sample. In particular, we find that the ECS sample has a higher proportion of passive galaxies with a lower star formation activity. This applies to both sub-samples (QSOs or radio galaxies) as well as to the entire sample. The star formation rates of the PSS-CSO host galaxies with available data are typically in the range ∼0 to 5 ℳ ⊙ yr −1 , and the stellar masses are in the range 3 × 10 11 to 10 12 ℳ ⊙ . Secondly, in agreement with previous results, we find a remarkably high proportion of PSS-CSO host galaxies with merger signatures. The merger fraction of the PSS-CSO sample is 0.61 ± 0.07, which is significantly higher than that of the comparison sample (0.15 ± 0.06). We suggest that this difference can be explained by assuming that the majority of the PSSs and CSOs cannot evolve to extended radio sources and are therefore not represented in our comparison sample.

Astronomy & Astrophysics↗
DOE OSTI2023

Morphological analysis of the polarized synchrotron emission with WMAP and Planck

The bright polarized synchrotron emission, away from the Galactic plane, originates mostly from filamentary structures. We implement a filament finder algorithm which allows the detection of bright elongated structures in polarized intensity maps. We analyse the sky at 23 and 30 GHz as observed respectively by WMAP and Planck. We identify 19 filaments, 13 of which have been previously observed. For each filament, we study the polarization fraction, finding values typically larger than for the areas outside the filaments, excluding the Galactic plane, and a fraction of about 30% is reached in two filaments. We study the polarization spectral indices of the filaments, and find a spectral index consistent with the values found in previous analysis (about -3.1) for more diffuse regions. Decomposing the polarization signals into the E and B families, we find that most of the filaments are detected in P E , but not in P B . We then focus on understanding the statistical properties of the diffuse regions of the synchrotron emission at 23 GHz. Using Minkowski functionals and tensors, we analyse the non-Gaussianity and statistical isotropy of the polarized intensity maps. For a sky coverage corresponding to 80% of the fainter emission, and on scales smaller than 6 degrees (ℓ > 30), the deviations from Gaussianity and isotropy are significantly higher than 3σ. The level of deviation decreases for smaller scales, however, it remains significantly high for the lowest analised scale (~ 1.5°). When 60% sky coverage is analysed, we find that the deviations never exceed 3σ. Finally, we present a simple data-driven model to generate non-Gaussian and anisotropic simulations of the synchrotron polarized emission. The simulations are fitted in order to match the spectral and statistical properties of the faintest 80% sky coverage of the data maps.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Constraining dark matter microphysics with the annihilation signal from subhalos

In the cold dark matter scenario, galactic dark matter halos are populated with a large number of smaller subhalos. Previous work has shown that dark matter annihilations in subhalos can generate a distinctive, non-Poisson signal in the gamma-ray photon counts probability distribution function (PDF). Here we show that the gamma-ray PDF also carries information about the velocity dependence of the dark matter annihilation cross section. After calculating the PDF assuming s-wave and Sommerfeld-enhanced annihilation, we perform a mock data analysis to illustrate how current and future observations can constrain the microphysics of the dark matter annihilation. We find that, with current Fermi data, and assuming a dark matter annihilation cross section roughly at the limit of current bounds from annihilation in dwarf spheroidal galaxies, one can potentially distinguish the non-Poissonian fluctuations expected from dark matter annihilation in subhalos from Poisson sources, as well as from dark matter models with an incorrect velocity-dependence. Here we explore how robust these results are to assumptions about the modeling of the galactic gamma-ray background, but further work is needed to determine the impact of realistic astrophysical source populations on our results. We also point out a four-parameter degeneracy between the velocity dependence of the dark matter annihilation, the minimum subhalo mass, the power law index of the subhalo mass function, and the normalization of the dark matter signal. This degeneracy can be broken with priors from N-body simulations or from observational constraints on the subhalo mass function.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Fuzzy dark matter and the Dark Energy Survey Year 1 data

ABSTRACT Gravitational weak lensing by dark matter haloes leads to a measurable imprint in the shear correlation function of galaxies. Fuzzy dark matter (FDM), composed of ultralight axion-like particles of mass m ∼ 10−22 eV, suppresses the matter power spectrum and shear correlation with respect to standard cold dark matter. We model the effect of FDM on cosmic shear using the optimized halo model HMCode, accounting for additional suppression of the mass function and halo concentration in FDM as observed in N-body simulations. We combine Dark Energy Survey Year 1 (DES-Y1) data with the Planck cosmic microwave background anisotropies to search for shear correlation suppression caused by FDM. We find no evidence of suppression compared to the preferred cold dark matter model, and thus set a new lower limit to the FDM particle mass. Using a log-flat prior and marginalizing over uncertainties related to the non-linear model of FDM, we find a new, independent 95 per cent C.L. lower limit log10m > −23 combining Planck and DES-Y1 shear, an improvement of almost two orders of magnitude on the mass bound relative to CMB-only constraints. Our analysis is largely independent of baryonic modelling, and of previous limits to FDM covering this mass range. Our analysis highlights the most important aspects of the FDM non-linear model for future investigation. The limit to FDM from weak lensing could be improved by up to three orders of magnitude with $\mathcal {O}(0.1)$ arcmin cosmic shear angular resolution, if FDM and baryonic feedback can be simultaneously modelled to high precision in the halo model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Core-collapse contamination in photometric samples of Type Ia Supernovae

This is an exciting time for cosmology with type Ia supernovae (SNe Ia). The recentlyconcluded Dark Energy Survey SN programme (DES-SN) has obtained the largest anddeepest high-redshift cosmological SN Ia sample, and the Vera Rubin Observatory isexpected to observe at least one order of magnitude more SNe Ia in the next decade.In both these experiments, only a limited fraction (.10 per cent) of the SNe can bespectroscopically classified. This leaves us with large ‘photometric’ SN samples, withthe potential for significant contamination by core-collapse SNe that may bias SN Iacosmological measurements. This thesis demonstrates how this contamination can bemodelled and accounted for in current and future cosmological analyses.First, I present state-of-the-art simulations of the SN universe. These are designedto accurately model the population of SNe Ia, peculiar SNe Ia and core-collapse SNe, aswell as their host galaxies. To improve the diversity and quality of the simulated corecollapseSNe, I build a new library of core-collapse SN templates using spectroscopicand photometric (optical and near-ultraviolet) data of 67 core-collapse SNe from theliterature. I account for our incomplete knowledge of core-collapse SN properties bygenerating a set of SN simulations (rather than a single one), each exploring differentmodelling choices and template libraries. I then characterise selection effects in theDES-SN survey and incorporate them in the simulations, thus obtaining a series ofDES-like simulated SN samples that can be compared to the observed DES-SN data.The agreement between the simulations and data is excellent across many observed SNproperties, including Hubble residuals. These simulations are the first to reproduce theobserved photometric SN and host galaxy properties in high-redshift surveys with no fine-tuning of the input parameters.I use my simulation framework to train and test the performance of SuperNNova,a photometric SN classifier based on recurrent neural networks. I explore differenttraining and validation strategies and show that, across all the DES-SN simulationstested, SuperNNova reduces core-collapse SN contamination to 0.8–3.5 per cent. Ithen show that biases due to contamination on the equation-of-state of dark energy,w, are < 0:008 when using our reference SuperNNova model. This compares to anexpected statistical uncertainty on w from the DES-SN sample of 0:039, thus showing that contamination is not a limiting systematic for the cosmological analysis of theDES-SN sample.The results presented in this thesis are the foundation of the DES SN Ia cosmologicalanalysis; they also provide important implications for the future of SN cosmology,as they demonstrate that contamination is not expected to significantly degrade thecosmological figure of merit of the Rubin SN Ia analysis.

79 ASTRONOMY AND ASTROPHYSICS↗