Source author record

Qiang Yuan

Qiang Yuan appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.HE hep-ph astro-ph.CO astro-ph.GA astro-ph.IM astro-ph gr-qc hep-ex astro-ph.SR hep-th physics.ins-det

Catalog footprint

What is connected

81works

11topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2026arXiv

Geomagnetic constraints on Millicharged Dark Matter

Millicharged particles are well-motivated dark matter candidates arising in many extensions of the Standard Model. We show that, despite their tiny coupling $e_m$ to photons, millicharged dark matter (mDM) in the Earth's geomagnetic field can generate a quasi-static, monochromatic magnetic signal with angular frequency twice the mDM mass. Using null results from the SuperMAG and SNIPE Hunt collaborations, we constrain the effective charge of bosonic mDM in the mass range $10^{-18}$--$10^{-14}\,\text{eV}$. The resulting upper bounds exceed stellar cooling constraints by over thirteen orders of magnitude, demonstrating the power of this method.

preprint2025arXiv

Ultrahigh-Energy Gamma-ray Emission Associated with Black Hole-Jet Systems

Black holes (BH), one of the most intriguing objects in the universe, can manifest themselves through electromagnetic radiation initiated by the accretion flow. Some stellar-mass BHs drive relativistic jets when accreting matter from their companion stars, forming microquasars. Non-thermal emission from the radio to tera-electronvolt (TeV) gamma-ray band has been observed from microquasars, indicating the acceleration of relativistic particles. Here we report detection of four microquasars (SS 433, V4641 Sgr, GRS 1915+105, MAXI J1820+070) of spectrum extending to the ultrahigh-energy (UHE; photon energy $E>100$ TeV) band and one microquasar (Cygnus X-1) of spectrum approaching 100 TeV, using the Large High Altitude Air Shower Observatory (LHAASO). Notably, the total emission associated with SS 433 cannot be interpreted with a single leptonic component. In the UHE band, its emission is in spatial coincidence with a giant atomic cloud, which is consistent with a hadronic origin. An elongated source is discovered from V4641 Sgr with the spectrum continuing up to 800 TeV. The detection of UHE gamma rays demonstrates that accreting BHs and their environments can operate as extremely efficient accelerators of particles out of 1 peta-electronvolt (PeV), suggesting microquasars to be important contributors to Galactic cosmic rays especially around the `knee' region.

preprint2022arXiv

Constraints on dark matter annihilation from the Event Horizon Telescope Observations of M87$^\star$

The fast developments of radio astronomy open a new window to explore the properties of Dark Matter (DM). The recent direct imaging of the supermassive black hole (SMBH) at the center of M87 radio galaxy by the Event Horizon Telescope (EHT) collaboration is expected to be very useful to search for possible new physics. In this work, we illustrate that such results can be used to detect the possible synchrotron radiation signature produced by DM annihilation from the innermost region of the SMBH. Assuming the existence of a spike DM density profile, we obtain the flux density due to DM annihilation induced electrons and positrons, and derive new limits on the DM annihilation cross section via the comparison with the EHT integral flux density at 230 GHz. Our results show that the parameter space can be probed by the EHT observations is largely complementary to other experiments. For DM with typical mass regions of being weakly interacting massive particles, the annihilation cross section several orders of magnitude below the thermal production level can be excluded by the EHT observations under the density spike assumption. Future EHT observations may further improve the sensitivity on the DM searches, and may also provide a unique opportunity to test the interplay between DM and the SMBH.

preprint2022arXiv

Divergence in mass ratio distributions between low-mass and high-mass coalescing binary black holes

Coalescing binary black hole (BBH) systems are likely formed via several channels, and it is challenging to understand their formation / evolutionary processes. Some features in the mass function of the primary components ($m_1$), such as the distinct Gaussian-like peak located at $\sim 34M_\odot$, have been previously found. In this work, we investigate the possible dependence of the mass ratio ($q=m_2/m_1$) distribution on the primary mass. We find a Bayesian odds ratio of 18.1 in favor of divergence in the mass ratio distributions between the low- and high-mass ranges over an invariable mass ratio distribution. The BBHs with $m_1\gtrsim29M_{\odot}$ have a stronger preference to be symmetric compared to those with $m_1\lesssim29M_{\odot}$ at a 97.6\% credible level. Additionally, we find mild evidence that the BBHs with $m_1$ located in the Gaussian-like peak have a mass ratio distribution different from that of other BBHs. Our findings may be in favor of some formation channels, such as the chemically homogeneous evolution and the dynamical assembly in globular clusters/nuclear star clusters, which are more likely to provide symmetric BBHs in the high-mass range.

preprint2022arXiv

Explanation of nearby SNRs for primary electron excess and proton spectral bump

Several groups have reported a possible excess of primary electrons at high energies with the joint fit of the positron fraction and total electron/positron spectra. With the latest release of high-precision electron/positron spectra measured by AMS-02, we further confirm this excess by fitting $ΔΦ$ $\rm(i.e., Φ_{e^-}-Φ_{e^+})$ data in this work. Then we investigate the contribution of a single nearby supernova remnant to the primary electron excess and find that Monogem can reasonably account for this excess. Moreover, we predict that the electron spectrum may harden again at a few TeVs due to Vela's contribution. DAMPE, which can accurately measure electrons at TeV scale, is expected to provide the robust test of this new spectral feature in the near future. Finally, we fit the proton spectrum data of DAMPE with Monogem or Loop I. We find that both the primary electron excess and the proton spectral bump could be mainly generated by Monogem.

preprint2022arXiv

Inverting cosmic ray propagation by Convolutional Neural Networks

We propose a machine learning method to investigate the propagation of cosmic rays based on the precisely measured spectra of the primary and secondary cosmic ray nuclei of Li, Be, B, C, and O from AMS-02, ACE, and Voyager-1. We train two convolutional neural networks. One network learns how to infer propagation and source parameters from the energy spectra of cosmic rays, and the other network, which is similar to the former, has the flexibility to learn from the data with added artificial fluctuations. Together with the simulated data generated by GALPROP, we find that both networks can properly invert the propagation process and infer the propagation and source parameters reasonably well. This approach can be much more efficient than the traditional Markov chain Monte Carlo fitting method for deriving the propagation parameters if users choose to update confidence intervals with new experimental data. Both of the trained networks are available at (https://github.com/alan200276/CR_ML).

preprint2022arXiv

Stringent axion constraints with Event Horizon Telescope polarimetric measurements of M87$^\star$

The hitherto unprecedented angular resolution of the Event Horizon Telescope (EHT) has created exciting opportunities in the search for new physics. Recently, the linear polarization of radiation emitted near the supermassive black hole M87$^\star$ was measured on four separate days, precisely enabling tests of the existence of a dense axion cloud produced by a spinning black hole. The presence of an axion cloud leads to a frequency-independent oscillation in the electric vector position angle (EVPA) of this linear polarization. For a nearly face-on M87$^\star$, this oscillation in the EVPA appears as a propagating wave along the photon ring. In this paper, we leverage the azimuthal distribution of EVPA measured by the EHT to study the axion-photon coupling. We propose a novel differential analysis procedure to reduce the astrophysical background, and derive stringent constraints on the existence of axions in the previously unexplored mass window $\sim (10^{-21}-10^{-20})$~eV.

preprint2022arXiv

The Large High Altitude Air Shower Observatory (LHAASO) Science Book (2021 Edition)

Since the science white paper of the Large High Altitude Air Shower Observatory (LHAASO) published on arXiv in 2019 [e-Print: 1905.02773 (astro-ph.HE)], LHAASO has completed the transition from a project to an operational gamma-ray astronomical observatory LHAASO is a new generation multi-component facility located in Daocheng, Sichuan province of China, at an altitude of 4410 meters. It aims at measuring with unprecedented sensitivity the spectrum, composition, and anisotropy of cosmic rays in the energy range between 10$^{12}$ and 10$^{18}$~eV, and acting simultaneously as a wide aperture (one stereoradiant) continuously operating gamma-ray telescope in the energy range between 10$^{11}$ and $10^{15}$~eV with the designed sensitivity of 1.3\% of the Crab Unit (CU) above 100 TeV. LHAASO's capability of measuring simultaneously different shower components (electrons, muons, and Cherenkov/fluorescence light), will allow it to investigate the origin, acceleration, and propagation of CR through measurement of the energy spectrum, elemental composition, and anisotropy with unprecedented resolution. The remarkable sensitivity of LHAASO will play a key role in CR physics and gamma-ray astronomy for a general and comprehensive exploration of the high energy universe and will allow important studies of fundamental physics (such as indirect dark matter search, Lorentz invariance violation, quantum gravity) and solar and heliospheric physics. The LHAASO Collaboration organized an editorial working group and finished all editorial work of this science book, to summarize the instrumental features and outline the prospects of scientific researches with the LHAASO experiment.

preprint2022arXiv

Ultrahigh-energy diffuse gamma-ray emission from cosmic-ray interactions with medium surrounding acceleration sources

The diffuse $γ$-ray spectrum at sub-PeV energy region has been measured for the first time by the Tibet-AS$γ$ experiment. It will shed new light on the understanding of origin and propagation of Galactic cosmic rays at very high energies. It has been pointed out that the traditional cosmic ray propagation model based on low energy measurements undershoot the new data, and modifications of the model with new ingredients or alternative propagation framework is required. In this work, we propose that the hadronic interactions between freshly accelerated cosmic rays and the medium surrounding the sources, which was neglected in the traditional model, can naturally account for the Tibet-AS$γ$ diffuse emission. We show that this scenario gives a consistent description of other secondary species such as the positron spectrum, the boron-to-carbon ratio, and the antiproton-to-proton ratio. As a result, the electron spectrum above 10 TeV will have a hardening due to this secondary component, which may be tested by future measurements.

preprint2022arXiv

Understanding the phase reversals of Galactic cosmic ray anisotropies

The energy spectra and anisotropies are very important probes of the origin of cosmic rays. Recent measurements show that complicated but very interesting structures exist, at similar energies, in both the spectra and energy-dependent anisotropies, indicating a common origin of these structures. Particularly interesting phenomenon is that there is a reversal of the phase of the dipole anisotropies, which challenges a theoretical modeling. In this work, for the first time, we identify that there might be an additional phase reversal at $\sim 100$ GeV energies of the dipole anisotropies as indicated by a few underground muon detectors and the first direct measurement by the Fermi satellite, coincident with the hundreds of GV hardenings of the spectra. We propose that these two phase reversals, together with the energy-evolution of the amplitudes and spectra, can be naturally explained with a nearby source overlapping onto the diffuse background. As a consequence, the spectra and anisotropies can be understood as the scalar and vector components of this model, and the two reversals of the phases characterize just the competition of the cosmic ray streamings between the nearby source and the background. The alignment of the cosmic ray streamings along the local large-scale magnetic field may play an important but sub-dominant role in regulating the cosmic ray propagation. More precise measurements of the anisotropy evolution at both low energies by space detectors and high energies by air shower experiments for individual species will be essential to further test this scenario.

preprint2021arXiv

A systematic study of $γ$-ray flares from the Crab Nebula with Fermi-LAT: I. flare detection

Significant flares of GeV $γ$-ray emission from the Crab Nebula have been found by AGILE and Fermi-LAT years ago, indicating that extreme particle acceleration and radiation occurs in young pulsar wind nebulae. To enlarge the flare sample and to investigate their statistical properties will be very useful in understanding the nature of the $γ$-ray flares. In this paper, we investigate the flaring emission from the Crab Nebula with eleven year observations of the Fermi-LAT. We identify 17 significant flares in the light curve of the low-energy (synchrotron) component of the $γ$-ray emission. The flare rate is about 1.5 per year, without any significant change or clustering during the 11 years of the observation. We detect a special flare with an extremely long duration of nearly one month, occurred in October, 2018, with synchrotron photons up to energies of about 1GeV. The synchrotron component could be fitted by a steady power-law background and a variable flare component with an expotentially cutoff power-law spectrum, not only for individual flare but also for the combined data, which may favor a similar emission mechanism for all flares. However, we do not find a universal relation between the cutoff energy and the energy fluxes of the flares, which may reflect the complicated acceleration and/or cooling processes of the involved particles.

preprint2021arXiv

Constraining $U(1)_{L_μ-L_τ}$ charged dark matter model for muon $g-2$ anomaly with AMS-02 electron and positron data

Very recently, the Fermi-Lab reported the new experimental combined results on the magnetic momentum of muon with a 4.2$σ$ discrepancy compared with the expectation of the Standard Model \cite{Fermi_Lab}. A new light gauge boson $X$ in the $L_μ-L_τ$ model provides a good explanation for the $g-2$ anomaly. A Dirac fermion dark matter with a large $L_μ-L_τ$ charge can explain both the $g-2$ anomaly and the dark matter relic density \cite{Asai_2021}. In this work, we focus on the case that the mass of the dark matter is larger than the mass of muon (i.e. $m_Ψ > m_μ$) for which the channel $ΨΨ\rightarrow μ^- μ^+$ opens. Although the cross section $(σv)_{μ^{-}μ^{+}}$ is smaller by a factor of $1/q_Ψ^2$ ($q_Ψ$ represents the $L_μ-L_τ$ charge of the dark matter) compared with the channel $ΨΨ\rightarrow XX \rightarrow νν\barν\barν$, the resulting secondary electrons and positrons could imprint on their spectra above GeV energies due to the reacceleration effect of cosmic ray propagation. We use the AMS-02 measurements of electrons and positrons to constrain the annihilation cross section of the channel $ΨΨ\rightarrow μ^{-}μ^{+}$, which rules out part of the parameter space of the large $L_μ-L_τ$ charged dark matter model to account for the muon $g-2$ anomaly.

preprint2021arXiv

Geminga SNR: Possible candidate of local cosmic-ray factory

The precise measurements of energy spectra and anisotropy could help us uncover the local cosmic-ray accelerators. Our recent works have shown that spectral hardening above $200$ GeV in the energy spectra and transition of large-scale anisotropy at $\sim 100$ TeV are of local source origin. Less than $100$ TeV, both spectral hardening and anisotropy explicitly indicate the dominant contribution from nearby sources. In this work, we further investigate the parameter space of sources allowed by the observational energy spectra and anisotropy amplitude. To obtain the best-fit source parameters, a numerical package to compute the parameter posterior distributions based on Bayesian inference, which is applied to perform an elaborate scan of parameter space. We find that by combining the energy spectra and anisotropy data, the permissible range of location and age of local source is considerably reduced. When comparing with the current local SNR catalog, only Geminga SNR could be the proper candidate of the local cosmic-ray source.

preprint2021arXiv

Klein-Nishina effect and the cosmic ray electron spectrum

Radiative energy losses are very important in regulating the cosmic ray electron and/or positron (CRE) spectrum during their propagation in the Milky Way. Particularly, the Klein-Nishina (KN) effect of the inverse Compton scattering (ICS) results in less efficient energy losses of high-energy electrons, which is expected to leave imprints on the propagated electron spectrum. It has been proposed that the hardening of CRE spectra around 50 GeV observed by Fermi-LAT, AMS-02, and DAMPE could be due to the KN effect. We show in this work that the transition from the Thomson regime to the KN regime of the ICS is actually quite smooth compared with the approximate treatment adopted in some previous works. As a result, the observed spectral hardening of CREs cannot be explained by the KN effect. It means that an additional hardening of the primary electrons spectrum is needed. We also provide a parameterized form for the accurate calculation of the ICS energy-loss rate in a wide energy range.

preprint2021arXiv

Lightcurve Evolution of the nearest Tidal Disruption Event: A late-time, radio-only flare

Tidal disruption events (TDEs) occur when a star passes close enough to a galaxy's supermassive black hole to be disrupted by tidal forces. We discuss new observations of IGRJ12580+0134, a TDE observed in NGC 4845 (d=17 Mpc) in November 2010, with the Karl G. Jansky Very Large Array (JVLA). We also discuss a reanalysis of 2010-2011 Swift and XMM-Newton observations, as well as new, late-time Swift observations. Our JVLA observations show a decay of the nuclear radio flux until 2015, when a plateau was seen, and then a significant (~factor 3) radio flare during 2016. The 2016 radio flare was also accompanied by radio spectral changes, but was not seen in the X-rays. We model the flare as resulting from the interaction of the nuclear jet with a cloud in the interstellar medium. This is distinct from late-time X-ray flares in a few other TDEs where changes in the accretion state and/or a fallback event were suggested, neither of which appears possible in this case. Our reanalysis of the Swift and XMM-Newton data from 2011 shows significant evidence for thermal emission from a disk, as well as a very soft power-law. This, in addition to the extreme X-ray flux increase seen in 2010 (a factor of >$100) bolsters the identification of IGRJ12580+0134 as a TDE, not an unusual AGN variability event.

preprint2020arXiv

A Comparative Study of Long and Short GRBs. II. A Multi-wavelength Method to distinguish Type II (massive star) and Type I (compact star) GRBs

Gamma Ray Burst (GRBs) are empirically classified as long-duration GRBs (LGRBs, $>$ 2s) and short-duration GRBs (SGRBs, $<$ 2s). Physically they can be grouped into two distinct progenitor categories: those originating from collapse of massive stars (also known as Type II) and those related to mergers of compact stars (also known as Type I). Even though most LGRBs are Type II and most SGRBs are Type I, the duration criterion is not always reliable to determine the physical category of a certain GRB. Based on our previous comprehensive study of the multi-wavelength properties of long and short GRBs, here we utilize the Naive Bayes method to physically classify GRBs as Type I and Type II GRBs based on multi-wavelength criteria. It results in 0.5\% training error rate and 1\% test error rate. Moreover, there is a gap [-1.2, -0.16] in the distribution of the posterior Odds, $\log O({\rm II:I})$, the Type II to Type I probability ratio. Therefore, we propose to use ${\cal O} = \log O({\rm II:I})+0.7$ as the parameter to classify GRBs into Type I ($<0$) or Type II ($>0$). The only confirmed Type I GRB, GRB 170817A, has log $O({\rm II:I})=-10$. According to this criterion, the supernova-less long GRBs 060614 and 060505 belong to Type I, and two controversial short GRBs 090426 and 060121 belong to Type II.

preprint2020arXiv

A serendipitous discovery of GeV gamma-ray emission from supernova 2004dj in a survey of nearby star-forming galaxies with Fermi-LAT

The interaction between a supernova ejecta and the circum-stellar medium drives a strong shock wave which accelerates particles (i.e., electrons and protons). The radio and X-ray emission observed after the supernova explosion constitutes the evidence of the electron acceleration. The accelerated protons are expected to produce GeV-TeV gamma-ray emission via $pp$ collisions, but the flux is usually low since only a small fraction of the supernova kinetic energy is converted into the shock energy at the very early time. The low gamma-ray flux of the nearest supernova explosion, SN 1987A, agrees with this picture. Here we report a serendipitous discovery of a fading GeV gamma-ray source in spatial coincidence with the second nearest supernova--SN 2004dj from our gamma-ray survey of nearby star-forming galaxies with Fermi-LAT. The total gamma-ray energy released by SN 2004dj is about $6\times10^{47}{\rm erg}$. We interpret this gamma-ray emission arising from the supernova ejecta interacting with a surrounding high-density shell, which decelerates the ejecta and converts ~1% of the ejecta's kinetic energy to relativistic protons. In addition, our gamma-ray survey of nearby star-forming galaxies discovers GeV emissions from two star-forming galaxies, i.e., Arp 299 and M33, for the first time.

preprint2020arXiv

Boresight Alignment of DArk Matter Particle Explorer

The DArk Matter Particle Explorer (DAMPE) can measure $γ$-rays in the energy range from a few GeV to about 10 TeV. The direction of each $γ$-ray is reconstructed with respect to the reference system of the DAMPE payload. In this paper, we adopt a maximum likelihood method and use the $γ$-ray data centered around several bright point-like sources to measure and correct the angular deviation from the real celestial coordinate system, the so called ``boresight alignment'' of the DAMPE payload. As a check, we also estimate the boresight alignment for some sets of simulation data with artificial orientation and obtain consistent results. The time-dependent boresight alignment analysis does not show evidence for significant variation of the parameters.

preprint2020arXiv

Correction Method for the Readout Saturation of the DAMPE Calorimeter

The DArk Matter Particle Explorer (DAMPE) is a space-borne high energy cosmic-ray and $γ$-ray detector which operates smoothly since the launch on December 17, 2015. The bismuth germanium oxide (BGO) calorimeter is one of the key sub-detectors of DAMPE used for energy measurement and electron proton identification. For events with total energy deposit higher than decades of TeV, the readouts of PMTs coupled on the BGO crystals would become saturated, which results in an underestimation of the energy measurement. Based on detailed simulations, we develop a correction method for the saturation effect according to the shower development topologies and energies measured by neighbouring BGO crystals. The verification with simulated and on-orbit events shows that this method can well reconstruct the energy deposit in the saturated BGO crystal.

preprint2020arXiv

Interpretation of the cosmic ray positron and electron excesses with an annihilating-decaying dark matter scenario

The precise measurements of energy spectra of cosmic ray positrons and/or electrons by recent experiments show clear excesses above 10 GeV. Moreover, a potential sharp spectral feature was suggested by the Dark Matter Particle Explorer (DAMPE) data. These results inspire quite a number of discussions on the connection with either the annihilation/decay of dark matter (DM) or the astrophysical origins. Here we discuss a DM scenario in which DM particles could annihilate and decay into standard model particle pairs simultaneously. In this model, the peak structure is due to the DM annihilation in a nearby subhalo and the broad positron/electron excesses are due to the decay of DM in the Milky Way. This model can reasonably explain the DAMPE and AMS-02 data of the total $e^+e^-$ spectra and the positron fraction, with model parameters being consistent with existing constraints. A simple realization of such a DM model is the spin-1 vector DM model.

preprint2020arXiv

Local interstellar spectra and solar modulation of cosmic ray electrons and positrons

Low energy cosmic rays are modulated by the solar activity when they propagation in the heliosphere, leading to ambiguities in understanding their acceleration at sources and propagation in the Milky Way. By means of the precise measurements of the $e^-$, $e^+$, $e^-+e^+$, and $e^+/(e^-+e^+)$ spectra by AMS-02 near the Earth, as well as the very low energy measurements of the $e^-+e^+$ fluxes by Voyager-1 far away from the Sun, we derive the local interstellar spectra (LIS) of $e^-$ and $e^+$ components individually. Our method is based on a non-parametric description of the LIS of $e^-$ and $e^+$ and a force-field solar modulation model. We then obtain the evolution of the solar modulation parameters based on the derived LIS and the monthly fluxes of cosmic ray $e^-$ and $e^+$ measured by AMS-02. {\bf To better fit the monthly data, additional renormalization factors for $e^-$ and $e^+$ have been multiplied to the modulated fluxes.} We find that the inferred solar modulation parameters of positrons are in good agreement with that of cosmic ray nuclei, and the time evolutions of the solar modulation parameters of electrons and positrons differ after the reversal of the heliosphere magnetic field polarity, which shows clearly the charge-sign dependent modulation effect.

preprint2020arXiv

Probing Axions with Event Horizon Telescope Polarimetric Measurements

With high spatial resolution, polarimetric imaging of a supermassive black hole, like M87$^\star$ or Sgr A$^\star$, by the Event Horizon Telescope can be used to probe the existence of ultralight bosonic particles, such as axions. Such particles can accumulate around a rotating black hole through the superradiance mechanism, forming an axion cloud. When linearly polarized photons are emitted from an accretion disk near the horizon, their position angles oscillate due to the birefringent effect when traveling through the axion background. In particular, the observations of supermassive black holes M87$^\star$ (Sgr A$^\star$) can probe the dimensionless axion-photon coupling $c = 2 πg_{a γ} f_a$ for axions with mass around $O(10^{-20})$~eV ($O( 10^{-17}$)~eV) and decay constant $f_a < O(10^{16})$ GeV, which is complimentary to other axion measurements.

preprint2020arXiv

Probing Flavor Structure of Cosmic Ray $e^\mp$ Spectrum and Implications for Dark Matter Indirect Searches

Measuring high energy cosmic ray electrons/positrons (CRE) provides important means for the dark matter (DM) indirect detection and for probing the nearby galactic sources. In this work, we perform a systematic analysis of the flavor structure of DM annihilations into charged leptons based on the cosmic ray CRE spectra measured by DAMPE, Fermi-LAT, AMS-02, and CALET experiments. We study the annihilations of possible TeV scale DM particles in a nearby subhalo, which is proposed to explain the possible peak-like structure of the DAMPE CRE data. We pay special attention to the possible non-resonant excess (besides the possible peak-like structure) and demonstrate that such non-resonant excess can mainly arise from the decay of muons produced by the DM annihilations in the subhalo. With these we study the flavor composition of the lepton final states from DM annihilations $χχ\to e^+e^-, μ^+μ^-, τ^+τ^-$ by fitting the CRE data. We demonstrate that decays of the final states $μ^+ μ^-$ and $τ^+ τ^-$ can provide the non-resonant excess, while the peak excess arises from the $e^+ e^-$ final state. We further analyze the constraints on the lepton flavor composition using the Fermi-LAT $γ$-ray measurements. We find that the flavor composition is consistent with the Fermi-LAT data at relatively low Galactic latitudes, while the fraction of the final state $τ^\pm$ is severely bounded.

preprint2019arXiv

Searching for the possible signal of the photon-axionlike particle oscillation in the combined GeV and TeV spectra of supernova remnants

The conversion between photons and axionlike particles (ALPs) in the Milky Way magnetic field could result in the detectable oscillation phenomena in $γ$-ray spectra of Galactic sources. In this work, the GeV (Fermi-LAT) and TeV (MAGIC/VERITAS/H.E.S.S.) data of three bright supernova remnants (SNRs, ie. IC443, W51C and W49B) have been adopted together to search such the oscillation effect. Different from our previous analysis of the sole Fermi-LAT data of IC443, we do not find any reliable signal for the photon-ALP oscillation in the joint broadband spectrum of each SNR. The reason for the inconsistence is that in this work we use the latest revision (P8R3) of Fermi-LAT data, updated diffuse emission templates and the new version of the source catalog (4FGL), which lead to some modification of the GeV spectrum of IC443. Then we set constraints on ALP parameters based on the combined analysis of all the three sources. Though these constraints are somewhat weaker than limits from the CAST experiment and globular clusters, they are supportive of and complementary to these other results.

preprint2016arXiv

A GeV source in the direction of Supernova Remnant CTB 37B

Supernova remnants (SNRs) are the most attractive candidates for the acceleration sites of Galactic cosmic rays. We report the detection of GeV $γ$-ray emission with the Pass 8 events recorded by Fermi Large Area Telescope (Fermi-LAT) in the vicinity of the shell type SNR CTB 37B that is likely associated with the TeV $γ-$ray source HESS J1713-381. The photon spectrum of CTB 37B is consistent with a power-law with an index of $1.89\pm0.08$ in the energy range of $0.5-500$ GeV, and the measured flux connects smoothly with that of HESS J1713-381 at a few hundred GeV. No significant spatial extension and time variation are detected. The multi-wavelength data can be well fitted with either a leptonic model or a hadronic one. However, parameters of both models suggest more efficient particle acceleration than typical SNRs. Meanwhile, the X-ray and $γ$-ray spectral properties of CTB 37B show that it is an interesting source bridging young SNRs dominated by non-thermal emission and old SNRs interacting with molecular clouds.

preprint2016arXiv

Catching jetted tidal disruption events early in millimetre

Relativistic jets can form from at least some tidal disruption events (TDEs) of (sub-)stellar objects around supermassive black holes. We detect the millimeter (MM) emission of IGR J12580+0134 --- the nearest TDE known in the galaxy NGC 4845 at the distance of only 17 Mpc, based on Planck all-sky survey data. The data show significant flux jumps after the event, followed by substantial declines, in all six high frequency Planck bands from 100 GHz to 857 GHz. We further show that the evolution of the MM flux densities are well consistent with our model prediction from an off-axis jet, as was initially suggested from radio and X-ray observations. This detection represents the second TDE with MM detections; the other is Sw J1644+57, an on-axis jetted TDE at redshift of 0.35. Using the on- and off-axis jet models developed for these two TDEs as templates, we estimate the detection potential of similar events with the Large Millimeter Telescope (LMT) and the Atacama Large Millimeter/submillimeter Array (ALMA). Assuming an exposure of one hour, we find that the LMT (ALMA) can detect jetted TDEs up to redshifts $z\sim1$ (2), for a typical disrupted star mass of $\sim1$ M$_\odot$. The detection rates of on- and off-axis TDEs can be as high as $\sim0.6$ (13) and 10 (220) per year, respectively, for the LMT (ALMA). We briefly discuss how such observations, together with follow-up radio monitoring, may lead to major advances in understanding the jetted TDEs themselves and the ambient environment of the CNM.

preprint2016arXiv

Constraints on the dark matter annihilation from Fermi-LAT observation of M31

Gamma-ray is a good probe of dark matter (DM) particles in the Universe. We search for the DM annihilation signals in the direction of the Andromeda galaxy (M31) using 7.5 year Fermi-LAT pass 8 data. Similar to Pshirkov et al. (2016), we find that there is residual excess emission from the direction of M31 if only the galactic disk as traced by the far infrared emission is considered. Adding a point-like source will improve the fitting effectively, although additional slight improvements can be found if an extended component such as a uniform disk or two bubbles is added instead. Taking the far infrared disk plus a point source as the background model, we search for the DM annihilation signals in the data. We find that there is strong degeneracy between the emission from the galaxy and that from 10s GeV mass DM annihilation in the main halo with quark final state. However, the required DM annihilation cross section is about $10^{-25}-10^{-24}$ cm$^3$s$^{-1}$, orders of magnitude larger than the constraints from observations of dwarf spheroidal galaxies, indicating a non-DM origin of the emission. If DM subhalos are taken into account, the degeneracy is broken. When considering the enhancement from DM subhalos, the constraints on DM model parameters are comparable to (or slightly weaker than) those from the population of dwarf spheroidal galaxies. We also discuss the inverse Compton scattering component from DM annihilation induced electrons/positrons. For the first time we include an energy dependent template of the inverse Compton emission (i.e., a template cube) in the data analysis to take into account the effect of diffusion of charged particles. We find a significant improvement of the constraints in the high mass range of DM particles after considering the inverse Compton emission.

preprint2016arXiv

Implications of the Tentative Association between GW150914 and a {\it Fermi}-GBM Transient

The merger-driven Gamma-ray Bursts (GRBs) and their associated gravitational wave (GW) radiation, if both successfully detected, have some far-reaching implications, including for instance: (i) The statistical comparison of the physical properties of the short/long-short GRBs with and without GW detection can test the general origin model; (ii) Revealing the physical processes taking place at the central engine; (iii) Measuring the velocity of the Gravitational wave directly/accurately. In this work we discuss these implications in the case of possible association of GW150914/ GBM transient 150914. We compared GBM transient 150914 with other SGRBs and found that such an event {may be} a distinct outlier in some statistical diagrams, possibly due to its specific binary-black-hole merger origin. However, the presence of a "new" group of SGRBs with "unusual" physical parameters is also possible. If the outflow of GBM transient 150914 was launched by the accretion onto the nascent black hole, the magnetic activity rather than the neutrino process is likely responsible for the energy extraction and the accretion disk mass is estimated to be $\sim 10^{-5}~M_\odot$. The GW150914/GBM transient 150914 association, {if confirmed, would} provide the first opportunity to directly measure the GW velocity and its departure from the speed of the light {should be within} a factor of $\sim 10^{-17}$.

preprint2016arXiv

Perspective of monochromatic gamma-ray line detection with the High Energy cosmic-Radiation Detection (HERD) facility onboard China's Space Station

HERD is the High Energy cosmic-Radiation Detection instrument proposed to operate onboard China's space station in the 2020s. It is designed to detect energetic cosmic ray nuclei, leptons and photons with a high energy resolution ($\sim1\%$ for electrons and photons and $20\%$ for nuclei) and a large geometry factor ($>3\,{ m^2\,sr}$ for electrons and diffuse photons and $>2\,{ m^2\,sr}$ for nuclei). In this work we discuss the capability of HERD to detect monochromatic $γ$-ray lines, based on simulations of the detector performance. It is shown that HERD will be one of the most sensitive instruments for monochromatic $γ$-ray searches at energies between $\sim10$ to a few hundred GeV. Above hundreds of GeV, Cherenkov telescopes will be more sensitive due to their large effective area. As a specific example, we show that a good portion of the parameter space of a supersymmetric dark matter model can be probed with HERD.

preprint2016arXiv

Search for gamma-ray line features from Milky Way satellites with Fermi-LAT Pass 8 data

With 91 months of the publicly available Fermi-LAT Pass 8 data, we analyze the gamma-ray emission from the Milky Way satellites to search for potential line signals due to the annihilation of dark matter particles into double photons. The searched targets include a sample of dwarf spheroidal galaxies, the Large Magellanic Cloud (LMC) and Small Magellanic Cloud (SMC). No significant line emission has been found neither in the stacked dwarf galaxy sample nor in the direction of LMC/SMC. The corresponding upper limits on the cross section of DM annihilation into two photons are derived. Compared with results of previous gamma-ray line searches with the Pass 7 data, the current constraints on the line emission from dwarf spheroidal galaxies has been significantly improved in a wide energy range. With the rapid increase of the sample of dwarf spheroidal galaxies (candidates), we expect that the sensitivity of gamma ray line searches will be significantly improved in the near future.

preprint2016arXiv

Testing the dark matter subhalo hypothesis of the gamma-ray source 3FGL J2212.5+0703

N-body simulations predict that galaxies at the Milky Way scale host a large number of dark matter (DM) subhalos. Some of these subhalos, if they are massive enough or close enough to the Earth, might be detectable in $γ$ rays due to the DM annihilation. 3FGL J2212.5+0703, an unidentified gamma-ray source, has been suggested to be the counterpart candidate of a DM subhalo by Bertoni et al. (2015, 2016). In this work we analyze the Fermi-LAT Pass 8 data of 3FGL J2212.5+0703 to independently test the DM subhalo hypothesis of this source. In order to suppress the possible contamination from two nearby very-bright blazars, we just take into account the front-converting gamma-rays which have better angular resolutions than that of the back-converting photons. In addition to the spatial distribution analysis, we have extended the spectrum analysis down to the energies of $\sim 100$ MeV, and thoroughly examined the variability of the emission during the past 8 years. We confirm that 3FGL J2212.5+0703 is a steady and spatially-extended gamma-ray emitter at a high confidence level. The spectrum is well consistent with that expected from DM annihilation into $b\bar{b}$. The introduction of a phenomenological LogParabola spectrum just improves the fit slightly. All these results suggest that 3FGL J2212.5+0703 could be indicative of a DM subhalo.

preprint2015arXiv

A Systematic Chandra study of Sgr A$^{\star}$: I. X-ray flare detection

Daily X-ray flaring represents an enigmatic phenomenon of Sgr A$^{\star}$ --- the supermassive black hole at the center of our Galaxy. We report initial results from a systematic X-ray study of this phenomenon, based on extensive {\it Chandra} observations obtained from 1999 to 2012, totaling about 4.5 Ms. We detect flares, using a combination of the maximum likelihood and Markov Chain Monte Carlo methods, which allow for a direct accounting for the pile-up effect in the modeling of the flare lightcurves and an optimal use of the data, as well as the measurements of flare parameters, including their uncertainties. A total of 82 flares are detected. About one third of them are relatively faint, which were not detected previously. The observation-to-observation variation of the quiescent emission has an average root-mean-square of $6\%-14\%$, including the Poisson statistical fluctuation of faint flares below our detection limits. We find no significant long-term variation in the quiescent emission and the flare rate over the 14 years. In particular, we see no evidence of changing quiescent emission and flare rate around the pericenter passage of the S2 star around 2002. We show clear evidence of a short-term clustering for the ACIS-S/HETG 0th-order flares on time scale of $20-70$ ks. We further conduct detailed simulations to characterize the detection incompleteness and bias, which is critical to a comprehensive follow-up statistical analysis of flare properties. These studies together will help to establish Sgr A$^{\star}$ as a unique laboratory to understand the astrophysics of prevailing low-luminosity black holes in the Universe.

preprint2015arXiv

Quantitative study of the AMS-02 electron/positron spectra: implications for the pulsar and dark matter properties

The AMS-02 has just published the unprecedentedly precise measurement of the cosmic electron and positron spectra. In this paper we try to give a quantitative study on the AMS-02 results by a global fitting to the electron and positron spectra, together with the updated positron fraction data. The Markov Chain Monte Carlo algorithm is adopted to do the fitting. The primary electron spectrum and the parameters for pulsars or dark matter which contribute extra positrons are determined simultaneously. We find that there is a hardening of the primary electron spectrum at $\sim 60$ GeV. With such a new feature at the background spectrum both the pulsars and dark matter can explain the AMS-02 results very well. The dark matter scenario shows a drop at positron fraction at $\sim 300$ GeV, however, suffers very strong constraints from Fermi $γ$-ray observations. The fitting results also suggest that the propagation model with convection may be more favored by the lepton data than the reacceleration model.

preprint2015arXiv

Statistics of X-ray flares of Sagittarius A*: evidence for solar-like self-organized criticality phenomenon

X-ray flares have routinely been observed from the supermassive black hole, Sagittarius A$^\star$ (Sgr A$^\star$), at our Galactic center. The nature of these flares remains largely unclear, despite of many theoretical models. In this paper, we study the statistical properties of the Sgr A$^\star$ X-ray flares, by fitting the count rate (CR) distribution and the structure function (SF) of the light curve with a Markov Chain Monte Carlo (MCMC) method. With the 3 million second \textit{Chandra} observations accumulated in the Sgr A$^\star$ X-ray Visionary Project, we construct the theoretical light curves through Monte Carlo simulations. We find that the $2-8$ keV X-ray light curve can be decomposed into a quiescent component with a constant count rate of $\sim6\times10^{-3}~$count s$^{-1}$ and a flare component with a power-law fluence distribution $dN/dE\propto E^{-α_{\rm E}}$ with $α_{\rm E}=1.65\pm0.17$. The duration-fluence correlation can also be modelled as a power-law $T\propto E^{α_{\rm ET}}$ with $α_{\rm ET} < 0.55$ ($95\%$ confidence). These statistical properties are consistent with the theoretical prediction of the self-organized criticality (SOC) system with the spatial dimension $S = 3$. We suggest that the X-ray flares represent plasmoid ejections driven by magnetic reconnection (similar to solar flares) in the accretion flow onto the black hole.

preprint2015arXiv

Systematic study of the uncertainties in fitting the cosmic positron data by AMS-02

The operation of AMS-02 opens a new era for the study of cosmic ray physics with unprecedentedly precise data which are comparable with the laboratory measurements. The high precision data allow a quantitative study on the cosmic ray physics and give strict constraints on the nature of cosmic ray sources. However, the intrinsic errors from the theoretical models to interpret the data become dominant over the errors in the data. In the present work we try to give a systematic study on the uncertainties of the models to explain the AMS-02 positron fraction data, which shows the cosmic ray $e^+e^-$ excesses together with the PAMELA and Fermi-LAT measurements. The excesses can be attributed to contributions from the extra $e^+e^-$ sources, such as pulsars or the dark matter annihilation. The possible systematic uncertainties of the theoretical models considered include the cosmic ray propagation, the treatment of the low energy data, the solar modulation, the $pp$ interaction models, the nuclei injection spectrum and so on. We find that in general a spectral hardening of the primary electron injection spectrum above $\sim50-100$ GeV is favored by the data. Furthermore, the present model uncertainties may lead to a factor of $\sim2$ enlargement in the determination of the parameter regions of the extra source, such as the dark matter mass, annihilation rate and so on.

preprint2015arXiv

Testing the millisecond pulsar scenario of the Galactic center gamma-ray excess with very high energy gamma-rays

The recent analyses of the Fermi Large Area Telescope data show an extended GeV $γ$-ray excess on top of the expected diffuse background in the Galactic center region, which can be explained with annihilating dark matter or a population of millisecond pulsars (MSPs). We propose to observe the very high energy $γ$-rays for distinguishing the MSP scenario from the dark matter scenario. The GeV $γ$-ray MSPs should release most energy to the relativistic $e^{\pm}$ wind, which will diffuse in the Galaxy and radiate TeV $γ$-rays through inverse Compton scattering and bremsstrahlung processes. By calculating the spectrum and spatial distribution, we show that such emission is detectable with the next generation very high energy $γ$-ray observatory, the Cherenkov Telescope Array (CTA), under reasonable model parameters. It is essential to search for the multi-wavelength counterparts to the GeV $γ$-ray excess for solving this mystery in the high energy universe.

preprint2014arXiv

An Updated Analysis of Inert Higgs Doublet Model in light of the Recent Results from LUX, PLANCK, AMS-02 and LHC

In light of the recent discovery by the ATLAS and CMS experiments at the Large Hadron Collider (LHC) of a Higgs-like particle with a narrow mass range of 125-126 GeV, we perform an updated analysis on one of the popular scalar dark matter models, the Inert Higgs Doublet Model (IHDM). We take into account in our likelihood analysis of various experimental constraints, including recent relic density measurement, dark matter direct and indirect detection constraints as well as the latest collider constraints on the invisible decay width of the Higgs boson and monojet search at the LHC. It is shown that if the invisible decay of the standard model Higgs boson is open, LHC as well as direct detection experiments like LUX and XENON100 could put stringent limits on the Higgs boson couplings to dark matter. We find that the most favoured parameter space for IHDM corresponds to dark matter with a mass less than 100 GeV or so. In particular, the best-fit points are at the dark matter mass around 70 GeV where the invisible Higgs decay to dark matter is closed. Scalar dark matter in the higher mass range of 0.5-4 TeV is also explored in our study. Projected sensitivities for the future experiments of monojet at LHC-14, XENON1T and AMS-02 one year antiproton flux are shown to put further constraints on the existing parameter space of IHDM.

preprint2014arXiv

Fermi Large Area Telescope Detection of Supernova Remnant RCW 86

Using $5.4$ year Fermi-LAT data, we report the detection of GeV $γ$-ray emission from the shell-type supernova remnant RCW 86 (G315.4-2.3) with a significance of $\sim5.1σ$. The data slightly favors an extended emission of this supernova remnant. The spectral index of RCW 86 is found to be very hard, $Γ\sim1.4$, in the $0.4$ to $300$ GeV range. A one zone leptonic model can well fit the multi-wavelength data from radio to very high energy $γ$-rays. The very hard GeV $γ$-ray spectrum and the inferred low gas density seem to disfavor the hadronic origin of the $γ$-rays. The $γ$-ray behavior of RCW 86 is very similar to several other TeV shell-type supernova remnants, e.g., RX J1713.7-3946, RX J0852.0-4622, SN 1006 and HESS J1731-347.

preprint2014arXiv

Fermi Large Area Telescope observations of the supernova remnant HESS J1731-347

Context: HESS J1731-347 has been identified as one of the few TeV-bright shell-type supernova remnants (SNRs). These remnants are dominated by nonthermal emission, and the nature of TeV emission has been continuously debated for nearly a decade. Aims: We carry out the detailed modeling of the radio to gamma-ray spectrum of HESS J1731-347 to constrain the magnetic field and energetic particles sources, which we compare with those of the other TeV-bright shell-type SNRs explored before. Methods: Four years of data from Fermi Large Area Telescope (LAT) observations for regions around this remnant are analyzed, leading to no detection correlated with the source discovered in the TeV band. The Markov Chain Monte Carlo method is used to constrain parameters of one-zone models for the overall emission spectrum. Results: Based on the 99.9% upper limits of fluxes in the GeV range, one-zone hadronic models with an energetic proton spectral slope greater than 1.8 can be ruled out, which favors a leptonic origin for the gamma-ray emission, making this remnant a sibling of the brightest TeV SNR RX J1713.7-3946, the Vela Junior SNR RX J0852.0-4622, and RCW 86. The best-fit leptonic model has an electron spectral slope of 1.8 and a magnetic field of about 30 muG, which is at least a factor of 2 higher than those of RX J1713.7-3946 and RX J0852.0-4622, posing a challenge to the distance estimate and/or the energy equipartition between energetic electrons and the magnetic field of this source. A measurement of the shock speed will address this challenge and has implications on the magnetic field evolution and electron acceleration driven by shocks of SNRs.

preprint2014arXiv

Hybrid model of GeV-TeV gamma ray emission from Galactic Center

The observations of high energy $γ$-ray emission from the Galactic center (GC) by HESS, and recently by Fermi, suggest the cosmic ray acceleration in the GC and possibly around the supermassive black hole. In this work we propose a lepton-hadron hybrid model to explain simultaneously the GeV-TeV $γ$-ray emission. Both electrons and hadronic cosmic rays were accelerated during the past activity of the GC. Then these particles would diffuse outwards and interact with the interstellar gas and background radiation field. The collisions between hadronic cosmic rays with gas is responsible to the TeV $γ$-ray emission detected by HESS. With fast cooling in the strong radiation field, the electrons would cool down and radiate GeV photons through inverse Compton scattering off the soft background photons. This scenario provides a natural explanation of the observed GeV-TeV spectral shape of $γ$-rays.

preprint2014arXiv

Implications of the AMS-02 positron fraction in cosmic rays

The AMS-02 collaboration has just released its first result of the cosmic positron fraction $e^+/(e^-+e^+)$ with high precision up to $\sim 350$ GeV. The AMS-02 result shows the same trend with the previous PAMELA result, which requires extra electron/positron sources on top of the conventional cosmic ray background, either from astrophysical sources or from dark matter annihilation/decay. In this paper we try to figure out the nature of the extra sources by fitting to the AMS-02 $e^+/(e^-+e^+)$ data, as well as the electron and proton spectra by PAMELA and the $(e^-+e^+)$ spectrum by Fermi and HESS. We adopt the GALPROP package to calculate the propagation of the Galactic cosmic rays and the Markov Chain Monte Carlo sampler to do the fit. We find that the AMS-02 data have implied essential difference from the PAMELA data. There is {\rm tension} between the AMS-02 $e^+/(e^-+e^+)$ data and the Fermi/HESS $(e^-+e^+)$ spectrum, that the AMS-02 data requires less contribution from the extra sources than Fermi/HESS. Then we redo the fit without including the Fermi/HESS data. In this case both the pulsars and dark matter annihilation/decay can explain the AMS-02 data. The pulsar scenario has a soft inject spectrum with the power-law index $\sim 2$, while the dark matter scenario needs $τ^+τ^-$ final state with mass $\sim 600$ GeV and a boost factor $\sim 200$.

preprint2014arXiv

Millisecond pulsar interpretation of the Galactic center gamma-ray excess

It was found in the Fermi-LAT data that there is an extended $γ$-ray excess in the Galactic center region. The proposed sources to be responsible for the excess include the dark matter annihilation or an astrophysical alternative from a population of millisecond pulsars (MSPs). Whether or not the MSP scenario can explain the data self-consistently has very important implications for the detection of particle dark matter, which is however, subject to debate in the literature. In this work we study the MSP scenario in detail, based on the detected properties of the MSPs by Fermi-LAT. We build a model of the Milky Way MSPs which can reproduce the $γ$-ray properties of the Fermi-LAT MSPs, and derive the intrinsic luminosity function of the MSPs. The model is then applied to a bulge population of MSPs. We find that the extended $γ$-ray excess can be well explained by the bulge MSPs without violating the detectable flux distribution of MSPs by Fermi-LAT. The spatial distribution of the bulge MSPs as implied by the distribution of low mass X-ray binaries follows a $r^{-2.4}$ profile, which is also consistent with the $γ$-ray excess data. We conclude that the MSP model can explain the Galactic center $γ$-ray excess self-consistently, satisfying all the current observational constraints.

preprint2014arXiv

Multi-wavelength study of the MGRO J2019+37

MGRO J2019+37 within the Cygnus region is a bright and extended source revealed by Milagro at 12-35 TeV. This source is almost as bright as Crab Nebula in northern sky, while it is not confirmed by ARGO-YBJ around TeV. Up to now, no obvious counterpart at low energy wavelengths has been found. Hence, MGRO J2019+37 becomes mysterious and its VHE gamma-ray emission mechanism is attractive. In this paper, a brief summary of the multi-wavelength observations from Radio to gamma-ray is presented. All the available data from the XMM-Newton and INTEGRAL at X-ray, and the Fermi-LAT at gamma-ray bands were used to get constraint on its emission flux at low energy wavelengths. Then, its possible counterparts and the VHE emission mechanism are discussed.

preprint2014arXiv

Reconcile the AMS-02 positron fraction and Fermi-LAT/HESS total $e^{\pm}$ spectra by the primary electron spectrum hardening

The recently reported positron fraction up to $\sim 350$ GeV by AMS-02 seems to have tension with the total electron/positron spectra detected by Fermi and HESS, for either pulsar or dark matter annihilation/decay scenario as the primary positron sources. In this work we will show that the tension will be removed by an adjustment of the primary electron spectrum. If the primary electron spectrum becomes harder above $\sim50$ GeV, similar as the cosmic ray nuclei spectrum, the AMS-02 positron fraction and Fermi/HESS data can be well fitted by both the pulsar and dark matter models. This result indicates that there should be a common origin of the cosmic ray nuclei and the primary electrons. Furthermore, this study also implies that the properties of the extra sources derived from the fitting to the AMS-02 data should depend on the form of background.

preprint2014arXiv

Status of Dark Matter Detection

The detection of dark matter has made great progresses in recent years. We give a brief review on the status and progress in dark matter detection, including the progresses in direct detection, collider detection at LHC and focus on the indirect detection. The results from PAMELA, ATIC, Fermi-LAT and relevant studies on these results are introduced. Then we give the progress on indirect detection of gamma rays from Fermi-LAT and ground based Cerenkov telescopes. Finally the detection of neutrinos and constraints on the nature of dark matter are reviewed briefly.

preprint2013arXiv

A generic method to constrain the dark matter model parameters from Fermi observations of dwarf spheroids

Observation of gamma-rays from dwarf galaxies is an effective way to search for particle dark matter. Using 4-year data of Fermi-LAT observations on a series of Milky Way satellites, we develop a general way to search for the signals from dark matter annihilation in such objects. Instead of giving prior information about the energy spectrum of dark matter annihilation, we bin the Fermi-LAT data into several energy bins and build a likelihood map in the "energy bin - flux" plane. The final likelihood of any spectrum can be easily derived through combining the likelihood of all the energy bins. It gives consistent result with that directly calculated using the Fermi Scientific Tool. This method is very efficient for the study of any specific dark matter models with gamma-rays. We use the new likelihood map with Fermi-LAT 4 year data to fit the parameter space in three representative dark matter models: i) toy dark matter model, ii) effective dark matter operators, and iii) supersymmetric neutralino dark matter.

preprint2013arXiv

Emitting electrons spectra and acceleration processes in the jet of Mrk 421: from low state to giant flare state

We investigate the electron energy distributions (EEDs) and the acceleration processes in the jet of Mrk 421 through fitting the spectral energy distributions (SEDs) in different active states in the frame of a one-zone synchrotron self-Compton (SSC) model. After assuming two possible EEDs formed in different acceleration models: the shock accelerated power-law with exponential cut-off (PLC) EED and the stochastic turbulence accelerated log-parabolic (LP) EED, we fit the observed SEDs of Mrk 421 in both low and giant flare states by using the Markov Chain Monte Carlo (MCMC) method which constrains the model parameters in a more efficient way. Our calculating results indicate that (1) the PLC and LP models give comparably good fits for the SED in low state, but the variations of model parameters from low state to flaring can be reasonably explained only in the case of the PLC in low state; and (2) the LP model gives better fits compared to the PLC model for the SED in flare state, and the intra-day/night variability observed at GeV-TeV bands can be accommodated only in the LP model. The giant flare may be attributed to the stochastic turbulence re-acceleration of the shock accelerated electrons in low state. Therefore, we may conclude that shock acceleration is dominant in low state, while stochastic turbulence acceleration is dominant in flare state. Moreover, our result shows that the extrapolated TeV spectra from the best-fit SEDs from optical through GeV with the two EEDs are different. It should be considered in caution when such extrapolated TeV spectra are used to constrain extragalactic background light (EBL) models.

preprint2013arXiv

On the Galactic Center Being the Main Source of Galactic Cosmic Rays as Evidenced by Recent Cosmic Ray and Gamma Ray Observations

We revisit the idea that the Galactic center (GC) is the dominant source of Galactic cosmic rays (GCRs), based on a series of new observational evidence. A unified model is proposed to explain the new phenomena of GCRs and $γ$-rays simultaneously. The GCRs are thought to be accelerated during past activities of the GC. The pair production process of GCRs in the strong radiation field due to the GC activity is responsible for the knee structure of the cosmic ray spectra. A f raction of $e^+e^-$ produced by pair production interactions, can be reaccelerated in the induced bipolar jets and be transported into the halo, leaving the Fermi $γ$-ray bubbles and WMAP microwave haze as the remnant signal. Finally, the CRs diffuse in the bulge could further interact with the interstellar medium to produce low energy $e^+e^-$. After cooling down, these positrons may annihilate to produce the 511 keV line emission as discovered by INTEGRAL.

preprint2013arXiv

PAMELA/Fermi-LAT electron cosmic ray spectrum at $\sim$100 GeV: implication for dark matter annihilation signal in accordance with the 130 GeV $γ$-ray line

Recently, a tentative 130 GeV $γ$-ray line signal was identified by quite a few groups. If correct it would constitute a ``smoking gun'' for dark matter annihilations. Interestingly, the spectra of the cosmic ray electrons detected by PAMELA and Fermi-LAT both show tiny wiggle-like structure at $\sim 100$ GeV, which might indicate a weak signal of the annihilation of $\sim 130$ GeV dark matter particles into electrons/positrons with a velocity-weighted cross section $\langleσv\rangle_{\rm χχ\rightarrow e^{+}e^{-}} \sim 4\times10^{-26}~{\rm cm^{3}~s^{-1}}$. The prospect of identifying such a potential weak dark-matter-annihilation electron and/or positron component by AMS-02, a mission to measure the high energy cosmic ray spectra with unprecedented accuracy, is investigated.

preprint2013arXiv

Pulsar interpretation for the AMS-02 result

The AMS-02 collaboration has just published a high precision measurement of the cosmic positron fraction $e^+/(e^- + e^+)$, which rises with energy from $\sim 5$ GeV to $\sim 350$ GeV. The result indicates the existence of primary electron/positron sources to account for the positron excess. In this work, we investigate the possibility that the nearby mature pulsars are the primary positron sources. By fitting the data we find that the positrons from a single nearby pulsar, such as Geminga or Monogem, with the spectral index $α\sim 2$ can interpret the AMS-02 result. We also investigate the possibility that high energy positrons are generated by multiple known pulsars in the ATNF catalogue. Such a scenario can also fit the AMS-02 data well. Future precise measurements of fine structures in the positron spectrum would be a support to the pulsar scenario.

preprint2012arXiv

Constraints on the dark matter annihilation scenario of Fermi 130 GeV $γ$-ray line emission by continuous gamma-rays, Milky Way halo, galaxy clusters and dwarf galaxies observations

It was recently reported that there may exist monochromatic $γ$-ray emission at $\sim 130$ GeV from the Galactic center in the Fermi Large Area Telescope data, which might be related with dark matter (DM) annihilation. In this work we carry out a comprehensive check of consistency of the results with the DM annihilation scenario, using the 3.7 yrs Fermi observation of the inner Galaxy, Galactic halo, clusters of galaxies and dwarf galaxies. The results found are as follows. 1) Very strong constraints on the DM annihilation into continuous $γ$-rays from the Galactic center are set, which are as stringent as the "natural" scale assuming thermal freeze-out of DM. Such limit sets strong constraint on the DM models to explain the line emission. 2) No line emission from the Galactic halo is found in the Fermi data, and the constraints on line emission is marginally consistent with the DM annihilation interpretation of the $\sim 130$ GeV line emission from the inner Galaxy. 3) No line emission from galaxy clusters and dwarf galaxies is detected, although possible concentration of photons from clusters in 120-140 GeV is revealed. The constraints from clusters and dwarf galaxies are weak and consistent with the DM annihilation scenario to explain the $\sim 130$ GeV line emission.

preprint2012arXiv

CosRayMC: a global fitting method in studying the properties of the new sources of cosmic e$^{\pm}$ excesses

Recently PAMELA collaboration published the cosmic nuclei and electron spectra with high precision, together with the cosmic antiproton data updated, and the Fermi-LAT collaboration also updated the measurement of the total $e^+e^-$ spectrum to lower energies. In this paper we develop a Markov Chain Monte Carlo (MCMC) package {\it CosRayMC}, based on the GALPROP cosmic ray propagation model to study the implications of these new data. It is found that if only the background electrons and secondary positrons are considered, the fit is very bad with $χ_{\rm red}^2 \approx 3.68$. Taking into account the extra $e^+e^-$ sources of pulsars or dark matter annihilation we can give much better fit to these data, with the minimum $χ_{\rm red}^2 \approx 0.83$. This means the extra sources are necessary with a very high significance in order to fit the data. However, the data show little difference between pulsar and dark matter scenarios. Both the background and extra source parameters are well constrained with this MCMC method. Including the antiproton data, we further constrain the branching ratio of dark matter annihilation into quarks $B_q<0.5%$ at $2σ$ confidence level. The possible systematical uncertainties of the present study are discussed.

preprint2012arXiv

Gamma-rays From Warm WIMP Dark Matter Annihilation

The weakly interacting massive particle (WIMP) often serves as a candidate for the cold dark matter, however when produced non-thermally it could behave like warm dark matter. In this paper we study the properties of the $γ$-ray emission from annihilation of WIMP dark matter in the halo of our own Milky-Way Galaxy with high resolution $N$-body simulations of a Milky-Way like dark matter halo, assuming different nature of WIMPs. Due to the large free-streaming length in the scenario of warm WIMPs, the substructure contend of the dark matter halo is significantly different from that of the cold WIMP counterpart, resulting in distinct predictions of the $γ$-ray signals from the dark matter annihilation. We illustrate these by comparing the predicted $γ$-ray signals from the warm WIMP annihilation to that of cold WIMPs. Pronounced differences from the subhalo skymap and statistical properties between two WIMP models are demonstrated. Due to the potentially enhanced cross section of the non-thermal production mechanism in warm WIMP scenario, the Galactic center might be prior for the indirect detection of warm WIMPs to dwarf galaxies, which might be different from the cold dark matter scenario. As a specific example we consider the non-thermally produced neutralino of supersymmetric model and discuss the detectability of warm WIMPs with Fermi $γ$-ray telescope.

preprint2012arXiv

Measuring the extragalactic background light from very high energy gamma-ray observations of blazars

The extragalactic background light (EBL) contains important information about stellar and galaxy evolution. It leaves imprint on the very high energy $γ$-ray spectra from sources at cosmological distances due to the process of pair production. In this work we propose to {\em measure} the EBL directly by extracting the collective attenuation effects in a number of $γ$-ray sources at different redshifts. Using a Markov Chain Monte Carlo fitting method, the EBL intensities and the intrinsic spectral parameters of $γ$-ray sources are derived simultaneously. No prior shape of EBL is assumed in the fit. With this method, we can for the first time to derive the spectral shape of the EBL model-independently. Our result shows the expected features predicted by the present EBL models and thus support the understanding of the EBL origin.

preprint2012arXiv

Neutrino emission from dark matter annihilation/decay in light of cosmic $e^{\pm}$ and $\bar{p}$ data

A self-consistent global fitting method based on the Markov Chain Monte Carlo technique to study the dark matter (DM) property associated with the cosmic ray electron/positron excesses was developed in our previous work. In this work we further improve the previous study to include the hadronic branching ratio of DM annihilation/decay. The PAMELA $\bar{p}/p$ data are employed to constrain the hadronic branching ratio. We find that the 95% ($2σ$) upper limits of the quark branching ratio allowed by the PAMELA $\bar{p}/p$ data is $\sim 0.032$ for DM annihilation and $\sim 0.044$ for DM decay respectively. This result shows that the DM coupling to pure leptons is indeed favored by the current data. Based on the global fitting results, we further study the neutrino emission from DM in the Galactic center. Our predicted neutrino flux is some smaller than previous works since the constraint from $γ$-rays is involved. However, it is still capable to be detected by the forth-coming neutrino detector such as IceCube. The improved points of the present study compared with previous works include: 1) the DM parameters, both the particle physical ones and astrophysical ones, are derived in a global fitting way, 2) constraints from various species of data sets, including $γ$-rays and antiprotons are included, and 3) the expectation of neutrino emission is fully self-consistent.

preprint2012arXiv

Search for dark matter signals with Fermi-LAT observation of globular clusters NGC 6388 and M 15

The globular clusters are probably good targets for dark matter (DM) searches in $γ$-rays due to the possible adiabatic contraction of DM by baryons. In this work we analyse the three-year data collected by {\it Fermi} Large Area Telescope of globular clusters NGC 6388 and M 15 to search for possible DM signals. For NGC 6388 the detection of $γ$-ray emission was reported by {\it Fermi} collaboration, which is consistent with the emission of a population of millisecond pulsars. The spectral shape of NGC 6388 is also shown to be consistent with a DM contribution if assuming the annihilation final state is $b\bar{b}$. No significant $γ$-ray emission from M 15 is observed. We give the upper limits of DM contribution to $γ$-ray emission in both NGC 6388 and M 15, for annihilation final states $b\bar{b}$, $W^+W^-$, $μ^+μ^-$, $τ^+τ^-$ and monochromatic line. The constraints are stronger than that derived from observation of dwarf galaxies by {\it Fermi}.

preprint2012arXiv

Statistical interpretation of the spatial distribution of current 130 GeV gamma-ray line signal within the dark matter annihilation scenario

Recently, several groups identified a tentative $γ$-ray line signal with energy $\sim 130$ GeV in the central Galaxy from the Fermi-LAT data. %The morphology study shows that the signal is consistent with dark matter %annihilation, but with an offset $\sim 220$ pc ($1.5^{\circ}$) of the %center from the Galactic center Sgr A$^{\star}$, Such a $γ-$ray line can be interpreted as the signal of dark matter annihilation. However, the offset $\sim 220$ pc ($1.5^{\circ}$) of the center of the most prominent signal region from the Galactic center Sgr A$^{\star}$ has been thought to challenge the dark matter annihilation interpretation. Considering the fact that such a 130 GeV $γ$-ray line signal consists of only $\sim14$ photons, we suggest that the "imperfect" consistency of these photons with the expected dark matter distribution is due to the limited statistics. The offset will be smaller as more signal photons have been collected in the near future. Our Monte Carlo simulation supports the above speculation.

preprint2012arXiv

Testing the 130 GeV gamma-ray line with high energy resolution detectors

Recently some hints of the existence of $γ$-ray line around 130 GeV are reported according to the analysis of Fermi-LAT data. If confirmed it would be the first direct evidence to show the existence of new physics beyond the standard model. Here we suggest that using the forthcoming high energy resolution $γ$-ray detectors, such as CALET and DAMPE, we may test whether it is real line structure or just the background effect. For DAMPE like detector with designed energy resolution $\sim1.5%$, a line significance will reach $11σ$ for the same statistics as Fermi-LAT. For about 1.4 yr survey observation, DAMPE may detect a $5σ$ signal of such a $γ$-ray line.

preprint2012arXiv

Unified model for the gamma-ray emission of supernova remnants

Shocks of supernova remnants (SNRs) are important (and perhaps the dominant) agents for production of the Galactic cosmic rays. Recent $γ$-ray observations of several SNRs have made this case more compelling. However, these broadband high-energy measurements also reveal a variety of spectral shape demanding more comprehensive modeling of emissions from SNRs. According to the locally observed fluxes of cosmic ray protons and electrons, the electron-to-proton number ratio is known to be about 1%. Assuming such a ratio is universal for all SNRs and identical spectral shape for all kinds of accelerated particles, we propose a unified model that ascribes the distinct $γ$-ray spectra of different SNRs to variations of the medium density and the spectral difference between cosmic ray electrons and protons observed at Earth to transport effects. For low density environments, the $γ$-ray emission is inverse-Compton dominated. For high density environments like systems of high-energy particles interacting with molecular clouds, the $γ$-ray emission is $π^0$-decay dominated. The model predicts a hadronic origin of $γ$-ray emission from very old remnants interacting mostly with molecular clouds and a leptonic origin for intermediate age remnants whose shocks propagate in a low density environment created by their progenitors via e.g., strong stellar winds. These results can be regarded as evidence in support of the SNR-origin of the Galactic cosmic rays.

preprint2011arXiv

A statistical model for the $γ$-ray variability of the Crab nebula

A statistical scenario is proposed to explain the $γ$-ray variability and flares of the Crab nebula, which were observed recently by the Fermi/LAT. In this scenario electrons are accelerated in a series of knots, whose sizes follow a power-law distribution. These knots presumably move outwards from the pulsar and have a distribution in the Doppler boost factor. The maximal electron energy is assumed to be proportional to the size of the knot. Fluctuations at the highest energy end of the overall electron distribution will result in variable $γ$-ray emission via the synchrotron process in the $\sim 100$ MeV range. Since highly boosted larger knots are rarer than smaller knots, the model predicts that the variability of the synchrotron emission increases with the photon energy. We realize such a scenario with a Monte-Carlo simulation and find that the model can reproduce both the two $γ$-ray flares over a period of $\sim$ year and the monthly scale $γ$-ray flux fluctuations as observed by the Fermi/LAT. The observed $γ$-ray spectra in both the steady and flaring states are also well reproduced.

preprint2011arXiv

A statistical model to explain the gamma-ray variability and flares of the Crab nebula

Recently the AGILE and Fermi/LAT detectors uncovered giant $γ$-ray flares from the Crab nebula. The duration of these flares is a few days. The Fermi/LAT data with monthly time binning further showed significant variability of the synchrotron tail of the emission, while the inverse Compton component was stable. The simultaneous or follow-up observations in X-ray, optical, infrared and radio bands did not find significant flux variation. Based on these observations, we propose that the $γ$-ray variability and flares are due to statistical fluctuations of knots that can accelerate electrons to $\sim$PeV energies. The maximum achievable energy of electrons is adopted to be proportional to the size of the knot, which is assumed to follow a power-law distribution. Thus the low energy electron flux will be stable due to the large number of small knots, while the high energy electron flux may experience large fluctuations. Monte Carlo realization of such a picture can reproduce the observational data quite well given proper model parameters.

preprint2011arXiv

Breit-Wigner Enhancement Considering the Dark Matter Kinetic Decoupling

In the paper we study the Breit-Wigner enhancement of dark matter (DM) annihilation considering the kinetic decoupling in the evolution of DM freeze-out at the early universe. Since the DM temperature decreases much faster (as $1/R^2$) after kinetic decoupling than that in kinetic equilibrium (as 1/R) we find the Breit-Wigner enhancement of DM annihilation rate after the kinetic decoupling will affect the DM relic density significantly. Focusing on the model parameters that trying to explain the anomalous cosmic positron/electron excesses observed by PAMELA/Fermi/ATIC we find the elastic scattering $Xf\to Xf$ is not efficient to keep dark matter in kinetic equilibrium, and the kinetic decoupling temperature $T_{kd}$ is comparable to the chemical decoupling temperature $T_f\sim O(10) GeV$. The reduction of the relic density after $T_{kd}$ is significant and leads to a limited enhancement factor $\sim O(10^2)$. Therefore it is difficult to explain the anomalous positron/electron excesses in cosmic rays by DM annihilation and give the correct DM relic density simultaneously in the minimal Breit-Wigner enhancement model.

preprint2011arXiv

Constraint on dark matter annihilation with dark star formation using Fermi extragalactic diffuse gamma-ray background data

It has been proposed that during the formation of the first generation stars there might be a "dark star" phase in which the power of the star comes from dark matter annihilation. The adiabatic contraction process to form the dark star would result in a highly concentrated density profile of the host halo at the same time, which may give enhanced indirect detection signals of dark matter. In this work we investigate the extragalactic $γ$-ray background from dark matter annihilation with such a dark star formation scenario, and employ the isotropic $γ$-ray data from Fermi-LAT to constrain the model parameters of dark matter. The results suffer from large uncertainties of both the formation rate of the first generation stars and the subsequent evolution effects of the host halos of the dark stars. We find, in the most optimistic case for $γ$-ray production via dark matter annihilation, the expected extragalactic $γ$-ray flux will be enhanced by 1-2 orders of magnitude. In such a case, the annihilation cross section of the supersymmetric dark matter can be constrained to the thermal production level, and the leptonic dark matter model which is proposed to explain the positron/electron excesses can be well excluded. Conversely, if the positron/electron excesses are of a dark matter annihilation origin, then the early Universe environment is such that no dark star can form.

preprint2011arXiv

Constraints and tests of the OPERA superluminal neutrinos

The superluminal neutrinos detected by OPERA indicates Lorentz invariance violation (LIV) of the neutrino sector at the order of $10^{-5}$. We study the implications of the result in this work. We find that such a large LIV implied by OPERA data will make the neutrino production process $π\to μ+ ν_μ$ kinematically forbidden for neutrino energy greater than about 5 GeV. The OPERA detection of neutrinos at 40 GeV can constrain the LIV parameter to be smaller than $3\times 10^{-7}$. Furthermore the neutrino decay in the LIV framework will modify the neutrino spectrum greatly. The atmospheric neutrino spectrum measured by IceCube can constrain the LIV parameter to the level of $10^{-12}$. The future detection of astrophysical neutrinos of Galactic sources is expected to be able to give even stronger constraint on the LIV parameter of neutrinos.

preprint2011arXiv

Cosmic ray spectral hardening due to dispersion in the source injection spectra

Recent cosmic ray (CR) experiments discovered that the CR spectra experience a remarkable hardening for rigidity above several hundred GV. We propose that this is caused by the superposition of the CR energy spectra of many sources that have a dispersion in the injection spectral indices. Adopting similar parameters as those of supernova remnants derived from the Fermi $γ$-ray observations, we can reproduce the observational CR spectra of different species well. This may be interpreted as evidence to support the supernova remnant origin of CRs below the knee. We further propose that the same mechanism may explain the "ankle" of the ultra high energy CR spectrum.

preprint2011arXiv

Cosmic ray spectral hardening due to dispersion of source injection spectra

The cosmic ray (CR) energy spectra measured with ATIC, CREAM and PAMELA showed that there is remarkable hardening for rigidity of several hundred GV. We propose that this hardening is due to the superposition of spectra from a population of sources, e.g., supernova remnants (SNRs), whose injection spectral indices have a dispersion. Adopting proper model parameters the observational data can be well explained. It is interesting that the injection source parameters are similar with that derived from gamma-ray observations of SNRs, which may support the SNR-origin of CRs. Furthermore this mechanism provides an alternative explanation of the "ankle-cutoff" structure of the ultra high energy CR spectra.

preprint2011arXiv

Modeling the Multi-Wavelength Emission of Shell-Type Supernova Remnant RX J1713.7-3946

Emission mechanisms of the shell-type supernova remnant (SNR) RX J1713.7-3946 are studied with multi-wavelength observational data from radio, X-ray, GeV $γ$-ray to TeV $γ$-ray band. A Markov Chain Monte Carlo method is employed to explore the high-dimensional model parameter space systematically. Three scenarios for the $γ$-ray emission are investigated: the leptonic, the hadronic and a hybrid one. Thermal emission from the background plasma is also included to constrain the gas density, assuming ionization equilibrium, and a 2$σ$ upper limit of about 0.03 cm$^{-3}$ is obtained as far as thermal energies account for a significant fraction of the dissipated kinetic energy of the SNR shock. Although systematic errors dominate the $χ^2$ of the spectral fit of all models, we find that 1) the leptonic model has the best constrained model parameters, whose values can be easily accommodated with a typical supernova, but gives relatively poor fit to the TeV $γ$-ray data; 2) The hybrid scenario has one more parameter than the leptonic one and improves the overall spectral fit significantly; 3) The hadronic one, which has three more parameters than the leptonic model, gives the best fit to the overall spectrum with relatively not-well-constrained model parameters and very hard spectra of accelerated particles. The uncertainties of the model parameters decrease significantly if the spectral indices of accelerated electrons and protons are the same. The hybrid and hadronic models also require an energy input into high-energy protons, which seems to be too high compared with typical values of a supernova explosion. Further investigations are required to reconcile these observations with SNR theories.

preprint2011arXiv

Neutrino emission of Fermi supernova remnants

The Fermi $γ$-ray space telescope reported the observation of several Galactic supernova remnants recently, with the $γ$-ray spectra well described by hadronic $pp$ collisions. The possible neutrino emissions from these Fermi detected supernova remnants are discussed in this work, assuming the hadronic origin of the $γ$-ray emission. The muon event rates induced by the neutrinos from these supernova remnants on typical km$^3$ neutrino telescopes, such as the IceCube and the KM3NeT, are calculated. The results show that for most of these supernova remnants the neutrino signals are too weak to be detected by the on-going or up-coming neutrino experiment. Only for the TeV bright sources RX J1713.7-3946 and possibly W28 the neutrino signals can be comparable with the atmospheric background in the TeV region, if the protons can be accelerated to very high energies. The northern hemisphere based neutrino telescope might detect the neutrinos from these two sources.

preprint2010arXiv

A Markov Chain Monte Carlo Study on Dark Matter Property Related to the Cosmic e$^{\pm}$ Excesses

In this paper we develop a Markov Chain Monte Carlo code to study the dark matter properties in interpreting the recent observations of cosmic ray electron/positron excesses. We assume that the dark matter particles couple dominantly to leptons and consider two cases separately with annihilating or decaying into lepton pairs. The constraint on the central density profile from diffuse $γ$-rays around the Galactic center is also included in the Markov Chain Monte Carlo code self-consistently. In the numerical study, 7 parameters are introduced. We fit two data sets independently and find that for the Data set I (PAMELA+ATIC), dark matter with $m_χ\approx0.7$ TeV for annihilation (or 1.4 TeV for decay) and pure $e^+e^-$ channel is favored, while for the Data set II (PAMELA+Fermi-LAT+H.E.S.S.) $m_χ\approx 2$ TeV for annihilation (or 4 TeV for decay) and the combination of $μ^+μ^-$ and $τ^+τ^-$ final states can best fit the data. The H.E.S.S. observation of the Galactic center $γ$-rays puts a strong constraint on the central density profile of the dark matter halo for the annihilation dark matter scenario. In this case the NFW profile which is regarded as the typical predication from the cold dark matter scenario, is excluded with a high significance ($>5σ$). For the decaying dark matter scenario, the constraint is much weaker.

preprint2010arXiv

A study on the sharp knee and fine structures of cosmic ray spectra

The paper investigates the overall and detailed features of cosmic ray (CR) spectra in the knee region using the scenario of nuclei-photon interactions around the acceleration sources. Young supernova remnants can be the physical realities of such kind of CR acceleration sites. The results show that the model can well explain the following problems simultaneously with one set of source parameters: the knee of CR spectra and the sharpness of the knee, the detailed irregular structures of CR spectra, the so-called "component B" of Galactic CRs, and the electron/positron excesses reported by recent observations. The coherent explanation serves as evidence that at least a portion of CRs might be accelerated at the sources similar to young supernova remnants, and one set of source parameters indicates that this portion mainly comes from standard sources or from a single source.

preprint2010arXiv

Dark matter annihilation and non-thermal Sunyaev-Zel'dovich effect: II. dwarf spheroidal galaxy

We calculate the CMB temperature distortion due to the energetic electrons and positrons produced by dark matter annihilation (Sunyaev-Zel'dovich effect), in dwarf spheroidal galaxies (dSphs). In the calculation we have included two important effects which were previously ignored. First we show that the electron-positron pairs with energy less than GeV, which were neglected in previous calculation, could contribute a significant fraction of the total signal. Secondly we also consider the full effects of diffusion loss, which could significantly reduce the density of electron-positron pairs at the center of cuspy halos. For neutralinos, we confirm that detecting such kind of SZ effect is beyond the capability of the current or even the next generation experiments. In the case of light dark matter (LDM) the signal is much larger, but even in this case it is only marginally detectable with the next generation of experiment such as ALMA. We conclude that similar to the case of galaxy clusters, in the dwarf galaxies the $SZ_2DM}$ effect is not a strong probe of DM annihilations.

preprint2010arXiv

Galactic diffuse gamma rays --- recalculation based on the new measurements of cosmic electron spectrum

In this work, we revisit the all-sky Galactic diffuse $γ$-ray emission taking into account the new measurements of cosmic ray electron/positron spectrum by PAMELA, ATIC and Fermi, which show excesses of cosmic electrons/positrons beyond the expected fluxes in the conventional model. Since the origins of the extra electrons/positrons are not clear, we consider three different scenarios to account for the excesses: the astrophysical sources such as the Galactic pulsars, dark matter decay and annihilation. Further, new results from Fermi-LAT of the (extra-)Galactic diffuse $γ$-ray are adopted. The background cosmic rays without the new sources give lower diffuse $γ$ rays compared to Fermi-LAT observation, which is consistent with previous analysis. The scenario with astrophysical sources predicts diffuse $γ$-rays with little difference with the background. The dark matter annihilation models with $τ^{\pm}$ final state are disfavored by the Fermi diffuse $γ$-ray data, while there are only few constraints on the decaying dark matter scenario. Furthermore, these is always a bump at higher energies ($\sim$ TeV) of the diffuse $γ$-ray spectra for the dark matter scenarios due to final state radiation. Finally we find that the Fermi-LAT diffuse $γ$-ray data can be explained by simply enlarging the normalization of the electron spectrum without introduce any new sources, which may indicate that the current constraints on the dark matter models can be much stronger given a precise background estimate.

preprint2010arXiv

Gamma rays and neutrinos from dark matter annihilation in galaxy clusters

The $γ$-ray and neutrino emissions from dark matter (DM) annihilation in galaxy clusters are studied. After about one year operation of Fermi-LAT, several nearby clusters are reported with stringent upper limits of GeV $γ$-ray emission. We use the Fermi-LAT upper limits of these clusters to constrain the DM model parameters. We find that the DM model distributed with substructures predicted in cold DM (CDM) scenario is strongly constrained by Fermi-LAT $γ$-ray data. Especially for the leptonic annihilation scenario which may account for the $e^{\pm}$ excesses discovered by PAMELA/Fermi-LAT/HESS, the constraint on the minimum mass of substructures is of the level $10^2-10^3$ M$_{\odot}$, which is much larger than that expected in CDM picture, but is consistent with a warm DM scenario. We further investigate the sensitivity of neutrino detections of the clusters by IceCube. It is found that neutrino detection is much more difficult than $γ$-rays. Only for very heavy DM ($\sim 10$ TeV) together with a considerable branching ratio to line neutrinos the neutrino sensitivity is comparable with that of $γ$-rays.

preprint2010arXiv

Leptonic dark matter annihilation in the evolving universe: constraints and implications

The cosmic electron and positron excesses have been explained as possible dark matter (DM) annihilation products. In this work we investigate the possible effects of such a DM annihilation scenario during the evolution history of the Universe. We first calculate the extragalactic $γ$-ray background (EGRB), which is produced through the final state radiation of DM annihilation to charged leptons and the inverse Compton scattering between electrons/positrons and the cosmic microwave background. The DM halo profile and the minimal halo mass, which are not yet well determined from the current N-body simulations, are constrained by the EGRB data from EGRET and Fermi telescopes. Then we discuss the impact of such leptonic DM models on cosmic evolution, such as the reionization and heating of intergalactic medium, neutral Hydrogen 21 cm signal and suppression of structure formation. We show that the impact on the Hydrogen 21 cm signal might show interesting signatures of DM annihilation, but the influence on star formation is not remarkable. Future observations of the 21 cm signals could be used to place new constraints on the properties of DM.

preprint2010arXiv

On the Detectability of Galactic Dark Matter Annihilation into Monochromatic Gamma-rays

Monochromatic gamma-rays are thought to be the smoking gun signal for identifying the dark matter annihilation. However, the flux of monochromatic gamma-rays is usually suppressed by the virtual quantum effects since dark matter should be neutral and does not couple with gamma-rays directly. In the work we study the detection strategy of the monochromatic gamma-rays in a future space-based detector. The monochromatic gamma-ray flux is calculated by assuming supersymmetric neutralino as a typical dark matter candidate. We discuss both the detection focusing on the Galactic center and in a scan mode which detects gamma-rays from the whole Galactic halo are compared. The detector performance for the purpose of monochromatic gamma-rays detection, with different energy and angular resolution, field of view, background rejection efficiencies, is carefully studied with both analytical and fast Monte-Carlo method.

preprint2009arXiv

Dark matter annihilation and non-thermal Sunyaev-Zel'dovich effect: I. galaxy cluster

In this work we calculate the Sunyaev-Zel'dovich (SZ) effect due to the $e^+e^-$ from dark matter (DM) annihilation in galaxy clusters. Two candidates of DM particle, (1) the weakly-interacting massive particle (WIMP) and (2) the light dark matter (LDM) are investigated. For each case, we also consider several DM profiles with and without central cusp. We generally find smaller signals than previously reported. Moreover, the diffusion of electrons and positrons in the galaxy clusters, which was generally thought to be negligible, is considered and found to have significant effect on the central electron/positron distribution for DM profile with large spatial gradient. We find that the SZ effect from WIMP is almost always non-observable, even for the highly cuspy DM profile, and using the next generation SZ interferometer such as ALMA. Although the signal of the LDM is much larger than that of the WIMP, the final SZ effect is still very small due to the smoothing effect of diffusion. Only for the configuration with large central cusp and extremely small diffusion effect, the LDM induced SZ effect might have a bit chance of being detected.

preprint2009arXiv

Non-Thermal Production of WIMPs, Cosmic $e^\pm$ Excesses and $γ$-rays from the Galactic Center

In this paper we propose a dark matter model and study aspects of its phenomenology. Our model is based on a new dark matter sector with a U(1)' gauge symmetry plus a discrete symmetry added to the Standard Model of particle physics. The new fields of the dark matter sector have no hadronic charges and couple only to leptons. Our model can not only give rise to the observed neutrino mass hierarchy, but can also generate the baryon number asymmetry via non-thermal leptogenesis. The breaking of the new U(1)' symmetry produces cosmic strings. The dark matter particles are produced non-thermally from cosmic string loop decay which allows one to obtain sufficiently large annihilation cross sections to explain the observed cosmic ray positron and electron fluxes recently measured by the PAMELA, ATIC, PPB-BETS, Fermi-LAT, and HESS experiments while maintaining the required overall dark matter energy density. The high velocity of the dark matter particles from cosmic string loop decay leads to a low phase space density and thus to a dark matter profile with a constant density core in contrast to what happens in a scenario with thermally produced cold dark matter where the density keeps rising towards the center. As a result, the flux of gamma rays radiated from the final leptonic states of dark matter annihilation from the Galactic center is suppressed and satisfies the constraints from the HESS gamma-ray observations.

preprint2009arXiv

Parameters in a Class of Leptophilic Dark Matter Models from PAMELA, ATIC and FERMI

In this work we study a class of leptophilic dark matter models, where the dark matter interacts with the standard model particles via the $U(1)_{L_i-L_j}$ gauge boson, to explain the $e^{\pm}$ excess in cosmic rays observed by ATIC and PAMELA experiments, and more recently by Fermi experiment. There are three types of $U(1)_{L_i-L_j}$ models: a) $U(1)_{L_e - L_μ}$, b) $U(1)_{L_e - L_τ}$, and c) $U(1)_ {L_e-L_τ}$. Although ATIC or Fermi data is consistent with PAMELA data separately, ATIC and Fermi data do not agree with each other. We therefore aim to identify which of the three models can explain which data set better. We find that models a) and b) can give correct dark matter relic density and explain the ATIC and PAMELA data simultaneously recur to the Breit-Wigner enhancement. Whereas model c) with a larger $Z^\prime$ mass can explain Fermi and PAMELA data simultaneously. In all cases the model parameters are restricted to narrow regions. Future improved data will decide which set of data are correct and also help to decide the correct dark matter model.

preprint2009arXiv

Perspective of Galactic dark matter subhalo detection on Fermi from the EGRET observation

The perspective of the detectability of Galactic dark matter subhaloes on the Fermi satellite is investigated in this work. Under the assumptions that dark matter annihilation accounts for the "GeV excess" of the Galactic diffuse $γ$-rays discovered by EGRET and the $γ$-ray flux is dominated by the contribution from subhaloes of dark matter, we calculate the expected number of dark matter subhaloes that Fermi may detect. We show that Fermi may detect a few tens to several hundred subhaloes in 1-year all sky survey. Since EGRET observation is taken as a normalization, this prediction is independent of the particle physics property of dark matter. The uncertainties of the prediction are discussed in detail. We find that the major uncertainty comes from the mass function of subhaloes, i.e., whether the subhaloes are "point like" (high-mass rich) or "diffuse like" (low-mass rich). Other uncertainties like the background estimation and the observational errors will contribute a factor of $2\sim 3$.

preprint2008arXiv

Cosmological Constraints from Hubble parameter $H(z)$ and SN Ia observations

In this paper, we use a set of observational $H(z)$ data (OHD) to constrain the $Λ$CDM cosmology. This data set can be derived from the differential ages of the passively evolving galaxies. Meanwhile, the $\mathcal {A}$-parameter, which describes the Baryonic Acoustic Oscillation (BAO) peak, and the newly measured value of the Cosmic Microwave Background (CMB) shift parameter $\mathcal {R}$ are used to present combinational constraints on the same cosmology. The combinational constraints favor an accelerating flat universe while the flat $Λ$CDM cosmology is also analyzed in the same way. We obtain a result compatible with that by many other independent cosmological observations. We find that the observational $H(z)$ data set is a complementarity to other cosmological probes.

preprint2007arXiv

The Galactic positron flux and dark matter substructures

In this paper we calculate the Galactic positron flux from dark matter annihilation in the frame of supersymmetry, taking the enhancement of the flux by existence of dark matter substructures into account. The propagation of positrons in the Galactic magnetic field is solved in a realistic numerical model GALPROP. The secondary positron flux is recalculated in the GLAPROP model. The total positron flux from secondary products and dark matter annihilation can fit the HEAT data well when taking a cuspy density profile of the substructures.

Qiang Yuan

What is connected

Connect this record

See the researcher in context

Building this map preview

81 published item(s)

Geomagnetic constraints on Millicharged Dark Matter

Ultrahigh-Energy Gamma-ray Emission Associated with Black Hole-Jet Systems

Constraints on dark matter annihilation from the Event Horizon Telescope Observations of M87$^\star$

Divergence in mass ratio distributions between low-mass and high-mass coalescing binary black holes

Explanation of nearby SNRs for primary electron excess and proton spectral bump

Inverting cosmic ray propagation by Convolutional Neural Networks

Stringent axion constraints with Event Horizon Telescope polarimetric measurements of M87$^\star$

The Large High Altitude Air Shower Observatory (LHAASO) Science Book (2021 Edition)

Ultrahigh-energy diffuse gamma-ray emission from cosmic-ray interactions with medium surrounding acceleration sources

Understanding the phase reversals of Galactic cosmic ray anisotropies

A systematic study of $γ$-ray flares from the Crab Nebula with Fermi-LAT: I. flare detection

Constraining $U(1)_{L_μ-L_τ}$ charged dark matter model for muon $g-2$ anomaly with AMS-02 electron and positron data

Geminga SNR: Possible candidate of local cosmic-ray factory

Klein-Nishina effect and the cosmic ray electron spectrum

Lightcurve Evolution of the nearest Tidal Disruption Event: A late-time, radio-only flare

A Comparative Study of Long and Short GRBs. II. A Multi-wavelength Method to distinguish Type II (massive star) and Type I (compact star) GRBs

A serendipitous discovery of GeV gamma-ray emission from supernova 2004dj in a survey of nearby star-forming galaxies with Fermi-LAT

Boresight Alignment of DArk Matter Particle Explorer

Correction Method for the Readout Saturation of the DAMPE Calorimeter

Interpretation of the cosmic ray positron and electron excesses with an annihilating-decaying dark matter scenario

Local interstellar spectra and solar modulation of cosmic ray electrons and positrons

Probing Axions with Event Horizon Telescope Polarimetric Measurements

Probing Flavor Structure of Cosmic Ray $e^\mp$ Spectrum and Implications for Dark Matter Indirect Searches

Searching for the possible signal of the photon-axionlike particle oscillation in the combined GeV and TeV spectra of supernova remnants

A GeV source in the direction of Supernova Remnant CTB 37B

Catching jetted tidal disruption events early in millimetre

Constraints on the dark matter annihilation from Fermi-LAT observation of M31

Implications of the Tentative Association between GW150914 and a {\it Fermi}-GBM Transient

Perspective of monochromatic gamma-ray line detection with the High Energy cosmic-Radiation Detection (HERD) facility onboard China's Space Station

Search for gamma-ray line features from Milky Way satellites with Fermi-LAT Pass 8 data

Testing the dark matter subhalo hypothesis of the gamma-ray source 3FGL J2212.5+0703

A Systematic Chandra study of Sgr A$^{\star}$: I. X-ray flare detection

Quantitative study of the AMS-02 electron/positron spectra: implications for the pulsar and dark matter properties

Statistics of X-ray flares of Sagittarius A*: evidence for solar-like self-organized criticality phenomenon

Systematic study of the uncertainties in fitting the cosmic positron data by AMS-02

Testing the millisecond pulsar scenario of the Galactic center gamma-ray excess with very high energy gamma-rays

An Updated Analysis of Inert Higgs Doublet Model in light of the Recent Results from LUX, PLANCK, AMS-02 and LHC

Fermi Large Area Telescope Detection of Supernova Remnant RCW 86

Fermi Large Area Telescope observations of the supernova remnant HESS J1731-347

Hybrid model of GeV-TeV gamma ray emission from Galactic Center

Implications of the AMS-02 positron fraction in cosmic rays

Millisecond pulsar interpretation of the Galactic center gamma-ray excess

Multi-wavelength study of the MGRO J2019+37

Reconcile the AMS-02 positron fraction and Fermi-LAT/HESS total $e^{\pm}$ spectra by the primary electron spectrum hardening

Status of Dark Matter Detection

A generic method to constrain the dark matter model parameters from Fermi observations of dwarf spheroids

Emitting electrons spectra and acceleration processes in the jet of Mrk 421: from low state to giant flare state

On the Galactic Center Being the Main Source of Galactic Cosmic Rays as Evidenced by Recent Cosmic Ray and Gamma Ray Observations

PAMELA/Fermi-LAT electron cosmic ray spectrum at $\sim$100 GeV: implication for dark matter annihilation signal in accordance with the 130 GeV $γ$-ray line

Pulsar interpretation for the AMS-02 result

Constraints on the dark matter annihilation scenario of Fermi 130 GeV $γ$-ray line emission by continuous gamma-rays, Milky Way halo, galaxy clusters and dwarf galaxies observations

CosRayMC: a global fitting method in studying the properties of the new sources of cosmic e$^{\pm}$ excesses

Gamma-rays From Warm WIMP Dark Matter Annihilation

Measuring the extragalactic background light from very high energy gamma-ray observations of blazars

Neutrino emission from dark matter annihilation/decay in light of cosmic $e^{\pm}$ and $\bar{p}$ data

Search for dark matter signals with Fermi-LAT observation of globular clusters NGC 6388 and M 15

Statistical interpretation of the spatial distribution of current 130 GeV gamma-ray line signal within the dark matter annihilation scenario

Testing the 130 GeV gamma-ray line with high energy resolution detectors

Unified model for the gamma-ray emission of supernova remnants

A statistical model for the $γ$-ray variability of the Crab nebula

A statistical model to explain the gamma-ray variability and flares of the Crab nebula

Breit-Wigner Enhancement Considering the Dark Matter Kinetic Decoupling

Constraint on dark matter annihilation with dark star formation using Fermi extragalactic diffuse gamma-ray background data

Constraints and tests of the OPERA superluminal neutrinos

Cosmic ray spectral hardening due to dispersion in the source injection spectra

Cosmic ray spectral hardening due to dispersion of source injection spectra

Modeling the Multi-Wavelength Emission of Shell-Type Supernova Remnant RX J1713.7-3946

Neutrino emission of Fermi supernova remnants

A Markov Chain Monte Carlo Study on Dark Matter Property Related to the Cosmic e$^{\pm}$ Excesses

A study on the sharp knee and fine structures of cosmic ray spectra

Dark matter annihilation and non-thermal Sunyaev-Zel'dovich effect: II. dwarf spheroidal galaxy

Galactic diffuse gamma rays --- recalculation based on the new measurements of cosmic electron spectrum

Gamma rays and neutrinos from dark matter annihilation in galaxy clusters

Leptonic dark matter annihilation in the evolving universe: constraints and implications