Source author record

K. Bechtol

K. Bechtol appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

81works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

81 published item(s)

preprint2026arXiv

Biasing from galaxy trough and peak profiles with the DES Y3 redMaGiC galaxies and the weak lensing mass map

We measure the correspondence between the distribution of galaxies and matter around troughs and peaks in the projected galaxy density, by comparing \texttt{redMaGiC} galaxies ($0.15<z<0.65$) to weak lensing mass maps from the Dark Energy Survey (DES) Y3 data release. We obtain stacked profiles, as a function of angle $θ$, of the galaxy density contrast $δ_{\rm g}$ and the weak lensing convergence $κ$, in the vicinity of these identified troughs and peaks, referred to as `void' and `cluster' superstructures. The ratio of the profiles depend mildly on $θ$, indicating good consistency between the profile shapes. We model the amplitude of this ratio using a function $F(η, θ)$ that depends on cosmological parameters $η$, scaled by the galaxy bias. We construct templates of $F(η, θ)$ using a suite of $N$-body (`Gower Street') simulations forward-modelled with DES Y3-like noise and systematics. We discuss and quantify the caveats of using a linear bias model to create galaxy maps from the simulation dark matter shells. We measure the galaxy bias in three lens tomographic bins (near to far): $2.32^{+0.86}_{-0.27}, 2.18^{+0.86}_{-0.23}, 1.86^{+0.82}_{-0.23}$ for voids, and $2.46^{+0.73}_{-0.27}, 3.55^{+0.96}_{-0.55}, 4.27^{+0.36}_{-1.14}$ for clusters, assuming the best-fit \textit{Planck} cosmology. Similar values with $\sim0.1σ$ shifts are obtained assuming the mean DES Y3 cosmology. The biases from troughs and peaks are broadly consistent, although a larger bias is derived for peaks, which is also larger than those measured from the DES Y3 $3\times2$-point analysis. This method shows an interesting avenue for measuring field-level bias that can be applied to future lensing surveys.

preprint2026arXiv

Constraining the Stellar-to-Halo Mass Relation with Galaxy Clustering and Weak Lensing from DES Year 3 Data

We develop a framework to study the relation between the stellar mass of a galaxy and the total mass of its host dark matter halo using galaxy clustering and galaxy-galaxy lensing measurements. We model a wide range of scales, roughly from $\sim 100 \; {\rm kpc}$ to $\sim 100 \; {\rm Mpc}$, using a theoretical framework based on the Halo Occupation Distribution and data from Year 3 of the Dark Energy Survey (DES) dataset. The new advances of this work include: 1) the generation and validation of a new stellar mass-selected galaxy sample in the range of $\log M_\star/M_\odot \sim 9.6$ to $\sim 11.5$; 2) the joint-modeling framework of galaxy clustering and galaxy-galaxy lensing that is able to describe our stellar mass-selected sample deep into the 1-halo regime; and 3) stellar-to-halo mass relation (SHMR) constraints from this dataset. In general, our SHMR constraints agree well with existing literature with various weak lensing measurements. We constrain the free parameters in the SHMR functional form $\log M_\star (M_h) = \log(εM_1) + f\left[ \log\left( M_h / M_1 \right) \right] - f(0)$, with $f(x) \equiv -\log(10^{αx}+1) + δ[\log(1+\exp(x))]^γ/ [1+\exp(10^{-x})]$, to be $\log M_1 = 11.506^{+0.325}_{-0.404}$, $\log ε= -1.632^{+0.306}_{-0.181}$, $α= -1.638^{+0.108}_{-0.099}$, $γ= 0.596^{+0.251}_{-0.210}$ and $δ= 3.810^{+2.045}_{-1.811}$. The inferred average satellite fraction is within $\sim 5-35\%$ for our fiducial results and we do not see any clear trends with redshift or stellar mass. Furthermore, we find that the inferred average galaxy bias values follow the generally expected trends with stellar mass and redshift. Our study is the first SHMR in DES in this mass range, and we expect the stellar mass sample to be of general interest for other science cases.

preprint2026arXiv

Dark Energy Survey Year 3 Results: Cosmological constraints from second and third-order shear statistics

We present a cosmological analysis of the third-order aperture mass statistic using Dark Energy Survey Year 3 (DES Y3) data. We perform a complete tomographic measurement of the three-point correlation function of the Y3 weak lensing shape catalog with the four fiducial source redshift bins. Building upon our companion methodology paper, we apply a pipeline that combines the two-point function $ξ_{\pm}$ with the mass aperture skewness statistic $\langle M_{\rm ap}^3\rangle$, which is an efficient compression of the full shear three-point function. We use a suite of simulated shear maps to obtain a joint covariance matrix. By jointly analyzing $ξ_\pm$ and $\langle M_{\rm ap}^3\rangle$ measured from DES Y3 data with a $Λ$CDM model, we find $S_8=0.780\pm0.015$ and $Ω_{\rm m}=0.266^{+0.039}_{-0.040}$, yielding 111% of figure-of-merit improvement in $Ω_m$-$S_8$ plane relative to $ξ_{\pm}$ alone, consistent with expectations from simulated likelihood analyses. With a $w$CDM model, we find $S_8=0.749^{+0.027}_{-0.026}$ and $w_0=-1.39\pm 0.31$, which gives an improvement of $22\%$ on the joint $S_8$-$w_0$ constraint. Our results are consistent with $w_0=-1$. Our new constraints are compared to CMB data from the Planck satellite, and we find that with the inclusion of $\langle M_{\rm ap}^3\rangle$ the existing tension between the data sets is at the level of $2.3σ$. We show that the third-order statistic enables us to self-calibrate the mean photometric redshift uncertainty parameter of the highest redshift bin with little degradation in the figure of merit. Our results demonstrate the constraining power of higher-order lensing statistics and establish $\langle M_{\rm ap}^3\rangle$ as a practical observable for joint analyses in current and future surveys.

preprint2022arXiv

A DECam Search for Explosive Optical Transients Associated with IceCube Neutrinos

In this work, we investigate the likelihood of association between realtime, TeV-PeV energy neutrino alerts from IceCube and optical counterparts in the form of core-collapse supernovae (CC SNe). The optical follow-up of IceCube alerts requires two main instrumental capabilities: (1) deep imaging, since 73\% of neutrinos would come from CC SNe at redshifts $z > 0.3$, and (2) a large field of view (FoV), since typical IceCube muon neutrino pointing accuracy is on the order of $\sim1$~deg. With Blanco/DECam ($gri$ to 24th magnitude and $2.2$~deg diameter FoV), we performed a triggered optical follow-up observation of two IceCube alerts, IC170922A and IC171106A on $\sim6$~nights during the $\sim3$~weeks following each alert. For the IC170922A (IC171106A) follow-up observations, we expect that 12.1\% (9.5\%) of coincident CC SNe at $z \lesssim 0.3$ are detectable, and that on average, 0.23 (0.07) unassociated SNe in the neutrino 90\% containment regions also pass our selection criteria. We find two candidate CC SNe that are temporally coincident with the neutrino alerts in the FoV, but none in the 90\% containment regions, which is statistically consistent with expected rates of background CC SNe for these observations. If CC SNe are the dominant source of TeV-PeV neutrinos, we would expect an excess of coincident CC SNe to be detectable at the $3σ$ confidence level using DECam observations similar to those of this work for $\sim60$ ($\sim200$) neutrino alerts with (without) redshift information for all candidates.

preprint2022arXiv

Constraints on the Physical Properties of GW190814 through Simulations based on DECam Follow-up Observations by the Dark Energy Survey

On 14 August 2019, the LIGO and Virgo Collaborations detected gravitational waves from a black hole and a 2.6 solar mass compact object, possibly the first neutron star -- black hole (NSBH) merger. In search of an optical counterpart, the Dark Energy Survey (DES) obtained deep imaging of the entire 90 percent confidence level localization area with Blanco/DECam 0, 1, 2, 3, 6, and 16 nights after the merger. Objects with varying brightness were detected by the DES Pipeline and we systematically reduced the candidate counterparts through catalog matching, light curve properties, host-galaxy photometric redshifts, SOAR spectroscopic follow-up observations, and machine-learning-based photometric classification. All candidates were rejected as counterparts to the merger. To quantify the sensitivity of our search, we applied our selection criteria to full light curve simulations of supernovae and kilonovae as they would appear in the DECam observations. Since the source class of the merger was uncertain, we utilized an agnostic, three-component kilonova model based on tidally-disrupted NS ejecta properties to quantify our detection efficiency of a counterpart if the merger included a NS. We find that if a kilonova occurred during this merger, configurations where the ejected matter is greater than 0.07 solar masses, has lanthanide abundance less than $10^{-8.56}$, and has a velocity between $0.18c$ and $0.21c$ are disfavored at the $2σ$ level. Furthermore, we estimate that our background reduction methods are capable of associating gravitational wave signals with a detected electromagnetic counterpart at the $4σ$ level in $95\%$ of future follow-up observations.

preprint2022arXiv

Dark Energy Survey Year 3 Results: Calibration of Lens Sample Redshift Distributions using Clustering Redshifts with BOSS/eBOSS

We present clustering redshift measurements for Dark Energy Survey (DES) lens sample galaxies to be used in weak gravitational lensing and galaxy clustering studies. To perform this measurement, we cross-correlate with spectroscopic galaxies from the Baryon Acoustic Oscillation Survey (BOSS) and its extension, eBOSS. We validate our methodology in simulations, including a new technique to calibrate systematic errors due to the galaxy clustering bias, finding our method to be generally unbiased in calibrating the mean redshift. We apply our method to the data, and estimate the redshift distribution for eleven different photometrically-selected bins. We find general agreement between clustering redshift and photometric redshift estimates, with differences on the inferred mean redshift to be below $|Δz|=0.01$ in most of the bins. We also test a method to calibrate a width parameter for redshift distributions, which we found necessary to use for some of our samples. Our typical uncertainties on the mean redshift ranged from 0.003 to 0.008, while our uncertainties on the width ranged from 4 to 9\%. We discuss how these results calibrate the photometric redshift distributions used in companion DES Year 3 Results papers.

preprint2022arXiv

Dark Energy Survey Year 3 Results: Cosmological Constraints from Galaxy Clustering and Weak Lensing

We present the first cosmology results from large-scale structure in the Dark Energy Survey (DES) spanning 5000 deg$^2$. We perform an analysis combining three two-point correlation functions (3$\times$2pt): (i) cosmic shear using 100 million source galaxies, (ii) galaxy clustering, and (iii) the cross-correlation of source galaxy shear with lens galaxy positions. The analysis was designed to mitigate confirmation or observer bias; we describe specific changes made to the lens galaxy sample following unblinding of the results. We model the data within the flat $Λ$CDM and $w$CDM cosmological models. We find consistent cosmological results between the three two-point correlation functions; their combination yields clustering amplitude $S_8=0.776^{+0.017}_{-0.017}$ and matter density $Ω_{\mathrm{m}} = 0.339^{+0.032}_{-0.031}$ in $Λ$CDM, mean with 68% confidence limits; $S_8=0.775^{+0.026}_{-0.024}$, $Ω_{\mathrm{m}} = 0.352^{+0.035}_{-0.041}$, and dark energy equation-of-state parameter $w=-0.98^{+0.32}_{-0.20}$ in $w$CDM. This combination of DES data is consistent with the prediction of the model favored by the Planck 2018 cosmic microwave background (CMB) primary anisotropy data, which is quantified with a probability-to-exceed $p=0.13$ to $0.48$. When combining DES 3$\times$2pt data with available baryon acoustic oscillation, redshift-space distortion, and type Ia supernovae data, we find $p=0.34$. Combining all of these data sets with Planck CMB lensing yields joint parameter constraints of $S_8 = 0.812^{+0.008}_{-0.008}$, $Ω_{\mathrm{m}} = 0.306^{+0.004}_{-0.005}$, $h=0.680^{+0.004}_{-0.003}$, and $\sum m_ν<0.13 \;\mathrm{eV\; (95\% \;CL)}$ in $Λ$CDM; $S_8 = 0.812^{+0.008}_{-0.008}$, $Ω_{\mathrm{m}} = 0.302^{+0.006}_{-0.006}$, $h=0.687^{+0.006}_{-0.007}$, and $w=-1.031^{+0.030}_{-0.027}$ in $w$CDM. (abridged)

preprint2022arXiv

Dark Energy Survey Year 3 results: cosmological constraints from the analysis of cosmic shear in harmonic space

We present cosmological constraints from the analysis of angular power spectra of cosmic shear maps based on data from the first three years of observations by the Dark Energy Survey (DES Y3). Our measurements are based on the pseudo-$C_\ell$ method and offer a view complementary to that of the two-point correlation functions in real space, as the two estimators are known to compress and select Gaussian information in different ways, due to scale cuts. They may also be differently affected by systematic effects and theoretical uncertainties, such as baryons and intrinsic alignments (IA), making this analysis an important cross-check. In the context of $Λ$CDM, and using the same fiducial model as in the DES Y3 real space analysis, we find ${S_8 \equiv σ_8 \sqrt{Ω_{\rm m}/0.3} = 0.793^{+0.038}_{-0.025}}$, which further improves to ${S_8 = 0.784\pm 0.026 }$ when including shear ratios. This constraint is within expected statistical fluctuations from the real space analysis, and in agreement with DES~Y3 analyses of non-Gaussian statistics, but favors a slightly higher value of $S_8$, which reduces the tension with the Planck cosmic microwave background 2018 results from $2.3σ$ in the real space analysis to $1.5σ$ in this work. We explore less conservative IA models than the one adopted in our fiducial analysis, finding no clear preference for a more complex model. We also include small scales, using an increased Fourier mode cut-off up to $k_{\rm max}={5}{h{\rm Mpc}^{-1}}$, which allows to constrain baryonic feedback while leaving cosmological constraints essentially unchanged. Finally, we present an approximate reconstruction of the linear matter power spectrum at present time, which is found to be about 20\% lower than predicted by Planck 2018, as reflected by the $1.5σ$ lower $S_8$ value.

preprint2022arXiv

Dark Energy Survey Year 3 Results: Cosmology from Cosmic Shear and Robustness to Modeling Uncertainty

This work and its companion paper, Amon et al. (2021), present cosmic shear measurements and cosmological constraints from over 100 million source galaxies in the Dark Energy Survey (DES) Year 3 data. We constrain the lensing amplitude parameter $S_8\equivσ_8\sqrt{Ω_\textrm{m}/0.3}$ at the 3% level in $Λ$CDM: $S_8=0.759^{+0.025}_{-0.023}$ (68% CL). Our constraint is at the 2% level when using angular scale cuts that are optimized for the $Λ$CDM analysis: $S_8=0.772^{+0.018}_{-0.017}$ (68% CL). With cosmic shear alone, we find no statistically significant constraint on the dark energy equation-of-state parameter at our present statistical power. We carry out our analysis blind, and compare our measurement with constraints from two other contemporary weak-lensing experiments: the Kilo-Degree Survey (KiDS) and Hyper-Suprime Camera Subaru Strategic Program (HSC). We additionally quantify the agreement between our data and external constraints from the Cosmic Microwave Background (CMB). Our DES Y3 result under the assumption of $Λ$CDM is found to be in statistical agreement with Planck 2018, although favors a lower $S_8$ than the CMB-inferred value by $2.3σ$ (a $p$-value of 0.02). This paper explores the robustness of these cosmic shear results to modeling of intrinsic alignments, the matter power spectrum and baryonic physics. We additionally explore the statistical preference of our data for intrinsic alignment models of different complexity. The fiducial cosmic shear model is tested using synthetic data, and we report no biases greater than 0.3$σ$ in the plane of $S_8\timesΩ_\textrm{m}$ caused by uncertainties in the theoretical models.

preprint2022arXiv

Dark Energy Survey Year 3 results: cosmology with moments of weak lensing mass maps

We present a cosmological analysis using the second and third moments of the weak lensing mass (convergence) maps from the first three years of data (Y3) data of the Dark Energy Survey (DES). The survey spans an effective area of 4139 square degrees and uses the images of over 100 million galaxies to reconstruct the convergence field. The second moment of the convergence as a function of smoothing scale contains information similar to standard shear 2-point statistics. The third moment, or the skewness, contains additional non-Gaussian information. The data is analysed in the context of the $Λ$CDM model, varying 5 cosmological parameters and 19 nuisance parameters modelling astrophysical and measurement systematics. Our modelling of the observables is completely analytical, and has been tested with simulations in our previous methodology study. We obtain a 1.7\% measurement of the amplitude of fluctuations parameter $S_8\equiv σ_8 (Ω_m/0.3)^{0.5} = 0.784\pm 0.013$. The measurements are shown to be internally consistent across redshift bins, angular scales, and between second and third moments. In particular, the measured third moment is consistent with the expectation of gravitational clustering under the $Λ$CDM model. The addition of the third moment improves the constraints on $S_8$ and $Ω_{\rm m}$ by $\sim$15\% and $\sim$25\% compared to an analysis that only uses second moments. We compare our results with {\it Planck} constraints from the Cosmic Microwave Background (CMB), finding a $2.2$ \textendash $2.8σ$ tension in the full parameter space, depending on the combination of moments considered. The third moment independently is in $2.8σ$ tension with {\it Planck}, and thus provides a cross-check on analyses of 2-point correlations.

preprint2022arXiv

Dark Energy Survey Year 3 Results: Deep Field Optical + Near-Infrared Images and Catalogue

We describe the Dark Energy Survey (DES) Deep Fields, a set of images and associated multi-wavelength catalogue ($ugrizJHKs$) built from Dark Energy Camera (DECam) and Visible and Infrared Survey Telescope for Astronomy (VISTA) data. The DES Deep Fields comprise 11 fields (10 DES supernova fields plus COSMOS), with a total area of $\sim30~$ square degrees in $ugriz$ bands and reaching a maximum $i$-band depth of 26.75 (AB, $10σ$, 2 arcsec). We present a catalogue for the DES 3-year cosmology analysis of those four fields with full 8-band coverage, totalling $5.88~$ sq. deg. after masking. Numbering $2.8~$million objects ($1.6~$million post masking), our catalogue is drawn from images coadded to consistent depths of $r=25.7, i=25, z=24.3$ mag. We use a new model-fitting code, built upon established methods, to deblend sources and ensure consistent colours across the $u$-band to $Ks$-band wavelength range. We further detail the tight control we maintain over the point-spread function modelling required for the model fitting, astrometry and consistency of photometry between the four fields. The catalogue allows us to perform a careful star-galaxy separation and produces excellent photometric redshift performance (${\rm NMAD} = 0.023$ at $i<23$). The Deep-Fields catalogue will be made available as part of the cosmology data products release, following the completion of the DES 3-year weak lensing and galaxy clustering cosmology work.

preprint2022arXiv

Dark Energy Survey Year 3 Results: Exploiting small-scale information with lensing shear ratios

Using the first three years of data from the Dark Energy Survey, we use ratios of small-scale galaxy-galaxy lensing measurements around the same lens sample to constrain source redshift uncertainties, intrinsic alignments and other nuisance parameters of our model. Instead of using a simple geometric approach for the ratios, we use the full modeling of the galaxy-galaxy lensing measurements, including the corresponding integration over the power spectrum and the contributions from intrinsic alignments and lens magnification. We perform extensive testing of the small-scale shear ratio (SR) modeling by studying the impact of different effects such as the inclusion of baryonic physics, non-linear biasing, halo occupation distribution descriptions and lens magnification, among others, and using realistic $N$-body simulations. We validate the robustness of our constraints in the data by using two independent lens samples, and by deriving constraints using the corresponding large-scale ratios for which the modeling is simpler. The DES Y3 results demonstrate how the ratios provide significant improvements in constraining power for several nuisance parameters in our model, especially on source redshift calibration and intrinsic alignments (IA). For source redshifts, SR improves the constraints from the prior by up to 38\% in some redshift bins. Such improvements, and especially the constraints it provides on IA, translate to tighter cosmological constraints when SR is combined with cosmic shear and other 2pt functions. In particular, for the DES Y3 data, SR improves $S_8$ constraints from cosmic shear by up to 31\%, and for the full combination of probes (3$\times$2pt) by up to 10\%. The shear ratios presented in this work are used as an additional likelihood for cosmic shear, 2$\times$2pt and the full 3$\times$2pt in the fiducial DES Y3 cosmological analysis.

preprint2022arXiv

Dark Energy Survey Year 3 results: Galaxy-halo connection from galaxy-galaxy lensing

Galaxy-galaxy lensing is a powerful probe of the connection between galaxies and their host dark matter halos, which is important both for galaxy evolution and cosmology. We extend the measurement and modeling of the galaxy-galaxy lensing signal in the recent Dark Energy Survey Year 3 cosmology analysis to the highly nonlinear scales ($\sim 100$ kpc). This extension enables us to study the galaxy-halo connection via a Halo Occupation Distribution (HOD) framework for the two lens samples used in the cosmology analysis: a luminous red galaxy sample (redMaGiC) and a magnitude-limited galaxy sample (MagLim). We find that redMaGiC (MagLim) galaxies typically live in dark matter halos of mass $\log_{10}(M_{h}/M_{\odot}) \approx 13.7$ which is roughly constant over redshift ($13.3-13.5$ depending on redshift). We constrain these masses to $\sim 15\%$, approximately $1.5$ times improvement over previous work. We also constrain the linear galaxy bias more than 5 times better than what is inferred by the cosmological scales only. We find the satellite fraction for redMaGiC (MagLim) to be $\sim 0.1-0.2$ ($0.1-0.3$) with no clear trend in redshift. Our constraints on these halo properties are broadly consistent with other available estimates from previous work, large-scale constraints and simulations. The framework built in this paper will be used for future HOD studies with other galaxy samples and extensions for cosmological analyses.

preprint2022arXiv

Dark Energy Survey Year 3 Results: High-precision measurement and modeling of galaxy-galaxy lensing

We present and characterize the galaxy-galaxy lensing signal measured using the first three years of data from the Dark Energy Survey (DES Y3) covering 4132 deg$^2$. These galaxy-galaxy measurements are used in the DES Y3 3$\times$2pt cosmological analysis, which combines weak lensing and galaxy clustering information. We use two lens samples: a magnitude-limited sample and the redMaGic sample, which span the redshift range $\sim 0.2-1$ with 10.7 M and 2.6 M galaxies respectively. For the source catalog, we use the Metacalibration shape sample, consisting of $\simeq$100 M galaxies separated into 4 tomographic bins. Our galaxy-galaxy lensing estimator is the mean tangential shear, for which we obtain a total S/N of $\sim$148 for MagLim ($\sim$120 for redMaGic), and $\sim$67 ($\sim$55) after applying the scale cuts of 6 Mpc/$h$. Thus we reach percent-level statistical precision, which requires that our modeling and systematic-error control be of comparable accuracy. The tangential shear model used in the 3$\times$2pt cosmological analysis includes lens magnification, a five-parameter intrinsic alignment model (TATT), marginalization over a point-mass to remove information from small scales and a linear galaxy bias model validated with higher-order terms. We explore the impact of these choices on the tangential shear observable and study the significance of effects not included in our model, such as reduced shear, source magnification and source clustering. We also test the robustness of our measurements to various observational and systematics effects, such as the impact of observing conditions, lens-source clustering, random-point subtraction, scale-dependent Metacalibration responses, PSF residuals, and B-modes.

preprint2022arXiv

Dark Energy Survey Year 3 results: imprints of cosmic voids and superclusters in the Planck CMB lensing map

The CMB lensing signal from cosmic voids and superclusters probes the growth of structure in the low-redshift cosmic web. In this analysis, we cross-correlated the Planck CMB lensing map with voids detected in the Dark Energy Survey Year 3 (Y3) data set ($\sim$5,000 deg$^{2}$), expanding on previous measurements that used Y1 catalogues ($\sim$1,300 deg$^{2}$). Given the increased statistical power compared to Y1 data, we report a $6.6σ$ detection of negative CMB convergence ($κ$) imprints using approximately 3,600 voids detected from a redMaGiC luminous red galaxy sample. However, the measured signal is lower than expected from the MICE N-body simulation that is based on the $Λ$CDM model (parameters $Ω_{\rm m} = 0.25$, $σ_8 = 0.8$), and the discrepancy is associated mostly with the void centre region. Considering the full void lensing profile, we fit an amplitude $A_κ=κ_{\rm DES}/κ_{\rm MICE}$ to a simulation-based template with fixed shape and found a moderate $2σ$ deviation in the signal with $A_κ\approx0.79\pm0.12$. We also examined the WebSky simulation that is based on a Planck 2018 $Λ$CDM cosmology, but the results were even less consistent given the slightly higher matter density fluctuations than in MICE. We then identified superclusters in the DES and the MICE catalogues, and detected their imprints at the $8.4σ$ level; again with a lower-than-expected $A_κ=0.84\pm0.10$ amplitude. The combination of voids and superclusters yields a $10.3σ$ detection with an $A_κ=0.82\pm0.08$ constraint on the CMB lensing amplitude, thus the overall signal is $2.3σ$ weaker than expected from MICE.

preprint2022arXiv

Dark Energy Survey Year 3 Results: Measuring the Survey Transfer Function with Balrog

We describe an updated calibration and diagnostic framework, Balrog, used to directly sample the selection and photometric biases of the Dark Energy Survey's (DES) Year 3 (Y3) dataset. We systematically inject onto the single-epoch images of a random 20% subset of the DES footprint an ensemble of nearly 30 million realistic galaxy models derived from DES Deep Field observations. These augmented images are analyzed in parallel with the original data to automatically inherit measurement systematics that are often too difficult to capture with traditional generative models. The resulting object catalog is a Monte Carlo sampling of the DES transfer function and is used as a powerful diagnostic and calibration tool for a variety of DES Y3 science, particularly for the calibration of the photometric redshifts of distant "source" galaxies and magnification biases of nearer "lens" galaxies. The recovered Balrog injections are shown to closely match the photometric property distributions of the Y3 GOLD catalog, particularly in color, and capture the number density fluctuations from observing conditions of the real data within 1% for a typical galaxy sample. We find that Y3 colors are extremely well calibrated, typically within ~1-8 millimagnitudes, but for a small subset of objects we detect significant magnitude biases correlated with large overestimates of the injected object size due to proximity effects and blending. We discuss approaches to extend the current methodology to capture more aspects of the transfer function and reach full coverage of the survey footprint for future analyses.

preprint2022arXiv

Dark Energy Survey Year 3 Results: Point-Spread Function Modeling

We introduce a new software package for modeling the point-spread function (PSF) of astronomical images, called Piff (PSFs In the Full FOV), which we apply to the first three years (known as Y3) of the Dark Energy Survey (DES) data. We describe the relevant details about the algorithms used by Piff to model the PSF, including how the PSF model varies across the field of view (FOV). Diagnostic results show that the systematic errors from the PSF modeling are very small over the range of scales that are important for the DES Y3 weak lensing analysis. In particular, the systematic errors from the PSF modeling are significantly smaller than the corresponding results from the DES year one (Y1) analysis. We also briefly describe some planned improvements to Piff that we expect to further reduce the modeling errors in future analyses.

preprint2022arXiv

Dark Energy Survey Year 3 Results: Three-Point Shear Correlations and Mass Aperture Moments

We present high signal-to-noise measurements of three-point shear correlations and the third moment of the mass aperture statistic using the first 3 years of data from the Dark Energy Survey. We additionally obtain the first measurements of the configuration and scale dependence of the four three-point shear correlations which carry cosmological information. With the third-order mass aperture statistic, we present tomographic measurements over angular scales of 4 to 60 arcminutes with a combined statistical significance of 15.0$σ$. Using the tomographic information and measuring also the second-order mass aperture, we additionally obtain a skewness parameter and its redshift evolution. We find that the amplitudes and scale-dependence of these shear 3pt functions are in qualitative agreement with measurements in a mock galaxy catalog based on N-body simulations, indicating promise for including them in future cosmological analyses. We validate our measurements by showing that B-modes, parity-violating contributions and PSF modeling uncertainties are negligible, and determine that the measured signals are likely to be of astrophysical and gravitational origin.

preprint2022arXiv

DeepZipper: A Novel Deep Learning Architecture for Lensed Supernovae Identification

Large-scale astronomical surveys have the potential to capture data on large numbers of strongly gravitationally lensed supernovae (LSNe). To facilitate timely analysis and spectroscopic follow-up before the supernova fades, an LSN needs to be identified soon after it begins. To quickly identify LSNe in optical survey datasets, we designed ZipperNet, a multi-branch deep neural network that combines convolutional layers (traditionally used for images) with long short-term memory (LSTM) layers (traditionally used for time series). We tested ZipperNet on the task of classifying objects from four categories -- no lens, galaxy-galaxy lens, lensed type Ia supernova, lensed core-collapse supernova -- within high-fidelity simulations of three cosmic survey data sets -- the Dark Energy Survey (DES), Rubin Observatory's Legacy Survey of Space and Time (LSST), and a Dark Energy Spectroscopic Instrument (DESI) imaging survey. Among our results, we find that for the LSST-like dataset, ZipperNet classifies LSNe with a receiver operating characteristic area under the curve of 0.97, predicts the spectroscopic type of the lensed supernovae with 79\% accuracy, and demonstrates similarly high performance for LSNe 1-2 epochs after first detection. We anticipate that a model like ZipperNet, which simultaneously incorporates spatial and temporal information, can play a significant role in the rapid identification of lensed transient systems in cosmic survey experiments.

preprint2022arXiv

Milky Way Satellite Census. IV. Constraints on Decaying Dark Matter from Observations of Milky Way Satellite Galaxies

We use a recent census of the Milky Way (MW) satellite galaxy population to constrain the lifetime of particle dark matter (DM). We consider two-body decaying dark matter (DDM) in which a heavy DM particle decays with lifetime $τ$ comparable to the age of the Universe to a lighter DM particle (with mass splitting $ε$) and to a dark radiation species. These decays impart a characteristic "kick velocity," $V_{\mathrm{kick}}=εc$, on the DM daughter particles, significantly depleting the DM content of low-mass subhalos and making them more susceptible to tidal disruption. We fit the suppression of the present-day DDM subhalo mass function (SHMF) as a function of $τ$ and $V_{\mathrm{kick}}$ using a suite of high-resolution zoom-in simulations of MW-mass halos, and we validate this model on new DDM simulations of systems specifically chosen to resemble the MW. We implement our DDM SHMF predictions in a forward model that incorporates inhomogeneities in the spatial distribution and detectability of MW satellites and uncertainties in the mapping between galaxies and DM halos, the properties of the MW system, and the disruption of subhalos by the MW disk using an empirical model for the galaxy--halo connection. By comparing to the observed MW satellite population, we conservatively exclude DDM models with $τ< 18\ \mathrm{Gyr}$ ($29\ \mathrm{Gyr}$) for $V_{\mathrm{kick}}=20\ \mathrm{km}\, \mathrm{s}^{-1}$ ($40\ \mathrm{km}\, \mathrm{s}^{-1}$) at $95\%$ confidence. These constraints are among the most stringent and robust small-scale structure limits on the DM particle lifetime and strongly disfavor DDM models that have been proposed to alleviate the Hubble and $S_8$ tensions.

preprint2022arXiv

The DECam Local Volume Exploration Survey Data Release 2

We present the second public data release (DR2) from the DECam Local Volume Exploration survey (DELVE). DELVE DR2 combines new DECam observations with archival DECam data from the Dark Energy Survey, the DECam Legacy Survey, and other DECam community programs. DELVE DR2 consists of ~160,000 exposures that cover >21,000 deg^2 of the high Galactic latitude (|b| > 10 deg) sky in four broadband optical/near-infrared filters (g, r, i, z). DELVE DR2 provides point-source and automatic aperture photometry for ~2.5 billion astronomical sources with a median 5σ point-source depth of g=24.3, r=23.9, i=23.5, and z=22.8 mag. A region of ~17,000 deg^2 has been imaged in all four filters, providing four-band photometric measurements for ~618 million astronomical sources. DELVE DR2 covers more than four times the area of the previous DELVE data release and contains roughly five times as many astronomical objects. DELVE DR2 is publicly available via the NOIRLab Astro Data Lab science platform.

preprint2022arXiv

The DES Bright Arcs Survey: Candidate Strongly Lensed Galaxy Systems from the Dark Energy Survey 5,000 Sq. Deg. Footprint

We report the combined results of eight searches for strong gravitational lens systems in the full 5,000 sq. deg. of Dark Energy Survey (DES) observations. The observations accumulated by the end of the third observing season fully covered the DES footprint in 5 filters (grizY), with an $i-$band limiting magnitude (at $10σ$) of 23.44. In four searches, a list of potential candidates was identified using a color and magnitude selection from the object catalogs created from the first three observing seasons. Three other searches were conducted at the locations of previously identified galaxy clusters. Cutout images of potential candidates were then visually scanned using an object viewer. An additional set of candidates came from a data-quality check of a subset of the color-coadd "tiles" created from the full DES six-season data set. A short list of the most promising strong lens candidates was then numerically ranked according to whether or not we judged them to be bona fide strong gravitational lens systems. These searches discovered a diverse set of 247 strong lens candidate systems, of which 81 are identified for the first time. We provide the coordinates, magnitudes, and photometric properties of the lens and source objects, and an estimate of the Einstein radius for 81 new systems and 166 previously reported. This catalog will be of use for selecting interesting systems for detailed follow-up, studies of galaxy cluster and group mass profiles, as well as a training/validation set for automated strong lens searches.

preprint2022arXiv

The SLAC T-510 experiment for radio emission from particle showers: detailed simulation study and interpretation

Over the last several decades, radio detection of air showers has been widely used to detect ultra-high-energy cosmic rays. We developed an experiment under controlled laboratory conditions at SLAC with which we measured the radio-frequency radiation from a charged particle shower produced by bunches of electrons as primaries with known energy. The shower took place in a target made of High Density Polyethylene located in a strong magnetic field. The experiment was designed so that Askaryan and magnetically-induced components of the radio emission could be measured independently. At the same time, we performed a detailed simulation of this experiment to predict the radio signal using two microscopic formalisms, endpoint and ZHS. In this paper, we present the simulation scheme and make a comparison with data characteristics such as linearity with magnetic field and amplitude. The simulations agree with the measurements within uncertainties and present a good description of the data. In particular, reflections within the target that accounted for the largest systematic uncertainties are addressed. The prediction of the amplitude of Askaryan emission agrees with measurements to within 5% for the endpoint formalism and 11% for the ZHS formalism. The amplitudes of magnetically-induced emission agree to within 5% for the endpoint formalism and less than 1% for the ZHS formalism. The agreement of the absolute scale of emission gives confidence in state-of-the-art air shower simulations which are based on the applied formalisms.

preprint2021arXiv

A machine learning approach to galaxy properties: joint redshift-stellar mass probability distributions with Random Forest

We demonstrate that highly accurate joint redshift-stellar mass probability distribution functions (PDFs) can be obtained using the Random Forest (RF) machine learning (ML) algorithm, even with few photometric bands available. As an example, we use the Dark Energy Survey (DES), combined with the COSMOS2015 catalogue for redshifts and stellar masses. We build two ML models: one containing deep photometry in the $griz$ bands, and the second reflecting the photometric scatter present in the main DES survey, with carefully constructed representative training data in each case. We validate our joint PDFs for $10,699$ test galaxies by utilizing the copula probability integral transform and the Kendall distribution function, and their univariate counterparts to validate the marginals. Benchmarked against a basic set-up of the template-fitting code BAGPIPES, our ML-based method outperforms template fitting on all of our predefined performance metrics. In addition to accuracy, the RF is extremely fast, able to compute joint PDFs for a million galaxies in just under $6$ min with consumer computer hardware. Such speed enables PDFs to be derived in real time within analysis codes, solving potential storage issues. As part of this work we have developed GALPRO, a highly intuitive and efficient Python package to rapidly generate multivariate PDFs on-the-fly. GALPRO is documented and available for researchers to use in their cosmology and galaxy evolution studies.

preprint2021arXiv

Dark Energy Survey Year 3 results: cosmology from combined galaxy clustering and lensing -- validation on cosmological simulations

We present a validation of the Dark Energy Survey Year 3 (DES Y3) $3\times2$-point analysis choices by testing them on Buzzard v2.0, a new suite of cosmological simulations that is tailored for the testing and validation of combined galaxy clustering and weak lensing analyses. We show that the Buzzard v2.0 simulations accurately reproduce many important aspects of the DES Y3 data, including photometric redshift and magnitude distributions, and the relevant set of two-point clustering and weak lensing statistics. We then show that our model for the $3\times2$-point data vector is accurate enough to recover the true cosmology in simulated surveys assuming the true redshift distributions for our source and lens samples, demonstrating robustness to uncertainties in the modeling of the non-linear matter power spectrum, non-linear galaxy bias and higher-order lensing corrections. Additionally, we demonstrate for the first time that our photometric redshift calibration methodology, including information from photometry, spectroscopy, clustering cross-correlations, and galaxy-galaxy lensing ratios, is accurate enough to recover the true cosmology in simulated surveys in the presence of realistic photometric redshift uncertainties.

preprint2021arXiv

Dark Energy Survey Year 3 results: Cosmology with peaks using an emulator approach

We constrain the matter density $Ω_{\mathrm{m}}$ and the amplitude of density fluctuations $σ_8$ within the $Λ$CDM cosmological model with shear peak statistics and angular convergence power spectra using mass maps constructed from the first three years of data of the Dark Energy Survey (DES Y3). We use tomographic shear peak statistics, including cross-peaks: peak counts calculated on maps created by taking a harmonic space product of the convergence of two tomographic redshift bins. Our analysis follows a forward-modelling scheme to create a likelihood of these statistics using N-body simulations, using a Gaussian process emulator. We include the following lensing systematics: multiplicative shear bias, photometric redshift uncertainty, and galaxy intrinsic alignment. Stringent scale cuts are applied to avoid biases from unmodelled baryonic physics. We find that the additional non-Gaussian information leads to a tightening of the constraints on the structure growth parameter yielding $S_8~\equiv~σ_8\sqrt{Ω_{\mathrm{m}}/0.3}~=~0.797_{-0.013}^{+0.015}$ (68% confidence limits), with a precision of 1.8%, an improvement of ~38% compared to the angular power spectra only case. The results obtained with the angular power spectra and peak counts are found to be in agreement with each other and no significant difference in $S_8$ is recorded. We find a mild tension of $1.5 \thinspace σ$ between our study and the results from Planck 2018, with our analysis yielding a lower $S_8$. Furthermore, we observe that the combination of angular power spectra and tomographic peak counts breaks the degeneracy between galaxy intrinsic alignment $A_{\mathrm{IA}}$ and $S_8$, improving cosmological constraints. We run a suite of tests concluding that our results are robust and consistent with the results from other studies using DES Y3 data.

preprint2021arXiv

Dark Energy Survey Year 3 Results: Galaxy clustering and systematics treatment for lens galaxy samples

In this work we present the galaxy clustering measurements of the two DES lens galaxy samples: a magnitude-limited sample optimized for the measurement of cosmological parameters, MagLim, and a sample of luminous red galaxies selected with the redMaGiC algorithm. MagLim / redMaGiC sample contains over 10 million / 2.5 million galaxies and is divided into six / five photometric redshift bins spanning the range $z\in[0.20,1.05]$ / $z\in[0.15,0.90]$. Both samples cover 4143 deg$^2$ over which we perform our analysis blind, measuring the angular correlation function with a S/N $\sim 63$ for both samples. In a companion paper (DES Collaboration et al. 2021)), these measurements of galaxy clustering are combined with the correlation functions of cosmic shear and galaxy-galaxy lensing of each sample to place cosmological constraints with a 3$\times$2pt analysis. We conduct a thorough study of the mitigation of systematic effects caused by the spatially varying survey properties and we correct the measurements to remove artificial clustering signals. We employ several decontamination methods with different configurations to ensure the robustness of our corrections and to determine the systematic uncertainty that needs to be considered for the final cosmology analyses. We validate our fiducial methodology using log-normal mocks, showing that our decontamination procedure induces biases no greater than $0.5σ$ in the $(Ω_m, b)$ plane, where $b$ is galaxy bias. We demonstrate that failure to remove the artificial clustering would introduce strong biases up to $\sim 7 σ$ in $Ω_m$ and of more than $4 σ$ in galaxy bias.

preprint2021arXiv

Milky Way Satellite Census. III. Constraints on Dark Matter Properties from Observations of Milky Way Satellite Galaxies

We perform a comprehensive study of Milky Way (MW) satellite galaxies to constrain the fundamental properties of dark matter (DM). This analysis fully incorporates inhomogeneities in the spatial distribution and detectability of MW satellites and marginalizes over uncertainties in the mapping between galaxies and DM halos, the properties of the MW system, and the disruption of subhalos by the MW disk. Our results are consistent with the cold, collisionless DM paradigm and yield the strongest cosmological constraints to date on particle models of warm, interacting, and fuzzy dark matter. At $95\%$ confidence, we report limits on (i) the mass of thermal relic warm DM, $m_{\rm WDM} > 6.5\ \mathrm{keV}$ (free-streaming length, $λ_{\rm{fs}} \lesssim 10\,h^{-1}\ \mathrm{kpc}$), (ii) the velocity-independent DM-proton scattering cross section, $σ_{0} < 8.8\times 10^{-29}\ \mathrm{cm}^{2}$ for a $100\ \mathrm{MeV}$ DM particle mass (DM-proton coupling, $c_p \lesssim (0.3\ \mathrm{GeV})^{-2}$), and (iii) the mass of fuzzy DM, $m_ϕ> 2.9 \times 10^{-21}\ \mathrm{eV}$ (de Broglie wavelength, $λ_{\rm{dB}} \lesssim 0.5\ \mathrm{kpc}$). These constraints are complementary to other observational and laboratory constraints on DM properties.

preprint2021arXiv

The DES view of the Eridanus supervoid and the CMB Cold Spot

The Cold Spot is a puzzling large-scale feature in the Cosmic Microwave Background temperature maps and its origin has been subject to active debate. As an important foreground structure at low redshift, the Eridanus supervoid was recently detected, but it was subsequently determined that, assuming the standard $Λ$CDM model, only about 10-20$\%$ of the observed temperature depression can be accounted for via its Integrated Sachs-Wolfe imprint. However, $R\gtrsim100~h^{-1}\mathrm{Mpc}$ supervoids elsewhere in the sky have shown ISW imprints $A_{\mathrm{ISW}}\approx5.2\pm1.6$ times stronger than expected from $Λ$CDM ($A_{\mathrm{ISW}}=1$), which warrants further inspection. Using the Year-3 redMaGiC catalogue of luminous red galaxies from the Dark Energy Survey, here we confirm the detection of the Eridanus supervoid as a significant under-density in the Cold Spot's direction at $z<0.2$. We also show, with $\mathrm{S/N}\gtrsim5$ significance, that the Eridanus supervoid appears as the most prominent large-scale under-density in the dark matter mass maps that we reconstructed from DES Year-3 gravitational lensing data. While we report no significant anomalies, an interesting aspect is that the amplitude of the lensing signal from the Eridanus supervoid at the Cold Spot centre is about $30\%$ lower than expected from similar peaks found in N-body simulations based on the standard $Λ$CDM model with parameters $Ω_{\rm m} = 0.279$ and $σ_8 = 0.82$. Overall, our results confirm the causal relation between these individually rare structures in the cosmic web and in the CMB, motivating more detailed future surveys in the Cold Spot region.

preprint2020arXiv

Birds of a Feather? Magellan/IMACS Spectroscopy of the Ultra-Faint Satellites Grus II, Tucana IV, and Tucana V

We present Magellan/IMACS spectroscopy of three recently discovered ultra-faint Milky Way satellites, Grus II, Tucana IV, and Tucana V. We measure systemic velocities of V_hel = -110.0 +/- 0.5 km/s, V_hel = 15.9 +/- 1.8 km/s, and V_hel = -36.2 +/-2.5 km/s for the three objects, respectively. Their large relative velocity differences demonstrate that the satellites are unrelated despite their close physical proximity. We determine a velocity dispersion for Tuc IV of sigma = 4.3^+1.7_-1.0 km/s, but we cannot resolve the velocity dispersions of the other two systems. For Gru II we place an upper limit (90% confidence) on the dispersion of sigma < 1.9 km/s, and for Tuc V we do not obtain any useful limits. All three satellites have metallicities below [Fe/H] = -2.1, but none has a detectable metallicity spread. We determine proper motions for each satellite based on Gaia astrometry and compute their orbits around the Milky Way. Gru II is on a tightly bound orbit with a pericenter of 25 kpc and orbital eccentricity of 0.45. Tuc V likely has an apocenter beyond 100 kpc, and could be approaching the Milky Way for the first time. The current orbit of Tuc IV is similar to that of Gru II, with a pericenter of 25 kpc and an eccentricity of 0.36. However, a backward integration of the position of Tuc IV demonstrates that it collided with the Large Magellanic Cloud at an impact parameter of 4 kpc ~120 Myr ago, deflecting its trajectory and possibly altering its internal kinematics. Based on their sizes, masses, and metallicities, we classify Gru II and Tuc IV as likely dwarf galaxies, but the nature of Tuc V remains uncertain.

preprint2020arXiv

Blinding multi-probe cosmological experiments

The goal of blinding is to hide an experiment's critical results -- here the inferred cosmological parameters -- until all decisions affecting its analysis have been finalised. This is especially important in the current era of precision cosmology, when the results of any new experiment are closely scrutinised for consistency or tension with previous results. In analyses that combine multiple observational probes, like the combination of galaxy clustering and weak lensing in the Dark Energy Survey (DES), it is challenging to blind the results while retaining the ability to check for (in)consistency between different parts of the data. We propose a simple new blinding transformation that works by modifying the summary statistics that are input to parameter estimation, such as two-point correlation functions. The transformation shifts the measured statistics to new values that are consistent with (blindly) shifted cosmological parameters, while preserving internal (in)consistency. We apply the blinding transformation to simulated data for the projected DES Year 3 galaxy clustering and weak lensing analysis, demonstrating that practical blinding is achieved without significant perturbation of internal-consistency checks, as measured here by degradation of the $χ^2$ between data and best-fitting model. Our blinding method conserves $χ^2$ more precisely as experiments evolve to higher precision.

preprint2020arXiv

Detection of cross-correlation between gravitational lensing and gamma rays

In recent years, many gamma-ray sources have been identified, yet the unresolved component hosts valuable information on the faintest emission. In order to extract it, a cross-correlation with gravitational tracers of matter in the Universe has been shown to be a promising tool. We report here the first identification of a cross-correlation signal between gamma rays and the distribution of mass in the Universe probed by weak gravitational lensing. We use the Dark Energy Survey Y1 weak lensing catalogue and the Fermi Large Area Telescope 9-year gamma-ray data, obtaining a signal-to-noise ratio of 5.3. The signal is mostly localised at small angular scales and high gamma-ray energies, with a hint of correlation at extended separation. Blazar emission is likely the origin of the small-scale effect. We investigate implications of the large-scale component in terms of astrophysical sources and particle dark matter emission.

preprint2020arXiv

Milky Way Satellite Census. I. The Observational Selection Function for Milky Way Satellites in DES Y3 and Pan-STARRS DR1

We report the results of a systematic search for ultra-faint Milky Way satellite galaxies using data from the Dark Energy Survey (DES) and Pan-STARRS1 (PS1). Together, DES and PS1 provide multi-band photometry in optical/near-infrared wavelengths over ~80% of the sky. Our search for satellite galaxies targets ~25,000 deg$^2$ of the high-Galactic-latitude sky reaching a 10$σ$ point-source depth of $\gtrsim$ 22.5 mag in the $g$ and $r$ bands. While satellite galaxy searches have been performed independently on DES and PS1 before, this is the first time that a self-consistent search is performed across both data sets. We do not detect any new high-significance satellite galaxy candidates, while recovering the majority of satellites previously detected in surveys of comparable depth. We characterize the sensitivity of our search using a large set of simulated satellites injected into the survey data. We use these simulations to derive both analytic and machine-learning models that accurately predict the detectability of Milky Way satellites as a function of their distance, size, luminosity, and location on the sky. To demonstrate the utility of this observational selection function, we calculate the luminosity function of Milky Way satellite galaxies, assuming that the known population of satellite galaxies is representative of the underlying distribution. We provide access to our observational selection function to facilitate comparisons with cosmological models of galaxy formation and evolution.

preprint2020arXiv

Milky Way Satellite Census. II. Galaxy--Halo Connection Constraints Including the Impact of the Large Magellanic Cloud

The population of Milky Way (MW) satellites contains the faintest known galaxies and thus provides essential insight into galaxy formation and dark matter microphysics. Here we combine a model of the galaxy--halo connection with newly derived observational selection functions based on searches for satellites in photometric surveys over nearly the entire high Galactic latitude sky. In particular, we use cosmological zoom-in simulations of MW-like halos that include realistic Large Magellanic Cloud (LMC) analogs to fit the position-dependent MW satellite luminosity function. We report decisive evidence for the statistical impact of the LMC on the MW satellite population due to an estimated $6\pm 2$ observed LMC-associated satellites, consistent with the number of LMC satellites inferred from Gaia proper-motion measurements, confirming the predictions of cold dark matter models for the existence of satellites within satellite halos. Moreover, we infer that the LMC fell into the MW within the last $2\ \rm{Gyr}$ at high confidence. Based on our detailed full-sky modeling, we find that the faintest observed satellites inhabit halos with peak virial masses below $3.2\times 10^{8}\ M_{\rm{\odot}}$ at $95\%$ confidence, and we place the first robust constraints on the fraction of halos that host galaxies in this regime. We predict that the faintest detectable satellites occupy halos with peak virial masses above $10^{6}\ M_{\rm{\odot}}$, highlighting the potential for powerful galaxy formation and dark matter constraints from future dwarf galaxy searches.

preprint2020arXiv

Modelling the Milky Way. I -- Method and first results fitting the thick disk and halo with DES-Y3 data

We present MWFitting, a method to fit the stellar components of the Galaxy by comparing Hess Diagrams (HDs) from TRILEGAL models to real data. We apply MWFitting to photometric data from the first three years of the Dark Energy Survey (DES). After removing regions containing known resolved stellar systems such as globular clusters, dwarf galaxies, nearby galaxies, the Large Magellanic Cloud and the Sagittarius Stream, our main sample spans a total area of $\sim$2,300 deg$^2$ distributed across the DES footprint. We further explore a smaller subset ($\sim$ 1,300 deg$^2$) that excludes all regions with known stellar streams and stellar overdensities. Validation tests on synthetic data possessing similar properties to the DES data show that the method is able to recover input parameters with a precision better than 3\%. Based on the best-fit models, we create simulated stellar catalogues covering the whole DES footprint down to $g = 24$ magnitude. Comparisons of data and simulations provide evidence for a break in the power law index describing the stellar density of the Milky Way (MW) halo. Several previously discovered stellar over-densities are recovered in the residual stellar density map, showing the reliability of MWFitting in determining the Galactic components. Simulations made with the best-fitting parameters are a promising way to predict MW star counts for surveys such as LSST and Euclid.

preprint2020arXiv

Noise from Undetected Sources in Dark Energy Survey Images

For ground-based optical imaging with current CCD technology, the Poisson fluctuations in source and sky background photon arrivals dominate the noise budget and are readily estimated. Another component of noise, however, is the signal from the undetected population of stars and galaxies. Using injection of artificial galaxies into images, we demonstrate that the measured variance of galaxy moments (used for weak gravitational lensing measurements) in Dark Energy Survey (DES) images is significantly in excess of the Poisson predictions, by up to 30\%, and that the background sky levels are overestimated by current software. By cross-correlating distinct images of "empty" sky regions, we establish that there is a significant image noise contribution from undetected static sources (US), which on average are mildly resolved at DES resolution. Treating these US as a stationary noise source, we compute a correction to the moment covariance matrix expected from Poisson noise. The corrected covariance matrix matches the moment variances measured on the injected DES images to within 5\%. Thus we have an empirical method to statistically account for US in weak lensing measurements, rather than requiring extremely deep sky simulations. We also find that local sky determinations can remove the bias in flux measurements, at a small penalty in additional, but quantifiable, noise.

preprint2020arXiv

STRIDES: a 3.9 per cent measurement of the Hubble constant from the strong lens system DES J0408-5354

We present a blind time-delay cosmographic analysis for the lens system DES J0408$-$5354. This system is extraordinary for the presence of two sets of multiple images at different redshifts, which provide the opportunity to obtain more information at the cost of increased modelling complexity with respect to previously analyzed systems. We perform detailed modelling of the mass distribution for this lens system using three band Hubble Space Telescope imaging. We combine the measured time delays, line-of-sight central velocity dispersion of the deflector, and statistically constrained external convergence with our lens models to estimate two cosmological distances. We measure the "effective" time-delay distance corresponding to the redshifts of the deflector and the lensed quasar $D_{Δt}^{\rm eff}=3382^{+146}_{-115}$ Mpc and the angular diameter distance to the deflector $D_{\rm d}=1711^{+376}_{-280}$ Mpc, with covariance between the two distances. From these constraints on the cosmological distances, we infer the Hubble constant $H_0 = 74.2^{+2.7}_{-3.0}$ km s$^{-1}$ Mpc$^{-1}$ assuming a flat $Λ$CDM cosmology and a uniform prior for $Ω_{\rm m}$ as $Ω_{\rm m} \sim \mathcal{U}(0.05, 0.5)$. This measurement gives the most precise constraint on $H_0$ to date from a single lens. Our measurement is consistent with that obtained from the previous sample of six lenses analyzed by the $H_0$ Lenses in COSMOGRAIL's Wellspring (H0LiCOW) collaboration. It is also consistent with measurements of $H_0$ based on the local distance ladder, reinforcing the tension with the inference from early Universe probes, for example, with 2.2$σ$ discrepancy from the cosmic microwave background measurement.

preprint2020arXiv

Supernova Siblings: Assessing the Consistency of Properties of Type Ia Supernovae that Share the Same Parent Galaxies

While many studies have shown a correlation between properties of the light curves of Type Ia SN (SNe Ia) and properties of their host galaxies, it remains unclear what is driving these correlations. We introduce a new direct method to study these correlations by analyzing `parent' galaxies that host multiple SNe Ia 'siblings'. Here, we search the Dark Energy Survey SN sample, one of the largest samples of discovered SNe, and find 8 galaxies that hosted two likely Type Ia SNe. Comparing the light-curve properties of these SNe and recovered distances from the light curves, we find no better agreement between properties of SNe in the same galaxy as any random pair of galaxies, with the exception of the SN light-curve stretch. We show at $2.8σ$ significance that at least 1/2 of the intrinsic scatter of SNe Ia distance modulus residuals is not from common host properties. We also discuss the robustness with which we could make this evaluation with LSST, which will find $100\times$ more pairs of galaxies, and pave a new line of study on the consistency of Type Ia supernovae in the same parent galaxies. Finally, we argue that it is unlikely some of these SNe are actually single, lensed SN with multiple images.

preprint2019arXiv

Observation and Confirmation of Nine Strong Lensing Systems in Dark Energy Survey Year 1 Data

We describe the observation and confirmation of \nbconfirmtext\ new strong gravitational lenses discovered in Year 1 data from the Dark Energy Survey (DES). We created candidate lists based on a) galaxy group and cluster samples and b) photometrically selected galaxy samples. We selected 46 candidates through visual inspection and then used the Gemini Multi-Object Spectrograph (GMOS) at the Gemini South telescope to acquire spectroscopic follow-up of 21 of these candidates. Through analysis of this spectroscopic follow-up data, we confirmed nine new lensing systems and rejected 2 candidates, but the analysis was inconclusive on 10 candidates. For each of the confirmed systems, we report measured spectroscopic properties, estimated \einsteinradiussub, and estimated enclosed masses. The sources that we targeted have an i-band surface brightness range of iSB ~ 22 - 24 mag arcsec^2 and a spectroscopic redshift range of zspec ~0.8 - 2.6. The lens galaxies have a photometric redshift range of zlens ~ 0.3 - 0.7. The lensing systems range in image-lens separation 2 - 9 arcsec and in enclosed mass 10^12 - 10^13 Msol.

preprint2019arXiv

The SPTpol Extended Cluster Survey

We describe the observations and resultant galaxy cluster catalog from the 2770 deg$^2$ SPTpol Extended Cluster Survey (SPT-ECS). Clusters are identified via the Sunyaev-Zel'dovich (SZ) effect, and confirmed with a combination of archival and targeted follow-up data, making particular use of data from the Dark Energy Survey (DES). With incomplete followup we have confirmed as clusters 244 of 266 candidates at a detection significance $ξ\ge 5$ and an additional 204 systems at $4<ξ<5$. The confirmed sample has a median mass of $M_{500c} \sim {4.4 \times 10^{14} M_\odot h_{70}^{-1}}$, a median redshift of $z=0.49$, and we have identified 44 strong gravitational lenses in the sample thus far. Radio data are used to characterize contamination to the SZ signal; the median contamination for confirmed clusters is predicted to be $\sim$1% of the SZ signal at the $ξ>4$ threshold, and $<4\%$ of clusters have a predicted contamination $>10\% $ of their measured SZ flux. We associate SZ-selected clusters, from both SPT-ECS and the SPT-SZ survey, with clusters from the DES redMaPPer sample, and find an offset distribution between the SZ center and central galaxy in general agreement with previous work, though with a larger fraction of clusters with significant offsets. Adopting a fixed Planck-like cosmology, we measure the optical richness-to-SZ-mass ($λ-M$) relation and find it to be 28% shallower than that from a weak-lensing analysis of the DES data---a difference significant at the 4 $σ$ level---with the relations intersecting at $λ=60$ . The SPT-ECS cluster sample will be particularly useful for studying the evolution of massive clusters and, in combination with DES lensing observations and the SPT-SZ cluster sample, will be an important component of future cosmological analyses.

preprint2016arXiv

A DECam Search for an Optical Counterpart to the LIGO Gravitational Wave Event GW151226

We report the results of a Dark Energy Camera (DECam) optical follow-up of the gravitational wave (GW) event GW151226, discovered by the Advanced LIGO detectors. Our observations cover 28.8 deg$^2$ of the localization region in the $i$ and $z$ bands (containing 3% of the BAYESTAR localization probability), starting 10 hours after the event was announced and spanning four epochs at $2-24$ days after the GW detection. We achieve $5σ$ point-source limiting magnitudes of $i\approx21.7$ and $z\approx21.5$, with a scatter of $0.4$ mag, in our difference images. Given the two day delay, we search this area for a rapidly declining optical counterpart with $\gtrsim 3σ$ significance steady decline between the first and final observations. We recover four sources that pass our selection criteria, of which three are cataloged AGN. The fourth source is offset by $5.8$ arcsec from the center of a galaxy at a distance of 187 Mpc, exhibits a rapid decline by $0.5$ mag over $4$ days, and has a red color of $i-z\approx 0.3$ mag. These properties roughly match the expectations for a kilonova. However, this source was detected several times, starting $94$ days prior to GW151226, in the Pan-STARRS Survey for Transients (dubbed as PS15cdi) and is therefore unrelated to the GW event. Given its long-term behavior, PS15cdi is likely a Type IIP supernova that transitioned out of its plateau phase during our observations, mimicking a kilonova-like behavior. We comment on the implications of this detection for contamination in future optical follow-up observations.

preprint2016arXiv

A Hard X-ray Study of the Normal Star-Forming Galaxy M83 with NuSTAR

We present results from sensitive, multi-epoch NuSTAR observations of the late-type star-forming galaxy M83 (d=4.6 Mpc), which is the first investigation to spatially resolve the hard (E>10 keV) X-ray emission of this galaxy. The nuclear region and ~ 20 off-nuclear point sources, including a previously discovered ultraluminous X-ray (ULX) source, are detected in our NuSTAR observations. The X-ray hardnesses and luminosities of the majority of the point sources are consistent with hard X-ray sources resolved in the starburst galaxy NGC 253. We infer that the hard X-ray emission is most likely dominated by intermediate accretion state black hole binaries and neutron star low-mass X-ray binaries (Z-sources). We construct the X-ray binary luminosity function (XLF) in the NuSTAR band for an extragalactic environment for the first time. The M83 XLF has a steeper XLF than the X-ray binary XLF in NGC 253, consistent with previous measurements by Chandra at softer X-ray energies. The NuSTAR integrated galaxy spectrum of M83 drops quickly above 10 keV, which is also seen in the starburst galaxies NGC253, NGC 3310 and NGC 3256. The NuSTAR observations constrain any AGN to be either highly obscured or to have an extremely low luminosity of $_{\sim}^<$10$^{38}$ erg/s (10-30 keV), implying it is emitting at a very low Eddington ratio. An X-ray point source consistent with the location of the nuclear star cluster with an X-ray luminosity of a few times 10$^{38}$ erg/s may be a low-luminosity AGN but is more consistent with being an X-ray binary.

preprint2016arXiv

A Technique for Detection of PeV Neutrinos Using a Phased Radio Array

The detection of high energy neutrinos ($10^{15}-10^{20}$ eV or $1-10^{5}$ PeV) is an important step toward understanding the most energetic cosmic accelerators and would enable tests of fundamental physics at energy scales that cannot easily be achieved on Earth. In this energy range, there are two expected populations of neutrinos: the astrophysical flux observed with IceCube at lower energies ($\sim1$ PeV) and the predicted cosmogenic flux at higher energies ($\sim10^{18}$ eV). Radio detector arrays such as RICE, ANITA, ARA, and ARIANNA exploit the Askaryan effect and the radio transparency of glacial ice, which together enable enormous volumes of ice to be monitored with sparse instrumentation. We describe here the design for a phased radio array that would lower the energy threshold of radio techniques to the PeV scale, allowing measurement of the astrophysical flux observed with IceCube over an extended energy range. Meaningful energy overlap with optical Cherenkov telescopes could be used for energy calibration. The phased radio array design would also provide more efficient coverage of the large effective volume required to discover cosmogenic neutrinos.

preprint2016arXiv

Accelerator measurements of magnetically-induced radio emission from particle cascades with applications to cosmic-ray air showers

For fifty years, cosmic-ray air showers have been detected by their radio emission. We present the first laboratory measurements that validate electrodynamics simulations used in air shower modeling. An experiment at SLAC provides a beam test of radio-frequency (RF) radiation from charged particle cascades in the presence of a magnetic field, a model system of a cosmic-ray air shower. This experiment provides a suite of controlled laboratory measurements to compare to particle-level simulations of RF emission, which are relied upon in ultra-high-energy cosmic-ray air shower detection. We compare simulations to data for intensity, linearity with magnetic field, angular distribution, polarization, and spectral content. In particular, we confirm modern predictions that the magnetically induced emission in a dielectric forms a cone that peaks at the Cherenkov angle and show that the simulations reproduce the data within systematic uncertainties.

preprint2016arXiv

An Ultra-Faint Galaxy Candidate Discovered in Early Data from the Magellanic Satellites Survey

We report a new ultra-faint stellar system found in Dark Energy Camera data from the first observing run of the Magellanic Satellites Survey (MagLiteS). MagLiteS J0644-5953 (Pictor II or Pic II) is a low surface brightness (μ = 28.5 mag arcsec$^{-2}$ within its half-light radius) resolved overdensity of old and metal-poor stars located at a heliocentric distance of 45 kpc. The physical size (r$_{1/2}$ = 46 pc) and low luminosity (Mv = -3.2 mag) of this satellite are consistent with the locus of spectroscopically confirmed ultra-faint galaxies. MagLiteS J0644-5953 (Pic II) is located 11.3 kpc from the Large Magellanic Cloud (LMC), and comparisons with simulation results in the literature suggest that this satellite was likely accreted with the LMC. The close proximity of MagLiteS J0644-5953 (Pic II) to the LMC also makes it the most likely ultra-faint galaxy candidate to still be gravitationally bound to the LMC.

preprint2016arXiv

Discovery of a Stellar Overdensity in Eridanus-Phoenix in the Dark Energy Survey

We report the discovery of an excess of main sequence turn-off stars in the direction of the constellations of Eridanus and Phoenix from the first year data of the Dark Energy Survey (DES). The Eridanus-Phoenix (EriPhe) overdensity is centered around l~285 deg and b~-60 deg and spans at least 30 deg in longitude and 10 deg in latitude. The Poisson significance of the detection is at least 9 sigma. The stellar population in the overdense region is similar in brightness and color to that of the nearby globular cluster NGC 1261, indicating that the heliocentric distance of EriPhe is about d~16 kpc. The extent of EriPhe in projection is therefore at least ~4 kpc by ~3 kpc. On the sky, this overdensity is located between NGC 1261 and a new stellar stream discovered by DES at a similar heliocentric distance, the so-called Phoenix Stream. Given their similar distance and proximity to each other, it is possible that these three structures may be kinematically associated. Alternatively, the EriPhe overdensity is morphologically similar to the Virgo overdensity and the Hercules-Aquila cloud, which also lie at a similar Galactocentric distance. These three overdensities lie along a polar plane separated by ~120 deg and may share a common origin. Spectroscopic follow-up observations of the stars in EriPhe are required to fully understand the nature of this overdensity.

preprint2016arXiv

Nearest Neighbor: The Low-Mass Milky Way Satellite Tucana III

We present Magellan/IMACS spectroscopy of the recently discovered Milky Way satellite Tucana III (Tuc III). We identify 26 member stars in Tuc III, from which we measure a mean radial velocity of v_hel = -102.3 +/- 0.4 (stat.) +/- 2.0 (sys.) km/s, a velocity dispersion of 0.1^+0.7_-0.1 km/s, and a mean metallicity of [Fe/H] = -2.42^+0.07_-0.08. The upper limit on the velocity dispersion is sigma < 1.5 km/s at 95.5% confidence, and the corresponding upper limit on the mass within the half-light radius of Tuc III is 9.0 x 10^4 Msun. We cannot rule out mass-to-light ratios as large as 240 Msun/Lsun for Tuc III, but much lower mass-to-light ratios that would leave the system baryon-dominated are also allowed. We measure an upper limit on the metallicity spread of the stars in Tuc III of 0.19 dex at 95.5% confidence. Tuc III has a smaller metallicity dispersion and likely a smaller velocity dispersion than any known dwarf galaxy, but a larger size and lower surface brightness than any known globular cluster. Its metallicity is also much lower than those of the clusters with similar luminosity. We therefore tentatively suggest that Tuc III is the tidally-stripped remnant of a dark matter-dominated dwarf galaxy, but additional precise velocity and metallicity measurements will be necessary for a definitive classification. If Tuc III is indeed a dwarf galaxy, it is one of the closest external galaxies to the Sun. Because of its proximity, the most luminous stars in Tuc III are quite bright, including one star at V=15.7 that is the brightest known member star of an ultra-faint satellite.

preprint2016arXiv

The Dark Energy Survey: more than dark energy - an overview

This overview article describes the legacy prospect and discovery potential of the Dark Energy Survey (DES) beyond cosmological studies, illustrating it with examples from the DES early data. DES is using a wide-field camera (DECam) on the 4m Blanco Telescope in Chile to image 5000 sq deg of the sky in five filters (grizY). By its completion the survey is expected to have generated a catalogue of 300 million galaxies with photometric redshifts and 100 million stars. In addition, a time-domain survey search over 27 sq deg is expected to yield a sample of thousands of Type Ia supernovae and other transients. The main goals of DES are to characterise dark energy and dark matter, and to test alternative models of gravity; these goals will be pursued by studying large scale structure, cluster counts, weak gravitational lensing and Type Ia supernovae. However, DES also provides a rich data set which allows us to study many other aspects of astrophysics. In this paper we focus on additional science with DES, emphasizing areas where the survey makes a difference with respect to other current surveys. The paper illustrates, using early data (from `Science Verification', and from the first, second and third seasons of observations), what DES can tell us about the solar system, the Milky Way, galaxy evolution, quasars, and other topics. In addition, we show that if the cosmological model is assumed to be Lambda+ Cold Dark Matter (LCDM) then important astrophysics can be deduced from the primary DES probes. Highlights from DES early data include the discovery of 34 Trans Neptunian Objects, 17 dwarf satellites of the Milky Way, one published z > 6 quasar (and more confirmed) and two published superluminous supernovae (and more confirmed).

preprint2015arXiv

Digging deeper into the Southern skies: a compact Milky-Way companion discovered in first-year Dark Energy Survey data

The Dark Energy Survey (DES) is a 5000 sq. degree survey in the southern hemisphere, which is rapidly reducing the existing north-south asymmetry in the census of MW satellites and other stellar substructure. We use the first-year DES data down to previously unprobed photometric depths to search for stellar systems in the Galactic halo, therefore complementing the previous analysis of the same data carried out by our group earlier this year. Our search is based on a matched filter algorithm that produces stellar density maps consistent with stellar population models of various ages, metallicities, and distances over the survey area. The most conspicuous density peaks in these maps have been identified automatically and ranked according to their significance and recurrence for different input models. We report the discovery of one additional stellar system besides those previously found by several authors using the same first-year DES data. The object is compact, and consistent with being dominated by an old and metal-poor population. DES J0034-4902 is found at high significance and appears in the DES images as a compact concentration of faint blue point sources at ~ 87 {kpc}. Its half-light radius of r_h = 9.88 +/- 4.31 {pc} and total luminosity of M_V ~ -3.05_{-0.42}^{+0.69} are consistent with it being a low mass halo cluster. It is also found to have a very elongated shape. In addition, our deeper probe of DES 1st year data confirms the recently reported satellite galaxy candidate Horologium II as a significant stellar overdensity. We also infer its structural properties and compare them to those reported in the literature.

preprint2015arXiv

Eight New Milky Way Companions Discovered in First-Year Dark Energy Survey Data

We report the discovery of eight new Milky Way companions in ~1,800 deg^2 of optical imaging data collected during the first year of the Dark Energy Survey (DES). Each system is identified as a statistically significant over-density of individual stars consistent with the expected isochrone and luminosity function of an old and metal-poor stellar population. The objects span a wide range of absolute magnitudes (M_V from -2.2 mag to -7.4 mag), physical sizes (10 pc to 170 pc), and heliocentric distances (30 kpc to 330 kpc). Based on the low surface brightnesses, large physical sizes, and/or large Galactocentric distances of these objects, several are likely to be new ultra-faint satellite galaxies of the Milky Way and/or Magellanic Clouds. We introduce a likelihood-based algorithm to search for and characterize stellar over-densities, as well as identify stars with high satellite membership probabilities. We also present completeness estimates for detecting ultra-faint galaxies of varying luminosities, sizes, and heliocentric distances in the first-year DES data.

preprint2015arXiv

Eight Ultra-faint Galaxy Candidates Discovered in Year Two of the Dark Energy Survey

We report the discovery of eight new ultra-faint dwarf galaxy candidates in the second year of optical imaging data from the Dark Energy Survey (DES). Six of these candidates are detected at high confidence, while two lower-confidence candidates are identified in regions of non-uniform survey coverage. The new stellar systems are found by three independent automated search techniques and are identified as overdensities of stars, consistent with the isochrone and luminosity function of an old and metal-poor simple stellar population. The new systems are faint (Mv > -4.7 mag) and span a range of physical sizes (17 pc < $r_{1/2}$ < 181 pc) and heliocentric distances (25 kpc < D < 214 kpc). All of the new systems have central surface brightnesses consistent with known ultra-faint dwarf galaxies (μ< 27.5 mag arcsec$^{-2}$). Roughly half of the DES candidates are more distant, less luminous, and/or have lower surface brightnesses than previously known Milky Way satellite galaxies. Most of the candidates are found in the southern part of the DES footprint close to the Magellanic Clouds. We find that the DES data alone exclude (p < 0.001) a spatially isotropic distribution of Milky Way satellites and that the observed distribution can be well, though not uniquely, described by an association between several of the DES satellites and the Magellanic system. Our model predicts that the full sky may hold ~100 ultra-faint galaxies with physical properties comparable to the DES satellites and that 20-30% of these would be spatially associated with the Magellanic Clouds.

preprint2015arXiv

Gamma-ray flaring activity from the gravitationally lensed blazar PKS 1830-211 observed by Fermi LAT

The Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope routinely detects the highly dust-absorbed, reddened, and MeV-peaked flat spectrum radio quasar PKS 1830-211 (z=2.507). Its apparent isotropic gamma-ray luminosity (E>100 MeV) averaged over $\sim$ 3 years of observations and peaking on 2010 October 14/15 at 2.9 X 10^{50} erg s^{-1}, makes it among the brightest high-redshift Fermi blazars. No published model with a single lens can account for all of the observed characteristics of this complex system. Based on radio observations, one expects time delayed variability to follow about 25 days after a primary flare, with flux about a factor 1.5 less. Two large gamma-ray flares of PKS 1830-211 have been detected by the LAT in the considered period and no substantial evidence for such a delayed activity was found. This allows us to place a lower limit of about 6 on the gamma rays flux ratio between the two lensed images. Swift XRT observations from a dedicated Target of Opportunity program indicate a hard spectrum and with no significant correlation of X-ray flux with the gamma-ray variability. The spectral energy distribution can be modeled with inverse Compton scattering of thermal photons from the dusty torus. The implications of the LAT data in terms of variability, the lack of evident delayed flare events, and different radio and gamma-ray flux ratios are discussed. Microlensing effects, absorption, size and location of the emitting regions, the complex mass distribution of the system, an energy-dependent inner structure of the source, and flux suppression by the lens galaxy for one image path may be considered as hypotheses for understanding our results.

preprint2015arXiv

Mass and galaxy distributions of four massive galaxy clusters from Dark Energy Survey Science Verification data

We measure the weak-lensing masses and galaxy distributions of four massive galaxy clusters observed during the Science Verification phase of the Dark Energy Survey. This pathfinder study is meant to 1) validate the DECam imager for the task of measuring weak-lensing shapes, and 2) utilize DECam's large field of view to map out the clusters and their environments over 90 arcmin. We conduct a series of rigorous tests on astrometry, photometry, image quality, PSF modeling, and shear measurement accuracy to single out flaws in the data and also to identify the optimal data processing steps and parameters. We find Science Verification data from DECam to be suitable for the lensing analysis described in this paper. The PSF is generally well-behaved, but the modeling is rendered difficult by a flux-dependent PSF width and ellipticity. We employ photometric redshifts to distinguish between foreground and background galaxies, and a red-sequence cluster finder to provide cluster richness estimates and cluster-galaxy distributions. By fitting NFW profiles to the clusters in this study, we determine weak-lensing masses that are in agreement with previous work. For Abell 3261, we provide the first estimates of redshift, weak-lensing mass, and richness. In addition, the cluster-galaxy distributions indicate the presence of filamentary structures attached to 1E 0657-56 and RXC J2248.7-4431, stretching out as far as 1 degree (approximately 20 Mpc), showcasing the potential of DECam and DES for detailed studies of degree-scale features on the sky.

preprint2015arXiv

Search for extended gamma-ray emission from the Virgo galaxy cluster with Fermi-LAT

Galaxy clusters are one of the prime sites to search for dark matter (DM) annihilation signals. Depending on the substructure of the DM halo of a galaxy cluster and the cross sections for DM annihilation channels, these signals might be detectable by the latest generation of $γ$-ray telescopes. Here we use three years of Fermi Large Area Telescope (LAT) data, which are the most suitable for searching for very extended emission in the vicinity of nearby Virgo galaxy cluster. Our analysis reveals statistically significant extended emission which can be well characterized by a uniformly emitting disk profile with a radius of 3°that moreover is offset from the cluster center. We demonstrate that the significance of this extended emission strongly depends on the adopted interstellar emission model (IEM) and is most likely an artifact of our incomplete description of the IEM in this region. We also search for and find new point source candidates in the region. We then derive conservative upper limits on the velocity-averaged DM pair annihilation cross section from Virgo. We take into account the potential $γ$-ray flux enhancement due to DM sub-halos and its complex morphology as a merging cluster. For DM annihilating into $b\overline{b}$, assuming a conservative sub-halo model setup, we find limits that are between 1 and 1.5 orders of magnitude above the expectation from the thermal cross section for $m_{\mathrm{DM}}\lesssim100\,\mathrm{GeV}$. In a more optimistic scenario, we exclude $\langle σv \rangle\sim3\times10^{-26}\,\mathrm{cm^{3}\,s^{-1}}$ for $m_{\mathrm{DM}}\lesssim40\,\mathrm{GeV}$ for the same channel. Finally, we derive upper limits on the $γ$-ray-flux produced by hadronic cosmic-ray interactions in the inter cluster medium. We find that the volume-averaged cosmic-ray-to-thermal pressure ratio is less than $\sim6\%$.

preprint2015arXiv

Search for Gamma-Ray Emission from DES Dwarf Spheroidal Galaxy Candidates with Fermi-LAT Data

Due to their proximity, high dark-matter content, and apparent absence of non-thermal processes, Milky Way dwarf spheroidal satellite galaxies (dSphs) are excellent targets for the indirect detection of dark matter. Recently, eight new dSph candidates were discovered using the first year of data from the Dark Energy Survey (DES). We searched for gamma-ray emission coincident with the positions of these new objects in six years of Fermi Large Area Telescope data. We found no significant excesses of gamma-ray emission. Under the assumption that the DES candidates are dSphs with dark matter halo properties similar to the known dSphs, we computed individual and combined limits on the velocity-averaged dark matter annihilation cross section for these new targets. If the estimated dark-matter content of these dSph candidates is confirmed, they will constrain the annihilation cross section to lie below the thermal relic cross section for dark matter particles with masses < 20 GeV annihilating via the b-bbar or tau+tau- channels.

preprint2015arXiv

Stellar Kinematics and Metallicities in the Ultra-Faint Dwarf Galaxy Reticulum II

We present Magellan/M2FS, VLT/GIRAFFE, and Gemini South/GMOS spectroscopy of the newly discovered Milky Way satellite Reticulum II. Based on the spectra of 25 Ret II member stars selected from Dark Energy Survey imaging, we measure a mean heliocentric velocity of 62.8 +/- 0.5 km/s and a velocity dispersion of 3.3 +/- 0.7 km/s. The mass-to-light ratio of Ret II within its half-light radius is 470 +/- 210 Msun/Lsun, demonstrating that it is a strongly dark matter-dominated system. Despite its spatial proximity to the Magellanic Clouds, the radial velocity of Ret II differs from that of the LMC and SMC by 199 and 83 km/s, respectively, suggesting that it is not gravitationally bound to the Magellanic system. The likely member stars of Ret II span 1.3 dex in metallicity, with a dispersion of 0.28 +/- 0.09 dex, and we identify several extremely metal-poor stars with [Fe/H] < -3. In combination with its luminosity, size, and ellipticity, these results confirm that Ret II is an ultra-faint dwarf galaxy. With a mean metallicity of [Fe/H] = -2.65 +/- 0.07, Ret II matches Segue~1 as the most metal-poor galaxy known. Although Ret II is the third-closest dwarf galaxy to the Milky Way, the line-of-sight integral of the dark matter density squared is log J = 18.8 +/- 0.6 Gev^2/cm^5 within 0.2 degrees, indicating that the predicted gamma-ray flux from dark matter annihilation in Ret II is lower than that of several other dwarf galaxies.

preprint2015arXiv

The Origin of the Extragalactic Gamma-Ray Background and Implications for Dark-Matter Annihilation

The origin of the extragalactic $γ$-ray background (EGB) has been debated for some time. { The EGB comprises the $γ$-ray emission from resolved and unresolved extragalactic sources, such as blazars, star-forming galaxies and radio galaxies, as well as radiation from truly diffuse processes.} This letter focuses on the blazar source class, the most numerous detected population, and presents an updated luminosity function and spectral energy distribution model consistent with the blazar observations performed by the {\it Fermi} Large Area Telescope (LAT). We show that blazars account for 50$^{+12}_{-11}$\,\% of the EGB photons ($>$0.1\,GeV), and that {\it Fermi}-LAT has already resolved $\sim$70\,\% of this contribution. Blazars, and in particular low-luminosity hard-spectrum nearby sources like BL Lacs, are responsible for most of the EGB emission above 100\,GeV. We find that the extragalactic background light, which attenuates blazars' high-energy emission, is responsible for the high-energy cut-off observed in the EGB spectrum. Finally, we show that blazars, star-forming galaxies and radio galaxies can naturally account for the amplitude and spectral shape of the background in the 0.1--820\,GeV range, leaving only modest room for other contributions. This allows us to set competitive constraints on the dark-matter annihilation cross section.

preprint2014arXiv

A Focused, Hard X-ray Look at Arp 299 with NuSTAR

We report on simultaneous observations of the local starburst system Arp 299 with NuSTAR and Chandra, which provides the first resolved images of this galaxy up to energies of ~ 45 keV. Fitting the 3-40 keV spectrum reveals a column density of $N_{\rm H}$ ~ 4 x10^{24} cm^{-2}, characteristic of a Compton-thick AGN, and a 10-30 keV luminosity of 1.2x 10^{43} ergs s^{-1}. The hard X-rays detected by NuSTAR above 10 keV are centered on the western nucleus, Arp 299-B, which previous X-ray observations have shown to be the primary source of neutral Fe-K emission. Other X-ray sources, including Arp 299-A, the eastern nucleus which is also thought to harbor an AGN, as well as X-ray binaries, contribute $\lesssim 10%$ to the 10-20 keV emission from the Arp 299 system. The lack of significant emission above 10 keV other than that attributed to Arp 299-B suggests that: a) any AGN in Arp 299-A must be heavily obscured ($N_{\rm H}$ > 10^{24} cm^{-2}) or have a much lower luminosity than Arp 299-B and b) the extranuclear X-ray binaries have spectra that cut-off above ~10 keV. Such soft spectra are characteristic of ultraluminous X-ray (ULX) sources observed to date by NuSTAR.

preprint2014arXiv

Dark Matter Constraints from Observations of 25 Milky Way Satellite Galaxies with the Fermi Large Area Telescope

The dwarf spheroidal satellite galaxies of the Milky Way are some of the most dark-matter-dominated objects known. Due to their proximity, high dark matter content, and lack of astrophysical backgrounds, dwarf spheroidal galaxies are widely considered to be among the most promising targets for the indirect detection of dark matter via gamma rays. Here we report on gamma-ray observations of 25 Milky Way dwarf spheroidal satellite galaxies based on 4 years of Fermi Large Area Telescope (LAT) data. None of the dwarf galaxies are significantly detected in gamma rays, and we present gamma-ray flux upper limits between 500 MeV and 500 GeV. We determine the dark matter content of 18 dwarf spheroidal galaxies from stellar kinematic data and combine LAT observations of 15 dwarf galaxies to constrain the dark matter annihilation cross section. We set some of the tightest constraints to date on the the annihilation of dark matter particles with masses between 2 GeV and 10 TeV into prototypical Standard Model channels. We find these results to be robust against systematic uncertainties in the LAT instrument performance, diffuse gamma-ray background modeling, and assumed dark matter density profile.

preprint2014arXiv

Inferred cosmic-ray spectrum from ${\it Fermi}$-LAT $γ$-ray observations of the Earth's limb

Recent accurate measurements of cosmic-ray (CR) species by ATIC-2, CREAM, and PAMELA reveal an unexpected hardening in the proton and He spectra above a few hundred GeV, a gradual softening of the spectra just below a few hundred GeV, and a harder spectrum of He compared to that of protons. These newly-discovered features may offer a clue to the origin of high-energy CRs. We use the ${\it Fermi}$ Large Area Telescope observations of the $γ$-ray emission from the Earth's limb for an indirect measurement of the local spectrum of CR protons in the energy range $\sim 90~$GeV-$6~$TeV (derived from a photon energy range $15~$GeV-$1~$TeV). Our analysis shows that single power law and broken power law spectra fit the data equally well and yield a proton spectrum with index $2.68 \pm 0.04$ and $2.61 \pm 0.08$ above $\sim 200~$GeV, respectively.

preprint2014arXiv

Search for cosmic-ray induced gamma-ray emission in Galaxy Clusters

Current theories predict relativistic hadronic particle populations in clusters of galaxies in addition to the already observed relativistic leptons. In these scenarios hadronic interactions give rise to neutral pions which decay into $γ$ rays, that are potentially observable with the Large Area Telescope (LAT) on board the Fermi space telescope. We present a joint likelihood analysis searching for spatially extended $γ$-ray emission at the locations of 50 galaxy clusters in 4 years of Fermi-LAT data under the assumption of the universal cosmic-ray model proposed by Pinzke & Pfrommer (2010). We find an excess at a significance of $2.7σ$ which upon closer inspection is however correlated to individual excess emission towards three galaxy clusters: Abell 400, Abell 1367 and Abell 3112. We discuss these cases in detail and conservatively attribute the emission to unmodeled background (for example, radio galaxies within the clusters). Through the combined analysis of 50 clusters we exclude hadronic injection efficiencies in simple hadronic models above 21% and establish limits on the cosmic-ray to thermal pressure ratio within the virial radius, $R_{200}$, to be below 1.2-1.4% depending on the morphological classification. In addition we derive new limits on the $γ$-ray flux from individual clusters in our sample.

preprint2014arXiv

The spectrum of isotropic diffuse gamma-ray emission between 100 MeV and 820 GeV

The γ-ray sky can be decomposed into individually detected sources, diffuse emission attributed to the interactions of Galactic cosmic rays with gas and radiation fields, and a residual all-sky emission component commonly called the isotropic diffuse γ-ray background (IGRB). The IGRB comprises all extragalactic emissions too faint or too diffuse to be resolved in a given survey, as well as any residual Galactic foregrounds that are approximately isotropic. The first IGRB measurement with the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope (Fermi) used 10 months of sky-survey data and considered an energy range between 200 MeV and 100 GeV. Improvements in event selection and characterization of cosmic-ray backgrounds, better understanding of the diffuse Galactic emission, and a longer data accumulation of 50 months, allow for a refinement and extension of the IGRB measurement with the LAT, now covering the energy range from 100 MeV to 820 GeV. The IGRB spectrum shows a significant high-energy cutoff feature, and can be well described over nearly four decades in energy by a power law with exponential cutoff having a spectral index of $2.32\pm0.02$ and a break energy of $(279\pm52)$ GeV using our baseline diffuse Galactic emission model. The total intensity attributed to the IGRB is $(7.2\pm0.6) \times 10^{-6}$ cm$^{-2}$ s$^{-1}$ sr$^{-1}$ above 100 MeV, with an additional $+15$%/$-30$% systematic uncertainty due to the Galactic diffuse foregrounds.

preprint2013arXiv

Constraints on the Galactic Halo Dark Matter from Fermi-LAT Diffuse Measurements

We have performed an analysis of the diffuse gamma-ray emission with the Fermi Large Area Telescope in the Milky Way Halo region searching for a signal from dark matter annihilation or decay. In the absence of a robust dark matter signal, constraints are presented. We consider both gamma rays produced directly in the dark matter annihilation/decay and produced by inverse Compton scattering of the e+e- produced in the annihilation/decay. Conservative limits are derived requiring that the dark matter signal does not exceed the observed diffuse gamma-ray emission. A second set of more stringent limits is derived based on modeling the foreground astrophysical diffuse emission using the GALPROP code. Uncertainties in the height of the diffusive cosmic-ray halo, the distribution of the cosmic-ray sources in the Galaxy, the index of the injection cosmic-ray electron spectrum and the column density of the interstellar gas are taken into account using a profile likelihood formalism, while the parameters governing the cosmic-ray propagation have been derived from fits to local cosmic-ray data. The resulting limits impact the range of particle masses over which dark matter thermal production in the early Universe is possible, and challenge the interpretation of the PAMELA/Fermi-LAT cosmic ray anomalies as annihilation of dark matter.

preprint2013arXiv

Constraints on the Galactic Population of TEV Pulsar Wind Nebulae Using Fermi Large Area Telescope Observations

Pulsar wind nebulae (PWNe) have been established as the most populous class of TeV gamma-ray emitters. Since launch, the Fermi Large Area Telescope (LAT)identified five high-energy (100MeV <E< 100 GeV) gamma-ray sources as PWNe, and detected a large number of PWNe candidates, all powered by young and energetic pulsars. The wealth of multi-wavelength data available and the new results provided by Fermi-LAT give us an opportunity to find new PWNe and to explore the radiative processes taking place in known ones. The TeV gamma-ray unidentifiedsources (UNIDs) are the best candidates for finding new PWNe. Using 45 months of Fermi-LAT data for energies above 10 GeV, an analysis was performed near the position of 58TeV PWNe and UNIDs within 5deg of the Galactic Plane to establish new constraints on PWNe properties and find new clues on the nature of UNIDs. Of the 58 sources, 30 were detected, and this work provides their gamma-rayfluxes for energies above 10 GeV. The spectral energy distributions (SED) andupper limits, in the multi-wavelength context, also provide new information on the source nature and can help distinguish between emission scenarios, i.e. between classification as a pulsar candidate or as a PWN candidate. Six new GeV PWNe candidates are described in detail and compared with existing models. A population study of GeV PWNe candidates as a function of the pulsar/PWN system characteristics is presented.

preprint2013arXiv

Detection of the Characteristic Pion-Decay Signature in Supernova Remnants

Cosmic rays are particles (mostly protons) accelerated to relativistic speeds. Despite wide agreement that supernova remnants (SNRs) are the sources of galactic cosmic rays, unequivocal evidence for the acceleration of protons in these objects is still lacking. When accelerated protons encounter interstellar material, they produce neutral pions, which in turn decay into gamma rays. This offers a compelling way to detect the acceleration sites of protons. The identification of pion-decay gamma rays has been difficult because high-energy electrons also produce gamma rays via bremsstrahlung and inverse Compton scattering. We detected the characteristic pion-decay feature in the gamma-ray spectra of two SNRs, IC 443 and W44, with the Fermi Large Area Telescope. This detection provides direct evidence that cosmic-ray protons are accelerated in SNRs.

preprint2012arXiv

Binary Millisecond Pulsar Discovery via Gamma-Ray Pulsations

Millisecond pulsars, old neutron stars spun-up by accreting matter from a companion star, can reach high rotation rates of hundreds of revolutions per second. Until now, all such "recycled" rotation-powered pulsars have been detected by their spin-modulated radio emission. In a computing-intensive blind search of gamma-ray data from the Fermi Large Area Telescope (with partial constraints from optical data), we detected a 2.5-millisecond pulsar, PSR J1311-3430. This unambiguously explains a formerly unidentified gamma-ray source that had been a decade-long enigma, confirming previous conjectures. The pulsar is in a circular orbit with an orbital period of only 93 minutes, the shortest of any spin-powered pulsar binary ever found.

preprint2012arXiv

Fermi LAT Search for Dark Matter in Gamma-ray Lines and the Inclusive Photon Spectrum

Dark matter particle annihilation or decay can produce monochromatic gamma-ray lines and contribute to the diffuse gamma-ray background. Flux upper limits are presented for gamma-ray spectral lines from 7 to 200 GeV and for the diffuse gamma-ray background from 4.8 GeV to 264 GeV obtained from two years of Fermi Large Area Telescope data integrated over most of the sky. We give cross section upper limits and decay lifetime lower limits for dark matter models that produce gamma-ray lines or contribute to the diffuse spectrum, including models proposed as explanations of the PAMELA and Fermi cosmic-ray data.

preprint2012arXiv

Gamma-ray observations of the Orion Molecular Clouds with the Fermi Large Area Telescope

We report on the gamma-ray observations of giant molecular clouds Orion A and B with the Large Area Telescope (LAT) on-board the Fermi Gamma-ray Space Telescope. The gamma-ray emission in the energy band between \sim100 MeV and \sim100 GeV is predicted to trace the gas mass distribution in the clouds through nuclear interactions between the Galactic cosmic rays (CRs) and interstellar gas. The gamma-ray production cross-section for the nuclear interaction is known to \sim10% precision which makes the LAT a powerful tool to measure the gas mass column density distribution of molecular clouds for a known CR intensity. We present here such distributions for Orion A and B, and correlate them with those of the velocity integrated CO intensity (WCO) at a 1° \times1° pixel level. The correlation is found to be linear over a WCO range of ~10 fold when divided in 3 regions, suggesting penetration of nuclear CRs to most of the cloud volumes. The Wco-to-mass conversion factor, Xco, is found to be \sim2.3\times10^20 cm-2(K km s-1)-1 for the high-longitude part of Orion A (l > 212°), \sim1.7 times higher than \sim1.3 \times 10^20 found for the rest of Orion A and B. We interpret the apparent high Xco in the high-longitude region of Orion A in the light of recent works proposing a non-linear relation between H2 and CO densities in the diffuse molecular gas. Wco decreases faster than the H2 column density in the region making the gas "darker" to Wco.

preprint2012arXiv

GeV Observations of Star-forming Galaxies with \textit{Fermi} LAT

Recent detections of the starburst galaxies M82 and NGC 253 by gamma-ray telescopes suggest that galaxies rapidly forming massive stars are more luminous at gamma-ray energies compared to their quiescent relatives. Building upon those results, we examine a sample of 69 dwarf, spiral, and luminous and ultraluminous infrared galaxies at photon energies 0.1-100 GeV using 3 years of data collected by the Large Area Telescope (LAT) on the \textit{Fermi Gamma-ray Space Telescope} (\textit{Fermi}). Measured fluxes from significantly detected sources and flux upper limits for the remaining galaxies are used to explore the physics of cosmic rays in galaxies. We find further evidence for quasi-linear scaling relations between gamma-ray luminosity and both radio continuum luminosity and total infrared luminosity which apply both to quiescent galaxies of the Local Group and low-redshift starburst galaxies (conservative $P$-values $\lesssim0.05$ accounting for statistical and systematic uncertainties). The normalizations of these scaling relations correspond to luminosity ratios of $\log(L_{0.1-100 \rm{GeV}}/L_{1.4 \rm{GHz}}) = 1.7 \pm 0.1_{\rm (statistical)} \pm 0.2_{\rm (dispersion)}$ and $\log(L_{0.1-100 \rm{GeV}}/L_{8-1000 μ\rm{m}}) = -4.3 \pm 0.1_{\rm (statistical)} \pm 0.2_{\rm (dispersion)}$ for a galaxy with a star formation rate of 1 $M_{\odot}$ yr$^{-1}$, assuming a Chabrier initial mass function. Using the relationship between infrared luminosity and gamma-ray luminosity, the collective intensity of unresolved star-forming galaxies at redshifts $0<z<2.5$ above 0.1 GeV is estimated to be 0.4-2.4 $\times 10^{-6}$ ph cm$^{-2}$ s$^{-1}$ sr$^{-1}$ (4-23% of the intensity of the isotropic diffuse component measured with the LAT). We anticipate that $\sim10$ galaxies could be detected by their cosmic-ray induced gamma-ray emission during a 10-year \textit{Fermi} mission.

preprint2012arXiv

Measurement of separate cosmic-ray electron and positron spectra with the Fermi Large Area Telescope

We measured separate cosmic-ray electron and positron spectra with the Fermi Large Area Telescope. Because the instrument does not have an onboard magnet, we distinguish the two species by exploiting the Earth's shadow, which is offset in opposite directions for opposite charges due to the Earth's magnetic field. We estimate and subtract the cosmic-ray proton background using two different methods that produce consistent results. We report the electron-only spectrum, the positron-only spectrum, and the positron fraction between 20 GeV and 200 GeV. We confirm that the fraction rises with energy in the 20-100 GeV range. The three new spectral points between 100 and 200 GeV are consistent with a fraction that is continuing to rise with energy.

preprint2012arXiv

Search for Dark Matter Satellites using the FERMI-LAT

Numerical simulations based on the Lambda-CDM model of cosmology predict a large number of as yet unobserved Galactic dark matter satellites. We report the results of a Large Area Telescope (LAT) search for these satellites via the gamma-ray emission expected from the annihilation of weakly interacting massive particle (WIMP) dark matter. Some dark matter satellites are expected to have hard gamma-ray spectra, finite angular extents, and a lack of counterparts at other wavelengths. We sought to identify LAT sources with these characteristics, focusing on gamma-ray spectra consistent with WIMP annihilation through the $b \bar b$ channel. We found no viable dark matter satellite candidates using one year of data, and we present a framework for interpreting this result in the context of numerical simulations to constrain the velocity-averaged annihilation cross section for a conventional 100 GeV WIMP annihilating through the $b \bar b$ channel.

preprint2012arXiv

TARGET: A multi-channel digitizer chip for very-high-energy gamma-ray telescopes

The next-generation very-high-energy (VHE) gamma-ray observatory, the Cherenkov Telescope Array, will feature dozens of imaging atmospheric Cherenkov telescopes (IACTs), each with thousands of pixels of photo-sensors. To be affordable and reliable, reading out such a mega-channel array requires event recording technology that is highly integrated and modular, with a low cost per channel. We present the design and performance of a chip targeted to this application: the TeV Array Readout with GSa/s sampling and Event Trigger (TARGET). This application-specific integrated circuit (ASIC) has 16 parallel input channels, a 4096-sample buffer for each channel, adjustable input termination, self-trigger functionality, and tight window-selected readout. We report the performance of TARGET in terms of sampling frequency, power consumption, dynamic range, current-mode gain, analog bandwidth, and cross talk. The large number of channels per chip allows a low cost per channel ($\10 to $\20 including front-end and back-end electronics but not including photosensors) to be achieved with a TARGET-based IACT readout system. In addition to basic performance parameters of the TARGET chip itself, we present a camera module prototype as well as a second-generation chip (TARGET~2), both of which have been produced.

preprint2011arXiv

Constraining Dark Matter Models from a Combined Analysis of Milky Way Satellites with the Fermi Large Area Telescope

Satellite galaxies of the Milky Way are among the most promising targets for dark matter searches in gamma rays. We present a search for dark matter consisting of weakly interacting massive particles, applying a joint likelihood analysis to 10 satellite galaxies with 24 months of data of the Fermi Large Area Telescope. No dark matter signal is detected. Including the uncertainty in the dark matter distribution, robust upper limits are placed on dark matter annihilation cross sections. The 95% confidence level upper limits range from about 1e-26 cm^3 s^-1 at 5 GeV to about 5e-23 cm^3 s^-1 at 1 TeV, depending on the dark matter annihilation final state. For the first time, using gamma rays, we are able to rule out models with the most generic cross section (~3e-26 cm^3 s^-1 for a purely s-wave cross section), without assuming additional boost factors.

preprint2011arXiv

Fermi-LAT Detection of the Young SuperNova Remnant Tycho

After almost three years of data taking in sky survey mode, the \emph{Fermi}-LAT has detected $γ$-ray emission toward the Tycho's Supernova Remnant (SNR). The Tycho SNR is among the youngest remnants in the Galaxy, originating from a Type Ia Supernova in AD 1572. The $γ$-ray integral flux from 400 MeV up to 100 GeV has been measured to be (3.5$\pm1.1_{stat}\pm0.7_{syst}$)$\times10^{-9}$ cm$^{-2}$s$^{-1}$ with a photon index of 2.3$\pm0.2_{stat}\pm0.1_{syst}$.

preprint2011arXiv

Simultaneous multi-wavelength campaign on PKS 2005-489 in a high state

The high-frequency peaked BL Lac object PKS 2005-489 was the target of a multi-wavelength campaign with simultaneous observations in the TeV gamma-ray (H.E.S.S.), GeV gamma-ray (Fermi/LAT), X-ray (RXTE, Swift), UV (Swift) and optical (ATOM, Swift) bands. This campaign was carried out during a high flux state in the synchrotron regime. The flux in the optical and X-ray bands reached the level of the historical maxima. The hard GeV spectrum observed with Fermi/LAT connects well to the very high energy (VHE, E>100GeV) spectrum measured with H.E.S.S. with a peak energy between ~5 and 500 GeV. Compared to observations with contemporaneous coverage in the VHE and X-ray bands in 2004, the X-ray flux was ~50 times higher during the 2009 campaign while the TeV gamma-ray flux shows marginal variation over the years. The spectral energy distribution during this multi-wavelength campaign was fit by a one zone synchrotron self-Compton model with a well determined cutoff in X-rays. The parameters of a one zone SSC model are inconsistent with variability time scales. The variability behaviour over years with the large changes in synchrotron emission and small changes in the inverse Compton emission does not warrant an interpretation within a one-zone SSC model despite an apparently satisfying fit to the broadband data in 2009.

preprint2010arXiv

Constraints on Cosmological Dark Matter Annihilation from the Fermi-LAT Isotropic Diffuse Gamma-Ray Measurement

The first published Fermi large area telescope (Fermi-LAT) measurement of the isotropic diffuse gamma-ray emission is in good agreement with a single power law, and is not showing any signature of a dominant contribution from dark matter sources in the energy range from 20 to 100 GeV. We use the absolute size and spectral shape of this measured flux to derive cross section limits on three types of generic dark matter candidates: annihilating into quarks, charged leptons and monochromatic photons. Predicted gamma-ray fluxes from annihilating dark matter are strongly affected by the underlying distribution of dark matter, and by using different available results of matter structure formation we assess these uncertainties. We also quantify how the dark matter constraints depend on the assumed conventional backgrounds and on the Universe's transparency to high-energy gamma-rays. In reasonable background and dark matter structure scenarios (but not in all scenarios we consider) it is possible to exclude models proposed to explain the excess of electrons and positrons measured by the Fermi-LAT and PAMELA experiments. Derived limits also start to probe cross sections expected from thermally produced relics (e.g. in minimal supersymmetry models) annihilating predominantly into quarks. For the monochromatic gamma-ray signature, the current measurement constrains only dark matter scenarios with very strong signals.

preprint2010arXiv

Discovery of a GeV blazar shining through the Galactic plane

The \emph{Fermi} Large Area Telescope (LAT) discovered a new gamma-ray source near the Galactic plane, \object{Fermi J0109+6134}, when it flared brightly in 2010 February. The low Galactic latitude (b =-1.2\degr) indicated that the source could be located within the Galaxy, which motivated rapid multi-wavelength follow-up including radio, optical, and X-ray observations. We report the results of analyzing all 19 months of LAT data for the source, and of X-ray observations with both \emph{Swift} and the \emph{Chandra X-ray Observatory}. We determined the source redshift, z =0.783, using a Keck LRIS observation. Finally, we compiled a broadband spectral energy distribution (SED) from both historical and new observations contemporaneous with the 2010 February flare. The redshift, SED, optical line width, X-ray absorption, and multi-band variability indicate that this new GeV source is a blazar seen through the Galactic plane. Because several of the optical emission lines have equivalent width >5Å, this blazar belongs in the flat-spectrum radio quasar category.

preprint2010arXiv

Fermi LAT Search for Photon Lines from 30 to 200 GeV and Dark Matter Implications

Dark matter (DM) particle annihilation or decay can produce monochromatic $γ$-rays readily distinguishable from astrophysical sources. $γ$-ray line limits from 30 GeV to 200 GeV obtained from 11 months of Fermi Large Area Space Telescope data from 20-300 GeV are presented using a selection based on requirements for a $γ$-ray line analysis, and integrated over most of the sky. We obtain $γ$-ray line flux upper limits in the range $0.6-4.5\times 10^{-9}\mathrm{cm}^{-2}\mathrm{s}^{-1}$, and give corresponding DM annihilation cross-section and decay lifetime limits. Theoretical implications are briefly discussed.

preprint2010arXiv

Radio and Gamma-Ray Constraints on the Emission Geometry and Birthplace of PSR J2043+2740

We report on the first year of Fermi gamma-ray observations of pulsed high-energy emission from the old PSR J2043+2740. The study of the gamma-ray efficiency of such old pulsars gives us an insight into the evolution of pulsars' ability to emit in gammma rays as they age. The gamma-ray lightcurve of this pulsar above 0.1 GeV is clearly defined by two sharp peaks, 0.353+/-0.035 periods apart. We have combined the gamma-ray profile characteristics of PSR J2043+2740 with the geometrical properties of the pulsar's radio emission, derived from radio polarization data, and constrained the pulsar-beam geometry in the framework of a Two Pole Caustic and an Outer Gap model. The ranges of magnetic inclination and viewing angle were determined to be {alpha,zeta}~{52-57,61-68} for the Two Pole Caustic model, and {alpha,zeta}~{62-73,74-81} and {alpha,zeta}~{72-83,60-75} for the Outer Gap model. Based on this geometry, we assess possible birth locations for this pulsar and derive a likely proper motion, sufficiently high to be measurable with VLBI. At a characteristic age of 1.2 Myr, PSR J2043+2740 is the third oldest of all discovered, non-recycled, gamma-ray pulsars: it is twice as old as the next oldest, PSR J0357+32, and younger only than the recently discovered PSR J1836+5925 and PSR J2055+25, both of which are at least 5 and 10 times less energetic, respectively.

preprint2010arXiv

The Discovery of gamma-Ray Emission From The Blazar RGB J0710+591

The high-frequency-peaked BL Lacertae object RGB J0710+591 was observed in the very high-energy (VHE; E > 100 GeV) wave band by the VERITAS array of atmospheric Cherenkov telescopes. The observations, taken between 2008 December and 2009 March and totaling 22.1 hr, yield the discovery of VHE gamma rays from the source. RGB J0710+591 is detected at a statistical significance of 5.5 standard deviations (5.5σ) above the background, corresponding to an integral flux of (3.9 +/- 0.8) x 10-12 cm-2 s-1 (3% of the Crab Nebula's flux) above 300 GeV. The observed spectrum can be fit by a power law from 0.31 to 4.6 TeV with a photon spectral index of 2.69 +/- 0.26stat +/- 0.20sys. These data are complemented by contemporaneous multiwavelength data from the Fermi Large Area Telescope, the Swift X-ray Telescope, the Swift Ultra-Violet and Optical Telescope, and the Michigan-Dartmouth-MIT observatory. Modeling the broadband spectral energy distribution (SED) with an equilibrium synchrotron self-Compton model yields a good statistical fit to the data. The addition of an external-Compton component to the model does not improve the fit nor brings the system closer to equipartition. The combined Fermi and VERITAS data constrain the properties of the high-energy emission component of the source over 4 orders of magnitude and give measurements of the rising and falling sections of the SED.

preprint2009arXiv

The Spectral Energy Distribution of Fermi bright blazars

(Abridged) We have conducted a detailed investigation of the broad-band spectral properties of the γ-ray selected blazars of the Fermi LAT Bright AGN Sample (LBAS). By combining our accurately estimated Fermi gamma-ray spectra with Swift, radio, infra-red, optical and other hard X-ray/gamma-ray data, collected within three months of the LBAS data taking period, we were able to assemble high-quality and quasi-simultaneous Spectral Energy Distributions (SED) for 48 LBAS blazars.The SED of these gamma-ray sources is similar to that of blazars discovered at other wavelengths, clearly showing, in the usual Log $ν$ - Log $ν$ F$_ν$ representation, the typical broad-band spectral signatures normally attributed to a combination of low-energy synchrotron radiation followed by inverse Compton emission of one or more components. We have used these SEDs to characterize the peak intensity of both the low and the high-energy components. The results have been used to derive empirical relationships that estimate the position of the two peaks from the broad-band colors (i.e. the radio to optical and optical to X-ray spectral slopes) and from the gamma-ray spectral index. Our data show that the synchrotron peak frequency $ν_p^S$ is positioned between 10$^{12.5}$ and 10$^{14.5}$ Hz in broad-lined FSRQs and between $10^{13}$ and $10^{17}$ Hz in featureless BL Lacertae objects.We find that the gamma-ray spectral slope is strongly correlated with the synchrotron peak energy and with the X-ray spectral index, as expected at first order in synchrotron - inverse Compton scenarios. However, simple homogeneous, one-zone, Synchrotron Self Compton (SSC) models cannot explain most of our SEDs, especially in the case of FSRQs and low energy peaked (LBL) BL Lacs. (...)