Showing 1–50 of 197 results for author: Cheng, K S

Influences of dynamical disruptions on the evolution of pulsars in globular clusters

Authors: Kwangmin Oh, C. Y. Hui, Jongsuk Hong, J. Takata, A. K. H. Kong, Pak-Hin Thomas Tam, Kwan-Lok Li, K. S. Cheng

Abstract: By comparing the physical properties of pulsars hosted by core-collapsed (CCed) and non-core-collapsed (Non-CCed) globular clusters (GCs), we find that pulsars in CCed GCs rotate significantly slower than their counterparts in Non-CCed GCs. Additionally, radio luminosities at 1.4 GHz in CCed GCs are higher. These findings are consistent with the scenario that dynamical interactions in GCs can inte… ▽ More By comparing the physical properties of pulsars hosted by core-collapsed (CCed) and non-core-collapsed (Non-CCed) globular clusters (GCs), we find that pulsars in CCed GCs rotate significantly slower than their counterparts in Non-CCed GCs. Additionally, radio luminosities at 1.4 GHz in CCed GCs are higher. These findings are consistent with the scenario that dynamical interactions in GCs can interrupt angular momentum transfer processes and surface magnetic field decay during the recycling phase. Our results suggest that such effects in CCed GCs are stronger due to more frequent disruptions of compact binaries. This is further supported by the observation that both estimated disruption rates and the fraction of isolated pulsars are predominantly higher in CCed GCs. △ Less

Submitted 9 August, 2023; originally announced August 2023.

Comments: 9 pages, 8 figures, 3 tables, Accepted in MNRAS

A comparison of millisecond pulsar populations between globular clusters and the Galactic field

Authors: Jongsu Lee, C. Y. Hui, J. Takata, A. K. H. Kong, Pak-Hin Thomas Tam, Kwan-Lok Li, K. S. Cheng

Abstract: We have performed a systematic study of the rotational, orbital and X-ray properties of millisecond pulsars (MSPs) in globular clusters (GCs) and compared their nature with those of the MSPs in the Galactic field (GF). We found that GC MSPs generally rotate slower than their counterparts in the GF. Different from the expectation of a simple recycling scenario, no evidence for the correlation betwe… ▽ More We have performed a systematic study of the rotational, orbital and X-ray properties of millisecond pulsars (MSPs) in globular clusters (GCs) and compared their nature with those of the MSPs in the Galactic field (GF). We found that GC MSPs generally rotate slower than their counterparts in the GF. Different from the expectation of a simple recycling scenario, no evidence for the correlation between the orbital period and the rotation period can be found from the MSP binaries in GCs. There is also an indication that the surface magnetic field of GC MSPs are stronger than those in the GF. All these suggest dynamical interactions in GCs can alter the evolution of MSPs/their progenitors which can leave an imprint on their X-ray emission properties. While the MSPs in both GF and GCs have similar distributions of X-ray luminosity and hardness, our sample supports the notion that these two populations follow different relation between the X-ray luminosity and spin-down power. We discuss this in terms of both pulsar emission model and the observational bias. △ Less

Submitted 17 February, 2023; originally announced February 2023.

Comments: 33 pages, 13 figures, 10 tables, Accepted for publication in ApJ

Repeating fast radio burst 20201124A originates from a magnetar/Be star binary

Authors: F. Y. Wang, G. Q. Zhang, Z. G. Dai, K. S. Cheng

Abstract: Fast radio bursts (FRBs) are cosmic sources emitting millisecond-duration radio bursts. Although several hundreds FRBs have been discovered, their physical nature and central engine remain unclear. The variations of Faraday rotation measure and dispersion measure, due to local environment, are crucial clues to understanding their physical nature. The recent observations on the rotation measure of… ▽ More Fast radio bursts (FRBs) are cosmic sources emitting millisecond-duration radio bursts. Although several hundreds FRBs have been discovered, their physical nature and central engine remain unclear. The variations of Faraday rotation measure and dispersion measure, due to local environment, are crucial clues to understanding their physical nature. The recent observations on the rotation measure of FRB 20201124A show a significant variation on a day time scale. Intriguingly, the oscillation of rotation measure supports that the local contribution can change sign, which indicates the magnetic field reversal along the line of sight. Here we present a physical model that explains observed characteristics of FRB 20201124A and proposes that repeating signal comes from a binary system containing a magnetar and a Be star with a decretion disk. When the magnetar approaches the periastron, the propagation of radio waves through the disk of the Be star naturally leads to the observed varying rotation measure, depolarization, large scattering timescale, and Faraday conversion. This study will prompt to search for FRB signals from Be/X-ray binaries. △ Less

Submitted 21 September, 2022; v1 submitted 17 April, 2022; originally announced April 2022.

Comments: 42 pages, 5 figures. Published in Nature Communications https://www.nature.com/articles/s41467-022-31923-y. Title changed, and new RM data of FRB 20190520B is fitted

Searches for Pulsar-like Candidates from Unidentified Objects in the Third Catalog of Hard Fermi-LAT (3FHL) sources with Machine Learning Techniques

Authors: C. Y. Hui, Jongsu Lee, K. L. Li, Sangin Kim, Kwangmin Oh, Shengda Luo, Alex P. Leung, A. K. H. Kong, J. Takata, K. S. Cheng

Abstract: We report the results of searching pulsar-like candidates from the unidentified objects in the $3^{\rm rd}$ Catalog of Hard Fermi-LAT sources (3FHL). Using a machine-learning based classification scheme with a nominal accuracy of $\sim98\%$, we have selected 27 pulsar-like objects from 200 unidentified 3FHL sources for an identification campaign. Using archival data, X-ray sources are found within… ▽ More We report the results of searching pulsar-like candidates from the unidentified objects in the $3^{\rm rd}$ Catalog of Hard Fermi-LAT sources (3FHL). Using a machine-learning based classification scheme with a nominal accuracy of $\sim98\%$, we have selected 27 pulsar-like objects from 200 unidentified 3FHL sources for an identification campaign. Using archival data, X-ray sources are found within the $γ-$ray error ellipses of 10 3FHL pulsar-like candidates. Within the error circles of the much better constrained X-ray positions, we have also searched for the optical/infrared counterparts and examined their spectral energy distributions. Among our short-listed candidates, the most secure identification is the association of 3FHL J1823.3-1339 and its X-ray counterpart with the globular cluster Mercer 5. The $γ-$rays from the source can be contributed by a population of millisecond pulsars residing in the cluster. This makes Mercer 5 as one of the slowly growing hard $γ-$ray population of globular clusters with emission $>10$ GeV. Very recently, another candidate picked by our classification scheme, 3FHL J1405.1-6118, has been identified as a new $γ-$ray binary with an orbital period of $13.7$ days. Our X-ray analysis with a short Chandra observation has found a possible periodic signal candidate of $\sim1.4$ hrs and a putative extended X-ray tail of $\sim20$ arcsec long. Spectral energy distribution of its optical/infrared counterpart conforms with a blackbody of $T_{\rm bb}\sim40000$ K and $R_{\rm bb}\sim12R_{\odot}$ at a distance of 7.7 kpc. This is consistent with its identification as an early O star as found by infrared spectroscopy. △ Less

Submitted 30 April, 2020; v1 submitted 22 April, 2020; originally announced April 2020.

Comments: MNRAS in press, 19 pages, 13 figures, 8 tables

Reply to "Rapid $^{14}$C excursion at 3372-3371 BCE not observed at two different locations"

Authors: F. Y. Wang, H. Yu, Y. C. Zou, Z. G. Dai, K. S. Cheng

Abstract: The nuclide $^{14}$C can be produced in the atmosphere by high energy particles and $γ$-rays from high-energy phenomena. Through the carbon cycle, some of $^{14}$CO$_2$ produced in the atmosphere can be retained in annual tree rings. Four events of rapid increase of the $^{14}$C content occurred in AD 775, AD 994, BC 660 and BC 3371 were found. Recently, the data of Jull et al. (2020) was inconsis… ▽ More The nuclide $^{14}$C can be produced in the atmosphere by high energy particles and $γ$-rays from high-energy phenomena. Through the carbon cycle, some of $^{14}$CO$_2$ produced in the atmosphere can be retained in annual tree rings. Four events of rapid increase of the $^{14}$C content occurred in AD 775, AD 994, BC 660 and BC 3371 were found. Recently, the data of Jull et al. (2020) was inconsistent with our records around BC 3371. We measured our sample again and found the $^{14}$C records are consistent with the value in Wang et al. (2017). Therefore, our $^{14}$C records are robust. The inconsistency may be caused by the difference of calendar ages for the wood samples, or the physical origin of the event. First, crossdating on ring width can be performed only between trees whose growth has the same environmental conditions. Because the master tree-ring for dendrochronology is lack for Chinese trees. The master tree-ring from California has to be used. Therefore, the calendar ages derived from dendrochronology may be not precise. Second, the $^{14}$C even may be not global. One evidence is the variation of $^{14}$C content around AD 1006. The $^{14}$C contents of Californian trees increase 12\textperthousand~ in two years, while Japanese trees show no $^{14}$C increase. △ Less

Submitted 25 March, 2020; originally announced March 2020.

Comments: 11 pages, 1 figure, 1 table

Repeated state change of variable gamma-ray pulsar, PSR J2021+4026

Authors: J. Takata, H. H. Wang, L. C. C. Lin, C. -P. Hu, C. Y. Hui, A. K. H. Kong, P. H. T. Tam, K. L. Li, K. S. Cheng

Abstract: PSR J2021+4026 is a radio-quiet gamma-ray pulsar and the first pulsar that shows state change of the gamma-ray emission and spin-down rate. The state change of PSR J2021+4026 was first observed at 2011 October, at which the pulsar changes the state from high gamma-ray flux/low spin-down rate state to low gamma-ray flux/high spin-down rate st\ ate. In December 2014, PSR J2021+4026 recovered the sta… ▽ More PSR J2021+4026 is a radio-quiet gamma-ray pulsar and the first pulsar that shows state change of the gamma-ray emission and spin-down rate. The state change of PSR J2021+4026 was first observed at 2011 October, at which the pulsar changes the state from high gamma-ray flux/low spin-down rate state to low gamma-ray flux/high spin-down rate st\ ate. In December 2014, PSR J2021+4026 recovered the state before the 2011 state change over a timescale of a few months. We report that the long term evolution of the gamma-ray flux and timing behavior suggests that PSR J2021+4026 changed the state near 2018 February 1st and entered a new low gamma-ray flux/high spin-down rate state. At the 2018 state change, the averaged flux dropped from $(1.29\pm 0.01)\times 10^{-6} {\rm cts~cm^{-2}s^{-1}}$ to $(1.12\pm 0.01)\times 10^{-6} {\rm cts~cm^{-2}s^{-1\ }}$, which has the similar behavior to the case of 2011 event. The spin-down rate has increased by $\sim 3\%$ in the new state since the 2018 state change. The shapes of pulse profile and spectrum in GeV bands also changed at the 2018 event, and they are consistent with behavior at the 2011 state change. Our results probably suggest that PSR J2021+4026 is switching between different states with a timescale of several years, like some radio pulsars (e.g. PSR~B1828-11). PSR J2021+4026 will provide a unique opportunity to study the mechanism of the state switching. △ Less

Submitted 24 January, 2020; v1 submitted 1 January, 2020; originally announced January 2020.

Comments: 20 pages, 6 figures. Accepted for publication in ApJ

Consequences of Energetic Magnetar-Like Outbursts of Nearby Neutron Stars: $^{14}$C Events and the Cosmic Electron Spectrum

Authors: F. Y. Wang, Xinyu Li, D. O. Chernyshov, C. Y. Hui, G. Q. Zhang, K. S. Cheng

Abstract: Four significant events of rapid $^{14}$C increase have taken place within the past several thousand years. The physical origin of these rapid increases is still a mystery but must be associated with extremely energetic cosmic processes. Pulsars are highly magnetized neutron stars that emit a beam of electromagnetic radiations. Any sudden release of the energy stored in the magnetic multipole fiel… ▽ More Four significant events of rapid $^{14}$C increase have taken place within the past several thousand years. The physical origin of these rapid increases is still a mystery but must be associated with extremely energetic cosmic processes. Pulsars are highly magnetized neutron stars that emit a beam of electromagnetic radiations. Any sudden release of the energy stored in the magnetic multipole field will trigger outbursts similar to the giant flares of magnetars. Here we show that the relativistic outflow from the outbursts of a nearby pulsar interacting with the interstellar medium generates a shock, which accelerates electrons to trillions of electron volts. The high-energy photons from synchrotron emission of the shock interact with Earth's atmosphere, producing the cosmogenic nuclide $^{14}$C, which can cause the rapid $^{14}$C increases discovered in tree rings. These same relativistic electrons can account for a significant fraction of the cosmic electron spectrum in the trillion electron volts energy range, as observed by space-borne satellites. Since these outburst events can significantly affect our environment, monitoring nearby pulsars for such outbursts may be important in the future. △ Less

Submitted 8 November, 2019; originally announced November 2019.

Comments: 24 pages, 2 figures, accepted for publication in ApJ

Where to find Electromagnetic Wave Counterparts of stellar-mass binary black hole mergers?

Authors: Shu-Xu Yi, K. S. Cheng

Abstract: Multi-messenger astronomy combining Gravitational Wave (GW) and Electromagnetic Wave (EM) observation brings huge impact on physics, astrophysics and cosmology. However, the majority of sources to be detected with currently running ground-based GW observatories are binary black hole (BBH) mergers, which are expected disappointedly to have no EM counterparts. In this letter, we propose that if the… ▽ More Multi-messenger astronomy combining Gravitational Wave (GW) and Electromagnetic Wave (EM) observation brings huge impact on physics, astrophysics and cosmology. However, the majority of sources to be detected with currently running ground-based GW observatories are binary black hole (BBH) mergers, which are expected disappointedly to have no EM counterparts. In this letter, we propose that if the BBH merger happens in a gaseous disk around a supermassive black hole, the merger can be accompanied by a transient radio flare alike a fast radio burst (FRB). We argue that the total mass and the effective spin derived from GW detection can be used to distinguish such a source from other channels of BBH mergers. If the prediction is confirmed with future observation, multi-messenger astronomy can be brought to a distance which is one order of magnitude farther than present. The mystery of the origin of FRBs can also be revealed partially. △ Less

Submitted 18 September, 2019; originally announced September 2019.

Comments: 7 pages, 4 figures; accepted for publication in ApJL

Modelling the multi-wavelength emissions from PSR B1259-63/LS 2883: the effects of the stellar disc on shock radiations

Authors: A. M. Chen, J. Takata, S. X. Yi, Y. W. Yu, K. S. Cheng

Abstract: PSR B1259-63/LS 2883 is an elliptical pulsar/Be star binary and emits broadband emissions from radio to TeV $γ$-rays. The massive star possesses an equatorial disc, which is inclined with the orbital plane of the pulsar. The non-thermal emission from the system is believed to be produced by the pulsar wind shock and the double-peak profiles in the X-ray and TeV $γ$-ray light curves are related to… ▽ More PSR B1259-63/LS 2883 is an elliptical pulsar/Be star binary and emits broadband emissions from radio to TeV $γ$-rays. The massive star possesses an equatorial disc, which is inclined with the orbital plane of the pulsar. The non-thermal emission from the system is believed to be produced by the pulsar wind shock and the double-peak profiles in the X-ray and TeV $γ$-ray light curves are related to the phases of the pulsar passing through the disc region of the star. In this paper, we investigate the interactions between the pulsar wind and stellar outflows, especially with the presence of the disc, and present a multi-wavelength modelling of the emission from this system. We show that the double-peak profiles of X-ray and TeV $γ$-ray light curves are caused by the enhancements of the magnetic field and the soft photons at the shock during the disc passages. As the pulsar is passing through the equatorial disc, the additional pressure of the disc pushes the shock surface closer to the pulsar, which causes the enhancement of magnetic field in the shock, and thus increases the synchrotron luminosity. The TeV $γ$-rays due to the inverse-Compton (IC) scattering of shocked electrons with seed photons from the star is expected to peak around periastron which is inconsistent with observations. However, the shock heating of the stellar disc could provide additional seed photons for IC scattering during the disc passages, and thus produces the double-peak profiles as observed in the TeV $γ$-ray light curve. Our model can possibly be examined and applied to other similar gamma-ray binaries, such as PSR J2032+4127/MT91 213, HESS J0632+057, and LS I+61$^{\circ}$303. △ Less

Submitted 5 July, 2019; v1 submitted 16 April, 2019; originally announced April 2019.

Comments: 14 pages, 6 figures

Journal ref: A&A 627, A87 (2019)

On the Polarized Nonthermal Emission from AR Scorpii

Authors: J. Takata, K. S. Cheng

Abstract: We study linear polarization of optical emission from white dwarf (WD) binary system AR~Scorpii. The optical emission from this binary is modulating with the beat frequency of the system, and it is highly polarized, with the degree of the polarization reaching $\sim 40$\%. The angle of the polarization monotonically increases with the spin phase, and the total swing angle can reach $360^{\circ}$ o… ▽ More We study linear polarization of optical emission from white dwarf (WD) binary system AR~Scorpii. The optical emission from this binary is modulating with the beat frequency of the system, and it is highly polarized, with the degree of the polarization reaching $\sim 40$\%. The angle of the polarization monotonically increases with the spin phase, and the total swing angle can reach $360^{\circ}$ over one spin phase. It is also observed that the morphology of the pulse profile and the degree of linear polarization evolve with the orbital phase. These polarization properties can constrain the scenario for nonthermal emission from AR Scorpii. In this paper, we study the polarization properties predicted by the emission model, in which (i) the pulsed optical emission is produced by the synchrotron emission from relativistic electrons trapped by magnetic field lines of the WD and (ii) the emission is mainly produced at magnetic mirror points of the electron motion. We find that this model can reproduce the large swing of the polarization angle, provided that the distribution of the initial pitch angle of the electrons that are leaving the M-type star is biased to a smaller angle rather than a uniform distribution. The observed direction of the swing suggests that the Earth viewing angle is less than $90^{\circ}$ measured from the WD spin axis. The current model prefers an Earth viewing angle of $50^{\circ}-60^{\circ}$ and a magnetic inclination angle of $50^{\circ}-60^{\circ}$ (or $120^{\circ}-130^{\circ}$). We discuss that the different contribution of the emission from M-type star to total emission causes a large variation in the pulsed fraction and the degree of the linear polarization along the orbital phase. △ Less

Submitted 11 April, 2019; originally announced April 2019.

Comments: 11 pages, 9 figures. accepted for publication in ApJ

Clumpy jets from black hole-massive star binaries as engines of Fast Radio Bursts

Authors: Shu-Xu Yi, K. S. Cheng, Rui Luo

Abstract: We propose a new model of Fast Radio Bursts (FRBs) based on stellar mass black hole-massive star binaries. We argue that the inhomogeneity of the circumstellar materials or/and the time varying wind activities of the stellar companion will cause the black hole to accrete at a transient super-Eddington rate. The collision among the clumpy ejecta in the resulted jet could trigger plasma instability.… ▽ More We propose a new model of Fast Radio Bursts (FRBs) based on stellar mass black hole-massive star binaries. We argue that the inhomogeneity of the circumstellar materials or/and the time varying wind activities of the stellar companion will cause the black hole to accrete at a transient super-Eddington rate. The collision among the clumpy ejecta in the resulted jet could trigger plasma instability. As a result, the plasma in the jet will emit coherent curvature radiation. When the jet cone aims toward the observer, the apparent luminosity can be $10^{41}-10^{42}$\,erg/s. The duration of the resulted flare is $\sim$ millisecond. The high event rate of the observed non-repeating FRBs can be explained. A similar scenario in the vicinity of a supermassive black hole can be used to explain the interval distribution of the repeating source FRB121102 qualitatively. △ Less

Submitted 27 November, 2018; originally announced November 2018.

Comments: Accepted for publication in MNRAS

Gamma-Ray Emission from Molecular Clouds Generated by Penetrating Cosmic Rays

Authors: V. A. Dogiel, D. O. Chernyshov, A. V. Ivlev, D. Malyshev, A. W. Strong, K. S. Cheng

Abstract: We analyze the processes governing cosmic-ray (CR) penetration into molecular clouds and the resulting generation of gamma-ray emission. The density of CRs inside a cloud is depleted at lower energies due to the self-excited MHD turbulence. The depletion depends on the effective gas column density ("size") of the cloud. We consider two different environments where the depletion effect is expected… ▽ More We analyze the processes governing cosmic-ray (CR) penetration into molecular clouds and the resulting generation of gamma-ray emission. The density of CRs inside a cloud is depleted at lower energies due to the self-excited MHD turbulence. The depletion depends on the effective gas column density ("size") of the cloud. We consider two different environments where the depletion effect is expected to be observed. For the Central Molecular Zone, the expected range of CR energy depletion is $E\lesssim 10$ GeV, leading to the depletion of gamma-ray flux below $E_γ\approx 2$ GeV. This effect can be important for the interpretation of the GeV gamma-ray excess in the Galactic Center, which has been revealed from the standard model of CR propagation (assuming the CR spectrum inside a cloud to be equal to the interstellar spectrum). Furthermore, recent observations of some local molecular clouds suggest the depletion of the gamma-ray emission, indicating possible self-modulation of the penetrating low-energy CRs. △ Less

Submitted 13 October, 2018; originally announced October 2018.

Comments: 10 pages, 5 figures, accepted for publication in ApJ

X-ray Census of Millisecond Pulsars in the Galactic field

Authors: Jongsu Lee, C. Y. Hui, J. Takata, A. K. H. Kong, P. H. T. Tam, K. S. Cheng

Abstract: We have conducted a systematic survey for the X-ray properties of millisecond pulsars (MSPs). Currently, there are 47 MSPs with confirmed X-ray detections. We have also placed the upper limits for the X-ray emission from the other 36 MSPs by using the archival data. We have normalized their X-ray luminosities $L_{x}$ and their effective photon indices $Γ$ into a homogeneous data set, which enable… ▽ More We have conducted a systematic survey for the X-ray properties of millisecond pulsars (MSPs). Currently, there are 47 MSPs with confirmed X-ray detections. We have also placed the upper limits for the X-ray emission from the other 36 MSPs by using the archival data. We have normalized their X-ray luminosities $L_{x}$ and their effective photon indices $Γ$ into a homogeneous data set, which enable us to carry out a detailed statistical analysis. Based on our censored sample, we report a relation of $L_{x}\simeq10^{31.05}\left(\dot{E}/10^{35}\right)^{1.31}$ erg/s (2-10 keV) for the MSPs. The inferred X-ray conversion efficiency is found to be lower than previously reported estimate that could be affected by selection bias. $L_{x}$ also correlates/anti-correlates with the magnetic field strength at the light cylinder $B_{LC}$/characteristic age $τ$. On the other hand, there is no correlation between $L_{x}$ and their surface magnetic field strength $B_{s}$. We have further divided the sample into four classes: (i) black-widows, (ii) redbacks, (iii) isolated MSPs and (iv) other MSP binaries, and compare the properties among them. We noted that while the rotational parameters and the orbital periods of redbacks and black-widow are similar, $L_{x}$ of redbacks are significantly higher than those of black-widows in the 2-10 keV band. Also the $Γ$ of redbacks are apparently smaller than those of black-widows, which indicates the X-ray emission of redbacks are harder than that of black-widows. This can be explained by the different contribution of intrabinary shocks in the X-ray emission of these two classes. △ Less

Submitted 17 July, 2018; originally announced July 2018.

Comments: 19 pages, 12 figures, 4 tables, accepted for publication in ApJ

Multi-wavelength observations of a new redback millisecond pulsar candidate: 3FGL J0954.8-3948

Authors: Kwan-Lok Li, Xian Hou, Jay Strader, Jumpei Takata, Albert K. H. Kong, Laura Chomiuk, Samuel J. Swihart, Chung Yue Hui, K. S. Cheng

Abstract: We present a multi-wavelength study of the unassociated Fermi-LAT source, 3FGL J0954.8-3948, which is likely the gamma-ray counterpart of a 9.3-hour binary in the field. With more than 9 years of Pass 8 LAT data, we updated the gamma-ray spectral properties and the LAT localization of the gamma-ray source. While the binary lies outside the cataloged 95% error ellipse, the optimized LAT ellipse is… ▽ More We present a multi-wavelength study of the unassociated Fermi-LAT source, 3FGL J0954.8-3948, which is likely the gamma-ray counterpart of a 9.3-hour binary in the field. With more than 9 years of Pass 8 LAT data, we updated the gamma-ray spectral properties and the LAT localization of the gamma-ray source. While the binary lies outside the cataloged 95% error ellipse, the optimized LAT ellipse is 0.1 degrees closer and encloses the binary. The system is likely spectrally hard in X-rays (photon index ~ 1.4) with orbital modulations detected in optical, UV, and possibly X-rays. A steep spectrum radio counterpart (spectral index ~ -1.6) is also found in the TIFR GMRT Sky Survey (TGSS), implying that it is a pulsar system. We obtained a series of SOAR and Gemini spectroscopic observations in 2017/2018, which show a low-mass secondary orbiting in a close circular orbit with K2 = 272 km/s under strong irradiation by the primary compact object. All the observations as well as the modelling of the X/gamma-ray high-energy emission suggest that 3FGL J0954.8-3948 is a redback millisecond pulsar in a rotation-powered state. △ Less

Submitted 6 July, 2018; originally announced July 2018.

Comments: 10 pages, including 5 figures and 3 tables. Accepted for publication in ApJ

Time-Variability of equivalent width of 6.4 keV line from the Arches Complex: reflected X-rays or charged particles?

Authors: D. O. Chernyshov, C. M. Ko, R. A. Krivonos, V. A. Dogiel, K. S. Cheng

Abstract: Molecular gas in Arches cloud located near the Arches cluster is one of the emitters of K-$α$ line of neutral iron and X-ray continuum in the Galactic center (GC). Similarly to the cloud Sgr B2, another well-known emitter of the iron line in the GC, the Arches cloud demonstrates temporal decline of the X-ray emission. The most natural origin of this emission is irradiation of primary photons of an… ▽ More Molecular gas in Arches cloud located near the Arches cluster is one of the emitters of K-$α$ line of neutral iron and X-ray continuum in the Galactic center (GC). Similarly to the cloud Sgr B2, another well-known emitter of the iron line in the GC, the Arches cloud demonstrates temporal decline of the X-ray emission. The most natural origin of this emission is irradiation of primary photons of an X-ray flare from a distant source, most likely Sgr A*. However, recent observations of the Arches cloud discovered variations of equivalent width of the 6.4 keV iron line, which indicated that the X-ray emission from the cloud is a combination of two components with different origin and different equivalent width, one of which is time-variable, while the other is stationary during the period of observations. We considered two different scenarios: a) this emission is formed by reflection from two clouds, which are at some distance from each other, when they are irradiated by two different flares; and b) the other scenario assumes a combination of X-ray fluxes produced in the same cloud by reflection of primary photons and by subrelativistic cosmic rays. We present restrictions for both model and conditions at which these scenarios can be realized. Although none of the models can be completely ruled out, we find that the X-ray reflection model requires less assumption and therefore is the most viable. △ Less

Submitted 2 July, 2018; originally announced July 2018.

Comments: 15 pages, 3 figures. accepted for publication in ApJ

The growth of stellar mass black hole binaries trapped in an accretion disk of active galactic nuclei

Authors: Shu-Xu Yi, K. S. Cheng, Ronald E. Taam

Abstract: Among the four black hole binary merger events detected by LIGO, six progenitor black holes have masses greater than 20\,$M_\odot$. The existence of such massive BHs calls for extreme metal-poor stars as the progenitors. An alternative possibility that a pair of stellar mass black holes each with mass $\sim7\,M_\odot$ increases to $>20\,M_\odot$ via accretion from a disk surrounding a super massiv… ▽ More Among the four black hole binary merger events detected by LIGO, six progenitor black holes have masses greater than 20\,$M_\odot$. The existence of such massive BHs calls for extreme metal-poor stars as the progenitors. An alternative possibility that a pair of stellar mass black holes each with mass $\sim7\,M_\odot$ increases to $>20\,M_\odot$ via accretion from a disk surrounding a super massive black hole in an active galactic nucleus is considered. The growth of mass of the binary and the transfer of orbital angular momentum to the disk accelerates the merger. Based on the recent numerical work of Tang et al. (2017), it is found that, in the disk of a low mass AGN with mass $\sim10^6\,M_\odot$ and Eddington ratio $>0.01$, the mass of an individual BH in the binary can grow to $>20\,M_\odot$ before coalescence provided that accretion takes place at a rate more than 10 times the Eddington value. The mechanism predicts a new class of gravitational wave sources involving the merger of two extreme Kerr black holes associated with active galactic nuclei and a possible electromagnetic wave counterpart. △ Less

Submitted 17 May, 2018; originally announced May 2018.

Comments: Accepted for publication in ApJL

Toward an inner connection of SNe Ic, SLSNe Ic, XRF connected SNe, SNe Ic-BL, and GRB connected SNe

Authors: Yuan-Chuan Zou, K. S. Cheng

Abstract: Type Ic supernovae (SNe) can be classified as: normal SNe Ic, type Ic super-luminous SNe (SLSNe Ic), X-ray flash (XRF) connected SNe, broad-line SNe Ic (SN Ic-BL), and gamma-ray burst (GRB) connected SNe. Here we suggest an inner connection for all kinds of SNe Ic based on a pair of jets being successfully launched: a normal SN Ic is a normal core collapsar without jets launched; a GRB associated… ▽ More Type Ic supernovae (SNe) can be classified as: normal SNe Ic, type Ic super-luminous SNe (SLSNe Ic), X-ray flash (XRF) connected SNe, broad-line SNe Ic (SN Ic-BL), and gamma-ray burst (GRB) connected SNe. Here we suggest an inner connection for all kinds of SNe Ic based on a pair of jets being successfully launched: a normal SN Ic is a normal core collapsar without jets launched; a GRB associated SN Ic is a core collapsar with relativistic jets launched and successfully breaking out the envelope of the progenitor; an XRF associated SN Ic is a core collapsar with jets launched but can only develop a relativistic shock breakout; an SN Ic-BL is an off-axis GRB or XRF associated supernova; and an SLSN Ic is close to the XRF-SN Ic but the shock breakout is not relativistic and most of the jet energy is deposited into the supernova component. Based on the luminosity-distance diagram, we derived the luminosity function of all different types of SNe Ic as a whole. We also show that the normal SNe Ic and GRB connected SNe Ic have similar accumulative distributions. △ Less

Submitted 30 March, 2018; originally announced March 2018.

The X-ray modulation of PSR J2032+4127/MT91 213 during the Periastron Passage in 2017

Authors: K. L. Li, J. Takata, C. W. Ng, A. K. H. Kong, P. H. T. Tam, C. Y. Hui, K. S. Cheng

Abstract: We present the Neil Gehrels Swift Observatory (Swift), Fermi Large Area Telescope (Fermi-LAT), and Karl G. Jansky Very Large Array (VLA) observations of the gamma-ray binary PSR J2032+4127/MT91 213, of which the periastron passage has just occurred in November 2017. In the Swift X-ray light curve, the flux was steadily increasing before mid-October 2017, however, a sharp X-ray dip on a weekly time… ▽ More We present the Neil Gehrels Swift Observatory (Swift), Fermi Large Area Telescope (Fermi-LAT), and Karl G. Jansky Very Large Array (VLA) observations of the gamma-ray binary PSR J2032+4127/MT91 213, of which the periastron passage has just occurred in November 2017. In the Swift X-ray light curve, the flux was steadily increasing before mid-October 2017, however, a sharp X-ray dip on a weekly time-scale is seen during the periastron passage, followed by a post-periastron X-ray flare lasting for ~20 days. We suggest that the X-ray dip is caused by (i) an increase of the magnetization parameter at the shock, and (ii) the suppression due to the Doppler boosting effect. The 20-day post-periastron flare could be a consequence of the Be stellar disk passage by the pulsar. An orbital GeV modulation is also expected in our model, however, no significant variability is seen in the Fermi-LAT light curve. We suspect that the GeV emission resulted from the interaction between the binary's members is hidden behind the bright magnetospheric emission of the pulsar. Pulsar gating technique would be useful to remove the magnetospheric emission and recover the predicted GeV modulation, if an accurate radio timing solution over the periastron passage is provided in the future. △ Less

Submitted 18 March, 2018; originally announced March 2018.

Comments: 6 pages, including 2 figures. Accepted for publication in ApJ

Penetration of cosmic rays into dense molecular clouds: role of diffuse envelope

Authors: A. V. Ivlev, V. A. Dogiel, D. O. Chernyshov, P. Caselli, C. -M. Ko, K. S. Cheng

Abstract: A flux of cosmic rays (CRs) propagating through a diffuse ionized gas can excite MHD waves, thus generating magnetic disturbances. We propose a generic model of CR penetration into molecular clouds through their diffuse envelopes, and identify the leading physical processes controlling their transport on the way from a highly ionized interstellar medium to a dense interior of the cloud. The model… ▽ More A flux of cosmic rays (CRs) propagating through a diffuse ionized gas can excite MHD waves, thus generating magnetic disturbances. We propose a generic model of CR penetration into molecular clouds through their diffuse envelopes, and identify the leading physical processes controlling their transport on the way from a highly ionized interstellar medium to a dense interior of the cloud. The model allows us to describe a transition between a free streaming of CRs and their diffusive propagation, determined by the scattering on the self-generated disturbances. A self consistent set of equations, governing the diffusive transport regime in an envelope and the MHD turbulence generated by the modulated CR flux, is essentially characterized by two dimensionless numbers. We demonstrate a remarkable mutual complementarity of different mechanisms leading to the onset of the diffusive regime, which results in a universal energy spectrum of the modulated CRs. In conclusion, we briefly discuss implications of our results for several fundamental astrophysical problems, such as the spatial distribution of CRs in the Galaxy as well as the ionization, heating, and chemistry in dense molecular clouds. △ Less

Submitted 7 February, 2018; originally announced February 2018.

Comments: The manuscript is accepted for publication in the Astrophysical Journal

A non-thermal pulsed X-ray emission of AR~Scorpii

Authors: J. Takata, C. -P. Hu, L. C. C. Lin, P. H. T. Tam, P. S. Pal, C. Y. Hui, A. K. H. Kong, K. S. Cheng

Abstract: We report the analysis result of UV/X-ray emission from AR~Scorpii, which is an intermediate polar (IP) composed of a magnetic white dwarf and a M-type star, with the XMM-Newton data. The X-ray/UV emission clearly shows a large variation over the orbit, and their intensity maximum (or minimum) is located at the superior conjunction (or inferior conjunction) of the M-type star orbit. The hardness r… ▽ More We report the analysis result of UV/X-ray emission from AR~Scorpii, which is an intermediate polar (IP) composed of a magnetic white dwarf and a M-type star, with the XMM-Newton data. The X-ray/UV emission clearly shows a large variation over the orbit, and their intensity maximum (or minimum) is located at the superior conjunction (or inferior conjunction) of the M-type star orbit. The hardness ratio of the X-ray emission shows a small variation over the orbital phase, and shows no indication of the absorption by an accretion column. These properties are naturally explained by the emission from the M-type star surface rather than from the accretion column on the WD's star similar to the usual IPs. Beside, the observed X-ray emission also modulates with WD's spin with a pulse fraction of $\sim 14\%$. The peak position is aligned in the optical/UV/X-ray band. This supports the hypothesis that the electrons in AR~Scorpii are accelerated to a relativistic speed, and emit non-thermal photons via the synchrotron radiation. In the X-ray bands, the evidence of the power-law spectrum is found in the pulsed component, although the observed emission is dominated by the optically thin thermal plasma emissions with several different temperatures. It is considered that the magnetic dissipation/reconnection process on the M-type star surface heats up the plasma to a temperature of several keV, and also accelerates the electrons to the relativistic speed. The relativistic electrons are trapped in the WD's closed magnetic field lines by the magnetic mirror effect. In this model, the observed pulsed component is explained by the emissions from the first magnetic mirror point. △ Less

Submitted 18 December, 2017; originally announced December 2017.

Comments: 8 pages, 11 figures, accepted for publication in ApJ

A Model for AR~Scorpii: Emission from Relativistic Electrons Trapped by Closed Magnetic Field Lines of Magnetic White Dwarfs

Authors: J. Takata, H. Yang, K. S. Cheng

Abstract: AR~Scorpii is an intermediate polar system composed of a magnetic white dwarf (WD) and an M-type star, and shows non-thermal, pulsed, and highly linearly polarized emission. The radio/optical emission modulates with the WD's spin and show the double peak structure in the light curves. In this paper, we discuss a possible scenario for the radiation mechanism of AR~Scorpii. The magnetic interaction… ▽ More AR~Scorpii is an intermediate polar system composed of a magnetic white dwarf (WD) and an M-type star, and shows non-thermal, pulsed, and highly linearly polarized emission. The radio/optical emission modulates with the WD's spin and show the double peak structure in the light curves. In this paper, we discuss a possible scenario for the radiation mechanism of AR~Scorpii. The magnetic interaction on the surface of the companion star produces an outflow from the companion star, the heating of the companion star surface, and the acceleration of electrons to a relativistic energy. The accelerated electrons, whose typical Lorentz factor is $\sim 50-100$, from the companion star move along the magnetic field lines toward the WD surface. The electrons injected with the pitch angle of $\sinθ_{p,0}>0.05$ are subject to the magnetic mirror effect and are trapped in the closed magnetic field line region.We find that the emission from the first magnetic mirror points mainly contributes to the observed pulsed emission and the formation of the double-peak structure in the light curve. For the inclined rotator, the pulse peak in the calculated light curve shifts the position in the spin phase, and a Fourier analysis exhibits a beat frequency feature, which are consistent with the optical/UV observations. The pulse profile also evolves with the orbital phase owing to the effect of the viewing geometry. The model also interprets the global features of the observed spectral energy distribution in radio to X-ray energy bands. We also discuss the curvature radiation and the inverse-Compton scattering process in the outer gap accelerator of the WD in AR Scorpii and discuss the possibility of the detection by future high-energy missions. △ Less

Submitted 10 December, 2017; originally announced December 2017.

Comments: 33 pages, 12 figures, accepted for publication in ApJ

A rapid cosmic-ray increase in BC 3372-3371 from ancient buried tree rings in China

Authors: F. Y. Wang, H. Yu, Y. C. Zou, Z. G. Dai, K. S. Cheng

Abstract: Cosmic rays interact with the Earth's atmosphere to produce $^{14}$C, which can be absorbed by trees. Therefore, rapid increases of $^{14}$C in tree rings can be used to probe previous cosmic-ray events. By this method, three $^{14}$C rapidly increasing events have been found. Plausible causes of these events include large solar proton events, supernovae or short gamma-ray bursts. However, due to… ▽ More Cosmic rays interact with the Earth's atmosphere to produce $^{14}$C, which can be absorbed by trees. Therefore, rapid increases of $^{14}$C in tree rings can be used to probe previous cosmic-ray events. By this method, three $^{14}$C rapidly increasing events have been found. Plausible causes of these events include large solar proton events, supernovae or short gamma-ray bursts. However, due to the lack of measurements of $^{14}$C by year, the occurrence frequency of such $^{14}$C rapidly increasing events is poorly known. In addition, rapid increases may be hidden in the IntCal13 data with five-year resolution. Here we report the result of $^{14}$C measurements using an ancient buried tree during the period between BC 3388 and 3358. We find a rapid increase of about 9\textperthousand~ in the $^{14}$C content from BC 3372 to BC 3371. We suggest that this event could originate from a large solar proton event. △ Less

Submitted 13 November, 2017; originally announced November 2017.

Comments: 23 pages, 3 figures, 2 tables, published in Nature Communications

Journal ref: Nature Communications, 2017, 8, 1487

Probing the properties of the pulsar wind via studying the dispersive effects in the pulses from the pulsar companion in a double neutron-star binary system

Authors: Shu-Xu Yi, K. S. Cheng

Abstract: The velocity and density distribution of $e^\pm$ in the pulsar wind are crucial distinction among magnetosphere models, and contains key parameters determining the high energy emission of pulsar binaries. In this work, a direct method is proposed, which might probe the properties of the wind from one pulsar in a double-pulsar binary. When the radio signals from the first-formed pulsar travel throu… ▽ More The velocity and density distribution of $e^\pm$ in the pulsar wind are crucial distinction among magnetosphere models, and contains key parameters determining the high energy emission of pulsar binaries. In this work, a direct method is proposed, which might probe the properties of the wind from one pulsar in a double-pulsar binary. When the radio signals from the first-formed pulsar travel through the relativistic $e^\pm$ flow in the pulsar wind from the younger companion, the components of different radio frequencies will be dispersed. It will introduce an additional frequency-dependent time-of-arrival delay of pulses, which is function of the orbital phase. In this paper, we formulate the above-mentioned dispersive delay with the properties of the pulsar wind. As examples, we apply the formula to the double pulsar system PSR J0737-3039A/B and the pulsar-neutron star binary PSR B1913+16. For PSR J0737-3039A/B, the time delay in 300\,MHz is $\lesssim10μ$s near the superior-conjunction, under the optimal pulsar wind parameters, which is $\sim$ half of the current timing accuracy. For PSR B1913+16, with the assumption that the neutron star companion has a typical spin down luminosity of $10^{33}$\,ergs/s, the time delay is as large as $10\sim20μ$s in 300\,MHz. The best timing precision of this pulsar is $\sim5μ$s in 1400\,MHz. Therefore, it is possible that we can find this signal in archival data. Otherwise, we can set an upper-limit on the spin down luminosity. Similar analysis can be apply to other eleven known pulsar-neutron star binaries △ Less

Submitted 30 August, 2017; originally announced August 2017.

Comments: 6 pages, 6 figures, accepted for publication in MNRAS

Rapid X-ray variations of the Geminga pulsar wind nebula

Authors: C. Y. Hui, Jongsu Lee, A. K. H. Kong, P. H. T. Tam, J. Takata, K. S. Cheng, Dongsu Ryu

Abstract: A recent study by Posselt et al. (2017) reported the deepest X-ray investigation of the Geminga pulsar wind nebula (PWN) by using \emph{Chandra X-ray Observatory}. In comparison with previous studies of this system, a number of new findings have been reported and we found these suggest the possible variabilities in various components of this PWN. This motivates us to carry out a dedicated search f… ▽ More A recent study by Posselt et al. (2017) reported the deepest X-ray investigation of the Geminga pulsar wind nebula (PWN) by using \emph{Chandra X-ray Observatory}. In comparison with previous studies of this system, a number of new findings have been reported and we found these suggest the possible variabilities in various components of this PWN. This motivates us to carry out a dedicated search for the morphological and spectral variations of this complex nebula. We have discovered variabilities on timescales from a few days to a few months from different components of the nebula. The fastest change occurred in the circumstellar environment at a rate of 80 per cent of the speed of light. One of the most spectacular results is the wiggling of a half light-year long tail as an extension of the jet, which is significantly bent by the ram pressure. The jet wiggling occurred at a rate of about 20 per cent of the speed of light. This twisted structure can possibly be a result of a propagating torsional Alfven wave. We have also found evidence of spectral hardening along this tail for a period of about nine months. △ Less

Submitted 14 August, 2017; originally announced August 2017.

Comments: 23 pages, 14 figures, 4 tables, Accepted for publication in ApJ

Modeling the GeV emission of HESS J0632+057

Authors: Shu-Xu Yi, Kwong Sang Cheng

Abstract: The binary system HESS J0632+057 was recently detected by {Fermi} to possess orbital modulated GeV emission. In this paper, we study the possibility that the compact companion of HESS J0632+057 is a pulsar. Under such a presumption, we focus on the high energy emission mechanism of this system, which is as follows. The pulsar companion travels through the circumstellar disc of the main sequence st… ▽ More The binary system HESS J0632+057 was recently detected by {Fermi} to possess orbital modulated GeV emission. In this paper, we study the possibility that the compact companion of HESS J0632+057 is a pulsar. Under such a presumption, we focus on the high energy emission mechanism of this system, which is as follows. The pulsar companion travels through the circumstellar disc of the main sequence star twice in each orbit, when some of the matter is gravity-captured. The captured matter develops an accretion disc around the pulsar, and the soft photons from which are inverse Compton scattered by the pulsar wind as the GeV emission from the system. With proper choice of parameters, SED and light curve which are in accordance with observations can be produced. We predict that the light curve of GeV emission has two peaks, the larger one is at around 0.4 after the periastron (or 0.1 after the X-ray maximum), while the smaller one is between phases 0 and 0.1, with integrated flux one forth of the larger one. △ Less

Submitted 25 July, 2017; originally announced July 2017.

Comments: 7pages, 7 figures. Accepted for publication in MNRAS

A new approach to the GeV flare of PSR B1259-63/LS2883

Authors: Shu-Xu Yi, K. S. Cheng

Abstract: PSR B1259-63/LS2883 is a binary system composed of a pulsar and a Be star. The Be star has an equatorial circumstellar disk (CD). The {\it Fermi} satellite discovered unexpected gamma-ray flares around 30 days after the last two periastron passages. The origin of the flares remain puzzling. In this work, we explore the possibility that, the GeV flares are consequences of inverse Compton-scattering… ▽ More PSR B1259-63/LS2883 is a binary system composed of a pulsar and a Be star. The Be star has an equatorial circumstellar disk (CD). The {\it Fermi} satellite discovered unexpected gamma-ray flares around 30 days after the last two periastron passages. The origin of the flares remain puzzling. In this work, we explore the possibility that, the GeV flares are consequences of inverse Compton-scattering of soft photons by the pulsar wind. The soft photons are from an accretion disk around the pulsar, which is composed by the matter from CD captured by the pulsar's gravity at disk-crossing before the periastron. At the other disk-crossing after the periastron, the density of the CD is not high enough so that accretion is prevented by the pulsar wind shock. This model can reproduce the observed SEDs and light curves satisfactorily. △ Less

Submitted 27 June, 2017; originally announced June 2017.

Comments: 14 pages, 8 figures, 1 table. Accepted for publication in ApJ

Mode change of a gamma-ray pulsar, PSR J2021+4026

Authors: J. Zhao, C. W. Ng, L. C. C. Lin, J. Takata, Y. Cai, C. P. Hu, D. C. C. Yen, P. H. T. Tam, C. Y. Hui, A. K. H. Kong, K. S. Cheng

Abstract: A glitch of a pulsar is known as a sudden increase in the spin frequency and spin-down rate (frequency time derivative), and it can be caused by a sudden rel\ ease of the stress built up in the solid crust of the star or pinned vortices in the superfluid interior. PSR J2021+4026 is the first pulsar that shows a significant change in the gamma-ray flux and pulse profile at the glitch that occurred… ▽ More A glitch of a pulsar is known as a sudden increase in the spin frequency and spin-down rate (frequency time derivative), and it can be caused by a sudden rel\ ease of the stress built up in the solid crust of the star or pinned vortices in the superfluid interior. PSR J2021+4026 is the first pulsar that shows a significant change in the gamma-ray flux and pulse profile at the glitch that occurred around 2011 October 16. We report the results of timing and spectral analysis of PSR~J2021+4026 using $\sim$ 8~yr Fermi-LAT data. We find that the pulsar stayed at a high spin-down rate ($\sim 4\%$ higher than the pre-glitch value) and a low gamma-ray state ($\sim 18\%$ lower) for about 3~yr after the glitch. Around 2014 December, the spin-down rate and gamma-ray flux gradually returned to pre-glitch values within a time scale of a few months. The phase-resolved spectra and pulse profiles after the relaxation are also consistent with those before the glitch. The observed long-term evolution of the spin-down rate and the gamma-ray flux indicates that the glitch triggered a mode change in the global magnetosphere. We speculate that the glitch changed the local magnetic field structure around the polar cap and/or the inclination angle of the dipole axis, leading to a change in the electric current circulating in the magnetosphere. △ Less

Submitted 11 June, 2017; v1 submitted 1 June, 2017; originally announced June 2017.

Comments: 23 pages, 6 figures. Accepted for publication in ApJ

High-Energy emissions from the Pulsar/Be binary system PSR J2032+4127/MT91 213

Authors: J. Takata, P. H. T. Tam, C. W. Ng, K. L. Li, A. K. H. Kong, C. Y. Hui, K. S. Cheng

Abstract: PSR J2032+4127 is a radio-loud gamma-ray-emitting pulsar; it is orbiting around a high-mass Be type star with a very long orbital period of 25-50years, and is approaching periastron, which will occur in late 2017/early 2018. This system comprises with a young pulsar and a Be type star, which is similar to the so-called gamma-ray binary PSR~B1259-63/LS2883. It is expected therefore that PSR J2032+4… ▽ More PSR J2032+4127 is a radio-loud gamma-ray-emitting pulsar; it is orbiting around a high-mass Be type star with a very long orbital period of 25-50years, and is approaching periastron, which will occur in late 2017/early 2018. This system comprises with a young pulsar and a Be type star, which is similar to the so-called gamma-ray binary PSR~B1259-63/LS2883. It is expected therefore that PSR J2032+4127 shows an enhancement of high-energy emission caused by the interaction between the pulsar wind and Be wind/disk around periastron. Ho et al. recently reported a rapid increase in the X-ray flux from this system. In this paper, we also confirm a rapid increase in the X-ray flux along the orbit, while the GeV flux shows no significant change. We discuss the high-energy emissions from the shock caused by the pulsar wind and stellar wind interaction and examine the properties of the pulsar wind in this binary system. We argue that the rate of increase of the X-ray flux observed by Swift indicates (1) a variation of the momentum ratio of the two-wind interaction region along the orbit, or (2) an evolution of the magnetization parameter of the pulsar wind with the radial distance from the pulsar. We also discuss the pulsar wind/Be disk interaction at the periastron passage, and propose the possibility of formation of an accretion disk around the pulsar. We model high-energy emissions through the inverse-Compton scattering process of the cold-relativistic pulsar wind off soft photons from the accretion disk. △ Less

Submitted 14 February, 2017; originally announced February 2017.

Comments: 18 pages, 23 figures, 1 Table, accepted for publication in ApJ

Gamma-ray Emission of the Kes 73/1E 1841-045 Region Observed with the Fermi Large Area Telescope

Authors: Paul K. H. Yeung, Albert K. H. Kong, P. H. Thomas Tam, C. Y. Hui, Jumpei Takata, K. S. Cheng

Abstract: The supernova remnant Kes 73 and/or the magnetar 1E 1841-045 at its center can deposit a large amount of energy to the surroundings and is potentially responsible for particle acceleration. Using the data taken with the Fermi Large Area Telescope (LAT), we confirmed the presence of an extended source whose centroid position is highly consistent with this magnetar/supernova-remnant pair. Its emissi… ▽ More The supernova remnant Kes 73 and/or the magnetar 1E 1841-045 at its center can deposit a large amount of energy to the surroundings and is potentially responsible for particle acceleration. Using the data taken with the Fermi Large Area Telescope (LAT), we confirmed the presence of an extended source whose centroid position is highly consistent with this magnetar/supernova-remnant pair. Its emission is intense from 100 MeV to $>$100 GeV. Its LAT spectrum can be decoupled into two components which are respectively governed by two different mechanisms. According to the young age of this system, the magnetar is seemingly a necessary and sufficient source for the downward-curved spectrum below 10 GeV, as the observed $<$10 GeV flux is too high for the supernova remnant to account for. On the other hand, the supernova remnant is reasonably responsible for the hard spectrum above 10 GeV. Further studies of this region in the TeV regime is required, so that we can perform physically meaningful comparisons of the $>$10 GeV spectrum and the TeV spectrum. △ Less

Submitted 21 February, 2017; v1 submitted 31 January, 2017; originally announced January 2017.

Comments: Accepted for publication in ApJ

Cosmic Ray (Stochastic) Acceleration from a Background Plasma

Authors: V. A. Dogiel, K. S. Cheng, D. O. Chernyshov, A. D. Erlykin, C. -M. Ko, A. W. Wolfendale

Abstract: We give a short review of processes of stochastic acceleration in the Galaxy. We discuss: how to estimate correctly the number of accelerated particles, and at which condition the stochastic mechanism is able to generate power-law nonthermal spectra. We present an analysis of stochastic acceleration in the Galactic halo and discuss whether this mechanism can be responsible for production of high e… ▽ More We give a short review of processes of stochastic acceleration in the Galaxy. We discuss: how to estimate correctly the number of accelerated particles, and at which condition the stochastic mechanism is able to generate power-law nonthermal spectra. We present an analysis of stochastic acceleration in the Galactic halo and discuss whether this mechanism can be responsible for production of high energy electrons there, which emit gamma-ray and microwave emission from the giant Fermi bubbles. Lastly, we discuss whether the effects of stochastic acceleration can explain the CR distribution in the Galactic disk (CR gradient). △ Less

Submitted 19 January, 2017; originally announced January 2017.

Comments: A talk presented at the conference "Cosmic Ray Origin - beyond the standard models" (San-Vito, Italy, 2016), submitted to Nuclear and Particle Physics Proceedings

X-ray/GeV emissions from Crab-like pulsars in LMC

Authors: J. Takata, K. S. Cheng

Abstract: We discuss X-ray and gamma-ray emissions from Crab-like pulsars, PSRs~J0537-6910 and~J0540-6919, in Large Magellanic Cloud. Fermi-LAT observations have resolved the gamma-ray emissions from these two pulsars and found the pulsed emissions from PSR~J0540-6919. The total pulsed radiation in the X-ray/gamma-ray energy bands of PSR~J0540-6919 is observed with the efficiency $η_{J0540}\sim 0.06$ (in 4… ▽ More We discuss X-ray and gamma-ray emissions from Crab-like pulsars, PSRs~J0537-6910 and~J0540-6919, in Large Magellanic Cloud. Fermi-LAT observations have resolved the gamma-ray emissions from these two pulsars and found the pulsed emissions from PSR~J0540-6919. The total pulsed radiation in the X-ray/gamma-ray energy bands of PSR~J0540-6919 is observed with the efficiency $η_{J0540}\sim 0.06$ (in 4$π$ sr), which is about a factor of ten larger than $η_{Crab}\sim 0.006$ of the Crab pulsar. Although PSR~J0537-6910 has the highest spin-down power among currently known pulsars, the efficiency of the observed X-ray emissions is about two orders of magnitude smaller than that of PSR~J0540-6919. This paper mainly discusses what causes the difference in the radiation efficiencies of these three energetic Crab-like pulsars. We discuss electron/positron acceleration and high-energy emission processes within the outer gap model. By solving the outer gap structure with the dipole magnetic field, we show that the radiation efficiency decreases as the inclination angle between the magnetic axis and the rotation axis increases. To explain the difference in the pulse profile and in the radiation efficiency, our model suggests that PSR~J0540-6919 has an inclination angle much smaller than the that of Crab pulsar (here we assume the inclination angles of both pulsars are $α<90^{\circ}$). On the other hand, we speculate that the difference in the radiation efficiencies between PSRs~J0537-6910 and J0549-6919 is mainly caused by the difference in the Earth viewing angle, and that we see PSR~J0537-6910 with an Earth viewing angle $ζ>>90^{\circ}$ (or $<<90^{\circ}$) measured from the spin axis, while we see PSR~J0540-6919 with $ζ\sim 90^{\circ}$. △ Less

Submitted 1 December, 2016; originally announced December 2016.

Comments: 23 pages, 7 figures, ApJ in press

Differences between radio-loud and radio-quiet gamma-ray pulsars as revealed by Fermi

Authors: C. Y. Hui, Jongsu Lee, J. Takata, C. W. Ng, K. S. Cheng

Abstract: By comparing the properties of non-recycled radio-loud $γ-$ray pulsars and radio-quiet $γ-$ray pulsars, we have searched for the differences between these two populations. We found that the $γ-$ray spectral curvature of radio-quiet pulsars can be larger than that of radio-loud pulsars. Based on the full sample of non-recycled $γ-$ray pulsars, their distributions of the magnetic field strength at t… ▽ More By comparing the properties of non-recycled radio-loud $γ-$ray pulsars and radio-quiet $γ-$ray pulsars, we have searched for the differences between these two populations. We found that the $γ-$ray spectral curvature of radio-quiet pulsars can be larger than that of radio-loud pulsars. Based on the full sample of non-recycled $γ-$ray pulsars, their distributions of the magnetic field strength at the light cylinder are also found to be different. We notice that this might be resulted from the observational bias. In re-examining the previously reported difference of $γ-$ray-to-X-ray flux ratios, we found the significance can be hampered by their statistical uncertainties. In the context of outer gap model, we discuss the expected properties of these two populations and compare with the possible differences identified in our analysis. △ Less

Submitted 22 November, 2016; originally announced November 2016.

Comments: 17 pages, 4 Figures, 4 Tables, accepted for publication in ApJ

Constraining warm dark matter mass with cosmic reionization and gravitational wave

Authors: W. W. Tan, F. Y. Wang, K. S. Cheng

Abstract: We constrain the warm dark matter (WDM) particle mass with the observations of cosmic reionization and CMB optical depth. We suggest that the GWs from stellar mass black holes (BHs) could give a further constraint on WDM particle mass for future observations. The star formation rates (SFRs) of Population I/II (Pop I/II) and Population III (Pop III) stars are also derived. If the metallicity of the… ▽ More We constrain the warm dark matter (WDM) particle mass with the observations of cosmic reionization and CMB optical depth. We suggest that the GWs from stellar mass black holes (BHs) could give a further constraint on WDM particle mass for future observations. The star formation rates (SFRs) of Population I/II (Pop I/II) and Population III (Pop III) stars are also derived. If the metallicity of the universe have been enriched beyond the critical value of $Z_{\rm crit}=10^{-3.5}Z_{\odot}$, the star formation shift from Pop III to Pop I/II stars. Our results show that the SFRs are quite dependent on the WDM particle mass, especially at high redshifts. Combing with the reionization history and CMB optical depth derived from the recent \emph{Planck} mission, we find that the current data requires the WDM particle mass in a narrow range of $1 < m_x < 3$ keV. Furthermore, we suggest that the stochastic gravitational wave background (SGWB) produced by stellar BHs could give a further constraint on the WDM particle mass for future observations. For $m_{x}=3$ keV with Salpeter (Chabrier) initial mass function (IMF), the SGWB from Pop I/II BHs has a peak amplitude of $Ω_{\rm GW}\approx2.8\times 10^{-9}~(5.0\times 10^{-9})$ at $f= 316 {\rm Hz}$, while the GW radiation at $f<10$Hz is seriously suppressed. For $m_{\rm x}=1$ keV, the SGWB peak amplitude is the same as that of $m_{\rm x}=1$ keV, but a little lower at low frequencies. Therefore, it is hard to constrain the WDM particle mass by the SGWB from Pop I/II BHs. To assess the detectability of GW signal, we also calculate the signal to noise ratio (SNR), which are $\rm SNR=37.7~ (66.5)$ and $27~(47.7)$ for $m_{\rm x}=3$ keV and $m_{\rm x}=1$ keV for Einstein Telescope (ET) with Salpeter (Chabrier) IMF, respectively. △ Less

Submitted 12 July, 2016; originally announced July 2016.

Comments: 11 pages, 5 figures, accepted for publication in ApJ

Studying the SGR 1806-20/Cl* 1806-20 region using the \emph{Fermi} Large Area Telescope

Authors: Paul K. H. Yeung, Albert K. H. Kong, P. H. Thomas Tam, Lupin C. C. Lin, C. Y. Hui, Chin-Ping Hu, K. S. Cheng

Abstract: The region around SGR 1806-20 and its host stellar cluster Cl* 1806-20 is a potentially important site of particle acceleration. The soft $γ-$ray repeater and Cl* 1806-20, which also contains several very massive stars including a luminous blue variable hypergiant LBV 1806-20, are capable of depositing a large amount of energy to the surroundings. Using the data taken with the \emph{Fermi} Large A… ▽ More The region around SGR 1806-20 and its host stellar cluster Cl* 1806-20 is a potentially important site of particle acceleration. The soft $γ-$ray repeater and Cl* 1806-20, which also contains several very massive stars including a luminous blue variable hypergiant LBV 1806-20, are capable of depositing a large amount of energy to the surroundings. Using the data taken with the \emph{Fermi} Large Area Telescope (LAT), we identified an extended LAT source to the south-west of Cl* 1806-20. The centroid of the 1-50~GeV emission is consistent with that of HESS J1808-204 (until now unidentified). The LAT spectrum is best-fit by a broken power-law with the break energy $E_\mathrm{b}=297\pm15$ MeV. The index above $E_\mathrm{b}$ is $2.60\pm0.04$, and is consistent with the flux and spectral index above 100 GeV for HESS J1808-204, suggesting an association between the two sources. Meanwhile, the interacting supernova remnant SNR G9.7-0.0 is also a potential contributor to the LAT flux. A tentative flux enhancement at the MeV band during a 45-day interval (2011 Jan 21 - 2011 Mar 7) is also reported. We discuss possible origins of the extended LAT source in the context of both leptonic and hadronic scenarios. △ Less

Submitted 6 June, 2016; originally announced June 2016.

Comments: Accepted for publication in ApJ

Observation and Simulation of the Variable Gamma-ray Emission from PSR~J2021+4026

Authors: C. W. Ng, J. Takata, K. S. Cheng

Abstract: Pulsars are rapidly spinning and highly magnetized neutron stars, with highly stable rotational period and gradual spin-down over a long timescale due to the loss of radiation. Glitches refer to the events that suddenly increase the rotational speed of a pulsar. The exact causes of glitches and the resulting processes are not fully understood. It is generally believed that couplings between the no… ▽ More Pulsars are rapidly spinning and highly magnetized neutron stars, with highly stable rotational period and gradual spin-down over a long timescale due to the loss of radiation. Glitches refer to the events that suddenly increase the rotational speed of a pulsar. The exact causes of glitches and the resulting processes are not fully understood. It is generally believed that couplings between the normal matter and the superfluid components, and the starquakes, are the common causes of glitches. In this study, one famous glitching pulsar, PSR~J2021+4026, is investigated. PSR~J2021+4026 is the first variable gamma-ray pulsar observed by Fermi. From the gamma-ray observations, it is found that the pulsar experienced a significant flux drop, an increase in the spin-down rate, a change in the pulse profile and a shift in the spectral cut-off to a lower energy, simultaneously around 2011 October 16. To explain these effects on the high-energy emissions by the glitch of PSR~J2021+4026, we hypothesized the glitch to be caused by the rearrangement of the surface magnetic field due to the crustal plate tectonic activities on the pulsar which is triggered by a starquake. In this glitch event, the inclination angle of the magnetic dipole axis is slightly shifted. This proposition is then tested by numerical modeling using a three-dimensional two-layer outer gap model. The simulation results indicate that a modification of the inclination angle can affect the pulse profile and the spectral properties, which can explain the observation changes after the glitch. △ Less

Submitted 29 May, 2016; originally announced May 2016.

Comments: 20 pages, 6 figures. Accepted for publication in ApJ

Observing two dark accelerators around the Galactic Centre with Fermi Large Area Telescope

Authors: C. Y. Hui, P. K. H. Yeung, C. W. Ng, L. C. C. Lin, P. H. T. Tam, K. S. Cheng, A. K. H. Kong, D. O. Chernyshov, V. A. Dogiel

Abstract: We report the results from a detailed $γ-$ray investigation in the field of two "dark accelerators", HESS J1745-303 and HESS J1741-302, with $6.9$ years of data obtained by the Fermi Large Area Telescope. For HESS J1745-303, we found that its MeV-GeV emission is mainly originated from the "Region A" of the TeV feature. Its $γ-$ray spectrum can be modeled with a single power-law with a photon index… ▽ More We report the results from a detailed $γ-$ray investigation in the field of two "dark accelerators", HESS J1745-303 and HESS J1741-302, with $6.9$ years of data obtained by the Fermi Large Area Telescope. For HESS J1745-303, we found that its MeV-GeV emission is mainly originated from the "Region A" of the TeV feature. Its $γ-$ray spectrum can be modeled with a single power-law with a photon index of $Γ\sim2.5$ from few hundreds MeV to TeV. Moreover, an elongated feature, which extends from "Region A" toward northwest for $\sim1.3^{\circ}$, is discovered for the first time. The orientation of this feature is similar to that of a large scale atomic/molecular gas distribution. For HESS J1741-302, our analysis does not yield any MeV-GeV counterpart for this unidentified TeV source. On the other hand, we have detected a new point source, Fermi J1740.1-3013, serendipitously. Its spectrum is apparently curved which resembles that of a $γ-$ray pulsar. This makes it possibly associated with PSR B1737-20 or PSR J1739-3023. △ Less

Submitted 25 January, 2016; originally announced January 2016.

Comments: 11 pages, 7 figures, 2 tables, accepted for publication in MNRAS

X-Ray afterglow of SWIFT J1644+57: a Compton echo?

Authors: K. S. Cheng, D. O. Chernyshov, V. A. Dogiel, Albert K. H. Kong, C. M. Ko

Abstract: Swift, Chandra and XMM have found a weak but nearly constant X-ray component from Swift J1644+57 that appeared at ~500 days and was visible at least until ~ 1400 days after the stellar capture, which cannot be explained by standard tidal disruption theories. We suggest that this X-ray afterglow component may result from Thomson scattering between the primary X-rays and its surrounding plasma, i.e.… ▽ More Swift, Chandra and XMM have found a weak but nearly constant X-ray component from Swift J1644+57 that appeared at ~500 days and was visible at least until ~ 1400 days after the stellar capture, which cannot be explained by standard tidal disruption theories. We suggest that this X-ray afterglow component may result from Thomson scattering between the primary X-rays and its surrounding plasma, i.e. the Compton echo effect. Similar phenomena has also been observed from molecular clouds in our Galactic Center, which were caused by the past activity of Srg A*. If this interpretation of Swift J1644+57 afterglow is correct, this is the first Compton Echo effect observed in the cosmological distances. △ Less

Submitted 18 December, 2015; v1 submitted 15 December, 2015; originally announced December 2015.

Comments: To be published in ApJL, accepted

Probing gamma-ray emissions of Fermi-LAT pulsars with a non-stationary outer gap model

Authors: J. Takata, C. W. Ng, K. S. Cheng

Abstract: We explore a non-stationary outer gap scenario for gamma-ray emission process in pulsar magnetosphere. Electrons/positrons that migrate along the magnetic field line and enter the outer gap from the outer/inner boundaries activate the pair-creation cascade and high-energy emission process. In our model, the rate of the particle injection at the gap boundaries is key physical quantity to control th… ▽ More We explore a non-stationary outer gap scenario for gamma-ray emission process in pulsar magnetosphere. Electrons/positrons that migrate along the magnetic field line and enter the outer gap from the outer/inner boundaries activate the pair-creation cascade and high-energy emission process. In our model, the rate of the particle injection at the gap boundaries is key physical quantity to control the gap structure and properties of the gamma-ray spectrum. Our model assumes that the injection rate is time variable and the observed gamma-ray spectrum are superposition of the emissions from different gap structures with different injection rates at the gap boundaries. The calculated spectrum superposed by assuming power law distribution of the particle injection rate can reproduce sub-exponential cut-off feature in the gamma-ray spectrum observed by Fermi-LAT. We fit the phase-averaged spectra for 43 young/middle-age pulsars and 14 millisecond pulsars with the model. Our results imply that (1) a larger particle injection at the gap boundaries is more frequent for the pulsar with a larger spin down power and (2) outer gap with an injection rate much smaller than the Goldreich-Julian value produces observe $>10$GeV emissions. Fermi-LAT gamma-ray pulsars show that (i) the observed gamma-ray spectrum below cut-off energy tends to be softer for the pulsar with a higher spin down rate and (ii) the second peak is more prominent in higher energy bands. Based on the results of the fitting, we describe possible theoretical interpretations for these observational properties. We also briefly discuss Crab-like millisecond pulsars that show phase-aligned radio and gamma-ray pulses. △ Less

Submitted 20 November, 2015; originally announced November 2015.

Comments: 38 pages, accepted for publication in MNRAS

Measuring dark energy with the $E_{\rm iso}-E_{\rm p}$ correlation of gamma-ray bursts using model-independent methods

Authors: J. S. Wang, F. Y. Wang, K. S. Cheng, Z. G. Dai

Abstract: In this paper, we use two model-independent methods to standardize long gamma-ray bursts (GRBs) using the $E_{\rm iso}-E_{\rm p}$ correlation, where $E_{\rm iso}$ is the isotropic-equivalent gamma-ray energy and $E_{\rm p}$ is the spectral peak energy. We update 42 long GRBs and try to make constraint on cosmological parameters. The full sample contains 151 long GRBs with redshifts from 0.0331 to… ▽ More In this paper, we use two model-independent methods to standardize long gamma-ray bursts (GRBs) using the $E_{\rm iso}-E_{\rm p}$ correlation, where $E_{\rm iso}$ is the isotropic-equivalent gamma-ray energy and $E_{\rm p}$ is the spectral peak energy. We update 42 long GRBs and try to make constraint on cosmological parameters. The full sample contains 151 long GRBs with redshifts from 0.0331 to 8.2. The first method is the simultaneous fitting method. The extrinsic scatter $σ_{\rm ext}$ is taken into account and assigned to the parameter $E_{\rm iso}$. The best-fitting values are $a=49.15\pm0.26$, $b=1.42\pm0.11$, $σ_{\rm ext}=0.34\pm0.03$ and $Ω_m=0.79$ in the flat $Λ$CDM model. The constraint on $Ω_m$ is $0.55<Ω_m<1$ at the 1$σ$ confidence level. If reduced $χ^2$ method is used, the best-fit results are $a=48.96\pm0.18$, $b=1.52\pm0.08$ and $Ω_m=0.50\pm0.12$. The second method is using type Ia supernovae (SNe Ia) to calibrate the $E_{\rm iso}-E_{\rm p}$ correlation. We calibrate 90 high-redshift GRBs in the redshift range from 1.44 to 8.1. The cosmological constraints from these 90 GRBs are $Ω_m=0.23^{+0.06}_{-0.04}$ for flat $Λ$CDM, and $Ω_m=0.18\pm0.11$ and $Ω_Λ=0.46\pm0.51$ for non-flat $Λ$CDM. For the combination of GRB and SNe Ia sample, we obtain $Ω_m=0.271\pm0.019$ and $h=0.701\pm0.002$ for the flat $Λ$CDM, and for the non-flat $Λ$CDM, the results are $Ω_m=0.225\pm0.044$, $Ω_Λ=0.640\pm0.082$ and $h=0.698\pm0.004$. These results from calibrated GRBs are consistent with that of SNe Ia. Meanwhile, the combined data can improve cosmological constraints significantly, comparing to SNe Ia alone. Our results show that the $E_{\rm iso}-E_{\rm p}$ correlation is promising to probe the high-redshift universe. △ Less

Submitted 22 October, 2015; v1 submitted 28 September, 2015; originally announced September 2015.

Comments: 10 pages, 6 figures, 4 table, accepted by A&A. Table 4 contains calibrated distance moduli of GRBs

Journal ref: A&A 585, A68 (2016)

Searches for millisecond pulsar candidates among the unidentified Fermi objects

Authors: C. Y. Hui, S. M. Park, C. P. Hu, L. C. C. Lin, K. L. Li, A. K. H. Kong, P. H. T. Tam, J. Takata, K. S. Cheng, Ruolan Jin, T. -C. Yen, Chunglee Kim

Abstract: Here we report the results of searching millisecond pulsar (MSP) candidates from the Fermi LAT second source catalog (2FGL). Seven unassociated $γ-$ray sources in this catalog are identified as promising MSP candidates based on their $γ$-ray properties. Through the X-ray analysis, we have detected possible X-ray counterparts, localized to an arcsecond accuracy. We have systematically estimated the… ▽ More Here we report the results of searching millisecond pulsar (MSP) candidates from the Fermi LAT second source catalog (2FGL). Seven unassociated $γ-$ray sources in this catalog are identified as promising MSP candidates based on their $γ$-ray properties. Through the X-ray analysis, we have detected possible X-ray counterparts, localized to an arcsecond accuracy. We have systematically estimated their X-ray fluxes and compared with the corresponding $γ$-ray fluxes. The X-ray to $γ$-ray flux ratios for 2FGL J1653.6-0159 and 2FGL J1946.4-5402 are comparable with the typical value for pulsars. For 2FGL J1625.2-0020, 2FGL J1653.6-0159 and 2FGL J1946.4-5402, their candidate X-ray counterparts are bright enough for performing a detailed spectral and temporal analysis to discriminate their thermal/non thermal nature and search for the periodic signal. We have also searched for possible optical/IR counterparts at the X-ray positions. For the optical/IR source coincident with the brightest X-ray object that associated with 2FGL J1120.0-2204, its spectral energy distribution is comparable with a late-type star. Evidence for the variability has also been found by examining its optical light curve. All the aforementioned 2FGL sources resemble a pulsar in one or more aspects, which make them as the promising targets for follow-up investigations. △ Less

Submitted 9 July, 2015; originally announced July 2015.

Comments: 37 pages, 12 figures, 6 tables, accepted for publication in ApJ

Spectrum of Relativistic and Subrelativistic Cosmic Rays in the 100 pc Central Region

Authors: V. A. Dogiel, D. O. Chernyshov, A. M. Kiselev, M. Nobukawa, K. S. Cheng, C. Y. Hui, C. M. Ko, K. K. Nobukawa, T. G. Tsuru

Abstract: From the rate of hydrogen ionization and the gamma ray flux, we derived the spectrum of relativistic and subrelativistic cosmic rays (CRs) nearby and inside the molecular cloud Sgr B2 near the Galactic Center (GC). We studied two cases of CR propagation in molecular clouds: free propagation and scattering of particles by magnetic fluctuations excited by the neutral gas turbulence. We showed that i… ▽ More From the rate of hydrogen ionization and the gamma ray flux, we derived the spectrum of relativistic and subrelativistic cosmic rays (CRs) nearby and inside the molecular cloud Sgr B2 near the Galactic Center (GC). We studied two cases of CR propagation in molecular clouds: free propagation and scattering of particles by magnetic fluctuations excited by the neutral gas turbulence. We showed that in the latter case CR propagation inside the cloud can be described as diffusion with the coefficient $\sim 3\times 10^{27}$ cm$^2$ s$^{-1}$. For the case of hydrogen ionization by subrelativistic protons, we showed that their spectrum outside the cloud is quite hard with the spectral index $δ>-1$. The energy density of subrelativistic protons ($>50$ eV cm$^{-3}$) is one order of magnitude higher than that of relativistic CRs. These protons generate the 6.4 keV emission from Sgr B2, which was about 30\% of the flux observed by Suzaku in 2013. Future observations for the period after 2013 may discover the background flux generated by subrelativistic CRs in Sgr B2. Alternatively hydrogen ionization of the molecular gas in Sgr B2 may be caused by high energy electrons. We showed that the spectrum of electron bremsstrahlung is harder than the observed continuum from Sgr B2, and in principle this X-ray component provided by electrons could be seen from the INTEGRAL data as a stationary high energy excess above the observed spectrum $E_x^{-2}$. △ Less

Submitted 9 July, 2015; originally announced July 2015.

Comments: 42 pages, 6 figures, accepted by ApJ

An unexpectedly low-redshift excess of Swift gamma-ray burst rate

Authors: H. Yu, F. Y. Wang, Z. G. Dai, K. S. Cheng

Abstract: Gamma-ray bursts (GRBs) are the most violent explosions in the Universe and can be used to explore the properties of high-redshift universe. It is believed that the long GRBs are associated with the deaths of massive stars. So it is possible to use GRBs to investigate the star formation rate (SFR). In this paper, we use Lynden-Bell's $c^-$ method to study the luminosity function and rate of \emph{… ▽ More Gamma-ray bursts (GRBs) are the most violent explosions in the Universe and can be used to explore the properties of high-redshift universe. It is believed that the long GRBs are associated with the deaths of massive stars. So it is possible to use GRBs to investigate the star formation rate (SFR). In this paper, we use Lynden-Bell's $c^-$ method to study the luminosity function and rate of \emph{Swift} long GRBs without any assumptions. We find that the luminosity of GRBs evolves with redshift as $L(z)\propto g(z)=(1+z)^k$ with $k=2.43_{-0.38}^{+0.41}$. After correcting the redshift evolution through $L_0(z)=L(z)/g(z)$, the luminosity function can be expressed as $ψ(L_0)\propto L_0^{-0.14\pm0.02}$ for dim GRBs and $ψ(L_0)\propto L_0^{-0.70\pm0.03}$ for bright GRBs, with the break point $L_{0}^{b}=1.43\times10^{51}~{\rm erg~s^{-1}}$. We also find that the formation rate of GRBs is almost constant at $z<1.0$ for the first time, which is remarkably different from the SFR. At $z>1.0$, the formation rate of GRB is consistent with the SFR. Our results are dramatically different from previous studies. Some possible reasons for this low-redshift excess are discussed. We also test the robustness of our results with Monte Carlo simulations. The distributions of mock data (i.e., luminosity-redshift distribution, luminosity function, cumulative distribution and $\log N-\log S$ distribution) are in good agreement with the observations. Besides, we also find that there are remarkable difference between the mock data and the observations if long GRB are unbiased tracers of SFR at $z<1.0$. △ Less

Submitted 10 April, 2015; v1 submitted 7 April, 2015; originally announced April 2015.

Comments: 33 pages, 10 figures, 1 table, accepted by ApJS

Multi-wavelength Emission from the Fermi Bubble III. Stochastic (Fermi) Re-Acceleration of Relativistic Electrons Emitted by SNRs

Authors: K. S. Cheng, D. O. Chernyshov, V. A. Dogiel, C. M. Ko

Abstract: We analyse the model of stochastic re-acceleration of electrons, which are emitted by supernova remnants (SNRs) in the Galactic Disk and propagate then into the Galactic halo, in order to explain the origin on nonthermal (radio and gamma-ray) emission from the Fermi Bubbles (FB). We assume that the energy for re-acceleration in the halo is supplied by shocks generated by processes of star accretio… ▽ More We analyse the model of stochastic re-acceleration of electrons, which are emitted by supernova remnants (SNRs) in the Galactic Disk and propagate then into the Galactic halo, in order to explain the origin on nonthermal (radio and gamma-ray) emission from the Fermi Bubbles (FB). We assume that the energy for re-acceleration in the halo is supplied by shocks generated by processes of star accretion onto the central black hole. Numerical simulations show that regions with strong turbulence (places for electron re-acceleration) are located high up in the Galactic Halo about several kpc above the disk. The energy of SNR electrons that reach these regions does not exceed several GeV because of synchrotron and inverse Compton energy losses. At appropriate parameters of re-acceleration these electrons can be re-accelerated up to the energy 10E12 eV which explains in this model the origin of the observed radio and gamma-ray emission from the FB. However although the model gamma-ray spectrum is consistent with the Fermi results, the model radio spectrum is steeper than the observed by WMAP and Planck. If adiabatic losses due to plasma outflow from the Galactic central regions are taken into account, then the re-acceleration model nicely reproduces the Planck datapoints. △ Less

Submitted 12 March, 2015; originally announced March 2015.

Comments: 33 pages, 8 figures, accepted by ApJ

Exploring the intrabinary shock from the redback millisecond pulsar PSR J2129-0429

Authors: C. Y. Hui, C. P. Hu, S. M. Park, J. Takata, K. L. Li, P. H. T. Tam, L. C. C. Lin, A. K. H. Kong, K. S. Cheng, Chunglee Kim

Abstract: We have investigated the intrabinary shock emission from the redback millisecond pulsar PSR J2129-0429 with XMM-Newton and Fermi. Orbital modulation in X-ray and UV can be clearly seen. Its X-ray modulation has a double-peak structure with a dip in between. The observed X-rays are non-thermal dominant which can be modeled by a power-law with a photon index of ~1.2. Intrabinary shock can be the ori… ▽ More We have investigated the intrabinary shock emission from the redback millisecond pulsar PSR J2129-0429 with XMM-Newton and Fermi. Orbital modulation in X-ray and UV can be clearly seen. Its X-ray modulation has a double-peak structure with a dip in between. The observed X-rays are non-thermal dominant which can be modeled by a power-law with a photon index of ~1.2. Intrabinary shock can be the origin of the observed X-rays. The UV light curve is resulted from the ellipsoidal modulation of the companion. Modeling the UV light curve prefers a large viewing angle. The heating effect of the UV light curve is found to be negligible which suggests the high energy radiation beam of PSR J2129-0429 does not direct toward its companion. On the other hand, no significant orbital modulation can be found in gamma-ray which suggests the majority of the gamma-rays come from the pulsar. △ Less

Submitted 13 February, 2015; originally announced February 2015.

Comments: 14 pages, 3 figures, 1 table, accepted for publication in ApJ Letter

Constraining the bulk Lorentz factor from the photosphere emission

Authors: Yuan-Chuan Zou, K. S. Cheng, F. Y. Wang

Abstract: We propose a direct and model-independent method to constrain the Lorentz factor of a relativistically expanding object, like gamma-ray bursts. Only the measurements, such as thermal component of the emission, the distance and the variable time scale of the light curve, are used. If the uncertainties are considered, we will obtain lower limits of the Lorentz factor instead. We apply this method to… ▽ More We propose a direct and model-independent method to constrain the Lorentz factor of a relativistically expanding object, like gamma-ray bursts. Only the measurements, such as thermal component of the emission, the distance and the variable time scale of the light curve, are used. If the uncertainties are considered, we will obtain lower limits of the Lorentz factor instead. We apply this method to GRB 090618 and get a lower limit of the Lorentz factor to be 22. The method can be used to any relativistically moving object, such as gamma-ray bursts, blazars, and soft gamma-ray repeaters, providing the thermal component of the emission being observed. △ Less

Submitted 25 January, 2015; originally announced January 2015.

Comments: 10 pages, 1 figure

The birth of radio millisecond pulsars and their high-energy signature

Authors: P. H. T. Tam, K. L. Li, A. K. H. Kong, J. Takata, G. C. K. Leung, K. S. Cheng, C. Y. Hui

Abstract: Millisecond pulsars (MSPs) are thought to born in low-mass X-ray binaries when the neutron star has gained enough angular momentum from the accreting materials of its companion star. It is generally believed that a radio MSP is born when the neutron star stops accreting and enters a rotation-powered state. Exactly what happens during the transition time was poorly understood until a year ago. In t… ▽ More Millisecond pulsars (MSPs) are thought to born in low-mass X-ray binaries when the neutron star has gained enough angular momentum from the accreting materials of its companion star. It is generally believed that a radio MSP is born when the neutron star stops accreting and enters a rotation-powered state. Exactly what happens during the transition time was poorly understood until a year ago. In the past year, observations have revealed a few objects that not only switched from one state to the other (as predicted in the above picture), but also have swung between the two states within weeks to years. In this work, we present observations of two of these transition objects (PSR J1023+0038 and XSS J12270-4859) and a theoretical framework that tries to explain their high-energy radiation. △ Less

Submitted 16 December, 2014; originally announced December 2014.

Comments: 7 pages, 6 figures; presented at the 10th INTEGRAL Workshop: 'A Synergistic View of the High Energy Sky' - Integral2014, held 15-19 September 2014, Annapolis, MD

Origin of Hydrogen Ionization in the 1 pc Galactic Central Region

Authors: D. O. Chernyshov, Gene C. K. Leung, K. S. Cheng, V. A. Dogiel, V. Tatischeff

Abstract: We study a possible connection between processes of gamma-ray emission and hydrogen ionization in a few pc of central region around Sgr A*. Previous investigations showed there is a discrepancy between interpretation of gamma-ray and ionization data if gamma-rays are generated by proton-proton collisions. Here we provided analysis of processes of ionization and emission basing on analytical and nu… ▽ More We study a possible connection between processes of gamma-ray emission and hydrogen ionization in a few pc of central region around Sgr A*. Previous investigations showed there is a discrepancy between interpretation of gamma-ray and ionization data if gamma-rays are generated by proton-proton collisions. Here we provided analysis of processes of ionization and emission basing on analytical and numerical calculations of kinetic equations which describe processes of particle propagation and their energy losses. The origin of gamma rays could be either due to collisions of relativistic protons with the dense gas of the surrounding circumnuclear disk (CND) or bremsstrahlung and inverse Compton scattering of relativistic electrons. The hydrogen ionization in this case is produced by a low energy component of the CR spectrum. We found that if ionization is produced by protons the expected ionization rate of hydrogen in the CND is of the same order as derived from IR observations. So we do not see any discrepancy between the gamma-ray and ionization data for the hadronic model. In the case of ionization by electrons we obtained the ionization rate one order of magnitude higher than follows from the IR data. In principle, a selection between the leptonic and hadronic interpretations can be performed basing on measurements of radio and X-ray fluxes from this region because the leptonic and hadronic models give different values of the fluxes from there. We do not exclude that gamma-ray production and hydrogen ionization in the CND are due to a past activity of Sgr A* which occurred about 100 year ago. Then we hypothesize that there may be connection between a past proton eruption and a flux of hard X-rays emitted by Sgr A* hundred years ago as follows from the observed time variability of the iron line seen in the direction of GC molecular clouds. △ Less

Submitted 8 December, 2014; originally announced December 2014.

Comments: 21 pages, 10 figures. Submitted

Fermi bubbles as a source of cosmic rays above 10^{15} eV

Authors: D. O. Chernyshov, K. S. Cheng, V. A. Dogiel, C. M. Ko

Abstract: Fermi bubbles are giant gamma-ray structures extended north and south of the Galactic center with characteristic sizes of order of 10 kpc recently discovered by Fermi Large Area Telescope. Good correlation between radio and gamma-ray emission in the region covered by Fermi bubbles implies the presence of high-energy electrons in this region. Since it is relatively difficult for relativistic electr… ▽ More Fermi bubbles are giant gamma-ray structures extended north and south of the Galactic center with characteristic sizes of order of 10 kpc recently discovered by Fermi Large Area Telescope. Good correlation between radio and gamma-ray emission in the region covered by Fermi bubbles implies the presence of high-energy electrons in this region. Since it is relatively difficult for relativistic electrons of this energy to travel all the way from the Galactic sources toward Fermi bubbles one can assume that they accelerated in-situ. The corresponding acceleration mechanism should also affect the distribution of the relativistic protons in the Galaxy. Since protons have much larger lifetimes the effect may even be observed near the Earth. In our model we suggest that Fermi bubbles are created by acceleration of electrons on series of shocks born due to periodic star accretions by supermassive black hole Sgr A*. We propose that hadronic CR within the "knee" of the observed CR spectrum are produced by Galactic supernova remnants distributed in the Galactic disk. Reacceleration of these particles in the Fermi Bubble produces CRs beyond the knee. This model provides a natural explanation of the observed CR flux, spectral indexes, and matching of spectra at the knee. △ Less

Submitted 8 December, 2014; originally announced December 2014.

Comments: 7 pages, 3 figures. Published in Nuclear Physics B, oral talk at the conference "Cosmic ray origin - beyond the standard models" San-Vito di Cadore, Italy, March 16-22, 2014

Journal ref: 2014, Nuclear Physics B (Proceedings Supplements), 2014, pp. 179-185

A detailed X-ray investigation of PSR J2021+4026 and $γ$-Cygni supernova remnant

Authors: C. Y. Hui, K. A. Seo, L. C. C. Lin, R. H. H. Huang, C. P. Hu, J. H. K. Wu, L. Trepl, J. Takata, Y. Wang, Y. Chou, K. S. Cheng, A. K. H. Kong

Abstract: We have investigated the field around the radio-quiet $γ$-ray pulsar, PSR J2021+4026, with a ~140 ks XMM-Newton observation and a ~56 ks archival Chandra data. Through analyzing the pulsed spectrum, we show that the X-ray pulsation is purely thermal in nature which suggests the pulsation is originated from a hot polar cap with $T\sim3\times10^{6}$ K on the surface of a rotating neutron star. On th… ▽ More We have investigated the field around the radio-quiet $γ$-ray pulsar, PSR J2021+4026, with a ~140 ks XMM-Newton observation and a ~56 ks archival Chandra data. Through analyzing the pulsed spectrum, we show that the X-ray pulsation is purely thermal in nature which suggests the pulsation is originated from a hot polar cap with $T\sim3\times10^{6}$ K on the surface of a rotating neutron star. On the other hand, the power-law component that dominates the pulsar emission in the hard band is originated from off-pulse phases, which possibly comes from a pulsar wind nebula. In re-analyzing the Chandra data, we have confirmed the presence of bow-shock nebula which extends from the pulsar to west by ~10 arcsec. The orientation of this nebular feature suggests that the pulsar is probably moving eastward which is consistent with the speculated proper motion by extrapolating from the nominal geometrical center of the supernova remnant (SNR) G78.2+2.1 to the current pulsar position. For G78.2+2.1, our deep XMM-Newton observation also enables a study of the central region and part of the southeastern region with superior photon statistics. The column absorption derived for the SNR is comparable with that for PSR J2021+4026, which supports their association. The remnant emission in both examined regions are in an non-equilibrium ionization state. Also, the elapsed time of both regions after shock-heating is apparently shorter than the Sedov age of G78.2+2.1. This might suggest the reverse shock has reached the center not long ago. Apart from PSR J2021+4026 and G78.2+2.1, we have also serendipitously detected an X-ray flash-like event XMM J202154.7+402855 from this XMM-Newton observation. △ Less

Submitted 9 November, 2014; originally announced November 2014.

Comments: 45 pages, 10 figures, 2 tables, accepted for publication in ApJ

Fermi-LAT Detection of Pulsed Gamma-rays Above 50 GeV from the Vela Pulsar

Authors: Gene C. K. Leung, J. Takata, C. W. Ng, A. K. H Kong, P. H. T. Tam, C. Y. Hui, K. S. Cheng

Abstract: The First \fermi-LAT Catalog of Sources Above 10 GeV reported evidence of pulsed emission above 25 GeV from 12 pulsars, including the Vela pulsar, which showed evidence of pulsation at $>37$ GeV energy bands. Using 62 months of \fermi-LAT data, we analyzed the gamma-ray emission from the Vela pulsar and searched for pulsed emission above 50 GeV. Having confirmed the significance of the pulsation i… ▽ More The First \fermi-LAT Catalog of Sources Above 10 GeV reported evidence of pulsed emission above 25 GeV from 12 pulsars, including the Vela pulsar, which showed evidence of pulsation at $>37$ GeV energy bands. Using 62 months of \fermi-LAT data, we analyzed the gamma-ray emission from the Vela pulsar and searched for pulsed emission above 50 GeV. Having confirmed the significance of the pulsation in 30-50 GeV with the H-test (p-value $\sim10^{-77}$), we extracted its pulse profile using the Bayesian block algorithm and compared it with the distribution of the 5 observed photons above 50 GeV using the likelihood ratio test. Pulsation was significantly detected for photons above 50 GeV with p-value $=3\times10^{-5}$ ($4.2σ$). The detection of pulsation is significant above $4σ$ at $>79$ GeV and above $3σ$ at $>90$ GeV energy bands, making this the highest energy pulsation significantly detected by the LAT. We explore non-stationary outer gap scenario of the very high-energy emissions from the Vela pulsar. △ Less

Submitted 21 October, 2014; v1 submitted 20 October, 2014; originally announced October 2014.

Comments: 14 pages, 4figures, accepted by ApJ Letter