Unresolved Astrophysical Sources Research Articles

Superflares on solar-like stars

Context. Over the past years, thousands of stellar flares have been detected by harvesting data from large photometric surveys. These detections, however, do not account for potential sources of contamination such as background stars or small Solar System objects appearing in the same aperture as the primary target. Aims. We present a new method for identifying the true flare sources in large photometric surveys using data from the Kepler mission as an illustrative example. The new method considers not only the brightness excess in the stellar light curves, but also the location of this excess in the pixel-level data. Methods. Potential flares are identified in two steps. First, we search the light curves for at least two subsequent data points exceeding a 5σ threshold above the running mean. For these two cadences, we subtract the “quiet” stellar flux from the Kepler pixel data to obtain new images where the potential flare is the main light source. In the second step, we use a Bayesian approach to fit the point spread function of the instrument to determine the most likely location of the flux excess on the detector. We match this location with the position of the primary target and other stars from the Gaia DR2 catalog within a radius of 10 arcsec around the primary Kepler target. When the location of the flux excess and the target star coincide, we associate the event with a flare on the target star. Results. We applied our method to 5862 main-sequence stars with near-solar effective temperatures. From the first step we found 2274 events exceeding the 5σ level in at least two consecutive points in the light curves. Applying the second step reduced this number to 342 superflares. Of these, 283 flares occurred on 178 target stars and 47 events are associated with fainter background stars; in 10 cases the flare location could not be distinguished between the target and a background star. We also present cases where flares were reported previously but our technique could not attribute them to the target star. Conclusions. We conclude that identifying outliers in the light curves alone is insufficient to attribute them to stellar flares and that flares can only be uniquely attributed to a certain star when the instrument pixel-level data together with the point spread function are taken into account. As a consequence, previous flare statistics are likely contaminated by instrumental effects and unresolved astrophysical sources.

Unified mapmaking for an anisotropic stochastic gravitational wave background

A stochastic gravitational wave background (SGWB), created by the superposition of signals from unresolved astrophysical sources, may be detected in the next few years. Estimating the variation of intensity across the sky can, therefore, play a key role in improving our understanding of astrophysical models. Skymaps have been produced for all the data-taking runs of the advanced ground-based interferometric detectors. While these maps are being produced in pixel and SpH bases, to probe, respectively, localized and diffuse astrophysical and cosmological sources, with algorithms that employ cross-correlation as the common strategy, the underlying algebra and numerical implementation remain different. As a consequence, there was a need for producing skymaps in both bases in those analyses. We show that these manifestly redundant methods could indeed be unified to a single analysis that can probe very different scales and demonstrate it by applying them on real data. We first develop the algebra to show that the results in two different bases are easily transformable. We then incorporate both the schemes in the now-standard analysis pipeline for anisotropic SGWB, PyStoch. This will enable SGWB anisotropy searches in SpH basis also to take full advantage of integrated \hpx tools and makes it computationally feasible to perform the search in every frequency bin. We, however, follow a different approach for direct estimation of the SpH moments. We show that the results obtained from these different methods match very well; the differences are less than $0.1$\% for the SpH moments and less than $0.01$\% for the Fisher information matrices. Thus we conclude that a single skymap will be sufficient to describe the anisotropies in a stochastic background. The multiple capabilities of PyStoch will be useful for estimating and constraining various measures that characterize an anisotropic background.

Foreground mismodeling and the point source explanation of the Fermi Galactic Center excess

The Fermi Large Area Telescope has observed an excess of ~GeV energy gamma rays from the center of the Milky Way, which may arise from near-thermal dark matter annihilation. Firmly establishing the dark matter origin for this excess is however complicated by challenges in modeling diffuse cosmic-ray foregrounds as well as unresolved astrophysical sources, such as millisecond pulsars. Non-Poissonian Template Fitting (NPTF) is one statistical technique that has previously been used to show that at least some fraction of the GeV excess is likely due to a population of dim point sources. These results were recently called into question by Leane and Slatyer (2019), who showed that a synthetic dark matter annihilation signal injected on top of the real Fermi data is not recovered by the NPTF procedure. In this work, we perform a dedicated study of the Fermi data and explicitly show that the central result of Leane and Slatyer (2019) is likely driven by the fact that their choice of model for the Galactic foreground emission does not provide a sufficiently good description of the data. We repeat the NPTF analyses using a state-of-the-art model for diffuse gamma-ray emission in the Milky Way and introduce a novel statistical procedure, based on spherical-harmonic marginalization, to provide an improved description of the Galactic diffuse emission in a data-driven fashion. With these improvements, we find that the NPTF results continue to robustly favor the interpretation that the Galactic Center excess is due, in part, to unresolved astrophysical point sources across the analysis variations that we have explored.

Characterizing the nature of the unresolved point sources in the Galactic Center: An assessment of systematic uncertainties

The Galactic Center Excess (GCE) of GeV gamma rays can be explained as a signal of annihilating dark matter or of emission from unresolved astrophysical sources, such as millisecond pulsars. Evidence for the latter is provided by a statistical procedure---referred to as Non-Poissonian Template Fitting (NPTF)---that distinguishes the smooth distribution of photons expected for dark matter annihilation from a "clumpy" photon distribution expected for point sources. In this paper, we perform an extensive study of the NPTF on simulated data, exploring its ability to recover the flux and luminosity function of unresolved sources at the Galactic Center. When astrophysical background emission is perfectly modeled, we find that the NPTF successfully distinguishes between the dark matter and point source hypotheses when either component makes up the entirety of the GCE. When the GCE is a mixture of dark matter and point sources, the NPTF may fail to reconstruct the correct contribution of each component. We further study the impact of mismodeling the Galactic diffuse backgrounds, finding that while a dark matter signal could be attributed to point sources in some outlying cases for the scenarios we consider, the significance of a true point source signal remains robust. Our work enables us to comment on a recent study by Leane and Slatyer (2019) that questions prior NPTF conclusions because the method does not recover an artificial dark matter signal injected on actual Fermi data. We demonstrate that the failure of the NPTF to extract an artificial dark matter signal can be natural when point sources are present in the data---with the effect further exacerbated by the presence of diffuse mismodeling---and does not on its own invalidate the conclusions of the NPTF analysis in the Inner Galaxy.

First Predictions of the Angular Power Spectrum of the Astrophysical Gravitational Wave Background.

We present the first predictions for the angular power spectrum of the astrophysical gravitational wave background constituted of the radiation emitted by all resolved and unresolved astrophysical sources. Its shape and amplitude depend on both the astrophysical properties on galactic scales and on cosmological properties. We show that the angular power spectrum behaves as C_{ℓ}∝1/ℓ on large scales and that relative fluctuations of the signal are of order 30% at 100Hz. We also present the correlations of the astrophysical gravitational wave background with weak lensing and galaxy distribution. These numerical results pave the way to the study of a new observable at the crossroad between general relativity, astrophysics, and cosmology.

Search for Tensor, Vector, and Scalar Polarizations in the Stochastic Gravitational-Wave Background.

The detection of gravitational waves with Advanced LIGO and Advanced Virgo has enabled novel tests of general relativity, including direct study of the polarization of gravitational waves. While general relativity allows for only two tensor gravitational-wave polarizations, general metric theories can additionally predict two vector and two scalar polarizations. The polarization of gravitational waves is encoded in the spectral shape of the stochastic gravitational-wave background, formed by the superposition of cosmological and individually unresolved astrophysical sources. Using data recorded by Advanced LIGO during its first observing run, we search for a stochastic background of generically polarized gravitational waves. We find no evidence for a background of any polarization, and place the first direct bounds on the contributions of vector and scalar polarizations to the stochastic background. Under log-uniform priors for the energy in each polarization, we limit the energy densities of tensor, vector, and scalar modes at 95% credibility to Ω_{0}^{T}<5.58×10^{-8}, Ω_{0}^{V}<6.35×10^{-8}, and Ω_{0}^{S}<1.08×10^{-7} at a reference frequency f_{0}=25 Hz.

Component separation of a isotropic Gravitational Wave Background

A Gravitational Wave Background (GWB) is expected in the universe from the superposition of a large number of unresolved astrophysical sources and phenomena in the early universe. Each component of the background (e.g., from primordial metric perturbations, binary neutron stars, milli-second pulsars etc.) has its own spectral shape. Many ongoing experiments aim to probe GWB at a variety of frequency bands. In the last two decades, using data from ground-based laser interferometric gravitational wave (GW) observatories, upper limits on GWB were placed in the frequency range of 0∼ 50−100 Hz, considering one spectral shape at a time. However, one strong component can significantly enhance the estimated strength of another component. Hence, estimation of the amplitudes of the components with different spectral shapes should be done jointly. Here we propose a method for “component separation” of a statistically isotropic background, that can, for the first time, jointly estimate the amplitudes of many components and place upper limits. The method is rather straightforward and needs negligible amount of computation. It utilises the linear relationship between the measurements and the amplitudes of the actual components, alleviating the need for a sampling based method, e.g., Markov Chain Monte Carlo (MCMC) or matched filtering, which are computationally intensive and cumbersome in a multi-dimensional parameter space. Using this formalism we could also study how many independent components can be separated using a given dataset from a network of current and upcoming ground based interferometric detectors.

Fast gravitational wave radiometry using data folding

Gravitational Waves (GWs) from the early universe and unresolved astrophysical sources are expected to create a stochastic GW background (SGWB). The GW radiometer algorithm is well suited to probe such a background using data from ground based laser interferometric detectors. Radiometer analysis can be performed in different bases, e.g., isotropic, pixel or spherical harmonic. Each of these analyses possesses a common temporal symmetry which we exploit here to fold the whole dataset for every detector pair, typically a few hundred to a thousand days of data, to only one sidereal day, without any compromise in precision. We develop the algebra and a software pipeline needed to fold data, accounting for the effect of overlapping windows and non-stationary noise. We implement this on LIGO's fifth science run data and validate it by performing a standard anisotropic SGWB search on both folded and unfolded data. Folded data not only leads to orders of magnitude reduction in computation cost, but it results in a conveniently small data volume of few gigabytes, making it possible to perform an actual analysis on a personal computer, as well as easy movement of data. A few important analyses, yet unaccomplished due to computational limitations, will now become feasible. Folded data, being independent of the radiometer basis, will also be useful in reducing processing redundancies in multiple searches and provide a common ground for mutual consistency checks. Most importantly, folded data will allow vast amount of experimentation with existing searches and provide substantial help in developing new strategies to find unknown sources.

The sensitivity of Cherenkov telescopes to dark matter and astrophysical anisotropies in the diffuse gamma-ray background

In this article, the capability of present (H.E.S.S., MAGIC, VERITAS) and planned (CTA) ground-based Cherenkov telescope systems for detecting angular anisotropies in the diffuse gamma-ray background is investigated. Following up on a study of the impact of instrumental characteristics (effective area, field of view, angular resolution, and background rejection efficiency), the first part examines the influence of different observational strategies, i.e. whether a single deep observation or a splitting over multiple shallow fields is preferred. In the second part, the sensitivity to anisotropies generated by self-annihilating dark matter is studied for different common dark matter models. We find that a relative contribution of ∼ 10% from dark matter annihilation to the extra-galactic diffuse gamma-ray background can be detected with planned configurations of CTA. In terms of the thermally-averaged self-annihilation cross section, the sensitivity of CTA corresponds to values below the thermal freeze-out expectation ⟨σv⟩ = 3 × 10−26 cm3 s−1 for dark matter particles lighter than ∼ 200 GeV. We stress the importance of constraining anisotropies from unresolved astrophysical sources with currently operating instruments already, as a novel and complementary method for investigating the properties of TeV sources.

Dark matter annihilation or unresolved astrophysical sources? Anisotropy probe of the origin of the cosmic gamma-ray background

The origin of the cosmic gamma-ray background (CGB) is a longstanding mystery in high-energy astrophysics. Possible candidates include ordinary astrophysical objects such as unresolved blazars, as well as more exotic processes such as dark matter annihilation. While it would be difficult to distinguish them from the mean intensity data alone, one can use anisotropy data instead. We investigate the CGB anisotropy both from unresolved blazars and dark matter annihilation (including contributions from dark matter substructures), and we find that the angular power spectra from these sources are very different. We then focus on detectability of dark matter annihilation signals using the anisotropy data by treating the unresolved blazar component as a known background. We find that the dark matter signature should be detectable in the angular power spectrum of the CGB from two-year all-sky observations with the Gamma Ray Large Area Space Telescope (GLAST), as long as the dark matter annihilation contributes to a reasonable fraction, e.g., $\ensuremath{\gtrsim}0.3$, of the CGB at around 10 GeV. We conclude that the anisotropy measurement of the CGB with GLAST should be a powerful tool for revealing the CGB origin, and potentially for the first detection of dark matter annihilation.

Relic neutrino masses and the highest energy cosmic rays

We consider the possibility that a large fraction of the ultrahigh energy cosmic rays are decay products of Z bosons which were produced in the scattering of ultrahigh energy cosmic neutrinos on cosmological relic neutrinos. We compare the observed ultrahigh energy cosmic ray spectrum with the one predicted in the above Z-burst scenario and determine the required mass of the heaviest relic neutrino as well as the necessary ultrahigh energy cosmic neutrino flux via a maximum likelihood analysis. We show that the value of the neutrino mass obtained in this way is fairly robust against variations in presently unknown quantities, like the amount of neutrino clustering, the universal radio background, and the extragalactic magnetic field, within their anticipated uncertainties. Much stronger systematics arises from different possible assumptions about the diffuse background of ordinary cosmic rays from unresolved astrophysical sources. In the most plausible case that these ordinary cosmic rays are protons of extragalactic origin, one is lead to a required neutrino mass in the range 0.08 eV - 1.3 eV at the 68 % confidence level. This range narrows down considerably if a particular universal radio background is assumed, e.g. to 0.08 eV - 0.40 eV for a large one. The required flux of ultrahigh energy cosmic neutrinos near the resonant energy should be detected in the near future by AMANDA, RICE, and the Pierre Auger Observatory, otherwise the Z-burst scenario will be ruled out.

Mapping thegravitational-wave background

The gravitational-wave sky is expected tohave isolated bright sources superimposedon a diffuse gravitational-wavebackground. The background radiation hastwo components: a confusion limitedbackground from unresolved astrophysicalsources; and a cosmological componentformed during the birth of the universe. Amap of the gravitational-wave backgroundcan be made by sweeping agravitational-wave detector across thesky. The detector output is a complicatedconvolution of the sky luminositydistribution, the detector responsefunction and the scan pattern. Here westudy the general deconvolution problem,and show how LIGO (laser interferometricgravitational observatory) and LISA (laserinterferometer space antenna) can be usedto detect anisotropies in thegravitational-wave background.