Showing papers by "David Maurin published in 2010"

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

Abstract: 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.

First detection of VHE gamma-rays from SN 1006 by HESS

Abstract: 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.

A Markov Chain Monte Carlo technique to sample transport and source parameters of Galactic cosmic rays - II. Results for the diffusion model combining B/C and radioactive nuclei

Abstract: Context. Ongoing measurements of the cosmic radiation (nuclear, electronic, and γ-ray) are providing additional insight into cosmicray physics. A comprehensive picture of these data relies on an accurate determination of the transport and source parameters of propagation models. Aims. A Markov Chain Monte Carlo method is used to obtain these parameters in a diffusion model. By measuring the B/C ratio and radioactive cosmic-ray clocks, we calculate their probability density functions, placing special emphasis on the halo size L of the Galaxy and the local underdense bubble of size rh. We also derive the mean, best-fit model parameters and 68% confidence level for the various parameters, and the envelopes of other quantities. Methods. The analysis relies on the USINE code for propagation and on a Markov Chain Monte Carlo technique previously developed by ourselves for the parameter determination. Results. The B/C analysis leads to a most probable diffusion slope δ = 0.86 +0.04

VHE gamma-ray emission of PKS 2155-304: spectral and temporal variability

Abstract: Observations of very high energy gamma-rays from blazars provide information about acceleration mechanisms occurring in their innermost regions. Studies of variability in these objects allow a better understanding of the mechanisms at play. To investigate the spectral and temporal variability of VHE (>100 GeV) gamma-rays of the well-known high-frequency-peaked BL Lac object PKS 2155-304 with the H.E.S.S. imaging atmospheric Cherenkov telescopes over a wide range of flux states. Data collected from 2005 to 2007 are analyzed. Spectra are derived on time scales ranging from 3 years to 4 minutes. Light curve variability is studied through doubling timescales and structure functions, and is compared with red noise process simulations. The source is found to be in a low state from 2005 to 2007, except for a set of exceptional flares which occurred in July 2006. The quiescent state of the source is characterized by an associated mean flux level of 4.32 +/-0.09 x 10^-11 cm^-2 s^-1 above 200 GeV, or approximately 15% of the Crab Nebula, and a power law photon index of 3.53 +/-0.06. During the flares of July 2006, doubling timescales of ~2 min are found. The spectral index variation is examined over two orders of magnitude in flux, yielding different behaviour at low and high fluxes,which is a new phenomenon in VHE gamma-ray emitting blazars. The variability amplitude characterized by the fractional r.m.s. is strongly energy-dependent and is proportional to E^(0.19 +/- 0.01). The light curve r.m.s. correlates with the flux. This is the signature of a multiplicative process which can be accounted for as a red noise with a Fourier index of ~2. This unique data set shows evidence for a low level gamma-ray emission state from PKS 2155-304, which possibly has a different origin than the outbursts. The discovery of the light curve lognormal behaviour might be an indicator ..

VHE gamma-ray emission of PKS 2155-304: spectral and temporal variability

Abstract: Observations of very high energy gamma-rays from blazars provide information about acceleration mechanisms occurring in their innermost regions. Studies of variability in these objects allow a better understanding of the mechanisms at play. To investigate the spectral and temporal variability of VHE (>100 GeV) gamma-rays of the well-known high-frequency-peaked BL Lac object PKS 2155-304 with the H.E.S.S. imaging atmospheric Cherenkov telescopes over a wide range of flux states. Data collected from 2005 to 2007 are analyzed. Spectra are derived on time scales ranging from 3 years to 4 minutes. Light curve variability is studied through doubling timescales and structure functions, and is compared with red noise process simulations. The source is found to be in a low state from 2005 to 2007, except for a set of exceptional flares which occurred in July 2006. The quiescent state of the source is characterized by an associated mean flux level of 4.32 +/-0.09 x 10^-11 cm^-2 s^-1 above 200 GeV, or approximately 15% of the Crab Nebula, and a power law photon index of 3.53 +/-0.06. During the flares of July 2006, doubling timescales of ~2 min are found. The spectral index variation is examined over two orders of magnitude in flux, yielding different behaviour at low and high fluxes,which is a new phenomenon in VHE gamma-ray emitting blazars. The variability amplitude characterized by the fractional r.m.s. is strongly energy-dependent and is proportional to E^(0.19 +/- 0.01). The light curve r.m.s. correlates with the flux. This is the signature of a multiplicative process which can be accounted for as a red noise with a Fourier index of ~2. This unique data set shows evidence for a low level gamma-ray emission state from PKS 2155-304, which possibly has a different origin than the outbursts. The discovery of the light curve lognormal behaviour might be an indicator ..

Systematic uncertainties on the cosmic-ray transport parameters - Is it possible to reconcile B/C data with δ = 1/3 or δ = 1/2?

Abstract: Context. The B/C ratio is used in cosmic-ray physics to constrain the transport parameters. However, from the same set of data, the various published values show a puzzling large scatter of these parameters. Aims. We investigate the impact of using different inputs (gas density and hydrogen fraction in the Galactic disc, source spectral shape, low-energy dependence of the diffusion coefficient, and nuclear fragmentation cross-sections) on the best-fit values of the transport parameters. We quantify the systematics produced when varying these inputs, and compare them to statistical uncertainties. We discuss the consequences for the slope of the diffusion coefficient δ. Methods. The analysis relies on the propagation code USINE interfaced with the Minuit minimisation routines. Results. We find the typical systematic uncertainties to be greater than the statistical ones. The several published values of δ (from 0.3 to 0.8) can be recovered when varying the low-energy shape of the diffusion coefficient and the convective wind strength. Models including a convective wind are characterised by δ 0.6, which cannot be reconciled with the expected theoretical values (1/ 3a nd 1/2). However, from a statistical point of view (χ 2 analysis), models with both reacceleration and convection – hence large δ –a re favoured. The next favoured models in line yield δ, which can be accommodated with 1/ 3a nd 1/2, but require a strong upturn of the diffusion coefficient at low energy (and no convection). Conclusions. To date, using the best statistical tools, the transport parameter determination is still plagued by many unknowns at low energy (∼GeV/n). To disentangle all these configurations, measurements of the B/C ratio at TeV/n energies and/or combination with other secondary-to-primary ratios is necessary.

Systematic uncertainties on the cosmic-ray transport parameters: Is it possible to reconcile B/C data with delta = 1/3 or delta = 1/2?

Abstract: The B/C ratio is used in cosmic-ray physics to constrain the transport parameters. However, from the same set of data, the various published values show a puzzling large scatter of these parameters. We investigate the impact of using different inputs (gas density and hydrogen fraction in the Galactic disc, source spectral shape, low-energy dependence of the diffusion coefficient, and nuclear fragmentation cross-sections) on the best-fit values of the transport parameters. We quantify the systematics produced when varying these inputs, and compare them to statistical uncertainties. We discuss the consequences for the slope of the diffusion coefficient delta. The analysis relies on the propagation code USINE interfaced with the Minuit minimisation routines. We find the typical systematic uncertainties to be larger than the statistical ones. The several published values of delta (from 0.3 to 0.8) can be recovered when varying the low-energy shape of the diffusion coefficient and the convective wind strength. Models including a convective wind are characterised by delta > 0.6, which cannot be reconcile with the expected theoretical values (1/3 and 1/2). However, from a statistical point of view (chi^2 analysis), models with both reacceleration and convection-hence large delta-are favoured. The next favoured models in line yield delta that can be accommodated with 1/3 and 1/2, but require a strong upturn of the diffusion coefficient at low energy (and no convection). To date, using the best statistical tools, the transport parameter determination is still plagued by many unknowns at low energy (~ GeV/n). To disentangle between all these configurations, measurements of the B/C ratio at TeV/n energies and/or combination with other secondary-to-primary ratios is necessary.

Localizing the VHE γ‐ray source at the Galactic Centre

Abstract: The inner 10 pc of our Galaxy contains many counterpart candidates of the very high energy (VHE; >100 GeV) γ-ray point source HESS J1745−290. Within the point spread function of the H.E.S.S. measurement, at least three objects are capable of accelerating particles to VHE and beyond and of providing the observed γ-ray flux. Previous attempts to address this source confusion were hampered by the fact that the projected distances between these objects were of the order of the error circle radius of the emission centroid (34 arcsec, dominated by the pointing uncertainty of the H.E.S.S. instrument). Here we present H.E.S.S. data of the Galactic Centre region, recorded with an improved control of the instrument pointing compared to H.E.S.S. standard pointing procedures. Stars observed during γ-ray observations by optical guiding cameras mounted on each H.E.S.S. telescope are used for off-line pointing calibration, thereby decreasing the systematic pointing uncertainties from 20 to 6 arcsec per axis. The position of HESS J1745−290 is obtained by fitting a multi-Gaussian profile to the background-subtracted γ-ray count map. A spatial comparison of the best-fitting position of HESS J1745−290 with the position and morphology of candidate counterparts is performed. The position is, within a total error circle radius of 13 arcsec, coincident with the position of the supermassive black hole Sgr A* and the recently discovered pulsar wind nebula candidate G359.95−0.04. It is significantly displaced from the centroid of the supernova remnant Sgr A East, excluding this object with high probability as the dominant source of the VHE γ-ray emission.

Revisiting the Westerlund 2 Field with the H.E.S.S. Telescope Array

Abstract: Aims. Previous observations with the H.E.S.S. telescope array revealed the existence of extended very-high-energy (VHE; E>100 GeV) {\gamma}-ray emission, HESS J1023-575, coincident with the young stellar cluster Westerlund 2. At the time of discovery, the origin of the observed emission was not unambiguously identified, and follow-up observations have been performed to further investigate the nature of this {\gamma}-ray source. Methods. The Carina region towards the open cluster Westerlund 2 has been re-observed, increasing the total exposure to 45.9 h. The combined dataset includes 33 h of new data and now permits a search for energy-dependent morphology and detailed spectroscopy. Results. A new, hard spectrum VHE {\gamma}-ray source, HESSJ1026-582, was discovered with a statistical significance of 7{\sigma}. It is positionally coincident with the Fermi LAT pulsar PSR J1028-5819. The positional coincidence and radio/{\gamma}-ray characteristics of the LAT pulsar favors a scenario where the TeV emission originates from a pulsar wind nebula. The nature of HESS J1023-575 is discussed in light of the deep H.E.S.S. observations and recent multi-wavelength discoveries, including the Fermi LAT pulsar PSRJ1022-5746 and giant molecular clouds in the region. Despite the improved VHE dataset, a clear identification of the object responsible for the VHE emission from HESS J1023-575 is not yet possible, and contribution from the nearby high-energy pulsar and/or the open cluster remains a possibility.

Multi-wavelength Observations of H 2356-309

Abstract: AIMS: The properties of the broad-band emission from the high-frequency peaked BL Lac H 2356-309 (z=0.165) are investigated. METHODS: Very High Energy (VHE; E > 100 GeV) observations of H 2356-309 were performed with the High Energy Stereoscopic System (HESS) from 2004 through 2007. Simultaneous optical/UV and X-ray observations were made with the XMM-Newton satellite on June 12/13 and June 14/15, 2005. NRT radio observations were also contemporaneously performed in 2005. ATOM optical monitoring observations were also made in 2007. RESULTS: A strong VHE signal, ~13 sigma total, was detected by HESS after the four years HESS observations (116.8 hrs live time). The integral flux above 240 GeV is I(>240 GeV) = (3.06 +- 0.26 {stat} +- 0.61 {syst}) x 10^{-12} cm^{-2} s^{-1}, corresponding to ~1.6% of the flux observed from the Crab Nebula. A time-averaged energy spectrum is measured from 200 GeV to 2 TeV and is characterized by a power law (photon index of Gamma = 3.06 +- 0.15 {stat} +- 0.10 {syst}). Significant small-amplitude variations in the VHE flux from H 2356-309 are seen on time scales of months and years, but not on shorter time scales. No evidence for any variations in the VHE spectral slope are found within these data. The XMM-Newton X-ray measurements show a historically low X-ray state, characterized by a hard, broken-power-law spectrum on both nights. CONCLUSIONS: The broad-band spectral energy distribution (SED) of the blazar can be adequately fit using a simple one-zone synchrotron self-Compton (SSC) model. In the SSC scenario, higher VHE fluxes could be expected in the future since the observed X-ray flux is at a historically low level.

Multi-wavelength Observations of H 2356-309

Abstract: AIMS: The properties of the broad-band emission from the high-frequency peaked BL Lac H 2356-309 (z=0.165) are investigated. METHODS: Very High Energy (VHE; E > 100 GeV) observations of H 2356-309 were performed with the High Energy Stereoscopic System (HESS) from 2004 through 2007. Simultaneous optical/UV and X-ray observations were made with the XMM-Newton satellite on June 12/13 and June 14/15, 2005. NRT radio observations were also contemporaneously performed in 2005. ATOM optical monitoring observations were also made in 2007. RESULTS: A strong VHE signal, ~13 sigma total, was detected by HESS after the four years HESS observations (116.8 hrs live time). The integral flux above 240 GeV is I(>240 GeV) = (3.06 +- 0.26 {stat} +- 0.61 {syst}) x 10^{-12} cm^{-2} s^{-1}, corresponding to ~1.6% of the flux observed from the Crab Nebula. A time-averaged energy spectrum is measured from 200 GeV to 2 TeV and is characterized by a power law (photon index of Gamma = 3.06 +- 0.15 {stat} +- 0.10 {syst}). Significant small-amplitude variations in the VHE flux from H 2356-309 are seen on time scales of months and years, but not on shorter time scales. No evidence for any variations in the VHE spectral slope are found within these data. The XMM-Newton X-ray measurements show a historically low X-ray state, characterized by a hard, broken-power-law spectrum on both nights. CONCLUSIONS: The broad-band spectral energy distribution (SED) of the blazar can be adequately fit using a simple one-zone synchrotron self-Compton (SSC) model. In the SSC scenario, higher VHE fluxes could be expected in the future since the observed X-ray flux is at a historically low level.

PKS2005-489 at VHE: four years of monitoring with HESS and simultaneous multi-wavelength observations

Abstract: VHE observations of PKS 2005-489 were made with HESS from 2004 through 2007, together with three simultaneous multi-wavelength campaigns performed with XMM-Newton and RXTE in 2004 and 2005. A strong VHE signal, ~17sigma total, is detected during the four years of HESS observations (90.3 hrs live time). The integral flux above the average analysis threshold of 400 GeV is ~3% of the Crab and varies weakly on time-scales from days to years. At X-ray energies the flux is observed to vary by more than an order of magnitude between 2004 and 2005. Strong changes in the X-ray spectrum (DeltaGamma~0.7) are also observed, which appear to be mirrored in the VHE band. The SED of PKS 2005-489, constructed for the first time with contemporaneous data on both humps, shows significant evolution. The large flux variations in the X-ray band, coupled with weak or no variations in the VHE band and a similar spectral behavior, suggest the emergence of a new, separate, harder emission component in September 2005.

Discovery of VHE gamma-rays from the BL Lac object PKS 0548-322

Abstract: PKS 0548-322 (z=0.069) is a ``high-frequency-peaked'' BL Lac object and a candidate very high energy (VHE, E>100 GeV) gamma-ray emitter, due to its high X-ray and radio flux. Observations at the VHE band provide insights into the origin of very energetic particles present in this source and the radiation processes at work. We report observations made between October 2004 and January 2008 with the H.E.S.S. array, a four imaging atmospheric-Cherenkov telescopes. Contemporaneous UV and X-ray observations with the Swift satellite in November 2006 are also reported. PKS 0548-322 is detected for the first time in the VHE band with H.E.S.S. We measure an excess of 216 gamma-rays corresponding to a significance of 5.6 standard deviations. The photon spectrum of the source is described by a power-law, with a photon index of Gamma=2.86 +/- 0.34 (stat) +/- 0.10 (sys). The integral flux above 200 GeV is 1.3 % of the flux of the Crab Nebula, and is consistent with being constant in time. Contemporaneous Swift/XRT observations reveal an X-ray flux between 2 and 10 keV of F_{2-10 keV}=2.3 +/- 0.2 x 10^{-11} erg.cm^{-2}. s^{-1}, an intermediate intensity state with respect to previous observations. The spectral energy distribution can be reproduced using a simple one-zone synchrotron self Compton model, with parameters similar those observed for other sources of this type.

Discovery of VHE γ-rays from the BL Lacertae object PKS 0548-322

Abstract: PKS 0548-322 (z=0.069) is a "high-frequency-peaked'' BL Lac object and a candidate very high energy (VHE, E>100 GeV) gamma-ray emitter, due to its high X-ray and radio flux. Observations at the VHE band provide insights into the origin of very energetic particles present in this source and the radiation processes at work. We report observations made between October 2004 and January 2008 with the H.E.S.S. array, a four imaging atmospheric-Cherenkov telescopes. Contemporaneous UV and X-ray observations with the Swift satellite in November 2006 are also reported. PKS 0548-322 is detected for the first time in the VHE band with H.E.S.S. We measure an excess of 216 gamma-rays corresponding to a significance of 5.6 standard deviations. The photon spectrum of the source is described by a power-law, with a photon index of Gamma=2.86 +/- 0.34 (stat) +/- 0.10 (sys). The integral flux above 200 GeV is 1.3 % of the flux of the Crab Nebula, and is consistent with being constant in time. Contemporaneous Swift/XRT observations reveal an X-ray flux between 2 and 10 keV of F_{2-10 keV}=2.3 +/- 0.2 x 10^{-11} erg.cm^{-2}. s^{-1}, an intermediate intensity state with respect to previous observations. The spectral energy distribution can be reproduced using a simple one-zone synchrotron self Compton model, with parameters similar those observed for other sources of this type.

Discovery and follow-up studies of the extended, off-plane, VHE gamma-ray source HESS J1507-622

Abstract: Context: The detection of gamma-rays in the very-high-energy (VHE) range (100 GeV-100 TeV) offers the possibility of studying the parent population of ultrarelativistic particles found in astrophysical sources, so it is useful for understanding the underlying astrophysical processes in nonthermal sources. Aim: The discovery of the VHE gamma-ray source HESS J1507-622 is reported and possibilities regarding its nature are investigated. Methods: The H.E.S.S. array of imaging atmospheric Cherenkov telescopes (IACTs) has a high sensitivity compared with previous instruments (~1% of the Crab flux in 25 hours observation time for a 5 sigma point-source detection) and has a large field of view (~5 deg in diameter). HESS J1507-622 was discovered within the ongoing H.E.S.S. survey of the inner Galaxy, and the source was also studied by means of dedicated multiwavelength observations. Results: A Galactic gamma-ray source, HESS J1507-622, located ~3.5 deg from the Galactic plane was detected with a statistical significance > 9 sigma. Its energy spectrum is well fitted by a power law with spectral index \Gamma = 2.24 +/- 0.16_{stat} +/- 0.20_{sys} and a flux above 1 TeV of (1.5 +/- 0.4_{stat} +/- 0.3_{sys}) X 10^{-12} cm^{-2} s^{-1}. Possible interpretations (considering both hadronic and leptonic models) of the VHE gamma-ray emission are discussed in the absence of an obvious counterpart.

A Markov Chain Monte Carlo technique to sample transport and source parameters of Galactic cosmic rays: II. Results for the diffusion model combining B/C and radioactive nuclei

Abstract: On-going measurements of the cosmic radiation (nuclear, electronic, and gamma-ray) are shedding new light on cosmic-ray physics. A comprehensive picture of these data relies on an accurate determination of the transport and source parameters of propagation models. A Markov Chain Monte Carlo is used to obtain these parameters in a diffusion model. From the measurement of the B/C ratio and radioactive cosmic-ray clocks, we calculate their probability density functions, with a special emphasis on the halo size L of the Galaxy and the local underdense bubble of size r_h. The analysis relies on the USINE code for propagation and on a Markov Chain Monte Carlo technique (Putze et al. 2009, paper I of this series) for the parameter determination. As found in previous studies, the B/C best-fit model favours diffusion/convection/reacceleration (Model III) over diffusion/reacceleration (Model II). A combined fit on B/C and the isotopic ratios (10Be/9Be, 26Al/27Al, 36Cl/Cl) leads to L ~ 8 kpc and r_h ~ 120 pc for the best-fit Model III. This value for r_h is consistent with direct measurements of the local interstallar medium. For Model II, L ~ 4 kpc and r_h is consistent with zero. We showed the potential and usefulness of the Markov Chain Monte Carlo technique in the analysis of cosmic-ray measurements in diffusion models. The size of the diffusive halo depends crucially on the value of the diffusion slope delta, and also on the presence/absence of the local underdensity damping effect on radioactive nuclei. More precise data from on-going experiments are expected to clarify this issue.

p, He, and C to Fe cosmic-ray primary fluxes in diffusion models: Source and transport signatures on fluxes and ratios

Abstract: The propagated fluxes of proton, helium, and heavier primary cosmic-ray species (up to Fe) are a means to indirectly access the source spectrum of cosmic rays. We check the compatibility of the primary fluxes with the transport parameters derived from the B/C analysis, but also if they bring further constraints. Proton data are well described in the simplest model defined by a power-law source spectrum and plain diffusion. They can also be accommodated by models with, e.g., convection and/or reacceleration. There is no need for breaks in the source spectral indices below $\sim 1$ TeV/n. Fits on the primary fluxes alone do not provide physical constraints on the transport parameters. If we let free the source spectrum $dQ/dE = q \beta^{\eta_S} {\cal R}^{-\alpha}$ and fix the diffusion coefficient $K(R)= K_0\beta^{\eta_T} {\cal R}^{\delta}$ such as to reproduce the B/C ratio, the MCMC analysis constrains the source spectral index $\alpha$ to be in the range $2.2-2.5$ for all primary species up to Fe, regardless of the value of the diffusion slope $\delta$. The $\eta_S$ low-energy shape of the source spectrum is degenerate with the low-energy shape $\eta_T$ of the diffusion coefficient: we find $\eta_S-\eta_T\approx 0$ for p and He data, but $\eta_S-\eta_T\approx 1$ for C to Fe primary species. This is consistent with the toy-model calculation in which the shape of the p/He and C/O to Fe/O data is reproduced if $\eta_S-\eta_T\approx 0-1$ (no need for different slopes $\alpha$). When plotted as a function of the kinetic energy per nucleon, the low-energy p/He ratio is shaped mostly by the modulation effect, whereas primary/O ratios are mostly shaped by their destruction rate.