Researcher profile

Yi-Zhong Fan

Yi-Zhong Fan contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
27works
0followers
6topics
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

27 published item(s)

preprint2026arXiv

Current Status of Inert Higgs Dark Matter with Dark Fermions

The precision measurements of the muon magnetic moment and the $W$ boson mass have sparked interest in the potential deviations from standard model (SM) predictions. While it may be premature to attribute any excesses in these precision measurements to new physics, they do offer a valuable indication of potential directions for physics beyond the SM. Additionally, the particle nature of dark matter (DM) remains a crucial enigma. Despite the absence of any definitive DM signal in direct detection and collider experiments, the Galactic Center GeV $γ$-ray excess and the AMS-02 antiproton ($\overline{p}$) excess could potentially offer hints related to the evidence of DM. Motivated by these observations, we propose a simple DM model that addresses all these issues. This model extends the SM by incorporating singlet and doublet Dirac fermion fields, along with a doublet complex scalar field. For the viable parameter regions in this model, we find that future upgrades of the Large Hadron Collider and DM direct detection experiments can only partially probe them, while future high-energy muon colliders hold promise for exploring the unexplored parameter space.

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.

preprint2025arXiv

Verification of the Black Hole Area Law with GW230814

We present an observational confirmation of Hawking's black-hole area theorem using the newly released gravitational-wave data from the GWTC-4.0. We analyze the high signal-to-noise ratio binary black hole (BBH) merger GW230814 and measure the (total) horizon area of the black holes before and after the merger. For preferred (and reasonable) choices of the post-truncation start time, the horizon area of the remnant black hole is found to be greater than the total horizon area of the two pre-merger black holes at a high possibility (at least $\gtrsim 99.5\%$). Importantly, our analysis accounts for sky-location uncertainty. These results provide a stringent observational confirmation of the black-hole area law, further bolstering the validity of classical general relativity in the dynamical, strong-field regime.

preprint2023arXiv

Potential Subpopulations and Assembling Tendency of the Merging Black Holes

The origins of coalescing binary black holes (BBHs) detected by the advanced LIGO/Virgo are still under debate, and clues may be present in the joint mass-spin distribution of these merger events. Here we construct phenomenological models containing two sub-populations to investigate the BBH population detected in gravitational wave observations. We find that our models can explain the GWTC-3 data rather well, and several constraints to our model are required by the data: first, the maximum mass for the component with a stellar-origin, $m_{\rm max}$, is $39.1^{+2.4}_{-2.7}M_{\odot}$ at 90\% credibility; second, about $15\%$ of the mergers happen in dynamical environments, in which $7-16\%$ of events are hierarchical mergers, and these BHs have an average spin magnitude significantly larger than the first-generation mergers, with ${\rm d}μ_{\rm a} > 0.4 $ at $99\%$ credibility; third, the dynamical component BHs tend to pair with each other with larger total mass and higher mass ratio. An independent analysis focusing on spins is also carried out, and we find that the spin amplitude of component BHs can be divided into two groups according to a division mass $m_{\rm d} = 46.1^{+5.6}_{-5.1}M_{\odot}$. These constraints can be naturally explained by current formation channels, and our results suggest that some of the observed events were likely from AGN disks.

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

GB6 J2113+1121: A multi-wavelength flaring gamma-ray blazar temporally and spatially coincident with the neutrino event IceCube-191001A

A radio-emitting tidal disruption event (AT2019dsg) is proposed as a likely counterpart of the IceCube neutrino event IC-191001A. In this work we have revisited the {\it Fermi}-LAT data in the direction of the neutrino and confirmed no signal at the site of AT2019dsg. Instead, at the edge of the 90\% confidential level error region of this neutrino there is a $γ$-ray transient source associated with a blazar GB6 J2113+1121. In May 2019, GB6 J2113+1121 was undergoing an unprecedented $γ$-ray flare since the start of the {\it Fermi}-LAT operation, with a variability amplitude about 20-fold. Similar violent flares of GB6 J2113+1121, unobserved before, have been also detected observed in optical bands. Moreover, the blazar remained in a high flux state in the infrared bands when IceCube-191001A arrived, though its $γ$-ray and optical activities has temporally ceased. Motivated by the spatial and temporal coincidence, we suggest that GB6 J2113+1121 is a candidate of the counterpart of IC-191001A. The jet properties of GB6 J2113+1121 are investigated, which are found to be comparable with that of the neutrino-emitting blazars (candidates). A specific analysis of archival IceCube data in this direction and future observations would put a further constraint on the origin of the neutrino.

preprint2022arXiv

Inert Higgs Dark Matter for CDF-II W-boson Mass and Detection Prospects

The $W$-boson mass, which was recently measured at FermiLab with an unprecedented precision, suggests the presence of new multiplets beyond the Standard Model (SM). One of the minimal extensions of the SM is to introduce an additional scalar doublet, in which the non-SM scalars can enhance $W$-boson mass via the loop corrections. On the other hand, with a proper discrete symmetry, the lightest new scalar in the doublet can be stable and play the role of dark matter particle. We show that the inert two Higgs doublet model can naturally handle the new $W$-boson mass without violating other constraints, and the preferred dark matter mass is between $54$ and $74$ GeV. We identify three feasible parameter regions for the thermal relic density: the $SA$ co-annihilation, the Higgs resonance, and the $SS \to WW^*$ annihilation. We find that the first region can be fully tested by the HL-LHC, the second region will be tightly constrained by direct detection experiments, and the third region could yield detectable GeV gamma-ray and antiproton signals in the Galaxy that may have been observed by Fermi-LAT and AMS-02.

preprint2022arXiv

Muonphilic Dark Matter explanation of gamma-ray galactic center excess: a comprehensive analysis

The Galactic center gamma-ray excess (GCE) is a long-standing unsolved problem. One of candidate solutions, the dark matter (DM) annihilation, has been recently tested with other astrophysical observations, such as AMS-02 electron-positron spectra, Fermi Dwarf spheroidal galaxies gamma-ray data, and so on. By assuming that the DM particles annihilate purely into a normal charged fermion pair, Di Mauro and Winkle (2021) claimed that only a muon-pair is compatible with the null detection of all the corresponding astrophysical measurements and can explain GCE simultaneously. On the other hand, a muonphilic DM model may also lead to a signal in the recent Fermilab muon $g-2$ measurement or be constrained by the latest PandaX-4T limit. In this work, we comprehensively study interactions between DM and muon, including various combinations of DM and mediator spins. In agreement with GCE (not only $2μ$ but also $4μ$ final states), we test these interactions against all the thermal DM constraints. Our results show that only the parameter space near the resonance region of mediator can explain GCE and relic density simultaneously, and larger parameter spaces are still allowed if other poorly-known systematic uncertainties are included. Regardless of the DM spin, only the interactions with the spin-0 mediator can explain the recent muon $g-2$ excess on top of GCE, relic density, and other DM and mediator constraints.

preprint2022arXiv

Prospects of calibrating afterglow modeling of short GRBs with gravitational wave inclination angle measurements and resolving the Hubble tension with a GW-GRB association event

In the numerical modeling of the GRB afterglow, some approximations have been made for simplicity, and different groups developed their codes. A robust test of these models/approaches is challenging because of the lack of directly measured physical parameters. Fortunately, the viewing angle inferred from the afterglow modeling is widely anticipated to be the same as the inclination angle of the binary neutron star (BNS) mergers that can be evaluated with the gravitational wave (GW) data. Therefore in the future, it is possible to calibrate the afterglow modeling with the GW inclination angle measurements. We take three methods, including both analytical estimations and direct simulations, to project the uncertainties of the inclination angle measurements. For some BNS mergers accompanied with electromagnetic counterparts detected in the O4/O5 runs of LIGO/Virgo/KAGRA/LIGO-India detectors, we show that the inclination angle can be determined within an uncertainty of $\leq 0.1$ rad, supposing that the Hubble constant is known with an accuracy of $\leq 3\%$ and the uncertainty of Hubble flow velocity is within $\sim 1 \%$. The off-axis GRB outflow will give rise to afterglow emission, and the most energetic ones may be detectable at the distance of $\sim$ 100-200 Mpc even for a viewing angle of $\geq 0.3$ rad. Such events can thus serve as a robust test of the afterglow modeling approach. We have also evaluated the prospect of resolving the so-called Hubble tension with a single GW/GRB association event. We find out that a $\sim 3\%$ precision Hubble constant is obtainable if the uncertainty of the viewing angle can be constrained to be within $\sim 0.1$ rad, which is expected to be the case for some nearby ($\leq 250$ Mpc) bright/on-axis GRBs with a well-behaved afterglow light curve displaying a clear achromatic break at early times.

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.

preprint2020arXiv

A kilonova associated with GRB 070809

For on-axis typical short gamma-ray bursts (sGRBs), the forward shock emission is usually so bright that renders the identification of kilonovae (also known as macronovae) in the early afterglow ($t<0.5$ d) phase rather challenging. This is why previously no thermal-like kilonova component has been identified at such early time except in the off-axis dim GRB 170817A associated with GW170817. Here we report the identification of an unusual optical radiation component in GRB 070809 at $t\sim 0.47$ d, thanks plausibly to the very-weak/subdominant forward shock emission. The optical emission with a very red spectrum is well in excess of the extrapolation of the X-ray emission that is distinguished by an unusually hard spectrum, which is at odds with the forward shock afterglow prediction but can be naturally interpreted as a kilonova. Our finding supports the speculation that kilonovae are ubiquitous , and demonstrates the possibility of revealing the neutron star merger origin with the early afterglow data of some typical sGRBs that take place well beyond the sensitive radius of the advanced gravitational wave detectors and hence the opportunity of organizing dedicated follow-up observations for events of interest.

preprint2020arXiv

Awakening of two $γ$-ray high redshift flat-spectrum radio quasars in the southern hemisphere

High-redshift blazars are valuable tools to study the early Universe. So far only a handful of $γ$-ray blazars have been found at redshifts above 3. $γ$-ray signals are detected in the direction of PMN J2219-2719 ($z=3.63$) and PMN J2321-0827 ($z=3.16$) by analyzing the 10-year $Fermi$-LAT Pass 8 data. PMN J2219-2719 is not distinguished from the background in the global analysis. During the 5-month epoch, the TS value is 47.8 and the flux is more than 10 times of the 10-year averaged flux. In addition, the angular distance between the $γ$-ray position and the radio position of PMN J2219-2719 is only ${0.04}^{\circ}$. Moreover, the $γ$-ray and infrared light curves of long time scale are very similar, which support the association between the $γ$-ray source and PMN J2219-2719. The global analysis of PMN J2321-0827 suggest a new $γ$-ray source, during the flare phase, the TS value is 61.4 and the $γ$-ray flux increased significantly. The association probability suggests that PMN J2321-0827 may be the counterpart of the new $γ$-ray source. In the future, the number of high-redshift $γ$-ray sources will increase by combining $Fermi$-LAT and the upcoming Large Synoptic Survey Telescope.

preprint2020arXiv

Black Hole Mass Function of Coalescing Neutron Star-Black Hole Binary Systems: The Prospect of Reconstruction with the Gravitational Wave Observations

The discovery of gravitational waves from compact objects coalescence opens a brand-new window to observe the universe. With more events being detected in the future, statistical examinations would be essential to better understand the underlying astrophysical processes. In this work we investigate the prospect of measuring the mass function of black holes that are merging with the neutron stars. Applying Bayesian parameter estimation for hundreds of simulated neutron star$-$black hole (NSBH) mergers, we find that the parameters for most of the injected events can be well recovered. We also take a Bayesian hierarchical model to reconstruct the population properties of the masses of black holes, in the presence of a low mass gap, both the mass gap and power-law index ($α$) of black hole mass function can be well measured, thus we can reveal where the $α$ is different for binary black hole (BBH) and NSBH systems. In the absence of a low mass gap, the gravitational wave data as well as the electromagnetic data can be used to pin down the nature of the merger event and then measure the mass of these very light black holes. However, as a result of the misclassification of BBH into NSBH, the measurement of $α$ is more challenging and further dedicated efforts are needed.

preprint2020arXiv

Detections of simultaneous brightening of gamma-ray and optical emissions of a distant blazar GB 1508+5714 at redshift 4.3

GB 1508+5714 is a high-redshift blazar ($z$ = 4.3) and a spectrally soft $γ$-ray source has been detected in its direction. By analyzing 11.4-yr {\it Fermi}-LAT data, significant long-term variability of the $γ$-ray source is confirmed. More importantly, a $γ$-ray emission flare appeared in an epoch of several tens of days in year 2018, when the flux is about four times of the value from the global fit. Meanwhile, optical flares displayed in both $r$- and $i$-bands from the {\it Zwicky} Transient Facility light curves. Detections of the simultaneous $γ$-ray and optical brightening provide a decisive evidence to pin down the association between the $γ$-ray source and GB 1508+5714, which makes it the {\it first} identified $γ$-ray blazar beyond redshift 4. A broadband spectral energy distribution in the high flux state is constructed and the origin of the multiwavelength brightening is also briefly discussed. Upcoming wide-deep-fast optical telescopes together with the $γ$-ray sky surveyors will shed lights on the role that the AGN jets play in the early cosmic time.

preprint2020arXiv

Estimating the maximum gravitational mass of nonrotating neutron stars from the GW170817/GRB 170817A/AT2017gfo observation

Assuming that the differential rotation of the massive neutron star (NS) formed in the double NS (DNS) mergers has been effectively terminated by the magnetic braking and a uniform rotation has been subsequently established (i.e., a supramassive NS is formed), we analytically derive in this work an approximated expression for the critical total gravitational mass ($M_{\rm tot,c}$) to form supramassive NS (SMNS) in the DNS mergers, benefited from some equation of state (EoS) insensitive relationships. The maximum gravitational mass of the nonrotating NSs ($M_{\rm TOV}$) as well as the dimensionless angular momentum of the remnant ($j$) play the dominant roles in modifying $M_{\rm tot,c}$, while the radius and mass differences of the premerger NSs do not. The GW170817/GRB 170817A/AT2017gfo observations have provided so far the best opportunity to quantitatively evaluate $M_{\rm TOV}$. Supposing the central engine for GRB 170817A is a black hole quickly formed in the collapse of an SMNS, we find $M_{\rm TOV}=2.13^{+0.09}_{-0.08}M_\odot$ (68.3% credibility interval, including also the uncertainties of the EoS insensitive relationships), which is consistent with the constraints set by current NS mass measurements.

preprint2020arXiv

GW190521 and the GWTC-1 Events: Implication on the Black Hole Mass Function of Coalescing Binary Black Hole Systems

With the black hole mass function (BHMF; assuming an exponential cutoff at a mass of $\sim 40\,M_\odot$) of coalescing binary black hole systems constructed with the events detected in the O1 run of the advanced LIGO/Virgo network, Liang et al.(2017) predicted that the birth of the lightest intermediate mass black holes (LIMBHs; with a final mass of $\gtrsim 100\,M_\odot$) is very likely to be caught by the advanced LIGO/Virgo detectors in their O3 run. The O1 and O2 observation run data, however, strongly favor a cutoff of the BHMF much sharper than the exponential one. In this work we show that a power-law function followed by a sudden drop at $\sim 40\,M_\odot$ by a factor of $\sim $a few tens and then a new power-law component extending to $\geq 100M_\odot$ are consistent with the O1 and O2 observation run data. With this new BHMF, quite a few LIMBH events can be detected in the O3 observation run of advanced LIGO/Virgo. The first LIMBH born in GW190521, an event detected in the early stage of the O3 run of advanced LIGO/Virgo network, provides additional motivation for our hypothesis.

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

Is GW190425 consistent with being a neutron star$-$black hole merger?

GW190425 is the second neutron star merger event detected by the Advanced LIGO/Virgo detectors. If interpreted as a double neutron star merger, the total gravitational mass is substantially larger than that of the binary systems identified in the Galaxy. In this work we analyze the gravitational-wave data within the neutron star$-$black hole merger scenario. For the black hole, we yield a mass of $2.40^{+0.36}_{-0.32}M_\odot$ and an aligned spin of $0.141^{+0.067}_{-0.064}$. As for the neutron star we find a mass of $1.15^{+0.15}_{-0.13}M_\odot$ and the dimensionless tidal deformability of $1.4^{+3.8}_{-1.2}\times 10^{3}$. These parameter ranges are for 90\% credibility. The inferred masses of the neutron star and the black hole are not in tension with current observations and we suggest that GW190425 is a viable candidate of a neutron star$-$black hole merger event. Benefitting from the continual enhancement of the sensitivities of the advanced gravitational detectors and the increase of the number of the observatories, similar events are anticipated to be much more precisely measured in the future and the presence of black holes below the so-called mass gap will be unambiguously clarified. If confirmed, the mergers of neutron stars with (quickly rotating) low-mass black holes are likely important production sites of the heaviest r-process elements.

preprint2020arXiv

Maximum mass cutoff in the neutron star mass distribution and the prospect of forming supramassive objects in the double neutron star mergers

The sample of neutron stars with a measured mass is growing quickly. With the latest sample, we adopt both a flexible Gaussian mixture model and a Gaussian plus Cauchy-Lorentz component model to infer the mass distribution of neutron stars and use the Bayesian model selection to explore evidence for multimodality and a sharp cutoff in the mass distribution. The two models yield rather similar results. Consistent with previous studies, we find evidence for a bimodal distribution together with a cutoff at a mass of $M_{\rm max}=2.26_{-0.05}^{+0.12}M_\odot$ (68% credible interval; for the Gaussian mixture model). If such a cutoff is interpreted as the maximum gravitational mass of nonrotating cold neutron stars, the prospect of forming supramassive remnants is found to be quite promising for the double neutron star mergers with a total gravitational mass less than or equal to 2.7$M_\odot$ unless the thermal pions could substantially soften the equation of state for the very hot neutron star matter. These supramassive remnants have a typical kinetic rotational energy of approximately $1-2\times 10^{53}$ ergs. Together with a high neutron star merger rate approximately $10^{3}~{\rm Gpc^{-3}~yr^{-3}}$, the neutron star mergers are expected to be significant sources of EeV($10^{18}$eV) cosmic-ray protons.

preprint2020arXiv

Plasma dark matter and electronic recoil events in XENON1T

Dark matter might be in the form of a dark plasma in the Milky Way halo. Specifically, we consider here a hidden sector consisting of a light `dark electron&#39; and a much heavier `dark proton&#39;, each charged under an unbroken $U(1)&#39;$ gauge symmetry. These self-interacting dark sector particles can also interact with ordinary matter via the kinetic mixing interaction, and lead to a signal in dark matter direct detection experiments. Indeed, keV electron recoils can arise quite naturally in such models from dark electron scattering off loosely bound atomic electrons. Here we examine the recently reported XENON1T excess in the context of such a plasma dark matter model. We find that the observed excess can be explained if kinetic mixing is in the approximate range: $10^{-12} \lesssim ε\lesssim 10^{-10}$. The allowed parameter space is consistent with astrophysical and cosmological constraints and consistent also with other direct detection experiments.

preprint2020arXiv

PSR J0030+0451, GW170817 and the nuclear data: joint constraints on equation of state and bulk properties of neutron stars

Very recently the NICER collaboration has published the first-ever accurate measurement of mass and radius together for PSR J0030+0451, a nearby isolated quickly-rotating neutron star (NS). In this work we set the joint constraints on the equation of state (EoS) and some bulk properties of NSs with the data of PSR J0030+0451, GW170817 and some nuclear experiments. The piecewise polytropic expansion method and the spectral decomposition method have been adopted to parameterize the EoS. The resulting constraints are consistent with each other. Assuming the maximal gravitational mass of non-rotating NS $M_{\rm TOV}$ lies between $2.04 \rm M_{\odot}$ and $2.4 \rm M_{\odot}$, with the piecewise method the pressure at twice nuclear saturation density is measured to be $3.38^{+2.43}_{-1.50}\times 10^{34}~{\rm dyn~cm^{-2}}$ at the $90\%$ level. For a NS with canonical mass of $1.4 \rm M_\odot$, we have the moment of inertia $I_{1.4} = {1.43}^{+0.28}_{-0.13} \times 10^{38}~{\rm kg \cdot m^2}$, tidal deformability $Λ_{1.4} = 390_{-140}^{+320}$, radius $R_{1.4} = 12.2_{-0.9}^{+1.0}~{\rm km}$, and binding energy $BE_{1.4} = {0.16}^{+0.01}_{-0.02} \rm M_{\odot}$ at the $90\%$ level, which are improved in comparison to the constraints with the sole data of GW170817.

preprint2020arXiv

The Masses of Isolated Neutron Stars Inferred from the Gravitational Redshift Measurements

For some neutron stars (NSs) in the binary systems, the masses have been accurately measured. While for the isolated neutron stars (INSs), no mass measurement has been reported yet. The situation will change soon thanks to the successful performance of the Neutron Star Interior Composition Explorer (NICER), with which the radius and mass of the isolated PSR J0030+0451 can be simultaneously measured. For most INSs, no mass measurements are possible for NICER because of observational limitations. Benefiting from recent significant progress made on constraining the equation of state of NSs, in this work we propose a way to estimate the masses of the INSs with the measured gravitational redshifts. We apply our method to RX J1856.5-3754, RX J0720.4-3125, and RBS 1223, three members of &#34;The Magnificent Seven&#34; (M7), and estimate their masses to be $1.24_{-0.29}^{+0.29}M_{\odot}$, $1.23_{-0.05}^{+0.10}M_{\odot}$, and $1.08_{-0.11}^{+0.20}M_{\odot}$, respectively. These masses are consistent with that of binary NS systems, suggesting no evidence for experiencing significant accretion of these isolated objects.

preprint2020arXiv

The roles of individual advanced LIGO detectors on revealing the neutron star properties

In this work we re-analyze the data of GW170817, the first binary neutron star (BNS) merger event, in two ways, including the parameterized equation of state (EoS) method and gravitational wave (GW) parameter estimation analysis. Three sets of data, including the combined LIGO/Virgo (HLV) data, the LIGO-Hanford (H1) data and the LIGO-Livingston (L1) data, have been analyzed, respectively. We confirm the bimodal probability distribution of the tidal deformability of GW170817 with the HLV data and resolve the different contributions of H1 and L1 detectors. Our simulation reveals a tendency of ``overestimating&#34; the tidal parameter with the decreasing signal-to-noise ratio (SNR). Such a tendency is likely caused by the dependence of the result on the prior (which is not well understood and has been assumed to be widely-distributed) in the low SNR case. Though this effect is interesting, the slightly higher SNR of H1 than L1 at frequencies above $\sim 100$ Hz is most-likely not the main reason for the lower tidal parameter region favored by H1. In view of the large fluctuation of the expected tide parameter in the case of low SNR, the different probability distributions favored by the L1 and H1 data of GW170817 are reasonable. We have also explored the dependence of expected difference of SNRs of H1 and L1 detectors on the source locations in the O3 and full sensitivity runs. Such maps are helpful in evaluating the relative powers of individual detectors on measuring the tidal deformabilities for the new double neutron star merger events.

preprint2019arXiv

Circinus Galaxy Revisted with 10 Years of Fermi-LAT Data

Circinus galaxy is a nearby composite starburst/AGN system. In this work we re-analyze the GeV emission from Circinus with 10 years of {\it Fermi}-LAT Pass 8 data. In the energy range of 1-500 GeV, the spectrum can be well fitted by a power-law model with a photon index of $Γ$ = $2.20\pm0.14$, and its photon flux is $(5.90\pm1.04) \times 10^{-10}$ photons cm$^{-2}$ s$^{-1}$. Our 0.1-500 GeV flux is several times lower than that reported in the previous literature, which is roughly in compliance with the empirical relation for star-forming and local group galaxies and might be reproduced by the interaction between cosmic rays and the interstellar medium. The ratio between the $γ$-ray luminosity and the total infrared luminosity is near the proton calorimetric limit, indicating that Circinus may be a proton calorimeter. However, marginal evidence for variability of the $γ$-ray emission is found in the timing analysis, which may indicate the activity of AGN jet. More {\it Fermi}-LAT data and future observation of CTA are required to fully reveal the origin of its $γ$-ray emission.

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.