Source author record

J. Vink

J. Vink 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

72works
5topics
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

72 published item(s)

preprint2026arXiv

Measuring the muon content of inclined air showers using AERA and the water-Cherenkov detectors of the Pierre Auger Observatory

We present a novel approach for assessing the muon content of air showers with large zenith angles on a combined analysis of their radio emission and particle footprint. We use the radiation energy reconstructed by the Auger Engineering Radio Array (AERA) as an energy estimator and determine the muon number independently with the water-Cherenkov detector array of the Pierre Auger Observatory, deployed on a 1500 m grid. We focus our analysis on air showers with primary energy above 4 EeV to ensure full detection efficiency. Over approximately ten years of accumulated data, we identify a set of 40 high-quality events that are used in the analysis. The estimated muon contents in data are compatible with those for iron primaries as predicted by current-generation hadronic interaction models. This result can be interpreted as a deficit of muons in simulations as a lighter mass composition has been established from Xmax measurements. This muon deficit was already observed in previous analyses of the Auger Collaboration and is confirmed using hybrid events that include radio measurements for the first time.

preprint2022arXiv

A search for photons with energies above $2{\times}10^{17}$ eV using hybrid data from the low-energy extensions of the Pierre Auger Observatory

Ultra-high-energy photons with energies exceeding $10^{17}$ eV offer a wealth of connections to different aspects of cosmic-ray astrophysics as well as to gamma-ray and neutrino astronomy. The recent observations of photons with energies in the $10^{15}$ eV range further motivate searches for even higher-energy photons. In this paper, we present a search for photons with energies exceeding $2{\times}10^{17}$ eV using about 5.5 years of hybrid data from the low-energy extensions of the Pierre Auger Observatory. The upper limits on the integral photon flux derived here are the most stringent ones to date in the energy region between $10^{17}$ and $10^{18}$ eV.

preprint2022arXiv

Accretion geometry of the neutron star low mass X-ray binary Cyg X-2 from X-ray polarization measurements

We report spectro-polarimetric results of an observational campaign of the bright neutron star low-mass X-ray binary Cyg X-2 simultaneously observed by IXPE, NICER and INTEGRAL. Consistently with previous results, the broad-band spectrum is characterized by a lower-energy component, attributed to the accretion disc with $kT_{\rm in} \approx$ 1 keV, plus unsaturated Comptonization in thermal plasma with temperature $kT_{\rm e} = 3$ keV and optical depth $τ\approx 4$, assuming a slab geometry. We measure the polarization degree in the 2-8 keV band $P=1.8 \pm 0.3$ per cent and polarization angle $ϕ= 140^{\circ} \pm 4^{\circ}$, consistent with the previous X-ray polarimetric measurements by OSO-8 as well as with the direction of the radio jet which was earlier observed from the source. While polarization of the disc spectral component is poorly constrained with the IXPE data, the Comptonized emission has a polarization degree $P =4.0 \pm 0.7$ per cent and a polarization angle aligned with the radio jet. Our results strongly favour a spreading layer at the neutron star surface as the main source of the polarization signal. However, we cannot exclude a significant contribution from reflection off the accretion disc, as indicated by the presence of the iron fluorescence line.

preprint2022arXiv

Completing the X-ray view of the recently discovered supernova remnant G53.41+0.03

Aims: We present a detailed X-ray study of the recently discovered supernova remnant (SNR) G53.41+0.03 that follows up and further expands on the previous, limited analysis of archival data covering a small portion of the SNR. Methods: With the new dedicated 70ks XMM-Newton observation we investigate the morphological structure of the SNR in X-rays, search for a presence of a young neutron star and characterise the plasma conditions in the selected regions by means of spectral fitting. Results: The first full view of SNR G53.41+0.03 shows an X-ray emission region well aligned with the reported half-shell radio morphology. We find three distinct regions of the remnant that vary in brightness and hardness of the spectra, and are all best characterised by a hot plasma model in a non-equilibrium ionisation state. Of the three regions, the brightest one contains the most mature plasma, with ionisation age $τ\approx 4\times10^{10}$s cm$^{-3}$ (where $τ= n_e t$), a lower electron temperature of kT$_\mathrm{e} \approx 1$ keV and the highest estimated gas density of n$_\mathrm{H}\approx 0.87$ cm$^{-3}$. The second, fainter but spectrally harder, region reveals a younger plasma ($τ\approx 1.7\times10^{10}$s cm$^{-3}$) with higher temperature (kT$_\mathrm{e} \approx 2$ keV) and two to three times lower density (n$_\mathrm{H}\approx 0.34$ cm$^{-3}$). The third region is very faint, but we identify spectroscopically the presence of a hot plasma.Employing several methods for age estimation, we find the remnant to be $t \approx 1000-5000$ yrs old, confirming the earlier reports of a relatively young age. The environment of the remnant also contains a number of point sources, of which most are expected to be positioned in the foreground. Of the two point sources in the geometrical centre of the remnant one is consistent with the characteristics of a young neutron star.

preprint2022arXiv

Investigating Hadronic Interactions at Ultra-High Energies with the Pierre Auger Observatory

The development of an extensive air shower depends not only on the nature of the primary ultra-high-energy cosmic ray but also on the properties of the hadronic interactions. For energies above those achievable in human-made accelerators, hadronic interactions are only accessible through the studies of extensive air showers, which can be measured at the Pierre Auger Observatory. With its hybrid detector design, the Pierre Auger Observatory measures both the longitudinal development of showers in the atmosphere and the lateral distribution of particles that arrive at the ground. This way, observables that are sensitive to hadronic interactions at ultra-high energies can be obtained. While the hadronic interaction cross-section can be assessed from the longitudinal profiles, the number of muons and their fluctuations measured with the ground detectors are linked to other physical properties. In addition to these direct studies, we discuss here how measurements of the atmospheric depth of the maximum of air-shower profiles and the characteristics of the muon signal at the ground can be used to test the self-consistency of the post-LHC hadronic models.

preprint2022arXiv

Search for dark matter annihilation signals in the H.E.S.S. Inner Galaxy Survey

The central region of the Milky Way is one of the foremost locations to look for dark matter (DM) signatures. We report the first results on a search for DM particle annihilation signals using new observations from an unprecedented gamma-ray survey of the Galactic Center (GC) region, ${\it i.e.}$, the Inner Galaxy Survey, at very high energies ($\gtrsim$ 100 GeV) performed with the H.E.S.S. array of five ground-based Cherenkov telescopes. No significant gamma-ray excess is found in the search region of the 2014-2020 dataset and a profile likelihood ratio analysis is carried out to set exclusion limits on the annihilation cross section $\langle σv\rangle$. Assuming Einasto and Navarro-Frenk-White (NFW) DM density profiles at the GC, these constraints are the strongest obtained so far in the TeV DM mass range. For the Einasto profile, the constraints reach $\langle σv\rangle$ values of $\rm 3.7\times10^{-26} cm^3s^{-1}$ for 1.5 TeV DM mass in the $W^+W^-$ annihilation channel, and $\rm 1.2 \times 10^{-26} cm^3s^{-1}$ for 0.7 TeV DM mass in the $τ^+τ^-$ annihilation channel. With the H.E.S.S. Inner Galaxy Survey, ground-based $γ$-ray observations thus probe $\langle σv\rangle$ values expected from thermal-relic annihilating TeV DM particles.

preprint2022arXiv

Search for Spatial Correlations of Neutrinos with Ultra-High-Energy Cosmic Rays

For several decades, the origin of ultra-high-energy cosmic rays (UHECRs) has been an unsolved question of high-energy astrophysics. One approach for solving this puzzle is to correlate UHECRs with high-energy neutrinos, since neutrinos are a direct probe of hadronic interactions of cosmic rays and are not deflected by magnetic fields. In this paper, we present three different approaches for correlating the arrival directions of neutrinos with the arrival directions of UHECRs. The neutrino data is provided by the IceCube Neutrino Observatory and ANTARES, while the UHECR data with energies above $\sim$50 EeV is provided by the Pierre Auger Observatory and the Telescope Array. All experiments provide increased statistics and improved reconstructions with respect to our previous results reported in 2015. The first analysis uses a high-statistics neutrino sample optimized for point-source searches to search for excesses of neutrinos clustering in the vicinity of UHECR directions. The second analysis searches for an excess of UHECRs in the direction of the highest-energy neutrinos. The third analysis searches for an excess of pairs of UHECRs and highest-energy neutrinos on different angular scales. None of the analyses has found a significant excess, and previously reported over-fluctuations are reduced in significance. Based on these results, we further constrain the neutrino flux spatially correlated with UHECRs.

preprint2022arXiv

Testing effects of Lorentz invariance violation in the propagation of astroparticles with the Pierre Auger Observatory

Lorentz invariance violation (LIV) is often described by dispersion relations of the form $E_i^2=m_i^2+p_i^2+δ_{i,n} E^{2+n}$ with delta different based on particle type $i$, with energy $E$, momentum $p$ and rest mass $m$. Kinematics and energy thresholds of interactions are modified once the LIV terms become comparable to the squared masses of the particles involved. Thus, the strongest constraints on the LIV coefficients $δ_{i,n}$ tend to come from the highest energies. At sufficiently high energies, photons produced by cosmic ray interactions as they propagate through the Universe could be subluminal and unattenuated over cosmological distances. Cosmic ray interactions can also be modified and lead to detectable fingerprints in the energy spectrum and mass composition observed on Earth. The data collected at the Pierre Auger Observatory are therefore possibly sensitive to both the electromagnetic and hadronic sectors of LIV. In this article, we explore these two sectors by comparing the energy spectrum and the composition of cosmic rays and the upper limits on the photon flux from the Pierre Auger Observatory with simulations including LIV. Constraints on LIV parameters depend strongly on the mass composition of cosmic rays at the highest energies. For the electromagnetic sector, while no constraints can be obtained in the absence of protons beyond $10^{19}$ eV, we obtain $δ_{γ,0} > -10^{-21}$, $δ_{γ,1} > -10^{-40}$ eV$^{-1}$ and $δ_{γ,2} > -10^{-58}$ eV$^{-2}$ in the case of a subdominant proton component up to $10^{20}$ eV. For the hadronic sector, we study the best description of the data as a function of LIV coefficients and we derive constraints in the hadronic sector such as $δ_{\mathrm{had},0} < 10^{-19}$, $δ_{\mathrm{had},1} < 10^{-38}$ eV$^{-1}$ and $δ_{\mathrm{had},2}< 10^{-57}$ eV$^{-2}$ at 5$σ$ CL.

preprint2022arXiv

The energy spectrum of cosmic rays beyond the turn-down around $10^{17}$ eV as measured with the surface detector of the Pierre Auger Observatory

We present a measurement of the cosmic-ray spectrum above 100\,PeV using the part of the surface detector of the Pierre Auger Observatory that has a spacing of 750~m. An inflection of the spectrum is observed, confirming the presence of the so-called \emph{second-knee} feature. The spectrum is then combined with that of the 1500\,m array to produce a single measurement of the flux, linking this spectral feature with the three additional breaks at the highest energies. The combined spectrum, with an energy scale set calorimetrically via fluorescence telescopes and using a single detector type, results in the most statistically and systematically precise measurement of spectral breaks yet obtained. These measurements are critical for furthering our understanding of the highest energy cosmic rays.

preprint2022arXiv

The LOFAR Two-metre Sky Survey -- V. Second data release

In this data release from the LOFAR Two-metre Sky Survey (LoTSS) we present 120-168MHz images covering 27% of the northern sky. Our coverage is split into two regions centred at approximately 12h45m +44$^\circ$30' and 1h00m +28$^\circ$00' and spanning 4178 and 1457 square degrees respectively. The images were derived from 3,451hrs (7.6PB) of LOFAR High Band Antenna data which were corrected for the direction-independent instrumental properties as well as direction-dependent ionospheric distortions during extensive, but fully automated, data processing. A catalogue of 4,396,228 radio sources is derived from our total intensity (Stokes I) maps, where the majority of these have never been detected at radio wavelengths before. At 6" resolution, our full bandwidth Stokes I continuum maps with a central frequency of 144MHz have: a median rms sensitivity of 83$μ$Jy/beam; a flux density scale accuracy of approximately 10%; an astrometric accuracy of 0.2"; and we estimate the point-source completeness to be 90% at a peak brightness of 0.8mJy/beam. By creating three 16MHz bandwidth images across the band we are able to measure the in-band spectral index of many sources, albeit with an error on the derived spectral index of +/-0.2 which is a consequence of our flux-density scale accuracy and small fractional bandwidth. Our circular polarisation (Stokes V) 20" resolution 120-168MHz continuum images have a median rms sensitivity of 95$μ$Jy/beam, and we estimate a Stokes I to Stokes V leakage of 0.056%. Our linear polarisation (Stokes Q and Stokes U) image cubes consist of 480 x 97.6 kHz wide planes and have a median rms sensitivity per plane of 10.8mJy/beam at 4' and 2.2mJy/beam at 20"; we estimate the Stokes I to Stokes Q/U leakage to be approximately 0.2%. Here we characterise and publicly release our Stokes I, Q, U and V images in addition to the calibrated uv-data.

preprint2022arXiv

Time-resolved hadronic particle acceleration in the recurrent Nova RS Ophiuchi

Recurrent Novae are repeating thermonuclear explosions in the outer layers of white dwarfs, due to the accretion of fresh material from a binary companion. The shock generated by ejected material slamming into the companion star's wind, accelerates particles to very-high-energies. We report very-high-energy (VHE, $\gtrsim100$\,GeV) gamma rays from the recurrent nova RS\,Ophiuchi up to a month after its 2021 outburst, using the High Energy Stereoscopic System. The VHE emission has a similar temporal profile to lower-energy GeV emission, indicating a common origin, with a two-day delay in peak flux. These observations constrain models of time-dependent particle energization, favouring a hadronic emission scenario over the leptonic alternative. This confirms that shocks in dense winds provide favourable environments for efficient cosmic-ray acceleration to very-high-energies.

preprint2021arXiv

A MeerKAT, e-MERLIN, H.E.S.S. and Swift search for persistent and transient emission associated with three localised FRBs

We report on a search for persistent radio emission from the one-off Fast Radio Burst (FRB) 20190714A, as well as from two repeating FRBs, 20190711A and 20171019A, using the MeerKAT radio telescope. For FRB 20171019A we also conducted simultaneous observations with the High Energy Stereoscopic System (H.E.S.S.) in very high energy gamma rays and searched for signals in the ultraviolet, optical, and X-ray bands. For this FRB, we obtain a UV flux upper limit of 1.39x10^-16 erg/cm^-2/s/Amstrong, X-ray limit of ~ 6.6x10^-14 erg/cm^-2/s and a limit on the very-high-energy gamma-ray flux (Phi) (E > 120 GeV) < 1.7 x 10^-12 erg/cm^-2/s. We obtain a radio upper limit of ~15 microJy/beam for persistent emission at the locations of both FRBs 20190711A and 20171019A, but detect diffuse radio emission with a peak brightness of ~53 microJy/beam associated with FRB 20190714A at z = 0.2365. This represents the first detection of the radio continuum emission potentially associated with the host (galaxy) of FRB 20190714A, and is only the third known FRB to have such an association. Given the possible association of a faint persistent source, FRB 20190714A may potentially be a repeating FRB whose age lies between that of FRB 20121102A and FRB 20180916A. A parallel search for repeat bursts from these FRBs revealed no new detections down to a fluence of 0.08 Jy ms for a 1 ms duration burst.

preprint2021arXiv

Design, upgrade and characterization of the silicon photomultiplier front-end for the AMIGA detector at the Pierre Auger Observatory

AMIGA (Auger Muons and Infill for the Ground Array) is an upgrade of the Pierre Auger Observatory to complement the study of ultra-high-energy cosmic rays (UHECR) by measuring the muon content of extensive air showers (EAS). It consists of an array of 61 water Cherenkov detectors on a denser spacing in combination with underground scintillation detectors used for muon density measurement. Each detector is composed of three scintillation modules, with 10 m$^2$ detection area per module, buried at 2.3 m depth, resulting in a total detection area of 30 m$^2$. Silicon photomultiplier sensors (SiPM) measure the amount of scintillation light generated by charged particles traversing the modules. In this paper, the design of the front-end electronics to process the signals of those SiPMs and test results from the laboratory and from the Pierre Auger Observatory are described. Compared to our previous prototype, the new electronics shows a higher performance, higher efficiency and lower power consumption, and it has a new acquisition system with increased dynamic range that allows measurements closer to the shower core. The new acquisition system is based on the measurement of the total charge signal that the muonic component of the cosmic ray shower generates in the detector.

preprint2021arXiv

Multi-messenger-Athena Synergy White Paper

In this paper we explore the scientific synergies between Athena and some of the key multi-messenger facilities that should be operative concurrently with Athena. These facilities include LIGO A+, Advanced Virgo+ and future detectors for ground-based observation of gravitational waves (GW), LISA for space-based observations of GW, IceCube and KM3NeT for neutrino observations, and CTA for very high energy observations. These science themes encompass pressing issues in astrophysics, cosmology and fundamental physics such as: the central engine and jet physics in compact binary mergers, accretion processes and jet physics in Super-Massive Binary Black Holes (SMBBHs) and in compact stellar binaries, the equation of state of neutron stars, cosmic accelerators and the origin of Cosmic Rays (CRs), the origin of intermediate and high-Z elements in the Universe, the Cosmic distance scale and tests of General Relativity and the Standard Model. Observational strategies for implementing the identified science topics are also discussed. A significant part of the sources targeted by multi-messenger facilities is of transient nature. We have thus also discussed the synergy of \textsl{Athena} with wide-field high-energy facilities, taking THESEUS as a case study for transient discovery. This discussion covers all the Athena science goals that rely on follow-up observations of high-energy transients identified by external observatories, and includes also topics that are not based on multi-messenger observations, such as the search for missing baryons or the observation of early star populations and metal enrichment at the cosmic dawn with Gamma-Ray Bursts (GRBs).

preprint2021arXiv

Observation of the gamma-ray binary HESS J0632+057 with the H.E.S.S., MAGIC, and VERITAS telescopes

The results of gamma-ray observations of the binary system HESS J0632+057 collected during 450 hours over 15 years, between 2004 and 2019, are presented. Data taken with the atmospheric Cherenkov telescopes H.E.S.S., MAGIC, and VERITAS at energies above 350 GeV were used together with observations at X-ray energies obtained with Swift-XRT, Chandra, XMM-Newton, NuSTAR, and Suzaku. Some of these observations were accompanied by measurements of the Hα emission line. A significant detection of the modulation of the VHE gamma-ray fluxes with a period of 316.7+-4.4 days is reported, consistent with the period of 317.3+-0.7 days obtained with a refined analysis of X-ray data. The analysis of data of four orbital cycles with dense observational coverage reveals short timescale variability, with flux-decay timescales of less than 20 days at very high energies. Flux variations observed over the time scale of several years indicate orbit-to-orbit variability. The analysis confirms the previously reported correlation of X-ray and gamma-ray emission from the system at very high significance, but can not find any correlation of optical Hα parameters with X-ray or gamma-ray energy fluxes in simultaneous observations. The key finding is that the emission of HESS J0632+057 in the X-ray and gamma-ray energy bands is highly variable on different time scales. The ratio of gamma-ray to X-ray flux shows the equality or even dominance of the gamma-ray energy range. This wealth of new data is interpreted taking into account the insufficient knowledge of the ephemeris of the system, and discussed in the context of results reported on other gamma-ray binary systems.

preprint2021arXiv

Sensitivity of the Cherenkov Telescope Array to a dark matter signal from the Galactic centre

We provide an updated assessment of the power of the Cherenkov Telescope Array (CTA) to search for thermally produced dark matter at the TeV scale, via the associated gamma-ray signal from pair-annihilating dark matter particles in the region around the Galactic centre. We find that CTA will open a new window of discovery potential, significantly extending the range of robustly testable models given a standard cuspy profile of the dark matter density distribution. Importantly, even for a cored profile, the projected sensitivity of CTA will be sufficient to probe various well-motivated models of thermally produced dark matter at the TeV scale. This is due to CTA's unprecedented sensitivity, angular and energy resolutions, and the planned observational strategy. The survey of the inner Galaxy will cover a much larger region than corresponding previous observational campaigns with imaging atmospheric Cherenkov telescopes. CTA will map with unprecedented precision the large-scale diffuse emission in high-energy gamma rays, constituting a background for dark matter searches for which we adopt state-of-the-art models based on current data. Throughout our analysis, we use up-to-date event reconstruction Monte Carlo tools developed by the CTA consortium, and pay special attention to quantifying the level of instrumental systematic uncertainties, as well as background template systematic errors, required to probe thermally produced dark matter at these energies. "Full likelihood tables complementing our analysis are provided here [ https://doi.org/10.5281/zenodo.4057987 ]"

preprint2021arXiv

Synthesis of radioactive elements in novae and supernovae and their use as a diagnostic tool

Novae and supernovae play a key role in many fields of Astrophysics and Cosmology. Despite their importance, an accurate description of which objects explode and why and how they explode is still lacking. One of the main characteristics of such explosions is that they are the main suppliers of newly synthesized chemical elements in the Galaxy. Since some of these isotopes are radioactive, it is possible to use the corresponding gamma-rays as a diagnostic tool of the explosion thanks to their independence on the thermal state of the debris. The drawback is the poor sensitivity of detectors in the MeV energy domain. As a consequence, the radioactive lines have only been detected in one core collapse supernova (SN 1987A), one Type Ia supernova (SN 2014J), and one supernova remnant (Cas A). Nevertheless these observations have provided and are providing important information about the explosion mechanisms. Unfortunately, novae are still eluding detection. These results emphasize the necessity to place as soon as possible a new instrument in orbit with enough sensitivity to noticeably enlarge the sample of detected events

preprint2020arXiv

An extreme particle accelerator in the Galactic plane: HESS J1826$-$130

The unidentified very-high-energy (VHE; E $>$ 0.1 TeV) $γ$-ray source, HESS J1826$-$130, was discovered with the High Energy Stereoscopic System (HESS) in the Galactic plane. The analysis of 215 h of HESS data has revealed a steady $γ$-ray flux from HESS J1826$-$130, which appears extended with a half-width of 0.21$^{\circ}$ $\pm$ 0.02$^{\circ}_{\text{stat}}$ $\pm$ 0.05$^{\circ}_{\text{sys}}$. The source spectrum is best fit with either a power-law function with a spectral index $Γ$ = 1.78 $\pm$ 0.10$_{\text{stat}}$ $\pm$ 0.20$_{\text{sys}}$ and an exponential cut-off at 15.2$^{+5.5}_{-3.2}$ TeV, or a broken power-law with $Γ_{1}$ = 1.96 $\pm$ 0.06$_{\text{stat}}$ $\pm$ 0.20$_{\text{sys}}$, $Γ_{2}$ = 3.59 $\pm$ 0.69$_{\text{stat}}$ $\pm$ 0.20$_{\text{sys}}$ for energies below and above $E_{\rm{br}}$ = 11.2 $\pm$ 2.7 TeV, respectively. The VHE flux from HESS J1826$-$130 is contaminated by the extended emission of the bright, nearby pulsar wind nebula (PWN), HESS J1825$-$137, particularly at the low end of the energy spectrum. Leptonic scenarios for the origin of HESS J1826$-$130 VHE emission related to PSR J1826$-$1256 are confronted by our spectral and morphological analysis. In a hadronic framework, taking into account the properties of dense gas regions surrounding HESS J1826$-$130, the source spectrum would imply an astrophysical object capable of accelerating the parent particle population up to $\gtrsim$200 TeV. Our results are also discussed in a multiwavelength context, accounting for both the presence of nearby supernova remnants (SNRs), molecular clouds, and counterparts detected in radio, X-rays, and TeV energies.

preprint2020arXiv

Cassiopeia A, Cygnus A, Taurus A, and Virgo A at ultra-low radio frequencies

The four persistent radio sources in the northern sky with the highest flux density at metre wavelengths are Cassiopeia A, Cygnus A, Taurus A, and Virgo A; collectively they are called the A-team. Their flux densities at ultra-low frequencies (<100 MHz) can reach several thousands of janskys, and they often contaminate observations of the low-frequency sky by interfering with image processing. Furthermore, these sources are foreground objects for all-sky observations hampering the study of faint signals, such as the cosmological 21 cm line from the epoch of reionisation. We aim to produce robust models for the surface brightness emission as a function of frequency for the A-team sources at ultra-low frequencies. These models are needed for the calibration and imaging of wide-area surveys of the sky with low-frequency interferometers. This requires obtaining images at an angular resolution better than 15 arcsec with a high dynamic range and good image fidelity. We observed the A-team with the Low Frequency Array (LOFAR) at frequencies between 30 MHz and 77 MHz using the Low Band Antenna (LBA) system. We reduced the datasets and obtained an image for each A-team source. The paper presents the best models to date for the sources Cassiopeia A, Cygnus A, Taurus A, and Virgo A between 30 MHz and 77 MHz. We were able to obtain the aimed resolution and dynamic range in all cases. Owing to its compactness and complexity, observations with the long baselines of the International LOFAR Telescope will be required to improve the source model for Cygnus A further.

preprint2020arXiv

Cosmic-ray anisotropies in right ascension measured by the Pierre Auger Observatory

We present measurements of the large-scale cosmic-ray anisotropies in right ascension, using data collected by the surface detector array of the Pierre Auger Observatory over more than 14 years. We determine the equatorial dipole component, $\vec{d}_\perp$, through a Fourier analysis in right ascension that includes weights for each event so as to account for the main detector-induced systematic effects. For the energies at which the trigger efficiency of the array is small, the ``East-West'' method is employed. Besides using the data from the array with detectors separated by 1500 m, we also include data from the smaller but denser sub-array of detectors with 750 m separation, which allows us to extend the analysis down to $\sim 0.03$ EeV. The most significant equatorial dipole amplitude obtained is that in the cumulative bin above 8~EeV, $d_\perp=6.0^{+1.0}_{-0.9}$%, which is inconsistent with isotropy at the 6$σ$ level. In the bins below 8 EeV, we obtain 99% CL upper-bounds on $d_\perp$ at the level of 1 to 3 percent. At energies below 1 EeV, even though the amplitudes are not significant, the phases determined in most of the bins are not far from the right ascension of the Galactic center, at $α_{\rm GC}=-94^\circ$, suggesting a predominantly Galactic origin for anisotropies at these energies. The reconstructed dipole phases in the energy bins above 4 EeV point instead to right ascensions that are almost opposite to the Galactic center one, indicative of an extragalactic cosmic ray origin.

preprint2020arXiv

Detection of very-high-energy γ-ray emission from the colliding wind binary η Car with H.E.S.S

Aims. Colliding wind binary systems have long been suspected to be high-energy (HE; 100 MeV < E < 100 GeV) γ-ray emitters. η Car is the most prominent member of this object class and is confirmed to emit phase-locked HE γ rays from hundreds of MeV to ~100 GeV energies. This work aims to search for and characterise the very-high-energy (VHE; E >100 GeV) γ-ray emission from η Car around the last periastron passage in 2014 with the ground-based High Energy Stereoscopic System (H.E.S.S.). Methods. The region around η Car was observed with H.E.S.S. between orbital phase p = 0.78 - 1.10, with a closer sampling at p {\approx} 0.95 and p {\approx} 1.10 (assuming a period of 2023 days). Optimised hardware settings as well as adjustments to the data reduction, reconstruction, and signal selection were needed to suppress and take into account the strong, extended, and inhomogeneous night sky background (NSB) in the η Car field of view. Tailored run-wise Monte-Carlo simulations (RWS) were required to accurately treat the additional noise from NSB photons in the instrument response functions. Results. H.E.S.S. detected VHE γ-ray emission from the direction of η Car shortly before and after the minimum in the X-ray light-curve close to periastron. Using the point spread function provided by RWS, the reconstructed signal is point-like and the spectrum is best described by a power law. The overall flux and spectral index in VHE γ rays agree within statistical and systematic errors before and after periastron. The γ-ray spectrum extends up to at least ~400 GeV. This implies a maximum magnetic field in a leptonic scenario in the emission region of 0.5 Gauss. No indication for phase-locked flux variations is detected in the H.E.S.S. data.

preprint2020arXiv

Mapping the spectral index of Cas A: evidence for flattening from radio to infrared

Synchrotron radiation from supernova remnants are caused by electrons accelerated through diffusive shock acceleration (DSA). The standard DSA theory predicts an electron spectral index of $p=2$, corresponding to a radio spectral index of $α=-0.5$. An extension of DSA predicts that the accelerated particles changes the shock structure, resulting in a spectrum that is steeper than $p>2$ ($α<-0.5$) at low energies and flattens with energy. For Cas A, a synchrotron spectral flattening was previously reported for a small part of the remnant in the mid-infrared regime. Here, we present new measurements for spectral flattening using archival radio (4.72~GHz) and mid-infrared (3.6~$μ$m) data, and produce a complete spectral index map to investigate the spatial variations within the remnant. We compare this to measurements of the radio spectral index from L-(1.285~GHz) and C-(4.64~GHz) band maps. Our result shows overall spectral flattening across the remnant ($α_\mathrm{R-IR} \sim -0.5$ to $-0.7$), to be compared to the radio spectral index of $α_{\rm R}=-0.77$. The flattest values coincide with the locations of most recent particle acceleration. In addition to overall flattening, we detect a relatively steeper region in the southeast of the remnant ($α_\mathrm{R-IR} \sim -0.67$). We explore whether these locally steeper spectra could be the result of synchrotron cooling, which provides constraints on the local magnetic-field strengths/age of the plasma, suggesting $B\lesssim 2$~mG for an age of 100~yr, and even $B\lesssim 1$~mG using the age of Cas A, in agreement with other estimates.

preprint2020arXiv

Probing the magnetic field in the GW170817 outflow using H.E.S.S. observations

The detection of the first electromagnetic counterpart to the binary neutron star (BNS) merger remnant GW170817 established the connection between short $γ$-ray bursts and BNS mergers. It also confirmed the forging of heavy elements in the ejecta (a so-called kilonova) via the r-process nucleosynthesis. The appearance of non-thermal radio and X-ray emission, as well as the brightening, which lasted more than 100 days, were somewhat unexpected. Current theoretical models attempt to explain this temporal behavior as either originating from a relativistic off-axis jet or a kilonova-like outflow. In either scenario, there is some ambiguity regarding how much energy is transported in the non-thermal electrons versus the magnetic field of the emission region. Combining the VLA (radio) and Chandra (X-ray) measurements with observations in the GeV-TeV domain can help break this ambiguity, almost independently of the assumed origin of the emission. Here we report for the first time on deep H.E.S.S. observations of GW170817 / GRB 170817A between 124 and 272 days after the BNS merger with the full H.E.S.S. array of telescopes, as well as on an updated analysis of the prompt (<5 days) observations with the upgraded H.E.S.S. phase-I telescopes. We discuss implications of the H.E.S.S. measurement for the magnetic field in the context of different source scenarios.

preprint2020arXiv

Resolving acceleration to very high energies along the Jet of Centaurus A

The nearby radio galaxy Centaurus A belongs to a class of Active Galaxies that are very luminous at radio wavelengths. The majority of these galaxies show collimated relativistic outflows known as jets, that extend over hundreds of thousands of parsecs for the most powerful sources. Accretion of matter onto the central super-massive black hole is believed to fuel these jets and power their emission, with the radio emission being related to the synchrotron radiation of relativistic electrons in magnetic fields. The origin of the extended X-ray emission seen in the kiloparsec-scale jets from these sources is still a matter of debate, although Cen A's X-ray emission has been suggested to originate in electron synchrotron processes. The other possible explanation is Inverse Compton (IC) scattering with CMB soft photons. Synchrotron radiation needs ultra-relativistic electrons ($\sim50$ TeV), and given their short cooling times, requires some continuous re-acceleration mechanism to be active. IC scattering, on the other hand, does not require very energetic electrons, but requires jets that stay highly relativistic on large scales ($\geq$1 Mpc) and that remain well-aligned with the line of sight. Some recent evidence disfavours inverse Compton-CMB models, although other evidence seems to be compatible with them. In principle, the detection of extended gamma-ray emission, directly probing the presence of ultra-relativistic electrons, could distinguish between these options, but instruments have hitherto been unable to resolve the relevant structures. At GeV energies there is also an unusual spectral hardening in Cen A, whose explanation is unclear. Here we report observations of Cen A at TeV energies that resolve its large-scale jet. We interpret the data as evidence for the acceleration of ultra-relativistic electrons in the jet, and favour the synchrotron explanation for the X-rays.

preprint2020arXiv

Search for dark matter signals towards a selection of recently-detected DES dwarf galaxy satellites of the Milky Way with H.E.S.S

Dwarf spheroidal galaxy satellites of the Milky Way are prime targets for indirect detection of dark matter with gamma rays due to their proximity, high dark matter content and absence of non-thermal emission processes. Recently, the Dark Energy Survey (DES) revealed the existence of new ultra-faint dwarf spheroidal galaxies in the southern-hemisphere sky, therefore ideally located for ground-based observations with the imaging atmospheric Cherenkov telescope array H.E.S.S. We present a search for very-high-energy ($E\gtrsim100$ GeV) gamma-ray emission using H.E.S.S. observations carried out recently towards Reticulum II, Tucana II, Tucana III, Tucana IV and Grus II satellites. No significant very-high-energy gamma-ray excess is found from the observations on any individual object nor in the combined analysis of all the datasets. Using the most recent modeling of the dark matter distribution in the dwarf galaxy halo, we compute for the first time on DES satellites individual and combined constraints from Cherenkov telescope observations on the annihilation cross section of dark matter particles in the form of Weakly Interacting Massive Particles. The combined 95% C.L. observed upper limits reach $\langle σv \rangle \simeq 1 \times 10^{-23}$ cm$^3$s$^{-1}$ in the $W^+W^-$ channel and $4 \times 10^{-26}$ cm$^3$s$^{-1}$ in the $γγ$ channels for a dark matter mass of 1.5 TeV. The H.E.S.S. constraints well complement the results from Fermi-LAT, HAWC, MAGIC and VERITAS and are currently the most stringent in the $γγ$ channels in the multi-GeV/multi-TeV mass range.

preprint2020arXiv

Search for magnetically-induced signatures in the arrival directions of ultra-high-energy cosmic rays measured at the Pierre Auger Observatory

We search for signals of magnetically-induced effects in the arrival directions of ultra-high-energy cosmic rays detected at the Pierre Auger Observatory. We apply two different methods. One is a search for sets of events that show a correlation between their arrival direction and the inverse of their energy, which would be expected if they come from the same point-like source, they have the same electric charge and their deflection is relatively small and coherent. We refer to these sets of events as "multiplets". The second method, called "thrust", is a principal axis analysis aimed to detect the elongated patterns in a region of interest. We study the sensitivity of both methods using a benchmark simulation and we apply them to data in two different searches. The first search is done assuming as source candidates a list of nearby active galactic nuclei and starburst galaxies. The second is an all-sky blind search. We report the results and we find no statistically significant features. We discuss the compatibility of these results with the indications on the mass composition inferred from data of the Pierre Auger Observatory.

preprint2020arXiv

Studies on the response of a water-Cherenkov detector of the Pierre Auger Observatory to atmospheric muons using an RPC hodoscope

Extensive air showers, originating from ultra-high energy cosmic rays, have been successfully measured through the use of arrays of water-Cherenkov detectors (WCDs). Sophisticated analyses exploiting WCD data have made it possible to demonstrate that shower simulations, based on different hadronic-interaction models, cannot reproduce the observed number of muons at the ground. The accurate knowledge of the WCD response to muons is paramount in establishing the exact level of this discrepancy. In this work, we report on a study of the response of a WCD of the Pierre Auger Observatory to atmospheric muons performed with a hodoscope made of resistive plate chambers (RPCs), enabling us to select and reconstruct nearly 600 thousand single muon trajectories with zenith angles ranging from 0$^\circ$ to 55$^\circ$. Comparison of distributions of key observables between the hodoscope data and the predictions of dedicated simulations allows us to demonstrate the accuracy of the latter at a level of 2%. As the WCD calibration is based on its response to atmospheric muons, the hodoscope data are also exploited to show the long-term stability of the procedure.

preprint2020arXiv

Very high energy $γ$-ray emission from two blazars of unknown redshift and upper limits on their distance

We report on the detection of very-high-energy (VHE; $E > 100$ GeV) $γ$-ray emission from the BL Lac objects KUV 00311-1938 and PKS 1440-389 with the High Energy Stereoscopic System (H.E.S.S.). H.E.S.S. observations were accompanied or preceded by multi-wavelength observations with Fermi/LAT, XRT and UVOT on board the Swift satellite, and ATOM. Based on an extrapolation of the Fermi/LAT spectrum towards the VHE $γ$-ray regime, we deduce a 95% confidence level upper limit on the unknown redshift of KUV 00311-1938 of z < 0.98, and of PKS 1440-389 of z < 0.53. When combined with previous spectroscopy results the redshift of KUV 00311-1938 is constrained to $0.51 \leq z < 0.98$ and for PKS 1440-389 to $0.14 \lessapprox z < 0.53$.

preprint2019arXiv

Constraining the cosmic ray spectrum in the vicinity of the supernova remnant W28: from sub-GeV to multi-TeV energies

Supernova remnants interacting with molecular clouds are ideal laboratories to study the acceleration of particles at shock waves and their transport and interactions in the surrounding interstellar medium. In this paper, we focus on the supernova remnant W28, which over the years has been observed in all energy domains from radio waves to very-high-energy gamma rays. The bright gamma-ray emission detected from molecular clouds located in its vicinity revealed the presence of accelerated GeV and TeV particles in the region. An enhanced ionization rate has also been measured by means of millimetre observations, but such observations alone cannot tell us whether the enhancement is due to low energy (MeV) cosmic rays (either protons or electrons) or the X-ray photons emitted by the shocked gas. The goal of this study is to determine the origin of the enhanced ionization rate and to infer from multiwavelength observations the spectrum of cosmic rays accelerated at the supernova remnant shock in the unprecedented range spanning from MeV to multi-TeV particle energies. We developed a model to describe the transport of X-ray photons into the molecular cloud, and we fitted the radio, millimeter, and gamma-ray data to derive the spectrum of the radiating particles. The contribution from X-ray photons to the enhanced ionization rate is negligible, and therefore the ionization must be due to cosmic rays. Even though we cannot exclude a contribution to the ionization rate coming from cosmic ray electrons, we show that a scenario where cosmic ray protons explain both the gamma-ray flux and the enhanced ionization rate provides the most natural fit to multiwavelength data. This strongly suggests that the intensity of CR protons is enhanced in the region for particle energies in a very broad range covering almost 6 orders of magnitude: from $\lesssim 100$ MeV up to several tens of TeV.

preprint2019arXiv

H.E.S.S. and Fermi-LAT observations of PSR B1259-63/LS 2883 during its 2014 and 2017 periastron passages

PSR B1259-63/LS 2883 is a gamma-ray binary system consisting of a pulsar in an eccentric orbit around a bright Oe stellar-type companion star that features a dense circumstellar disc. The high- and very-high-energy (HE, VHE) gamma-ray emission from PSR B1259-63/LS 2883 around the times of its periastron passage are characterised, in particular, at the time of the HE gamma-ray flares reported to have occurred in 2011, 2014, and 2017. Spectra and light curves were derived from observations conducted with the H.E.S.S.-II array in 2014 and 2017. A local double-peak profile with asymmetric peaks in the VHE light curve is measured, with a flux minimum at the time of periastron $t_p$ and two peaks coinciding with the times at which the neutron star crosses the companion's circumstellar disc ($\sim t_p \pm 16$ d). A high VHE gamma-ray flux is also observed at the times of the HE gamma-ray flares ($\sim t_p + 30$ d) and at phases before the first disc crossing ($\sim t_p - 35$ d). PSR B1259-63/LS 2883 displays periodic flux variability at VHE gamma-rays without clear signatures of super-orbital modulation in the time span covered by H.E.S.S. observations. In contrast, the photon index of the measured power-law spectra remains unchanged within uncertainties for about 200 d around periastron. Lower limits on exponential cut-off energies up to $\sim 40$ TeV are placed. At HE gamma-rays, PSR B1259-63/LS 2883 has now been detected also before and after periastron, close to the disc crossing times. Repetitive flares with distinct variability patterns are detected in this energy range. Such outbursts are not observed at VHEs, although a relatively high emission level is measured. The spectra obtained in both energy regimes displays a similar slope, although a common physical origin either in terms of a related particle population, emission mechanism, or emitter location is ruled out.

preprint2019arXiv

H.E.S.S. detection of very-high-energy gamma-ray emission from the quasar PKS 0736+017

Flat-spectrum radio-quasars (FSRQs) are rarely detected at very-high-energies (VHE; E>100 GeV) due to their low-frequency-peaked SEDs. At present, only 6 FSRQs are known to emit VHE photons, representing only 7% of the VHE extragalactic catalog. Following the detection of MeV-GeV gamma-ray flaring activity from the FSRQ PKS 0736+017 (z=0.189) with Fermi, the H.E.S.S. array of Cherenkov telescopes triggered ToO observations on February 18, 2015, with the goal of studying the gamma-ray emission in the VHE band. H.E.S.S. ToO observations were carried out during the nights of February 18, 19, 21, and 24, 2015. Together with Fermi-LAT, the multi-wavelength coverage of the flare includes Swift observations in soft-X-rays and optical/UV, and optical monitoring (photometry and spectro-polarimetry) by the Steward Observatory, the ATOM, the KAIT and the ASAS-SN telescope. VHE emission from PKS 0736+017 was detected with H.E.S.S. during the night of February 19, 2015, only. Fermi data indicate the presence of a gamma-ray flare, peaking at the time of the H.E.S.S. detection, with a flux doubling time-scale of around six hours. The gamma-ray flare was accompanied by at least a 1 mag brightening of the non-thermal optical continuum. No simultaneous observations at longer wavelengths are available for the night of the H.E.S.S. detection. The gamma-ray observations with H.E.S.S. and Fermi are used to put constraints on the location of the gamma-ray emitting region during the flare: it is constrained to be just outside the radius of the broad-line-region with a bulk Lorentz factor $\simeq 20$, or at the level of the radius of the dusty torus with Gamma > 60. PKS 0736+017 is the seventh FSRQ known to emit VHE photons and, at z=0.189, is the nearest so far. The location of the gamma-ray emitting region during the flare can be tightly constrained thanks to opacity, variability, and collimation arguments.

preprint2016arXiv

Acceleration of petaelectronvolt protons in the Galactic Centre

Galactic cosmic rays reach energies of at least a few Peta-electronvolts (1 PeV =$10^\mathbf{15}$ electron volts). This implies our Galaxy contains PeV accelerators (PeVatrons), but all proposed models of Galactic cosmic-ray accelerators encounter non-trivial difficulties at exactly these energies. Tens of Galactic accelerators capable of accelerating particle to tens of TeV (1 TeV =$10^\mathbf{12}$ electron volts) energies were inferred from recent gamma-ray observations. None of the currently known accelerators, however, not even the handful of shell-type supernova remnants commonly believed to supply most Galactic cosmic rays, have shown the characteristic tracers of PeV particles: power-law spectra of gamma rays extending without a cutoff or a spectral break to tens of TeV. Here we report deep gamma-ray observations with arcminute angular resolution of the Galactic Centre regions, which show the expected tracer of the presence of PeV particles within the central 10~parsec of the Galaxy. We argue that the supermassive black hole Sagittarius A* is linked to this PeVatron. Sagittarius A* went through active phases in the past, as demonstrated by X-ray outbursts and an outflow from the Galactic Centre. Although its current rate of particle acceleration is not sufficient to provide a substantial contribution to Galactic cosmic rays, Sagittarius A* could have plausibly been more active over the last $\gtrsim 10^{6-7}$ years, and therefore should be considered as a viable alternative to supernova remnants as a source of PeV Galactic cosmic rays.

preprint2016arXiv

Deep morphological and spectral study of the SNR RCW 86 with Fermi-LAT

RCW 86 is a young supernova remnant (SNR) showing a shell-type structure at several wavelengths and is thought to be an efficient cosmic-ray (CR) accelerator. Earlier \textit{Fermi} Large Area Telescope results reported the detection of $γ$-ray emission coincident with the position of RCW 86 but its origin (leptonic or hadronic) remained unclear due to the poor statistics. Thanks to 6.5 years of data acquired by the \textit{Fermi}-LAT and the new event reconstruction Pass 8, we report the significant detection of spatially extended emission coming from RCW 86. The spectrum is described by a power-law function with a very hard photon index ($Γ= 1.42 \pm 0.1_{\rm stat} \pm 0.06_{\rm syst}$) in the 0.1--500 GeV range and an energy flux above 100 MeV of ($2.91$ $\pm$ $0.8_{\rm stat}$ $\pm$ $0.12_{\rm syst}$) $\times$ $10^{-11}$ erg cm$^{-2}$ s$^{-1}$. Gathering all the available multiwavelength (MWL) data, we perform a broadband modeling of the nonthermal emission of RCW 86 to constrain parameters of the nearby medium and bring new hints about the origin of the $γ$-ray emission. For the whole SNR, the modeling favors a leptonic scenario in the framework of a two-zone model with an average magnetic field of 10.2 $\pm$ 0.7 $μ$G and a limit on the maximum energy injected into protons of 2 $\times$ 10$^{49}$ erg for a density of 1 cm$^{-3}$. In addition, parameter values are derived for the North-East (NE) and South-West (SW) regions of RCW 86, providing the first indication of a higher magnetic field in the SW region.

preprint2016arXiv

Detailed spectral and morphological analysis of the shell type SNR RCW 86

Aims: We aim for an understanding of the morphological and spectral properties of the supernova remnant RCW~86 and for insights into the production mechanism leading to the RCW~86 very high-energy gamma-ray emission. Methods: We analyzed High Energy Spectroscopic System data that had increased sensitivity compared to the observations presented in the RCW~86 H.E.S.S. discovery publication. Studies of the morphological correlation between the 0.5-1~keV X-ray band, the 2-5~keV X-ray band, radio, and gamma-ray emissions have been performed as well as broadband modeling of the spectral energy distribution with two different emission models. Results:We present the first conclusive evidence that the TeV gamma-ray emission region is shell-like based on our morphological studies. The comparison with 2-5~keV X-ray data reveals a correlation with the 0.4-50~TeV gamma-ray emission.The spectrum of RCW~86 is best described by a power law with an exponential cutoff at $E_{cut}=(3.5\pm 1.2_{stat})$ TeV and a spectral index of $Γ$~$1.6\pm 0.2$. A static leptonic one-zone model adequately describes the measured spectral energy distribution of RCW~86, with the resultant total kinetic energy of the electrons above 1 GeV being equivalent to $\sim$0.1\% of the initial kinetic energy of a Type I a supernova explosion. When using a hadronic model, a magnetic field of $B$~100$μ$G is needed to represent the measured data. Although this is comparable to formerly published estimates, a standard E$^{-2}$ spectrum for the proton distribution cannot describe the gamma-ray data. Instead, a spectral index of $Γ_p$~1.7 would be required, which implies that ~$7\times 10^{49}/n_{cm^{-3}}$erg has been transferred into high-energy protons with the effective density $n_{cm^{-3}}=n/ 1$ cm^-3. This is about 10\% of the kinetic energy of a typical Type Ia supernova under the assumption of a density of 1~cm^-3.

preprint2016arXiv

H.E.S.S. limits on line-like dark matter signatures in the 100 GeV to 2 TeV energy range close to the Galactic Centre

A search for dark matter line-like signals was performed in the vicinity of the Galactic Centre by the H.E.S.S. experiment on observational data taken in 2014. An unbinned likelihood analysis was developed to improve the sensitivity to line-like signals. The upgraded analysis along with newer data extend the energy coverage of the previous measurement down to 100 GeV. The 18 h of data collected with the H.E.S.S. array allow one to rule out at 95% CL the presence of a 130 GeV line (at $l = -1.5^{\circ}, b = 0^{\circ}$ and for a dark matter profile centered at this location) previously reported in Fermi-LAT data. This new analysis overlaps significantly in energy with previous Fermi-LAT and H.E.S.S. results. No significant excess associated with dark matter annihilations was found in the energy range 100 GeV to 2 TeV and upper limits on the gamma-ray flux and the velocity weighted annihilation cross-section are derived adopting an Einasto dark matter halo profile. Expected limits for present and future large statistics H.E.S.S. observations are also given.

preprint2016arXiv

Interstellar gas towards the TeV gamma-ray sources HESS J1640-465 and HESS J1641-463

We present a detailed analysis of the interstellar medium towards the TeV $γ$-ray sources HESS J1640$-$465 and HESS J1641$-$463 using results from the Mopra Southern Galactic Plane CO Survey and from a Mopra 7 mm-wavelength study. The $γ$-ray sources are positionally coincident with two supernova remnants G338.3$-$0.0 and G338.5+0.1 respectively. A bright complex of HII regions connect the two SNRs and TeV objects. Observations in the CO(1-0) transition lines reveal substantial amounts of diffuse gas positionally coincident with the $γ$-ray sources at multiple velocities along the line of sight, while 7 mm observations in CS, SiO, HC$_{3}$N and CH$_{3}$OH transition lines reveal regions of dense, shocked gas. Archival HI data from the Southern Galactic Plane Survey was used to account for the diffuse atomic gas. Physical parameters of the gas towards the TeV sources were calculated from the data. We find that for a hadronic origin for the $γ$-ray emission, the cosmic-ray enhancement rates are $\sim 10^{3}$ and $10^{2}$ times the local solar value for HESS J1640$-$465 and HESS J1641$-$463 respectively.

preprint2016arXiv

Search for dark matter annihilations towards the inner Galactic halo from 10 years of observations with H.E.S.S

The inner region of the Milky Way halo harbors a large amount of dark matter (DM). Given its proximity, it is one of the most promising targets to look for DM. We report on a search for the annihilations of DM particles using $γ$-ray observations towards the inner 300 parsecs of the Milky Way, with the H.E.S.S. array of ground-based Cherenkov telescopes. The analysis is based on a 2D maximum likelihood method using Galactic center (GC) data accumulated by H.E.S.S. over the last 10 years (2004-2014), and does not show any significant $γ$-ray signal above background. Assuming Einasto and Navarro-Frenk-White DM density profiles at the GC, we derive upper limits on the annihilation cross section $\langle σv\rangle$. These constraints are the strongest obtained so far in the TeV DM mass range and improve upon previous limits by a factor 5. For the Einasto profile, the constraints reach $\langle σv\rangle$ values of $\rm 6\times10^{-26} cm^3s^{-1}$ in the $W^+W^-$ channel for a DM particle mass of 1.5 TeV, and $\rm 2\times10^{-26} cm^3s^{-1}$ in the $τ^+τ^-$ channel for 1 TeV mass. For the first time, ground-based $γ$-ray observations have reached sufficient sensitivity to probe $\langle σv\rangle$ values expected from the thermal relic density for TeV DM particles.

preprint2016arXiv

The GCT camera for the Cherenkov Telescope Array

The Gamma-ray Cherenkov Telescope (GCT) is proposed for the Small-Sized Telescope component of the Cherenkov Telescope Array (CTA). GCT's dual-mirror Schwarzschild-Couder (SC) optical system allows the use of a compact camera with small form-factor photosensors. The GCT camera is ~0.4 m in diameter and has 2048 pixels; each pixel has a ~0.2 degree angular size, resulting in a wide field-of-view. The design of the GCT camera is high performance at low cost, with the camera housing 32 front-end electronics modules providing full waveform information for all of the camera's 2048 pixels. The first GCT camera prototype, CHEC-M, was commissioned during 2015, culminating in the first Cherenkov images recorded by a SC telescope and the first light of a CTA prototype. In this contribution we give a detailed description of the GCT camera and present preliminary results from CHEC-M's commissioning.

preprint2015arXiv

Constraints on an Annihilation Signal from a Core of Constant Dark Matter Density around the Milky Way Center with H.E.S.S

An annihilation signal of dark matter is searched for from the central region of the Milky Way. Data acquired in dedicated ON/OFF observations of the Galactic center region with H.E.S.S. are analyzed for this purpose. No significant signal is found in a total of $\sim 9$ h of ON/OFF observations. Upper limits on the velocity averaged cross section, $<σv >$, for the annihilation of dark matter particles with masses in the range of $\sim 300$ GeV to $\sim 10$ TeV are derived. In contrast to previous constraints derived from observations of the Galactic center region, the constraints that are derived here apply also under the assumption of a central core of constant dark matter density around the center of the Galaxy. Values of $<σv >$ that are larger than $3\cdot 10^{-24}\:\mathrm{cm^3/s}$ are excluded for dark matter particles with masses between $\sim 1$ and $\sim 4$ TeV at 95% CL if the radius of the central dark matter density core does not exceed $500$ pc. This is the strongest constraint that is derived on $<σv>$ for annihilating TeV mass dark matter without the assumption of a centrally cusped dark matter density distribution in the search region.

preprint2015arXiv

Discovery of variable VHE gamma-ray emission from the binary system 1FGL J1018.6-5856

Re-observations with the H.E.S.S. telescope array of the very-high-energy (VHE) source HESS J1018-589 A coincident with the Fermi-LAT $γ$-ray binary 1FGL J1018.6-5856 have resulted in a source detection significance of more than 9$σ$, and the detection of variability ($χ^2$/$ν$ of 238.3/155) in the emitted $γ$-ray flux. This variability confirms the association of HESS J1018-589 A with the high-energy $γ$-ray binary detected by Fermi-LAT, and also confirms the point-like source as a new very-high-energy binary system. The spectrum of HESS J1018-589 A is best fit with a power-law function with photon index $Γ= 2.20 \pm 0.14_{\rm stat} \pm 0.2_{\rm sys}$. Emission is detected up to ~20 TeV. The mean differential flux level is $(2.9 \pm 0.4)\times10^{-13}$ TeV$^{-1}$ cm$^{-2}$ s$^{-1}$ at 1 TeV, equivalent to ~1% of the flux from the Crab Nebula at the same energy. Variability is clearly detected in the night-by-night lightcurve. When folded on the orbital period of 16.58 days, the rebinned lightcurve peaks in phase with the observed X-ray and high-energy phaseograms. The fit of the H.E.S.S. phaseogram to a constant flux provides evidence of periodicity at the level of 3$σ$. The shape of the VHE phaseogram and measured spectrum suggest a low inclination, low eccentricity system with a modest impact from VHE $γ$-ray absorption due to pair production ($τ$ < 1 at 300 GeV).

preprint2015arXiv

H.E.S.S. reveals a lack of TeV emission from the supernova remnant Puppis A

Puppis A is an interesting ~4 kyr-old supernova remnant (SNR) that shows strong evidence of interaction between the forward shock and a molecular cloud. It has been studied in detail from radio frequencies to high-energy (HE, 0.1-100 GeV) gamma-rays. An analysis of the Fermi-LAT data has shown an extended HE gamma-ray emission with a 0.2-100 GeV spectrum exhibiting no significant deviation from a power law, unlike most of the GeV-emitting SNRs known to be interacting with molecular clouds. This makes it a promising target for imaging atmospheric Cherenkov telescopes (IACTs) to probe the gamma-ray emission above 100 GeV. Very-high-energy (VHE, E >= 0.1 TeV) gamma-ray emission from Puppis A is for the first time searched for with the High Energy Stereoscopic System (H.E.S.S.). The analysis of the H.E.S.S. data does not reveal any significant emission towards Puppis A. The derived upper limits on the differential photon flux imply that its broadband gamma-ray spectrum must exhibit a spectral break or cutoff. By combining Fermi-LAT and H.E.S.S. measurements, the 99% confidence level upper limits on such a cutoff are found to be 450 and 280 GeV, assuming a power law with a simple exponential and a sub-exponential cutoff, respectively. It is concluded that none of the standard limitations (age, size, radiative losses) on the particle acceleration mechanism, assumed to be still on-going at present, can explain the lack of VHE signal. The scenario in which particle acceleration has ceased some time ago is considered as an alternative explanation. The HE/VHE spectrum of Puppis A could then exhibit a break of non-radiative origin, (as observed in several other interacting SNRs, albeit at somewhat higher energies) owing to the interaction with dense and neutral material in particular towards the northeastern region.

preprint2015arXiv

Revealing the Large-Scale Structures of Interstellar Gas Associated with the Magellanic SNR N132D

We report preliminary results of large-scale distribution toward the Magellanic supernova remnant N132D using Mopra and Chandra archival datasets. We identified a cavity-like CO structure along the X-ray shell toward the southern half of it. The total mass of associating molecular gas is $\sim10^4 M_\odot$, which is smaller than the previous study by an order of magnitude. Further observations using ALMA, ASTE, and Mopra will reveal the detailed spatial structures and its physical conditions.

preprint2015arXiv

The 2012 flare of PG 1553+113 seen with H.E.S.S. and Fermi-LAT

Very high energy (VHE, $E>$100 GeV) $γ$-ray flaring activity of the high-frequency peaked BL Lac object \pg\ has been detected by the \hess\ telescopes. The flux of the source increased by a factor of 3 during the nights of 2012 April 26 and 27 with respect to the archival measurements with hint of intra-night variability. No counterpart of this event has been detected in the \fla\ data. This pattern is consistent with VHE $γ$ ray flaring being caused by the injection of ultrarelativistic particles, emitting $γ$ rays at the highest energies. The dataset offers a unique opportunity to constrain the redshift of this source at \bestz\ using a novel method based on Bayesian statistics. The indication of intra-night variability is used to introduce a novel method to probe for a possible Lorentz Invariance Violation (LIV), and to set limits on the energy scale at which Quantum Gravity (QG) effects causing LIV may arise. For the subluminal case, the derived limits are $\textrm{E}_{\rm QG,1}>4.10\times 10^{17}$ GeV and $\textrm{E}_{\rm QG,2}>2.10\times 10^{10}$ GeV for linear and quadratic LIV effects, respectively.

preprint2015arXiv

The exceptionally powerful TeV gamma-ray emitters in the Large Magellanic Cloud

The Large Magellanic Cloud, a satellite galaxy of the Milky Way, has been observed with the High Energy Stereoscopic System (H.E.S.S.) above an energy of 100 billion electron volts for a deep exposure of 210 hours. Three sources of different types were detected: the pulsar wind nebula of the most energetic pulsar known N 157B, the radio-loud supernova remnant N 132D and the largest non-thermal X-ray shell - the superbubble 30 Dor C. The unique object SN 1987A is, surprisingly, not detected, which constrains the theoretical framework of particle acceleration in very young supernova remnants. These detections reveal the most energetic tip of a gamma-ray source population in an external galaxy, and provide via 30 Dor C the unambiguous detection of gamma-ray emission from a superbubble.

preprint2014arXiv

Diffuse Galactic gamma-ray emission with H.E.S.S

Diffuse $γ$-ray emission is the most prominent observable signature of celestial cosmic-ray interactions at high energies. While already being investigated at GeV energies over several decades, assessments of diffuse $γ$-ray emission at TeV energies remain sparse. After completion of the systematic survey of the inner Galaxy, the H.E.S.S. experiment is in a prime position to observe large-scale diffuse emission at TeV energies. Data of the H.E.S.S. Galactic Plane Survey are investigated in regions off known $γ$-ray sources. Corresponding $γ$-ray flux measurements were made over an extensive grid of celestial locations. Longitudinal and latitudinal profiles of the observed $γ$-ray fluxes show characteristic excess emission not attributable to known $γ$-ray sources. For the first time large-scale $γ$-ray emission along the Galactic Plane using imaging atmospheric Cherenkov telescopes has been observed. While the background subtraction technique limits the ability to recover modest variation on the scale of the H.E.S.S. field of view or larger, which is characteristic of the inverse Compton scatter-induced Galactic diffuse emission, contributions of neutral pion decay as well as emission from unresolved $γ$-ray sources can be recovered in the observed signal to a large fraction. Calculations show that the minimum $γ$-ray emission from $π^0$-decay represents a significant contribution to the total signal. This detection is interpreted as a mix of diffuse Galactic $γ$-ray emission and unresolved sources.

preprint2014arXiv

Discovery of the hard spectrum VHE gamma-ray source HESS J1641-463

This letter reports the discovery of a remarkably hard spectrum source, HESS J1641-463, by the High Energy Stereoscopic System (H.E.S.S.) in the very-high energy (VHE) domain. HESS J1641-463 remained unnoticed by the usual analysis techniques due to confusion with the bright nearby source HESS J1640-465. It emerged at a significance level of 8.5 standard deviations after restricting the analysis to events with energies above 4 TeV. It shows a moderate flux level of F(E > 1 TeV) = (3.64 +/- 0.44_stat +/- 0.73_sys) x 10^-13 cm^-2s-1, corresponding to 1.8% of the Crab Nebula flux above the same energy, and a hard spectrum with a photon index of Gamma = 2.07 +/- 0.11_stat +/- 0.20_sys. It is a point-like source, although an extension up to Gaussian width of sigma = 3 arcmin cannot be discounted due to uncertainties in the H.E.S.S. PSF. The VHE gamma-ray flux of HESS J1641-463 is found to be constant over the observed period when checking time binnings from year-by-year to the 28 min exposures timescales. HESS J1641-463 is positionally coincident with the radio supernova remnant SNR G338.5+0.1. No X-ray candidate stands out as a clear association, however Chandra and XMM-Newton data reveal some potential weak counterparts. Various VHE gamma-ray production scenarios are discussed. If the emission from HESS J1641-463 is produced by cosmic ray protons colliding with the ambient gas, then their spectrum must extend close to 1 PeV. This object may represent a source population contributing significantly to the galactic cosmic ray flux around the knee.

preprint2014arXiv

H.E.S.S. detection of TeV emission from the interaction region between the supernova remnant G349.7+0.2 and a molecular cloud

G349.7+0.2 is a young Galactic supernova remnant (SNR) located at the distance of 11.5 kpc and observed across the entire electromagnetic spectrum from radio to high energy (HE) Gamma-rays. Radio and infrared observations indicate that the remnant is interacting with a molecular cloud. In this paper, the detection of very high energy (VHE) Gamma-ray emission coincident with this SNR with the High Energy Stereoscopic System (H.E.S.S.) is reported. An integral flux F(E>400GeV)=(6.5 +-1.1stat +-1.3syst) x 10^{-13} ph/cm/s corresponding to 0.7% of that of the Crab Nebula and to a luminosity of 10^34 erg/s above the same energy threshold, and a steep photon index Gamma_VHE = 2.8 +-0.27stat +-0.20syst are measured. The analysis of more than 5 yr of Fermi-LAT data towards this source shows a power-law like spectrum with a best-fit photon index Gamma_HE = 2.2 +-0.04stat +0.13-0.31syst. The combined Gamma-ray spectrum of G349.7+0.2 can be described by either a broken power-law (BPL) or a power-law with exponential (or sub-exponential) cutoff (PLC). In the former case, the photon break energy is found at E_br,gamma = 55 +70-30 GeV, slightly higher than what is usually observed in the HE/VHE Gamma-ray emitting middle-aged SNRs known to be interacting with molecular clouds. In the latter case, the exponential (respectively sub-exponential) cutoff energy is measured at E_cut,gamma = 1.4 +1.6-0.55 (respectively 0.35 +0.75-0.21) TeV. A pion-decay process resulting from the interaction of the accelerated protons and nuclei with the dense surrounding medium is clearly the preferred scenario to explain the Gamma-ray emission. The BPL with a spectral steepening of 0.5-1 and the PLC provide equally good fits to the data. The product of the average gas density and the total energy content of accelerated protons and nuclei amounts to nH Wp ~ 5 x 10^51 erg/cm3.

preprint2014arXiv

HESS J1640-465 - an exceptionally luminous TeV gamma-ray supernova remnant

The results of follow-up observations of the TeV gamma-ray source HESSJ 1640-465 from 2004 to 2011 with the High Energy Stereoscopic System (H.E.S.S.) are reported in this work. The spectrum is well described by an exponential cut-off power law with photon index Gamma=2.11 +/- 0.09_stat +/- 0.10_sys, and a cut-off energy of E_c = (6.0 +2.0 -1.2) TeV. The TeV emission is significantly extended and overlaps with the north-western part of the shell of the SNR G338.3-0.0. The new H.E.S.S. results, a re-analysis of archival XMM-Newton data, and multi-wavelength observations suggest that a significant part of the gamma-ray emission from HESS J1640-465 originates in the SNR shell. In a hadronic scenario, as suggested by the smooth connection of the GeV and TeV spectra, the product of total proton energy and mean target density could be as high as W_p n_H ~ 4 x 10^52 (d/10kpc)^2 erg cm^-3.

preprint2014arXiv

Long-term monitoring of PKS 2155$-$304 with ATOM and H.E.S.S.: investigation of optical/$γ$-ray correlations in different spectral states

In this paper we report on the analysis of all the available optical and very high-energy $γ$-ray ($>$200 GeV) data for the BL Lac object PKS 2155$-$304, collected simultaneously with the ATOM and H.E.S.S. telescopes from 2007 until 2009. This study also includes X-ray (RXTE, Swift) and high-energy $γ$-ray (Fermi-LAT) data. During the period analysed, the source was transitioning from its flaring to quiescent optical states,and was characterized by only moderate flux changes at different wavelengths on the timescales of days and months. A flattening of the optical continuum with an increasing optical flux can be noted in the collected dataset, but only occasionally and only at higher flux levels. We did not find any universal relation between the very high-energy $γ$-ray and optical flux changes on the timescales from days and weeks up to several years. On the other hand, we noted that at higher flux levels the source can follow two distinct tracks in the optical flux-colour diagrams, which seem to be related to distinct $γ$-ray states of the blazar. The obtained results therefore indicate a complex scaling between the optical and $γ$-ray emission of PKS 2155$-$304, with different correlation patterns holding at different epochs, and a $γ$-ray flux depending on the combination of an optical flux and colour rather than a flux alone.

preprint2014arXiv

Search for dark matter annihilation signatures in H.E.S.S. observations of Dwarf Spheroidal Galaxies

Dwarf spheroidal galaxies of the Local Group are close satellites of the Milky Way characterized by a large mass-to-light ratio and are not expected to be the site of non-thermal high-energy gamma-ray emission or intense star formation. Therefore they are amongst the most promising candidates for indirect dark matter searches. During the last years the High Energy Stereoscopic System (H.E.S.S.) of imaging atmospheric Cherenkov telescopes observed five of these dwarf galaxies for more than 140 hours in total, searching for TeV gamma-ray emission from annihilation of dark matter particles. The new results of the deep exposure of the Sagittarius dwarf spheroidal galaxy, the first observations of the Coma Berenices and Fornax dwarves and the re-analysis of two more dwarf spheroidal galaxies already published by the H.E.S.S. Collaboration, Carina and Sculptor, are presented. In the absence of a significant signal new constraints on the annihilation cross-section applicable to Weakly Interacting Massive Particles (WIMPs) are derived by combining the observations of the five dwarf galaxies. The combined exclusion limit depends on the WIMP mass and the best constraint is reached at 1-2 TeV masses with a cross-section upper bound of ~3.9x10-24 cm^3 s-1 at a 95% confidence level.

preprint2014arXiv

Theoretical uncertainties of the type Ia supernova rate

It is thought that type Ia supernovae (SNe Ia) are explosions of carbon-oxygen white dwarfs (CO WDs). Two main evolutionary channels are proposed for the WD to reach the critical density required for a thermonuclear explosion: the single degenerate scenario (SD), in which a CO WD accretes from a non-degenerate companion, and the double degenerate scenario (DD), in which two CO WDs merge. However, it remains difficult to reproduce the observed SN Ia rate with these two scenarios. With a binary population synthesis code we study the main evolutionary channels that lead to SNe Ia and we calculate the SN Ia rates and the associated delay time distributions. We find that the DD channel is the dominant formation channel for the longest delay times. The SD channel with helium-rich donors is the dominant channel at the shortest delay times. Our standard model rate is a factor five lower than the observed rate in galaxy clusters. We investigate the influence of ill-constrained aspects of single- and binary-star evolution and uncertain initial binary distributions on the rate of type Ia SNe. These distributions, as well as uncertainties in both helium star evolution and common envelope evolution, have the greatest influence on our calculated rates. Inefficient common envelope evolution increases the relative number of SD explosions such that for $α_{\rm ce} = 0.2$ they dominate the SN Ia rate. Our highest rate is a factor three less than the galaxy-cluster SN Ia rate, but compatible with the rate determined in a field-galaxy dominated sample. If we assume unlimited accretion onto WDs, to maximize the number of SD explosions, our rate is compatible with the observed galaxy-cluster rate.

preprint2013arXiv

Modeling the interaction of thermonuclear supernova remnants with circumstellar structures: The case of Tycho's supernova remnant

The well-established Type Ia remnant of Tycho's supernova (SN 1572) reveals discrepant ambient medium density estimates based on either the measured dynamics or on the X-ray emission properties. This discrepancy can potentially be solved by assuming that the supernova remnant (SNR) shock initially moved through a stellar wind bubble, but is currently evolving in the uniform interstellar medium with a relatively low density. We investigate this scenario by combining hydrodynamical simulations of the wind-loss phase and the supernova remnant evolution with a coupled X-ray emission model, which includes non-equilibrium ionization. For the explosion models we use the well-known W7 deflagration model and the delayed detonation model that was previously shown to provide good fits to the X-ray emission of Tycho's SNR. Our simulations confirm that a uniform ambient density cannot simultaneously reproduce the dynamical and X-ray emission properties of Tycho. In contrast, models that considered that the remnant was evolving in a dense, but small, wind bubble reproduce reasonably well both the measured X-ray emission spectrum and the expansion parameter of Tycho's SNR. Finally, we discuss possible mass loss scenarios in the context of single- and double-degenerate models which possible could form such a small dense wind bubble.

preprint2013arXiv

Proper Motions of H-alpha filaments in the Supernova Remnant RCW 86

We present a proper motion study of the eastern shock-region of the supernova remnant RCW 86 (MSH 14-63, G315.4-2.3), based on optical observations carried out with VLT/FORS2 in 2007 and 2010. For both the northeastern and southeastern regions, we measure an average proper motion of H-alpha filaments of 0.10 +/- 0.02 arcsec/yr, corresponding to 1200 +/- 200 km/s at 2.5kpc. There is substantial variation in the derived proper motions, indicating shock velocities ranging from just below 700 km/s to above 2200 km/s. The optical proper motion is lower than the previously measured X-ray proper motion of northeastern region. The new measurements are consistent with the previously measured proton temperature of 2.3 +/- 0.3 keV, assuming no cosmic-ray acceleration. However, within the uncertainties, moderately efficient (< 27 per cent) shock acceleration is still possible. The combination of optical proper motion and proton temperature rule out the possibility that RCW 86 has a distance less than 1.5kpc. The similarity of the proper motions in the northeast and southeast is peculiar, given the different densities and X-ray emission properties of the regions. The northeastern region has lower densities and the X-ray emission is synchrotron dominated, suggesting that the shock velocities should be higher than in the southeastern, thermal X-ray dominated, region. A possible solution is that the H-alpha emitting filaments are biased toward denser regions, with lower shock velocities. Alternatively, in the northeast the shock velocity may have decreased rapidly during the past 200yr, and the X-ray synchrotron emission is an afterglow from a period when the shock velocity was higher.

preprint2013arXiv

The Hot and Energetic Universe: The astrophysics of supernova remnants and the interstellar medium

The study of both supernova remnants and the hot and cold phases of the interstellar medium are essential for understanding the final stages of stellar evolution and their feedback on the evolution of galaxies through injection of energy and heavy elements. These studies are also crucial for understanding the physics of supernovae, their cosmological implication, and the origin of galactic cosmic rays. The unique capabilities of Athena+ will allow us to explore a new parameter space. Spatially-resolved high-resolution spectroscopy using Athena+ X-IFU of young remnants will allow to characterize individual parcels of ejected material in the line of sight in terms of kinematics, ionization and composition, providing access to the three dimensional geometry of the explosion. Athena+ will also allow studying shock physics and particle acceleration in supernova remnants, as well as their interaction with their environment. Athena+ X-IFU will also characterize the ionization mechanisms competing in forming the complex structures of the hot interstellar medium, likely to keep the echo of past star formation activity, both in our Galaxy and nearby ones. For the first time the dust and gas of the densest cold medium, like in the Galactic Centre environment, will be studied. Athena+ X-IFU will observe, along with the Mg K and Si K edges, which are the main tracers of the silicates content of the ISM, the Fe K edge with unprecedented sensitivity and energy-resolution. This will allow us to study for the first time the nature of Fe-bearing dust in such regions.

preprint2013arXiv

The loss-limited electron energy in SN 1006: effects of the shock velocity and of the diffusion process

The spectral shape of the synchrotron X-ray emission from SN 1006 reveals the fundamental role played by radiative losses in shaping the high-energy tail of the electron spectrum. We analyze data from the XMM-Newton SN 1006 Large Program and confirm that in both nonthermal limbs the loss-limited model correctly describes the observed spectra. We study the physical origin of the observed variations of the synchrotron cutoff energy across the shell. We investigate the role played by the shock velocity and by the electron gyrofactor. We found that the cutoff energy of the syncrotron X-ray emission reaches its maximum value in regions where the shock has experienced its highest average speed. This result is consistent with the loss-limited framework. We also find that the electron acceleration in both nonthermal limbs of SN 1006 proceeds close to the Bohm diffusion limit, the gyrofactor being in the range 1.5-4. We finally investigate possible explanations for the low values of cutoff energy measured in thermal limbs.

preprint2013arXiv

The shape of the cutoff in the synchrotron emission of SN 1006 observed with XMM-Newton

Synchrotron X-ray emission from the rims of young supernova remnants allows us to study the high-energy tail of the electrons accelerated at the shock front. The analysis of X-ray spectra can provide information on the physical mechanisms that limit the energy achieved by the electrons in the acceleration process. We aim at verifying whether the maximum electron energy in SN 1006 is limited by synchrotron losses and at obtaining information on the shape of the cutoff in the X-ray synchrotron emission. We analyzed the deep observations of the XMM-Newton SN 1006 Large Program. We performed spatially resolved spectral analysis of a set of small regions in the nonthermal limbs and studied the X-ray spectra by adopting models that assume different electron spectra. We found out that a loss-limited model provides the best fit to all the spectra and this indicates that the shape of the cutoff in the electron momentum (p) distribution has the form exp[-(p/p_cut)^2]. We also detected residual thermal emission from shocked ambient medium and confirmed the reliability of previous estimates of the post-shock density. Our results indicate that radiative losses play a fundamental role in shaping the electron spectrum in SN 1006.

preprint2012arXiv

Constraints on cosmic-ray efficiency in the supernova remnant RCW 86 using multi-wavelength observations

Several young supernova remnants (SNRs) have recently been detected in the high-energy and very-high-energy gamma-ray domains. As exemplified by RX J1713.7-3946, the nature of this emission has been hotly debated, and direct evidence for the efficient acceleration of cosmic-ray protons at the SNR shocks still remains elusive. We analyzed more than 40 months of data acquired by the Large Area Telescope (LAT) on-board the Fermi Gamma-Ray Space Telescope in the HE domain, and gathered all of the relevant multi-wavelength (from radio to VHE gamma-rays) information about the broadband nonthermal emission from RCW 86. For this purpose, we re-analyzed the archival X-ray data from the ASCA/Gas Imaging Spectrometer (GIS), the XMM-Newton/EPIC-MOS, and the RXTE/Proportional Counter Array (PCA). Beyond the expected Galactic diffuse background, no significant gamma-ray emission in the direction of RCW 86 is detected in any of the 0.1-1, 1-10 and 10-100 GeV Fermi-LAT maps. In the hadronic scenario, the derived HE upper limits together with the HESS measurements in the VHE domain can only be accommodated by a spectral index Gamma <= 1.8, i.e. a value in-between the standard (test-particle) index and the asymptotic limit of theoretical particle spectra in the case of strongly modified shocks. The interpretation of the gamma-ray emission by inverse Compton scattering of high energy electrons reproduces the multi-wavelength data using a reasonable value for the average magnetic field of 15-25 muG. For these two scenarios, we assessed the level of acceleration efficiency. We discuss these results in the light of existing estimates of the magnetic field strength, the effective density and the acceleration efficiency in RCW 86.

preprint2012arXiv

Observational signatures of particle acceleration in supernova remnants

We evaluate the current status of supernova remnants as the sources of Galactic cosmic rays. We summarize observations of supernova remnants, covering the whole electromagnetic spectrum and describe what these obser- vations tell us about the acceleration processes by high Mach number shock fronts. We discuss the shock modification by cosmic rays, the shape and maximum energy of the cosmic-ray spectrum and the total energy budget of cosmic rays in and surrounding supernova remnants. Additionally, we discuss problems with supernova remnants as main sources of Galactic cosmic rays, as well as alternative sources.

preprint2012arXiv

XMM-Newton evidence of shocked ISM in SN 1006: indications of hadronic acceleration

Shock fronts in young supernova remnants are the best candidates for being sites of cosmic ray acceleration up to a few PeV, though conclusive experimental evidence is still lacking. Hadron acceleration is expected to increase the shock compression ratio, providing higher postshock densities, but X-ray emission from shocked ambient medium has not firmly been detected yet in remnants where particle acceleration is at work. We exploited the deep observations of the XMM-Newton Large Program on SN 1006 to verify this prediction. We performed spatially resolved spectral analysis of a set of regions covering the southeastern rim of SN 1006. We studied the spatial distribution of the thermodynamic properties of the ambient medium and carefully verified the robustness of the result with respect to the analysis method. We detected the contribution of the shocked ambient medium. We also found that the postshock density of the interstellar medium significantly increases in regions where particle acceleration is efficient. Under the assumption of uniform preshock density, we found that the shock compression ratio reaches a value of ~6 in regions near the nonthermal limbs. Our results support the predictions of shock modification theory and indicate that effects of acceleration of cosmic ray hadrons on the postshock plasma can be observed in supernova remnants.

preprint2011arXiv

Effects of Neutral Hydrogen on Cosmic Ray Precursors in Supernova Remnant Shock Waves

Many fast supernova remnant shocks show spectra dominated by Balmer lines. The H$α$ profiles have a narrow component explained by direct excitations and a thermally Doppler broadened component due to atoms that undergo charge exchange in the post-shock region. However, the standard model does not take into account the cosmic-ray shock precursor, which compresses and accelerates plasma ahead of the shock. In strong precursors with sufficiently high densities, the processes of charge exchange, excitation and ionization will affect the widths of both narrow and broad line components. Moreover, the difference in velocity between the neutrals and the precursor plasma gives rise to frictional heating due to charge exchange and ionization in the precursor. In extreme cases, all neutrals can be ionized by the precursor. In this paper we compute the ion and electron heating for a wide range of shock parameters, along with the velocity distribution of the neutrals that reach the shock. Our calculations predict very large narrow component widths for some shocks with efficient acceleration, along with changes in the broad- to-narrow intensity ratio used as a diagnostic for the electron-ion temperature ratio. Balmer lines may therefore provide a unique diagnostic of precursor properties. We show that heating by neutrals in the precursor can account for the observed H$α$ narrow component widths, and that the acceleration efficiency is modest in most Balmer line shocks observed thus far.

preprint2011arXiv

Modeling supernova remnants: effects of diffusive cosmic-ray acceleration on the evolution, application to observations

We present numerical models for supernova remnant evolution, using a new version of the hydrodynamical code SUPREMNA. We added cosmic ray diffusion equation to the code scheme, employing two-fluid approximation. We investigate the dynamics of the simulated supernova remnants with different values of cosmic ray acceleration efficiency and diffusion coefficient. We compare the numerical models with observational data of Tycho's and SN1006 supernova remnants. We find models which reproduce the observed locations of the blast wave, contact discontinuity, and reverse shock for the both remnants, thus allowing us to estimate the contribution of cosmic ray particles into total pressure and cosmic-ray energy losses in these supernova remnants. We derive that the energy losses due to cosmic rays escape in Tycho's supernova remnant are 10-20% of the kinetic energy flux and 20-50% in SN1006.

preprint2011arXiv

Temperature Equilibration Behind the Shock Front: an Optical and X-ray Study of RCW 86

We study the electron-proton temperature equilibration behind several shocks of the RCW 86 supernova remnant. To measure the proton temperature, we use published and new optical spectra, all from different locations on the remnant. For each location, we determine the electron temperature from X-ray spectra, and correct for temperature equilibration between the shock front and the location of the X-ray spectrum. We confirm the result of previous studies that the electron and proton temperatures behind shock fronts are consistent with equilibration for slow shocks and deviate for faster shocks. However, we can not confirm the previously reported trend of the electron temperature to proton temperature ratio of 1/v^2.

preprint2011arXiv

The high resolution X-ray spectrum of SNR 0506-68 using XMM-Newton

Aims: We study the supernova remnant 0506-68 in order to obtain detailed information about, among other things, the ionisation state and age of the ionised plasma. Methods: Using the Reflection Grating Spectrometer (RGS) onboard the XMM-Newton satellite we are able to take detailed spectra of the remnant. In addition, we use the MOS data to obtain spectral information at higher energies. Results: The spectrum shows signs of recombination and we derive the conditions for which the remnant and SNR in general are able to cool rapidly enough to become over-ionised. The elemental abundances found are mostly in agreement with the mean LMC abundances. Our models and calculations favour the lower age estimate mentioned in the literature of $\sim4000$ year.

preprint2010arXiv

First detection of VHE gamma-rays from SN 1006 by H.E.S.S

Recent theoretical predictions of the lowest very high energy (VHE) luminosity of SN 1006 are only a factor 5 below the previously published H.E.S.S. upper limit, thus motivating further in-depth observations of this source. Deep observations at VHE energies (above 100 GeV) were carried out with the High Energy Stereoscopic System (H.E.S.S.) of Cherenkov Telescopes from 2003 to 2008. More than 100 hours of data have been collected and subjected to an improved analysis procedure. Observations resulted in the detection of VHE gamma-rays from SN 1006. The measured gamma-ray spectrum is compatible with a power-law, the flux is of the order of 1% of that detected from the Crab Nebula, and is thus consistent with the previously established H.E.S.S. upper limit. The source exhibits a bipolar morphology, which is strongly correlated with non-thermal X-rays. Because the thickness of the VHE-shell is compatible with emission from a thin rim, particle acceleration in shock waves is likely to be the origin of the gamma-ray signal. The measured flux level can be accounted for by inverse Compton emission, but a mixed scenario that includes leptonic and hadronic components and takes into account the ambient matter density inferred from observations also leads to a satisfactory description of the multi-wavelength spectrum.

preprint2009arXiv

Spectral and temporal variations of the isolated neutron star RX J0720.4-3125: new XMM-Newton observations

In the past, the isolated, radio-quiet neutron star RX J0720.4-3125 showed variations in the spectral parameters (apparent radius, temperature of the emitting area and equivalent width of the absorption feature) seen in the X-ray spectra, not only during the spin period of 8.39s, but also over time scales of years. New X-ray observations of RX J0720.4-3125 with XMM Newton extend the coverage to about 7.5 years with the latest pointing performed in November 2007. Out of a total of fourteen available EPIC-pn datasets, eleven have been obtained with an identical instrumental setup (full frame read-out mode with thin filter), and are best suited for a comparative investigations of the spectral and timing properties of this enigmatic X-ray pulsar. We analysed the new XMM Newton observations together with archival data in order to follow the spectral and temporal evolution of RX J0720.4-3125 All XMM-Newton data were reduced with the standard XMM-SAS software package. A systematic and consistent data reduction of all these observations was emphasised in order to reduce systematic errors as far as possible. We investigate the phase residuals derived from data from different energy bands using different timing solutions for the spin period evolution and confirm the phase lag between hard and soft photons. The phase shift in the X-ray pulses between hard and soft photons varies with time and changes sign around MJD=52800 days, regardless of the chosen timing solution. The phase residuals[abridge]

preprint2009arXiv

Suzaku X-Ray Imaging and Spectroscopy of Cassiopeia A

Suzaku X-ray observations of a young supernova remnant, Cassiopeia A, were carried out. K-shell transition lines from highly ionized ions of various elements were detected, including Chromium (Cr-Kalpha at 5.61 keV). The X-ray continuum spectra were modeled in the 3.4--40 keV band, summed over the entire remnant, and were fitted with a simplest combination of the thermal bremsstrahlung and the non-thermal cut-off power-law models. The spectral fits with this assumption indicate that the continuum emission is likely to be dominated by the non-thermal emission with a cut-off energy at > 1 keV. The thermal-to-nonthermal fraction of the continuum flux in the 4-10 keV band is best estimated as ~0.1. Non-thermal-dominated continuum images in the 4--14 keV band were made. The peak of the non-thermal X-rays appears at the western part. The peak position of the TeV gamma-rays measured with HEGRA and MAGIC is also shifted at the western part with the 1-sigma confidence. Since the location of the X-ray continuum emission was known to be presumably identified with the reverse shock region, the possible keV-TeV correlations give a hint that the accelerated multi-TeV hadrons in Cassiopeia A are dominated by heavy elements in the reverse shock region.

preprint2009arXiv

The hydrodynamics of the supernova remnant Cas A: The influence of the progenitor evolution on the velocity structure and clumping

There are large differences in the proposed progenitor models for the Cas A SNR. One of these differences is the presence or absence of a Wolf-Rayet (WR) phase of the progenitor star. The mass loss history of the progenitor star strongly affects the shape of the Supernova remnant (SNR). In this paper we investigate whether the progenitor star of Cas A had a WR phase or not and how long it may have lasted. We performed two-dimensional multi-species hydrodynamical simulations of the CSM around the progenitor star for several WR life times, each followed by the interaction of the supernova ejecta with the CSM. We then looked at the influence of the length of the WR phase and compared the results of the simulations with the observations of Cas A. The difference in the structure of the CSM, for models with different WR life times, has a strong impact on the resulting SNR. With an increasing WR life time the reverse shock velocity of the SNR decreases and the range of observed velocities in the shocked material increases. Furthermore, if a WR phase occurs, the remainders of the WR shell will be visible in the resulting SNR. Comparing our results with the observations suggests that the progenitor star of Cas A did not have a WR phase. We also find that the quasi-stationary flocculi (QSF) in Cas A are not consistent with the clumps from a WR shell that have been shocked and accelerated by the interaction with the SN ejecta. We can also conclude that for a SN explosion taking place in a CSM that is shaped by the wind during a short < 15000 yr WR phase, the clumps from the WR shell will be visible inside the SNR.