• Probing X-ray burst -- accretion disk interaction in low mass X-ray binaries through kilohertz quasiperiodic oscillations(1403.5471)

Aug. 20, 2014 astro-ph.HE
The intense radiation flux of Type I X-ray bursts is expected to interact with the accretion flow around neutron stars. High frequency quasiperiodic oscillations (kHz QPOs), observed at frequencies matching orbital frequencies at tens of gravitational radii, offer a unique probe of the innermost disk regions. In this paper, we follow the lower kHz QPOs, in response to Type I X-ray bursts, in two prototypical QPO sources, namely 4U 1636-536 and 4U 1608-522, as observed by the Proportional Counter Array of the Rossi X-ray Timing Explorer. We have selected a sample of 15 bursts for which the kHz QPO frequency can be tracked on timescales commensurable with the burst durations (tens of seconds). We find evidence that the QPOs are affected for over ~200 s during one exceptionally long burst and ~100 s during two others (although at a less significant level), while the burst emission has already decayed to a level that would enable the pre-burst QPO to be detected. On the other hand, for most of our burst-kHz QPO sample, we show that the QPO is detected as soon as the statistics allow and in the best cases, we are able to set an upper limit of ~20 s on the recovery time of the QPO. This diversity of behavior cannot be related to differences in burst peak luminosity. We discuss these results in the framework of recent findings that accretion onto the neutron star may be enhanced during Type I X-ray bursts. The subsequent disk depletion could explain the disappearance of the QPO for ~100 s, as possibly observed in two events. However, alternative scenarios would have to be invoked for explaining the short recovery timescales inferred from most bursts. Clearly the combination of fast timing and spectral information of Type I X-ray bursts holds great potential in the study of the dynamics of the inner accretion flow around neutron stars.
• The high-frequency peaked BL Lac object PKS 2005-489 was the target of a multi-wavelength campaign with simultaneous observations in the TeV gamma-ray (H.E.S.S.), GeV gamma-ray (Fermi/LAT), X-ray (RXTE, Swift), UV (Swift) and optical (ATOM, Swift) bands. This campaign was carried out during a high flux state in the synchrotron regime. The flux in the optical and X-ray bands reached the level of the historical maxima. The hard GeV spectrum observed with Fermi/LAT connects well to the very high energy (VHE, E>100GeV) spectrum measured with H.E.S.S. with a peak energy between ~5 and 500 GeV. Compared to observations with contemporaneous coverage in the VHE and X-ray bands in 2004, the X-ray flux was ~50 times higher during the 2009 campaign while the TeV gamma-ray flux shows marginal variation over the years. The spectral energy distribution during this multi-wavelength campaign was fit by a one zone synchrotron self-Compton model with a well determined cutoff in X-rays. The parameters of a one zone SSC model are inconsistent with variability time scales. The variability behaviour over years with the large changes in synchrotron emission and small changes in the inverse Compton emission does not warrant an interpretation within a one-zone SSC model despite an apparently satisfying fit to the broadband data in 2009.
• Search for a Dark Matter annihilation signal from the Galactic Center halo with H.E.S.S(1103.3266)

March 16, 2011 hep-ph, astro-ph.HE
A search for a very-high-energy (VHE; >= 100 GeV) gamma-ray signal from self-annihilating particle Dark Matter (DM) is performed towards a region of projected distance r ~ 45-150 pc from the Galactic Center. The background-subtracted gamma-ray spectrum measured with the High Energy Stereoscopic System (H.E.S.S.) gamma-ray instrument in the energy range between 300 GeV and 30 TeV shows no hint of a residual gamma-ray flux. Assuming conventional Navarro-Frenk-White (NFW) and Einasto density profiles, limits are derived on the velocity-weighted annihilation cross section < \sigma v> as a function of the DM particle mass. These are among the best reported so far for this energy range. In particular, for the DM particle mass of ~1 TeV, values for <\sigma v> above 3 * 10^(-25) cm^3 s^(-1) are excluded for the Einasto density profile. The limits derived here differ much less for the chosen density profile parametrizations, as opposed to limits from gamma-ray observations of dwarf galaxies or the very center of the Milky Way, where the discrepancy is significantly larger.
• We report on a newly detected point-like source, HESS J1943+213 located in the Galactic plane. This source coincides with an unidentified hard X-ray source IGR J19443+2117, which was proposed to have radio and infrared counterparts. HESS J1943+213 is detected at the significance level of 7.9 \sigma (post-trials) at RA(J2000)=19h 43m 55s +- 1s (stat) +- 1s (sys), DEC(J2000) = +21deg 18' 8" +- 17" (stat) +- 20" (sys). The source has a soft spectrum with photon index Gamma = 3.1 +- 0.3 (stat) +- 0.2 (sys) and a flux above 470 GeV of 1.3 +- 0.2 (stat) +- 0.3 (sys) x 10^{-12} cm^{-2} s^{-1}. There is no Fermi/LAT counterpart down to a flux limit of 6 x 10^{-9} cm^{-2} s^{-1} in the 0.1-100 GeV energy range (95% confidence upper limit calculated for an assumed power-law model with a photon index Gamma=2.0). The data from radio to VHE gamma-rays do not show any significant variability. We combine new H.E.S.S., Fermi/LAT and Nancay Radio Telescope observations with pre-existing non-simultaneous multi-wavelength observations of IGR J19443+2117 and discuss the likely source associations as well as the interpretation as an active galactic nucleus, a gamma-ray binary or a pulsar wind nebula. The lack of a massive stellar counterpart disfavors the binary hypothesis, while the soft VHE spectrum would be very unusual in case of a pulsar wind nebula. In addition, the distance estimates for Galactic counterparts places them outside of the Milky Way. All available observations favor an interpretation as an extreme, high-frequency peaked BL Lac object with a redshift z>0.14. This would be the first time a blazar is detected serendipitously from ground-based VHE observations, and the first VHE AGN detected in the Galactic Plane.
• Detection of very-high-energy gamma-ray emission from the vicinity of PSR B1706-44 and G343.1-2.3 with H.E.S.S(1102.0773)

Feb. 10, 2011 astro-ph.HE
The gamma-ray pulsar PSR B1706-44 and the adjacent supernova remnant (SNR) candidate G343.1-2.3 were observed by H.E.S.S. during a dedicated observation campaign in 2007. As a result of this observation campaign, a new source of very-high-energy (VHE; E > 100 GeV) gamma-ray emission, HESS J1708-443, was detected with a statistical significance of 7 sigma, although no significant point-like emission was detected at the position of the energetic pulsar itself. In this paper, the morphological and spectral analyses of the newly-discovered TeV source are presented. The centroid of HESS J1708-443 is considerably offset from the pulsar and located near the apparent center of the SNR, at RA(J2000) = 17h08m11s +/- 17s and Dec(J2000) = -44d20' +/- 4'. The source is found to be significantly more extended than the H.E.S.S. point spread function (~0.1 deg), with an intrinsic Gaussian width of 0.29 deg +/- 0.04 deg. Its integral flux between 1 and 10 TeV is ~ 3.8 x 10^-12 ph cm^-2 s^-1, equivalent to 17% of the Crab Nebula flux in the same energy range. The measured energy spectrum is well-fit by a power law with a relatively hard photon index Gamma = 2.0 +/- 0.1 (stat) +/- 0.2 (sys). Additional multi-wavelength data, including 330 MHz VLA observations, were used to investigate the VHE gamma-ray source's possible associations with the pulsar wind nebula of PSR B1706-44 and/or with the complex radio structure of the partial shell-type SNR G343.1-2.3.
• Search for Lorentz Invariance breaking with a likelihood fit of the PKS 2155-304 flare data taken on MJD 53944(1101.3650)

Feb. 2, 2011 gr-qc, astro-ph.HE
Several models of Quantum Gravity predict Lorentz Symmetry breaking at energy scales approaching the Planck scale (10^{19} GeV). With present photon data from the observations of distant astrophysical sources, it is possible to constrain the Lorentz Symmetry breaking linear term in the standard photon dispersion relations. Gamma-ray Bursts (GRB) and flaring Active Galactic Nuclei (AGN) are complementary to each other for this purpose, since they are observed at different distances in different energy ranges and with different levels of variability. Following a previous publication of the High Energy Stereoscopic System (H.E.S.S.) collaboration, a more sensitive event-by-event method consisting of a likelihood fit is applied to PKS 2155-304 flare data of MJD 53944 (July 28, 2006) as used in the previous publication. The previous limit on the linear term is improved by a factor of ~3 up to M^{l}_{QG} > 2.1x10^{18} GeV and is currently the best result obtained with blazars. The sensitivity to the quadratic term is lower and provides a limit of M^{q}_{QG} > 6.4x10^10 GeV, which is the best value obtained so far with an AGN and similar to the best limits obtained with GRB.
• H.E.S.S. constraints on Dark Matter annihilations towards the Sculptor and Carina Dwarf Galaxies(1012.5602)

Dec. 27, 2010 astro-ph.HE
The Sculptor and Carina Dwarf spheroidal galaxies were observed with the H.E.S.S. Cherenkov telescope array between January 2008 and December 2009. The data sets consist of a total of 11.8 and 14.8 hours of high quality data, respectively. No gamma-ray signal was detected at the nominal positions of these galaxies above 220 GeV and 320 GeV, respectively. Upper limits on the gamma-ray fluxes at 95% C.L. assuming two forms for the spectral energy distribution (a power law shape and one derived from dark matter annihilation) are obtained at the level of 10^-13 to 10^-12 cm^-2s^-1 in the TeV range. Constraints on the velocity weighted dark matter particle annihilation cross section for both Sculptor and Carina dwarf galaxies range from <\sigma v> ~ 10^-21 cm^3s^-1 down to <\sigma v> ~ 10^-22 cm^3s^-1 depending on the dark matter halo model used. Possible enhancements of the gamma-ray flux are studied: the Sommerfeld effect, which is found to exclude some dark matter particle masses, the internal Bremsstrahlung and clumps in the dark-matter halo distributions.
• Discovery and follow-up studies of the extended, off-plane, VHE gamma-ray source HESS J1507-622(1010.4907)

Oct. 23, 2010 astro-ph.HE
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.
• 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.
• 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.
• 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 ..
• 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.
• 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.
• High Energy Observations of XRF 030723: Evidence for an Off-axis Gamma-Ray Burst?(astro-ph/0408453)

Jan. 13, 2005 astro-ph
We report High Energy Transient Explorer 2 (HETE-2) Wide Field X-ray Monitor/French Gamma Telescope observations of XRF030723 along with observations of the XRF afterglow made using the 6.5m Magellan Clay telescope and the Chandra X-ray Observatory. The observed peak energy E_pk_obs of the nu F_nu burst spectrum is found to lie within (or below) the WXM 2-25 keV passband at 98.5% confidence, and no counts are detected above 30 keV. Our best fit value is E_pk_obs=8.4+3.5/-3.4 keV. The ratio of X-ray to Gamma-ray flux for the burst follows a correlation found for GRBs observed with HETE-2, and the duration of the burst is similar to that typical of long-duration GRBs. If we require that the burst isotropic equivalent energy E_iso and E_pk_rest satisfy the relation discovered by Amati et al. (2002), a redshift of z=0.38+0.36/-0.18 can be determined, in agreement with constraints determined from optical observations. We are able to fit the X-ray afterglow spectrum and to measure its temporal fade. Although the best-fit fade is shallower than the concurrent fade in the optical, the spectral similarity between the two bands indicates that the X-ray fade may actually trace the optical fade. If this is the case, the late time rebrightening observed in the optical cannot be due to a supernova bump. We interpret the prompt and afterglow X-ray emission as arising from a jetted GRB observed off-axis and possibly viewed through a complex circumburst medium due to a progenitor wind.
• Global Characteristics of X-Ray Flashes and X-Ray-Rich GRBs Observed by HETE-2(astro-ph/0409128)

Sept. 6, 2004 astro-ph
We describe and discuss the global properties of 45 gamma-ray bursts (GRBs) observed by HETE-2 during the first three years of its mission, focusing on the properties of X-Ray Flashes (XRFs) and X-ray-rich GRBs (XRRs). We find that the numbers of XRFs, XRRs, and GRBs are comparable. We find that the durations and the sky distributions of XRFs and XRRs are similar to those of GRBs. We also find that the spectral properties of XRFs and XRRs are similar to those of GRBs, except that the values of the peak energy $E^{\rm obs}_{\rm peak}$ of the burst spectrum in $\nu F_\nu$, the peak energy flux $\Fp$, and the energy fluence $S_E$ of XRFs are much smaller -- and those of XRRs are smaller -- than those of GRBs. Finally, we find that the distributions of all three kinds of bursts form a continuum in the [$S_E$(2-30 keV),$S_E$(30-400) keV]-plane, the [$S_E$(2-400 keV), $E_{\rm peak}$]-plane, and the [$F_{\rm peak}$(50-300 keV), $E_{\rm peak}$]-plane. These results provide strong evidence that all three kinds of bursts arise from the same phenomenon.
• HETE-2 Observations of the Extremely Soft X-Ray Flash XRF 020903(astro-ph/0309455)

Nov. 10, 2003 astro-ph
We report HETE-2 WXM/FREGATE observations of the X-ray flash, XRF 020903. This event was extremely soft: the ratio log(S_X/S_gamma) = 0.7, where S_X and S_gamma are the fluences in the 2-30 and 30-400 keV energy bands, is the most extreme value observed so far by HETE-2. In addition, the spectrum has an observed peak energy E^{obs}_{peak} < 5.0 keV (99.7 % probability upper limit) and no photons were detected above ~10 keV. The burst is shorter at higher energies, which is similar to the behavior of long GRBs. We consider the possibility that the burst lies at very high redshift and that the low value of E^{obs}_{peak} is due to the cosmological redshift, and show that this is very unlikely. We find that the properties of XRF 020903 are consistent with the relation between the fluences S(7-30 keV) and S(30-400 keV) found by Barraud et al. for GRBs and X-ray-rich GRBs, and are consistent with the extension by a decade of the hardness-intensity correlation (Mallozzi et al. 1995) found by the same authors. Assuming that XRF 020903 lies at a redshift z = 0.25 as implied by the host galaxy of the candidate optical and radio afterglows of this burst, we find that the properties of XRF 020903 are consistent with an extension by a factor ~300 of the relation between the isotropic-equivalent energy E_iso and the peak E_peak of the nu F_nu spectrum (in the source frame of the burst) found by Amati et al. for GRBs. The results presented in this paper therefore provide evidence that XRFs, X-ray-rich GRBs, and GRBs form a continuum and are a single phenomenon. The results also impose strong constraints on models of XRFs and X-ray-rich GRBs.
• Spectral analysis of 35 GRBs/XRFs observed with HETE-2/FREGATE(astro-ph/0206380)

Jan. 13, 2003 astro-ph
We present a spectral analysis of 35 GRBs detected with the HETE-2 gamma-ray detectors (the FREGATE instrument) in the energy range 7-400 keV. The GRB sample analyzed is made of GRBs localized with the Wide Field X-ray Monitor onboard HETE-2 or with the GRB Interplanetary Network. We derive the spectral parameters of the time-integrated spectra, and present the distribution of the low-energy photon index, alpha, and of the peak energy, e_peak . We then discuss the existence and nature of the recently discovered X-Ray Flashes and their relationship with classical GRBs.
• In-flight verification of the FREGATE spectral response(astro-ph/0203079)

March 6, 2002 astro-ph
We present the first results of the in-flight validation of the spectral response of the FREGATE X/gamma detectors on-board the HETE-2 satellite. This validation uses the Crab pulsar and nebula as reference spectra.
• FREGATE observation of a strong burst from SGR1900+14(astro-ph/0203080)

March 6, 2002 astro-ph
After a long period of quiescence, the soft gamma repeater SGR1900+14 was suddenly reactivated on April 2001. On July 2, 2001, a bright flare emitted by this source triggerred the WXM and FREGATE instruments onboard the HETE-2 satellite. Unlike typical short (0.1 s) and spiky SGRs recurrent bursts, this event features a 4.1 s long main peak, with a sharp rise (50 ms) and a slower cutoff (250 ms). This main peak is followed by a 2 sec decreasing tail. We found no evidence of any precursor or any extended `afterglow' tail to this burst. We present the preliminary spectral fits of the total emission of this flare as observed by the FREGATE instrument between 7 and 150 keV. The best fit is obtained with a model consisting of two blackbody components of temperatures 4.15 keV and 10.4 keV. A thermal bremsstrahlung can not be fitted to this spectrum. We compare these features and the burst energetics with the other strong or giant flares from SGR1900+14.
• Measurement of the Crab Flux Above 60 GeV with the CELESTE Cherenkov Telescope(astro-ph/0107301)

Dec. 5, 2001 astro-ph
We have converted the former solar electrical plant THEMIS (French Pyrenees) into an atmospheric Cherenkov detector called CELESTE, which records gamma rays above 30 GeV (7E24 Hz). Here we present the first sub-100 GeV detection by a ground based telescope of a gamma ray source, the Crab nebula, in the energy region between satellite measurements and imaging atmospheric Cherenkov telescopes. At our analysis threshold energy of 60 +/- 20 GeV we measure a gamma ray rate of 6.1 +/- 0.8 per minute. Allowing for 30% systematic uncertainties and a 30% error on the energy scale yields an integral gamma ray flux of I(E>60 GeV) = 6.2^{+5.3}_{-2.3} E-6 photons m^-2 s^-1. The analysis methods used to obtain the gamma ray signal from the raw data are detailed. In addition, we determine the upper limit for pulsed emission to be <12% of the Crab flux at the 99% confidence level, in the same energy range. Our result indicates that if the power law observed by EGRET is attenuated by a cutoff of form e^{-E/E_0} then E_0 < 26 GeV. This is the lowest energy probed by a Cherenkov detector and leaves only a narrow range unexplored beyond the energy range studied by EGRET.