Source author record

X. F. Wu

X. F. Wu 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

40works
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

40 published item(s)

preprint2026arXiv

An Enigmatic PeVatron in an Area around HII Region G35.6$-$0.5

Identifying Galactic PeVatrons (PeV particle accelerators) from the ultra-high-energy (UHE, >100 TeV) $γ$-ray sources plays a crucial role in revealing the origin of Galactic cosmic rays. The UHE source 1LHAASO J1857+0203u is suggested to be associated with HESS J1858+020, which may be attributed to the possible PeVatron candidate supernova remnant (SNR) G35.6$-$0.4 or HII region G35.6$-$0.5. We perform detailed analysis on the very-high-energy and UHE $γ$-ray emissions towards this region with data from the Large High Altitude Air Shower Observatory (LHAASO). 1LHAASO J1857+0203u is detected with a significance of 11.6$σ$ above 100 TeV, indicating the presence of a PeVatron. It has an extension of $\sim 0.18^\circ$ with a power-law (PL) spectral index of $\sim$2.5 in 1-25 TeV and a point-like emission with a PL spectral index of $\sim$3.2 above 25 TeV. Using the archival CO and HI data, we identify some molecular and atomic clouds that may be associated with the TeV $γ$-ray emissions. Our modelling indicates that the TeV $γ$-ray emissions are unlikely to arise from the clouds illuminated by the protons that escaped from SNR G35.6$-$0.4. In the scenario that HII region G35.6$-$0.5 could accelerate particles to the UHE band, the observed GeV-TeV $γ$-ray emission could be well explained by a hadronic model with a PL spectral index of $\sim$2.0 and cutoff energy of $\sim$450 TeV. However, an evolved pulsar wind nebula origin cannot be ruled out.

preprint2026arXiv

An Ultrahigh-energy $γ$-ray Bubble Powered by a Super PeVatron

We report the detection of a $γ$-ray bubble spanning at least 100$\rm deg^2$ in ultra high energy (UHE) up to a few PeV in the direction of the star-forming region Cygnus X, implying the presence Super PeVatron(s) accelerating protons to at least 10 PeV. A log-parabola form with the photon index $Γ(E) = (2.71 \pm 0.02) + (0.11 \pm 0.02) \times \log_{10} (E/10 \ {\rm TeV})$ is found fitting the gamma-ray energy spectrum of the bubble well. UHE sources, `hot spots' correlated with very massive molecular clouds, and a quasi-spherical amorphous $γ$-ray emitter with a sharp central brightening are observed in the bubble. In the core of $\sim 0.5^{\circ}$, spatially associating with a region containing massive OB association (Cygnus OB2) and a microquasar (Cygnus X-3), as well as previously reported multi-TeV sources, an enhanced concentration of UHE $γ$-rays are observed with 2 photons at energies above 1 PeV. The general feature of the bubble, the morphology and the energy spectrum, are reasonably reproduced by the assumption of a particle accelerator in the core, continuously injecting protons into the ambient medium.

preprint2026arXiv

Constraining the Cosmic-ray Energy Based on Observations of Nearby Galaxy Clusters by LHAASO

Galaxy clusters act as reservoirs of high-energy cosmic rays (CRs). As CRs propagate through the intracluster medium, they generate diffuse $γ$-rays detectable by arrays such as LHAASO. These $γ$-rays result from proton-proton ($pp$) collisions of very high-energy cosmic rays (VHECRs) or inverse Compton (IC) scattering of positron-electron pairs created by $pγ$ interactions of ultra-high-energy cosmic rays (UHECRs). We analyzed diffuse $γ$-ray emission from the Coma, Perseus, and Virgo clusters using LHAASO data. Diffuse emission was modeled as a disk of radius $R_{500}$ for each cluster while accounting for point sources. No significant diffuse emission was detected, yielding 95\% confidence level (C.L.) upper limits on the $γ$-ray flux: for WCDA (1-25~TeV) and KM2A ($>25$~TeV), less than $(49.4, 13.7, 54.0)$ and $(1.34, 1.14, 0.40) \times 10^{-14}$~ph~cm$^{-2}$~s$^{-1}$ for Coma, Perseus, and Virgo, respectively. The $γ$-ray upper limits can be used to derive model-independent constraints on the integral energy of CRp above 10~TeV (corresponding to the LHAASO observational range $>1$~TeV under the $pp$ scenario) to be less than $(1.96, 0.59, 0.08) \times 10^{61}$~erg. The absence of detectable annuli/ring-like structures, indicative of cluster accretion or merging shocks, imposes further constraints on models in which the UHECRs are accelerated in the merging shocks of galaxy clusters.

preprint2026arXiv

Constraints on heavy decaying dark matter from 570 days of LHAASO observations

The Kilometer Square Array~(KM2A) of the Large High Altitude Air Shower Observatory (LHAASO) aims at surveying the northern gamma-ray sky at energies above 10 TeV with unprecedented sensitivity. Gamma-ray observations have long been one of the most powerful tools for dark matter searches, as e.g., high-energy gamma-rays could be produced by the decays of heavy dark matter particles. In this letter, we present the first dark matter analysis with LHAASO-KM2A, using the first 340~days of data from 1/2-KM2A and 230~days of data from 3/4-KM2A. Several regions of interest are used to search for a signal and account for the residual cosmic-ray background after gamma/hadron separation. We find no excess of dark matter signals, and thus place some of the strongest gamma-ray constraints on the lifetime of heavy dark matter particles with mass between 10^5 and 10^9~GeV. Our results with LHAASO are robust, and have important implications for dark matter interpretations of the diffuse astrophysical high-energy neutrino emission.

preprint2026arXiv

Deep view of Composite SNR CTA1 with LHAASO in $γ$-rays up to 300 TeV

The ultra-high-energy (UHE) gamma-ray source 1LHAASO J0007+7303u is positionally associated with the composite SNR CTA1 that is located at high Galactic Latitude $b\approx 10.5^\circ$. This provides a rare opportunity to spatially resolve the component of the pulsar wind nebula (PWN) and supernova remnant (SNR) at UHE. This paper conducted a dedicated data analysis of 1LHAASO J0007+7303u using the data collected from December 2019 to July 2023. This source is well detected with significances of 21$σ$ and 17$σ$ at 8$-$100 TeV and $>$100 TeV, respectively. The corresponding extensions are determined to be 0.23$^{\circ}\pm$0.03$^{\circ}$ and 0.17$^{\circ}\pm$0.03$^{\circ}$. The emission is proposed to originate from the relativistic electrons and positrons accelerated within the PWN of PSR J0007+7303. The energy spectrum is well described by a power-law with an exponential cutoff function $dN/dE = (42.4\pm4.1)(\frac{E}{20\rm\ TeV})^{-2.31\pm0.11}\exp(-\frac{E}{110\pm25\rm\ TeV})$ $\rm\ TeV^{-1}\ cm^{-2}\ s^{-1}$in the energy range from 8 TeV to 300 TeV, implying a steady-state parent electron spectrum $dN_e/dE_e\propto (\frac{E_e}{100\rm\ TeV})^{-3.13\pm0.16}\exp[(\frac{-E_e}{373\pm70\rm\ TeV})^2]$ at energies above $\approx 50 \rm\ TeV$. The cutoff energy of the electron spectrum is roughly equal to the expected current maximum energy of particles accelerated at the PWN terminal shock. Combining the X-ray and gamma-ray emission, the current space-averaged magnetic field can be limited to $\approx 4.5\rm\ μG$. To satisfy the multi-wavelength spectrum and the $γ$-ray extensions, the transport of relativistic particles within the PWN is likely dominated by the advection process under the free-expansion phase assumption.

preprint2026arXiv

Discovery of a new $γ$-ray source LHAASO J0341+5258 with emission up to 200TeV

We report the discovery of a new unidentified extended $γ$-ray source in the Galactic plane named LHAASO J0341+5258 with a pre-trial significance of 8.2 standard deviations above 25 TeV. The best fit position is R.A.$=55.34^{\circ}\pm0.11^{\circ}$ and Dec$=52.97^{\circ}\pm0.07^{\circ}$. The angular size of LHAASO J0341+5258 is $0.29^\circ \pm 0.06^\circ_{stat} \pm0.02^\circ_{sys}$. The flux above 25 TeV is about $20\%$ of the flux of Crab Nebula. Although a power-law fit of the spectrum from 10 TeV to 200 TeV with the photon index $α=2.98 \pm 0.19_{stat} \pm 0.02_{sys}$ is not excluded, the LHAASO data together with the flux upper limit at 10 GeV set by the Fermi LAT observation, indicate a noticeable steepening of an initially hard power-law spectrum %($α\leq 1.75$) spectrum with a cutoff at $\approx 50$ TeV. We briefly discuss the origin of UHE gamma-rays. The lack of an energetic pulsar and a young SNR inside or in the vicinity of LHAASO J0341+5258 challenge, but do not exclude both the leptonic and hadronic scenarios of gamma-ray production.

preprint2026arXiv

Discovery of the Ultra-high energy gamma-ray source LHAASO J2108+5157

We report the discovery of a UHE gamma-ray source, LHAASO J2108+5157, by analyzing the LHAASO-KM2A data of 308.33 live days. Significant excess of gamma-ray induced showers is observed in both energy bands of 25-100 TeV and $\gt$100 TeV with 9.5 sigma and 8.5 sigma, respectively. This source is not significantly favored as an extensive source with the angular extension smaller than the point-spread function of KM2A. The measured energy spectrum from 20 to 200 TeV can be approximately described by a power-law function with an index of -2.83$\pm$ 0.18stat. A harder spectrum is demanded at lower energies considering the flux upper limit set by Fermi-LAT observations. The position of the gamma-ray emission is correlated with a giant molecular cloud, which favors a hadronic origin. No obvious counterparts have been found, deeper multiwavelength observations will help to shed new light on this intriguing UHE source.

preprint2026arXiv

Energy calibration of LHAASO-KM2A using the cosmic ray Moon shadow

We present a precise measurement of the westward, rigidity-dependent shift of the Moon's shadow using three and a half years of cosmic-ray data collected by the Kilometer Square Array (KM2A) of the Large High Altitude Air Shower Observatory (LHAASO). These measurements enable us to calibrate the detector energy response in the range 20-260 TeV, with results showing excellent agreement with the response derived from Monte Carlo (MC) simulations of the KM2A detector. We also measure a best-fit parameter $ε= 0.015 \pm 0.08$, corresponding to a 95% confidence interval of [-14%, +17%] for the energy-scale estimation. This result establishes the exceptional accuracy of the KM2A-MC in simulating the detector's response within this energy range.

preprint2026arXiv

Energy-Dependent Shifts of Medium-Scale Anisotropies in Very-High-Energy Cosmic Rays Observed by LHAASO-KM2A

Small deviations from isotropy in the arrival directions of Galactic cosmic rays serve as a unique probe of the local magnetic environment. In this Letter, we report observations of medium-scale anisotropies (MSA) at energies above 10 TeV using the LHAASO-KM2A array. Our analysis identifies four regions of excess and four regions of deficit, each spanning angular scales of approximately ten degrees. Crucially, we detect significant energy-dependent shifts in the centroids of two excess regions: Region B and the newly identified Region $\mathrm{\widetilde{D}}$. We also characterize the energy evolution of the fractional relative intensity across both excess and deficit regions. These findings imply that the observed anisotropies are shaped by the specific realization of the local turbulent magnetic field within the cosmic ray scattering length. Such energy-dependent behaviors impose strict constraints on local turbulence models and cosmic ray propagation theories.

preprint2026arXiv

Evidence for particle acceleration approaching PeV energies in the W51 complex

The $γ$-ray emission from the W51 complex is widely acknowledged to be attributed to the interaction between the cosmic rays (CRs) accelerated by the shock of supernova remnant (SNR) W51C and the dense molecular clouds in the adjacent star-forming region, W51B. However, the maximum acceleration capability of W51C for CRs remains elusive. Based on observations conducted with the Large High Altitude Air Shower Observatory (LHAASO), we report a significant detection of $γ$ rays emanating from the W51 complex, with energies from 2 TeV to 200 TeV. The LHAASO measurements, for the first time, extend the $γ$-ray emission from the W51 complex beyond 100 TeV and reveal a significant spectrum bending at tens of TeV. By combining the ``$π^0$-decay bump" featured data from Fermi-LAT, the broadband $γ$-ray spectrum of the W51 region can be well-characterized by a simple pp-collision model. The observed spectral bending feature suggests an exponential cutoff at $\sim400$~TeV or a power-law break at $\sim200$~TeV in the CR proton spectrum, most likely providing the first evidence of SNRs serving as CR accelerators approaching the PeV regime. Additionally, two young star clusters within W51B could also be theoretically viable to produce the most energetic $γ$ rays observed by LHAASO. Our findings strongly support the presence of extreme CR accelerators within the W51 complex and provide new insights into the origin of Galactic CRs.

preprint2026arXiv

Exploring Lorentz Invariance Violation from Ultra-high-energy Gamma Rays Observed by LHAASO

Recently the LHAASO Collaboration published the detection of 12 ultra-high-energy gamma-ray sources above 100 TeV, with the highest energy photon reaching 1.4 PeV. The first detection of PeV gamma rays from astrophysical sources may provide a very sensitive probe of the effect of the Lorentz invariance violation (LIV), which results in decay of high-energy gamma rays in the superluminal scenario and hence a sharp cutoff of the energy spectrum. Two highest energy sources are studied in this work. No signature of the existence of LIV is found in their energy spectra, and the lower limits on the LIV energy scale are derived. Our results show that the first-order LIV energy scale should be higher than about 10^5 times the Planck scale M_{pl} and that the second-order LIV scale is >10^{-3}M_{pl}. Both limits improve by at least one order of magnitude the previous results.

preprint2026arXiv

Extended Very-High-Energy Gamma-Ray Emission Surrounding PSR J0622 + 3749 Observed by LHAASO-KM2A

We report the discovery of an extended very-high-energy (VHE) gamma-ray source around the location of the middle-aged (207.8 kyr) pulsar PSR J0622+3749 with the Large High Altitude Air Shower Observatory (LHAASO). The source is detected with a significance of $8.2σ$ for $E>25$~TeV assuming a Gaussian template. The best-fit location is (R.A., Dec.)$=(95^{\circ}\!.47\pm0^{\circ}\!.11,\,37^{\circ}\!.92 \pm0^{\circ}\!.09)$, and the extension is $0^{\circ}\!.40\pm0^{\circ}\!.07$. The energy spectrum can be described by a power-law spectrum with an index of ${-2.92 \pm 0.17_{\rm stat} \pm 0.02_{\rm sys} }$. No clear extended multi-wavelength counterpart of the LHAASO source has been found from the radio to sub-TeV bands. The LHAASO observations are consistent with the scenario that VHE electrons escaped from the pulsar, diffused in the interstellar medium, and scattered the interstellar radiation field. If interpreted as the pulsar halo scenario, the diffusion coefficient, inferred for electrons with median energies of $\sim160$~TeV, is consistent with those obtained from the extended halos around Geminga and Monogem and much smaller than that derived from cosmic ray secondaries. The LHAASO discovery of this source thus likely enriches the class of so-called pulsar halos and confirms that high-energy particles generally diffuse very slowly in the disturbed medium around pulsars.

preprint2026arXiv

LHAASO Detection of Ultra-High-Energy Gamma-Ray Emission toward the Giant Molecular Clouds

The $γ$-ray from Giant molecular clouds (GMCs) is regarded as the most ideal tool to perform in-situ measurement of cosmic ray (CR) density and spectra in our Galaxy. We report the first detection of $γ$-ray emissions in the very-high-energy (VHE) domain from the five nearby GMCs with a stacking analysis based on a 4.5-year $γ$-ray observation with the Large High Altitude Air Shower Observatory (LHAASO) experiment. The spectral energy distributions derived from the GMCs are consistent with the expected $γ$-ray flux produced via CR interacting with the ISM in the energy interval 1 - 100 $~\rm$ TeV. In addition, we investigate the presence of the CR spectral `knee' by introducing a spectral break in the $γ$-ray data. While no significant evidence for the CR knee is found, the current KM2A measurements from GMCs strongly favor a proton CR knee located above 0.9$~\rm$ PeV, which is consistent with the latest measurement of the CR spectrum by ground-based experiments.

preprint2026arXiv

Measurement of attenuation length of the muon content in extensive air showers from 0.3 to 30 PeV with LHAASO

The attenuation length of the muon content in extensive air showers provides important information regarding the generation and development of air showers. This information can be used not only to improve the description of such showers but also to test fundamental models of hadronic interactions. Using data from the LHAASO-KM2A experiment, the development of the muon content in high-energy air showers was studied. The attenuation length of muon content in the air showers was measured from experimental data in the energy range from 0.3 to 30 PeV using the constant intensity cut method. By comparing the attenuation length of the muon content with predictions from high-energy hadronic interaction models (QGSJET-II-04, SIBYLL 2.3d, and EPOS-LHC), it is evident that LHAASO results are significantly shorter than those predicted by the first two models (QGSJET-II-04 and SIBYLL 2.3d) but relatively close to those predicted by the third model (EPOS-LHC). Thus, the LHAASO data favor the EPOS-LHC model over the other two models. The three interaction models confirmed an increasing trend in the attenuation length as the cosmic-ray energy increases.

preprint2026arXiv

Measurement of Very-high-energy Diffuse Gamma-ray Emissions from the Galactic Plane with LHAASO-WCDA

The diffuse Galactic gamma-ray emission is a very important tool used to study the propagation and interaction of cosmic rays in the Milky Way. In this work, we report the measurements of the diffuse emission from the Galactic plane, covering Galactic longitudes from $15^{\circ}$ to $235^{\circ}$ and latitudes from $-5^{\circ}$ to $+5^{\circ}$, in an energy range of 1 TeV to 25 TeV, with the Water Cherenkov Detector Array (WCDA) of the Large High Altitude Air Shower Observatory (LHAASO). After masking the sky regions of known sources, the diffuse emission is detected with $24.6σ$ and $9.1σ$ significance in the inner Galactic plane and outer Galactic plane, respectively. The WCDA spectra in both regions can be well described by a power-law function, with spectral indices of $-2.67\pm0.05_{\rm stat}$ in the inner region and $-2.83\pm0.19_{\rm stat}$ in the outer region, respectively. Combined with the Square Kilometer Array (KM2A) measurements at higher energies, a clear softening of the spectrum is found in the inner region, with change of spectral indices by $\sim0.5$ at a break energy around $30$ TeV. The fluxes of the diffuse emission are higher by a factor of $1.5-2.7$ than the model prediction assuming local CR spectra and the gas column density, which are consistent with those measured by the KM2A. Along Galactic longitude, the spatial distribution of the diffuse emission shows deviation from that of the gas column density. The spectral shape of the diffuse emission are possibly variation in different longitude region. The WCDA measurements bridge the gap between the low-energy measurements by space detectors and the ultra-high-energy observations by LHAASO-KM2A and other experiments. These results suggest that improved modeling of the wide-band diffuse emission is required.

preprint2026arXiv

Transient Large-Scale Anisotropy in TeV Cosmic Rays due to an Interplanetary Coronal Mass Ejection

Large- or medium-scale cosmic ray anisotropy at TeV energies has not previously been confirmed to vary with time. Transient anisotropy changes have been observed below 150 GeV, especially near the passage of an interplanetary shock and coronal mass ejection containing a magnetic flux rope ejected by a solar storm, which can trigger a geomagnetic storm with practical consequences. In such events, cosmic rays provide remote sensing of the magnetic field properties. Here we report the observation of transient large-scale anisotropy in TeV cosmic ray ions using data from the Large High Altitude Air Shower Observatory (LHAASO). We analyze hourly skymaps of the transient cosmic ray intensity excess or deficit, the gradient of which indicates the direction and magnitude of transient large-scale anisotropy across the field of view. We observe enhanced anisotropy above typical hourly fluctuations with $>$5$σ$ significance during some hours of November 4, 2021, in separate data sets for four primary cosmic ray energy ranges of median energy from $E$=0.7 to 3.1 TeV. The gradient varies with energy as $E^γ$, where $γ\approx-0.5$. At a median energy $\leq$1.0 TeV, this gradient corresponds to a dipole anisotropy of at least 1\%, or possibly a weaker anisotropy of higher order. This new type of observation opens the opportunity to study interplanetary magnetic structures using air shower arrays around the world, complementing existing in situ and remote measurements of plasma properties.

preprint2025arXiv

Observation of the $γ$-ray Emission from W43 with LHAASO

In this paper, we report the detection of the very-high-energy (VHE, $ 100{\rm\ GeV} < E < 100{\rm\ TeV} $) and ultra-high-energy (UHE, $E > 100\rm\ TeV$) $γ$-ray emissions from the direction of the young star-forming region W43, observed by the Large High Altitude Air Shower Observation (LHAASO). The extended $γ$-ray source was detected with a significance of ${\sim}16\,σ$ by KM2A and ${\sim}17\,σ$ by WCDA, respectively. The angular extension of this $γ$-ray source is about 0.5 degrees, corresponding to a physical size of about 50 pc. We discuss the origin of the $γ$-ray emission and possible cosmic ray acceleration in the W43 region using multi-wavelength data. Our findings suggest that W43 is likely another young star cluster capable of accelerating cosmic rays (CRs) to at least several hundred TeV.

preprint2025arXiv

Study of Ultra-High-Energy Gamma-Ray Source 1LHAASO J0056+6346u and Its Possible Origins

We report a dedicated study of the newly discovered extended UHE $γ$-ray source 1LHAASO J0056+6346u. Analyzing 979 days of LHAASO-WCDA data and 1389 days of LHAASO-KM2A data, we observed a significant excess of $γ$-ray events with both WCDA and KM2A. Assuming a point power-law source with a fixed spectral index, the significance maps reveal excesses of ${\sim}12.65\,σ$, ${\sim}22.18\,σ$, and ${\sim}10.24\,σ$ in the energy ranges of 1--25 TeV, 25--100 TeV, and $> 100$ TeV, respectively. We use a 3D likelihood algorithm to derive the morphological and spectral parameters, and the source is detected with significances of $12.65\,σ$ by WCDA and $25.27\,σ$ by KM2A. The best-fit positions derived from WCDA and KM2A data are (R.A. = $13.96^\circ\pm0.09^\circ$, Decl. = $63.92^\circ\pm0.05^\circ$) and (R.A. = $14.00^\circ\pm0.05^\circ$, Decl. = $63.79^\circ\pm0.02^\circ$), respectively. The angular size ($r_{39}$) of 1LHAASO J0056+6346u is $0.34^\circ\pm0.04^\circ$ at 1--25 TeV and $0.24^\circ\pm0.02^\circ$ at $> 25$ TeV. The differential flux of this UHE $γ$-ray source can be described by an exponential cutoff power-law function: $(2.67\pm0.25) \times 10^{-15} (E/20\,\text{TeV})^{-1.97\pm0.10} e^{-E/(55.1\pm7.2)\,\text{TeV}} \,\text{TeV}^{-1}\,\text{cm}^{-2}\,\text{s}^{-1}$. To explore potential sources of $γ$-ray emission, we investigated the gas distribution around 1LHAASO J0056+6346u. 1LHAASO J0056+6346u is likely to be a TeV PWN powered by an unknown pulsar, which would naturally explain both its spatial and spectral properties. Another explanation is that this UHE $γ$-ray source might be associated with gas content illuminated by a nearby CR accelerator, possibly the SNR candidate G124.0+1.4.

preprint2020arXiv

The observation of the Crab Nebula with LHAASO-KM2A for the performance study

As a sub-array of the Large High Altitude Air Shower Observatory (LHAASO), KM2A is mainly designed to cover a large fraction of the northern sky to hunt for gamma-ray sources at energies above 10 TeV. Even though the detector construction is still underway, a half of the KM2A array has been operating stably since the end of 2019. In this paper, we present the pipeline of KM2A data analysis and the first observation on the Crab Nebula, a standard candle in very high energy gamma-ray astronomy. We detect gamma-ray signals from the Crab Nebula in both energy ranges of 10$-$100 TeV and $>$100 TeV with high significance, by analyzing the KM2A data of 136 live days between December 2019 and May 2020. With the observations, we test the detector performance including angular resolution, pointing accuracy and cosmic ray background rejection power. The energy spectrum of the Crab Nebula in the energy range 10-250 TeV fits well with a single power-law function dN/dE =(1.13$\pm$0.05$_{stat}$$\pm$0.08$_{sys}$)$\times$10$^{-14}$$\cdot$(E/20TeV)$^{-3.09\pm0.06_{stat}\pm0.02_{sys}}$ cm$^{-2}$ s$^{-1}$ TeV$^{-1}$. It is consistent with previous measurements by other experiments. This opens a new window of gamma-ray astronomy above 0.1 PeV through which ultrahigh-energy gamma-ray new phenomena, such as cosmic PeVatrons, might be discovered.

preprint2016arXiv

A Triple-Energy-Source Model for Superluminous Supernova iPTF13ehe

Almost all superluminous supernovae (SLSNe) whose peak magnitudes are $\lesssim -21$ mag can be explained by the $^{56}$Ni-powered model, magnetar-powered (highly magnetized pulsar) model or ejecta-circumstellar medium (CSM) interaction model. Recently, iPTF13ehe challenges these energy-source models, because the spectral analysis shows that $\sim 2.5M_\odot$ of $^{56}$Ni have been synthesized but are inadequate to power the peak bolometric emission of iPTF13ehe, while the rebrightening of the late-time light-curve (LC) and the H$α$ emission lines indicate that the ejecta-CSM interaction must play a key role in powering the late-time LC. Here we propose a triple-energy-source model, in which a magnetar together with some amount ($\lesssim 2.5M_\odot$) of $^{56}$Ni may power the early LC of iPTF13ehe while the late-time rebrightening can be quantitatively explained by an ejecta-CSM interaction. Furthermore, we suggest that iPTF13ehe is a genuine core-collapse supernova rather than a pulsational pair-instability supernova candidate. Further studies on similar SLSNe in the future would eventually shed light on their explosion and energy-source mechanisms.

preprint2016arXiv

Could FRB 131104 Originate from the Merger of Binary Neutron Stars?

Recently, DeLaunay et al. (2016) discovered a gamma-ray transient, Swift J0644.5-5111, associated with the fast radio burst (FRB) 131104. They also reported follow-up broadband observations beginning two days after the FRB and provided upper limits on a putative afterglow of this transient. In this paper, we show that if such a transient drives a relativistic shock as in a cosmological gamma-ray burst (GRB), these upper limits are consistent with an environment of which density is much less than that of an interstellar medium but typical for the outskirts' density of a galaxy when the typical values of three microphysical parameters of the shock are taken. This appears to be inconsistent with the catastrophic event models in which the central engine of Swift J0644.5-5111 is surrounded by an interstellar medium, but together with the properties of the gamma-ray transient, favors the binary neutron star merger origin. We further constrain the physical parameters of the postmerger object by assuming that Swift J0644.5-5111 results from internal dissipation of a spinning-down pulsar wind, and we find that the postmerger object is an ultra-strongly magnetized, very rapidly rotating pulsar. This merger event should have given birth to a gravitational wave burst, an FRB, and a short GRB or an extended X-ray/gamma-ray emission if a relativistic jet of the GRB is missed. Such "triplets" would be testable in the near future with the advanced LIGO and Virgo gravitational-wave observatories.

preprint2016arXiv

eXTP -- enhanced X-ray Timing and Polarimetry Mission

eXTP is a science mission designed to study the state of matter under extreme conditions of density, gravity and magnetism. Primary targets include isolated and binary neutron stars, strong magnetic field systems like magnetars, and stellar-mass and supermassive black holes. The mission carries a unique and unprecedented suite of state-of-the-art scientific instruments enabling for the first time ever the simultaneous spectral-timing-polarimetry studies of cosmic sources in the energy range from 0.5-30 keV (and beyond). Key elements of the payload are: the Spectroscopic Focusing Array (SFA) - a set of 11 X-ray optics for a total effective area of about 0.9 m^2 and 0.6 m^2 at 2 keV and 6 keV respectively, equipped with Silicon Drift Detectors offering <180 eV spectral resolution; the Large Area Detector (LAD) - a deployable set of 640 Silicon Drift Detectors, for a total effective area of about 3.4 m^2, between 6 and 10 keV, and spectral resolution <250 eV; the Polarimetry Focusing Array (PFA) - a set of 2 X-ray telescope, for a total effective area of 250 cm^2 at 2 keV, equipped with imaging gas pixel photoelectric polarimeters; the Wide Field Monitor (WFM) - a set of 3 coded mask wide field units, equipped with position-sensitive Silicon Drift Detectors, each covering a 90 degrees x 90 degrees FoV. The eXTP international consortium includes mostly major institutions of the Chinese Academy of Sciences and Universities in China, as well as major institutions in several European countries and the United States. The predecessor of eXTP, the XTP mission concept, has been selected and funded as one of the so-called background missions in the Strategic Priority Space Science Program of the Chinese Academy of Sciences since 2011. The strong European participation has significantly enhanced the scientific capabilities of eXTP. The planned launch date of the mission is earlier than 2025.

preprint2016arXiv

GRBs and fundamental physics

Gamma-ray bursts (GRBs) are short and intense flashes at the cosmological distances, which are the most luminous explosions in the Universe. The high luminosities of GRBs make them detectable out to the edge of the visible universe. So, they are unique tools to probe the properties of high-redshift universe: including the cosmic expansion and dark energy, star formation rate, the reionization epoch and the metal evolution of the Universe. First, they can be used to constrain the history of cosmic acceleration and the evolution of dark energy in a redshift range hardly achievable by other cosmological probes. Second, long GRBs are believed to be formed by collapse of massive stars. So they can be used to derive the high-redshift star formation rate, which can not be probed by current observations. Moreover, the use of GRBs as cosmological tools could unveil the reionization history and metal evolution of the Universe, the intergalactic medium (IGM) properties and the nature of first stars in the early universe. But beyond that, the GRB high-energy photons can be applied to constrain Lorentz invariance violation (LIV) and to test Einstein's Equivalence Principle (EEP). In this paper, we review the progress on the GRB cosmology and fundamental physics probed by GRBs.

preprint2016arXiv

Optical Transients Powered by Magnetars: Dynamics, Light Curves, and Transition to the Nebular Phase

Millisecond magnetars can be formed via several channels: core-collapse of massive stars, accretion-induced collapse of white dwarfs (WDs), double WD mergers, double neutron star (NS) mergers, and WD-NS mergers. Because the mass of ejecta from these channels could be quite different, their light curves are also expected to be diverse. We evaluate the dynamic evolution of optical transients powered by millisecond magnetars. We find that the magnetar with short spin-down timescale converts its rotational energy mostly into the kinetic energy of the transient, while the energy of a magnetar with long spin-down timescale goes into radiation of the transient. This leads us to speculate that hypernovae could be powered by magnetars with short spin-down timescales. At late times the optical transients will gradually evolve into a nebular phase because of the photospheric recession. We treat the photosphere and nebula separately because their radiation mechanisms are different. In some cases the ejecta could be light enough that the magnetar can accelerate it to a relativistic speed. It is well known that the peak luminosity of a supernova (SN) occurs when the luminosity is equal to the instantaneous energy input rate, as shown by Arnett (1979). We show that photospheric recession and relativistic motion can modify this law. The photospheric recession always leads to a delay of the peak time $t_{\mathrm{pk}}$ relative to the time $t_{\times }$ at which the SN luminosity equals the instantaneous energy input rate. Relativistic motion, however, may change this result significantly.

preprint2016arXiv

Probing the Birth of Post-merger Millisecond Magnetars by X-ray and Gamma-ray Emission

There is growing evidence that a stable magnetar could be formed from the coalescence of double neutron stars. In previous papers, we investigated the signature of formation of stable millisecond magnetars in radio and optical/ultraviolet bands by assuming that the central rapidly rotating magnetar deposits its rotational energy in the form of a relativistic leptonized wind. We found that the optical transient PTF11agg could be the first evidence for the formation of post-merger millisecond magnetars. To enhance the probability of finding more evidence for the post-merger magnetar formation, it is better to extend the observational channel to other photon energy bands. In this paper we propose to search the signature of post-merger magnetar formation in X-ray and especially gamma-ray bands. We calculate the SSC emission of the reverse shock powered by post-merger millisecond magnetars. We find that the SSC component peaks at $1\,{\rm GeV}$ in the spectral energy distribution and extends to $\gtrsim 10\,{\rm TeV}$ for typical parameters. These energy bands are quite suitable for Fermi/LAT and CTA, which, with their current observational sensitivities, can detect the SSC emission powered by post-merger magnetars up to $1\,{\rm Gpc}$. NuSTAR, sensible in X-ray bands, can detect the formation of post-merger millisecond magnetars at redshift $z\sim 1$. Future improvement in sensitivity of CTA can also probe the birth of post-merger millisecond magnetars at redshift $z\sim 1$. However, because of the $γ$-$γ$ collisions, strong high-energy emission is clearly predicted only for ejecta masses lower than $10^{-3}M_\odot$.

preprint2016arXiv

Repeating Fast Radio Bursts from Highly Magnetized Pulsars Travelling through Asteroid Belts

Very recently Spitler et al. (2016) and Scholz et al. (2016) reported their detections of sixteen additional bright bursts from the direction of the fast radio burst (FRB) 121102. This repeating FRB is inconsistent with all the catastrophic event models put forward previously for hypothetically non-repeating FRBs. Here we propose a different model, in which highly magnetized pulsars travel through asteroid belts of other stars. We show that a repeating FRB could originate from such a pulsar encountering lots of asteroids in the belt. During each pulsar-asteroid impact, an electric field induced outside the asteroid has such a large component parallel to the stellar magnetic field that electrons are torn off the asteroidal surface and accelerated to ultra-relativistic energies instantaneously. Subsequent movement of these electrons along magnetic field lines will cause coherent curvature radiation, which can account for all the properties of an FRB. In addition, this model can self-consistently explain the typical duration, luminosity, and repetitive rate of the seventeen bursts of FRB 121102. The predicted occurrence rate of repeating FRB sources may imply that our model would be testable in the next few years.

preprint2015arXiv

A Unified Energy-Reservoir Model Containing Contributions from $^{56}$Ni and Neutron Stars and Its Implication to Luminous Type Ic Supernovae

Most type-Ic core-collapse supernovae (CCSNe) produce $^{56}$Ni and neutron stars (NSs) or black holes (BHs). The dipole radiation of nascent NSs has usually been neglected in explaining supernovae (SNe) with peak absolute magnitude $M_{\rm peak}$ in any band are $\gtrsim -19.5$~mag, while the $^{56}$Ni can be neglected in fitting most type-Ic superluminous supernovae (SLSNe Ic) whose $M_{\rm peak}$ in any band are $\lesssim -21$~mag, since the luminosity from a magnetar (highly magnetized NS) can outshine that from a moderate amount of $^{56}$Ni. For luminous SNe Ic with $-21 \lesssim M_{\rm peak}\lesssim -19.5$~mag, however, both contributions from $^{56}$Ni and NSs cannot be neglected without serious modeling, since they are not SLSNe and the $^{56}$Ni mass could be up to $\sim 0.5 M_{\odot}$. In this paper we propose a unified model that contain contributions from both $^{56}$Ni and a nascent NS. We select three luminous SNe Ic-BL, SN~2010ay, SN~2006nx, and SN~14475, and show that, if these SNe are powered by $^{56}$Ni, the ratio of $M_{\rm Ni}$ to $M_{\rm ej}$ are unrealistic. Alternatively, we invoke the magnetar model and the hybrid ($^{56}$Ni + NS) model and find that they can fit the observations, indicating that our models are valid and necessary for luminous SNe Ic. Owing to the lack of late-time photometric data, we cannot break the parameter degeneracy and thus distinguish among the model parameters, but we can expect that future multi-epoch observations of luminous SNe can provide stringent constraints on $^{56}$Ni yields and the parameters of putative magnetars.

preprint2015arXiv

Fall back accretion and energy injections in gamma-ray bursts

Intense flares that occur at late times relative to the prompt phase have been observed by the $Swift$ satellite in the X-ray afterglows of gamma-ray bursts (GRBs). Here, we present a detailed analysis on the fall back accretion process to explain the intense flare phase in the very early X-ray afterglow light curves. To reproduce the afterglow at late times, we resort to the external shock by engaging energy injections. By applying our model to GRBs 080810, 081028 and 091029, we show that their X-ray afterglow light curves can be reproduced well. We then apply our model to the ultra-long $Swift$ GRB 111209A, which is the longest burst ever observed. The very early X-ray afterglow of GRB 111209A showed many interesting features, such as a significant bump observed at around 2000 s after the $Swift$/BAT trigger. We assume two constant energy injection processes in our model. These can explain the observed plateau at X-ray wavelength in the relatively early stage ($8.0\times10^{3}$ s) and a second X-ray plateau and optical rebrightening at about $10^{5}$ s. Our analysis supports the scenario that a significant amount of material may fall back toward the central engine after the prompt phase, causing an enhanced and long lived mass accretion rate powering a Poynting-flux-dominated outflow.

preprint2015arXiv

Second-order correlation function from asymmetric to symmetric transitions due to spectrally indistinguishable biexciton cascade emission

We report the observed photon bunching statistics of biexciton cascade emission at zero time delay in single quantum dots by second-order correlation function measurements under continuous wave excitation. It is found that the bunching phenomenon is independent of the biexciton binding energy when it varies from 0.59 meV to nearly zero. The photon bunching takes place when the exciton photon is not spectrally distinguishable from biexciton photon, and either of them can trigger the start in a Hanbury-Brown and Twiss setup. However, if the exciton energy is spectrally distinguishable from the biexciton the photon statistics becomes asymmetric and a cross-bunching lineshape is obtained. The theoretical calculations based on a model of three-level rate-equation analysis are consistent with the result of second-order correlation function measurements.

preprint2015arXiv

Signature of a spin-up magnetar from multi-band afterglow rebrightening of GRB 100814A

In recent years, more and more gamma-ray bursts with late rebrightenings in multi-band afterglows unveil the late-time activities of the central engines. GRB 100814A is a special one among the well-sampled events, with complex temporal and spectral evolution. The single power-law shallow decay index of the optical light curve observed by GROND between 640 s and 10 ks is $α_{\rm opt} = 0.57 \pm 0.02$, which apparently conflicts with the simple external shock model expectation. Especially, there is a remarkable rebrightening in the optical to near infrared bands at late time, challenging the external shock model with synchrotron emission coming from the interaction of the blast wave with the surrounding interstellar medium. In this paper, we invoke a magnetar with spin evolution to explain the complex multi-band afterglow emission of GRB 100814A. The initial shallow decay phase in optical bands and the plateau in X-ray can be explained as due to energy injection from a spin-down magnetar. At late time, with the falling of materials from the fall-back disk onto the central object of the burster, angular momentum of the accreted materials is transferred to the magnetar, which leads to a spin-up process. As a result, the magnetic dipole radiation luminosity will increase, resulting in the significant rebrightening of the optical afterglow. It is shown that the observed multi-band afterglow emission can be well reproduced by the model.

preprint2014arXiv

Superluminous Supernovae Powered by Magnetars: Late-time Light Curves and Hard Emission Leakage

Recently, researches performed by two groups have revealed that the magnetar spin-down energy injection model with full energy trapping can explain the early-time light curves of SN 2010gx, SN 2013dg, LSQ12dlf, SSS120810 and CSS121015, but fails to fit the late-time light curves of these Superluminous Supernovae (SLSNe). These results imply that the original magnetar-powered model is challenged in explaining these SLSNe. Our paper aims to simultaneously explain both the early- and late-time data/upper limits by considering the leakage of hard emissions. We incorporate quantitatively the leakage effect into the original magnetar-powered model and derive a new semi-analytical equation. Comparing the light curves reproduced by our revised magnetar-powered model to the observed data and/or upper limits of these five SLSNe, we found that the late-time light curves reproduced by our semi-analytical equation are in good agreement with the late-time observed data and/or upper limits of SN 2010gx, CSS121015, SN 2013dg and LSQ12dlf and the late-time excess of SSS120810, indicating that the magnetar-powered model might be responsible for these SLSNe and that the gamma ray and X-ray leakage are unavoidable when the hard photons were down-Comptonized to softer photons. To determine the details of the leakage effect and unveil the nature of SLSNe, more high quality bolometric light curves and spectra of SLSNe are required.

preprint2014arXiv

The Optical Luminosity Function of Gamma-ray Bursts deduced from ROTSE-III Observations

We present the optical luminosity function (LF) of gamma-ray bursts (GRBs) estimated from a uniform sample of 58 GRBs from observations with the Robotic Optical Transient Search Experiment III (ROTSE-III). Our GRB sample is divided into two sub-samples: detected afterglows (18 GRBs), and those with upper limits (40 GRBs). The $R$ band fluxes 100s after the onset of the burst for these two sub-samples are derived. The optical LFs at 100s are fitted by assuming that the co-moving GRB rate traces the star-formation rate. The detection function of ROTSE-III is taken into account during the fitting of the optical LFs by using Monte Carlo simulations. We find that the cumulative distribution of optical emission at 100s is well-described with an exponential rise and power-law decay (ERPLD), broken power-law (BPL), and Schechter LFs. A single power-law (SPL) LF, on the other hand, is ruled out with high confidence.

preprint2014arXiv

The origin of the plateau and late rebrightening in the afterglow of GRB 120326A

GRB 120326A is an unusual gamma-ray burst (GRB) which has a quite long plateau and a very late rebrightening both in X-ray and optical bands. The similar behavior of the optical and X-ray light curves suggests that they maybe have a common origin. The long plateau starts from several hundred seconds and ends at tens of thousands seconds. The peak time of the late rebrightening is about 30000 s. We analyze the energy injection model by means of numerical and analytical solutions, considering both the wind environment and ISM environment for GRB afterglows. We especially study the influence of the injection starting time, ending time, stellar wind density (or density of the circumburst environment), and injection luminosity on the shape of the afterglow light curves, respectively. We find that the light curve is largely affected by the parameters in the wind model. There is a "bump" at the late time only in the wind model too. In the wind case, it is interesting that the longer the energy injected, the more obvious the rebrightening will be. We also find the peak time of bump is determined by the stellar wind density. We use the late continuous injection model to interpret the unusual afterglow of GRB 120326A. The model can well fit the observational data, however, we find that the time scale of the injection must be larger than ten thousands seconds. This implies that the time scale of the central engine activity must be more than ten thousands seconds. This can give useful constraints on the central engine of GRBs. We consider a new born millisecond pulsar with strong magnetic field as the central engine. On the other hand, our results suggest that the circumburst environment of GRB 120326A is very likely a stellar wind.

preprint2014arXiv

X-Ray Afterglow Plateaus of Long Gamma-Ray Bursts: Further Evidence for Millisecond Magnetars

Many long-duration gamma-ray bursts (GRBs) were observed by {\it Swift}/XRT to have plateaus in their X-ray afterglow light curves. This plateau phase has been argued to be evidence for long-lasting activity of magnetar (ultra-strongly magnetized neutron stars) central engines. However, the emission efficiency of such magnetars in X-rays is still unknown. Here we collect 24 long GRB X-ray afterglows showing plateaus followed by steep decays. We extend the well-known relationship between the X-ray luminosity ${L_{\mathrm{X}}}$ and spin-down luminosity $L_{\mathrm{sd}}$ of pulsars to magnetar central engines, and find that the initial rotation period $P_{0}$ ranges from 1 ms to 10 ms and that the dipole magnetic field $B$ is centered around $10^{15}$ G. These constraints not only favor the suggestion that the central engines of some long GRBs are very likely to be rapidly rotating magnetars but also indicate that the magnetar plateau emission efficiency in X-rays is close to 100%.

preprint2012arXiv

Panchromatic observations of the textbook GRB 110205A: constraining physical mechanisms of prompt emission and afterglow

We present a comprehensive analysis of a bright, long duration (T90 ~ 257 s) GRB 110205A at redshift z= 2.22. The optical prompt emission was detected by Swift/UVOT, ROTSE-IIIb and BOOTES telescopes when the GRB was still radiating in the gamma-ray band. Nearly 200 s of observations were obtained simultaneously from optical, X-ray to gamma-ray, which makes it one of the exceptional cases to study the broadband spectral energy distribution across 6 orders of magnitude in energy during the prompt emission phase. By fitting the time resolved prompt spectra, we clearly identify, for the first time, an interesting two-break energy spectrum, roughly consistent with the standard GRB synchrotron emission model in the fast cooling regime. Although the prompt optical emission is brighter than the extrapolation of the best fit X/gamma-ray spectra, it traces the gamma-ray light curve shape, suggesting a relation to the prompt high energy emission. The synchrotron + SSC scenario is disfavored by the data, but the models invoking a pair of internal shocks or having two emission regions can interpret the data well. Shortly after prompt emission (~ 1100 s), a bright (R = 14.0) optical emission hump with very steep rise (alpha ~ 5.5) was observed which we interpret as the emission from the reverse shock. It is the first time that the rising phase of a reverse shock component has been closely observed. The full optical and X-ray afterglow lightcurves can be interpreted within the standard reverse shock (RS) + forward shock (FS) model. In general, the high quality prompt emission and afterglow data allow us to apply the standard fireball shock model to extract valuable information about the GRB including the radiation mechanism, radius of prompt emission R, initial Lorentz factor of the outflow, the composition of the ejecta, as well as the collimation angle and the total energy budget.

preprint2011arXiv

A Photometric Redshift of z ~ 9.4 for GRB 090429B

Gamma-ray bursts (GRBs) serve as powerful probes of the early Universe, with their luminous afterglows revealing the locations and physical properties of star forming galaxies at the highest redshifts, and potentially locating first generation (Population III) stars. Since GRB afterglows have intrinsically very simple spectra, they allow robust redshifts from low signal to noise spectroscopy, or photometry. Here we present a photometric redshift of z~9.4 for the Swift-detected GRB 090429B based on deep observations with Gemini-North, the Very Large Telescope, and the GRB Optical and Near-infrared Detector. Assuming a Small Magellanic Cloud dust law (which has been found in a majority of GRB sight-lines), the 90% likelihood range for the redshift is 9.06 < z < 9.52, although there is a low-probability tail to somewhat lower redshifts. Adopting Milky Way or Large Magellanic Cloud dust laws leads to very similar conclusions, while a Maiolino law does allow somewhat lower redshift solutions, but in all cases the most likely redshift is found to be z>7. The non-detection of the host galaxy to deep limits (Y_AB >~ 28 mag, which would correspond roughly to 0.001 L* at z=1) in our late time optical and infrared observations with the Hubble Space Telescope strongly supports the extreme redshift origin of GRB 090429B, since we would expect to have detected any low-z galaxy, even if it were highly dusty. Finally, the energetics of GRB 090429B are comparable to those of other GRBs, and suggest that its progenitor is not greatly different to those of lower redshift bursts.

preprint2010arXiv

The afterglow and host galaxy of GRB 090205: evidence for a Ly-alpha emitter at z=4.65

Gamma-ray bursts have been proved to be detectable up to distances much larger than any other astrophysical object, providing the most effective way, complementary to ordinary surveys, to study the high redshift universe. To this end, we present here the results of an observational campaign devoted to the study of the high-z GRB 090205. We carried out optical/NIR spectroscopy and imaging of GRB 090205 with the ESO-VLT starting from hours after the event up to several days later to detect the host galaxy. We compared the results obtained from our optical/NIR observations with the available Swift high-energy data of this burst. Our observational campaign led to the detection of the optical afterglow and host galaxy of GRB 090205 and to the first measure of its redshift, z=4.65. Similar to other, recent high-z GRBs, GRB 090205 has a short duration in the rest-frame with T_{90,rf}=1.6 s, which suggests the possibility that it might belong to the short GRBs class. The X-ray afterglow of GRB 090205 shows a complex and interesting behaviour with a possible rebrightening at 500-1000s from the trigger time and late flaring activity. Photometric observations of the GRB 090205 host galaxy argue in favor of a starburst galaxy with a stellar population younger than ~ 150 Myr. Moreover, the metallicity of Z > 0.27 Z_Sun derived from the GRB afterglow spectrum is among the highest derived from GRB afterglow measurement at high-z, suggesting that the burst occurred in a rather enriched envirorment. Finally, a detailed analysis of the afterglow spectrum shows the existence of a line corresponding to Lyman-alpha emission at the redshift of the burst. GRB 090205 is thus hosted in a typical Lyman-alpha emitter (LAE) at z=4.65. This makes the GRB 090205 host the farthest GRB host galaxy, spectroscopically confirmed, detected to date.

preprint2008arXiv

GRB 060206: hints of precession of the central engine?

Aims. The high-redshift (z = 4:048) gamma-ray burst GRB 060206 showed unusual behavior, with a significant rebrightening by a factor of ~ 4 at about 3000 s after the burst. We argue that this rebrightening implies that the central engine became active again after the main burst produced by the first ejecta, then drove another more collimated jet-like ejecta with a larger viewing angle. The two ejecta both interacted with the ambient medium, giving rise to forward shocks that propagated into the ambient medium and reverse shocks that penetrated into the ejecta. The total emission was a combination of the emissions from the reverse- and forward- shocked regions. We discuss how this combined emission accounts for the observed rebrightening. Methods. We apply numerical models to calculate the light curves from the shocked regions, which include a forward shock originating in the first ejecta and a forward-reverse shock for the second ejecta. Results. We find evidence that the central engine became active again 2000 s after the main burst. The combined emission produced by interactions of these two ejecta with the ambient medium can describe the properties of the afterglow of this burst. We argue that the rapid rise in brightness at ~ 3000 s in the afterglow is due to the off-axis emission from the second ejecta. The precession of the torus or accretion disk of the central engine is a natural explanation for the departure of the second ejecta from the line of sight.

preprint2008arXiv

GRB 080913 at redshift 6.7

We report on the detection by Swift of GRB 080913, and subsequent optical/near-infrared follow-up observations by GROND which led to the discovery of its optical/NIR afterglow and the recognition of its high-z nature via the detection of a spectral break between the i' and z' bands. Spectroscopy obtained at the ESO-VLT revealed a continuum extending down to lambda = 9400 A, and zero flux for 7500 A < lambda<9400 A, which we interpret as the onset of a Gunn-Peterson trough at z=6.695+-0.025 (95.5% conf. level), making GRB 080913 the highest redshift GRB to date, and more distant than the highest-redshift QSO. We note that many redshift indicators which are based on promptly available burst or afterglow properties have failed for GRB 080913. We report on our follow-up campaign and compare the properties of GRB 080913 with bursts at lower redshift. In particular, since the afterglow of this burst is fainter than typical for GRBs, we show that 2m-class telescopes can identify most high-redshift GRBs.

preprint2004arXiv

Afterglow Light Curves of Jetted Gamma-ray Burst Ejecta in Stellar Winds

Optical and radio afterglows arising from the shocks by relativistic conical ejecta running into pre-burst massive stellar winds are revisited. Under the homogeneous thin-shell approximation and a realistic treatment for the lateral expansion of jets, our results show that a notable break exists in the optical light curve in most cases we calculated in which the physical parameters are varied within reasonable ranges. For a relatively tenuous wind which cannot decelerate the relativistic jet to cause a light curve break within days, the wind termination shock due to the ram pressure of the surrounding medium occurs at a small radius, namely, a few times 10^{17} cm. In such a structured wind environment, the jet will pass through the wind within several hours and run into the outer uniform dense medium. The resulting optical light curve flattens with a shallower drop after the jet encounters the uniform medium, and then declines deeply, triggered by runaway lateral expansion.