Researcher profile

Qiang Yuan

Qiang Yuan contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
24works
0followers
9topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

24 published item(s)

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&#39;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&#39;&#39; 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.