Source author record

J. Takata

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

61works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

61 published item(s)

preprint2026arXiv

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

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

preprint2026arXiv

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

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

preprint2026arXiv

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

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

preprint2026arXiv

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

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

preprint2026arXiv

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

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

preprint2026arXiv

Evidence for particle acceleration approaching PeV energies in the W51 complex

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

preprint2026arXiv

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

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

preprint2026arXiv

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

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

preprint2026arXiv

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

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

preprint2026arXiv

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

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

preprint2026arXiv

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

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

preprint2025arXiv

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

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

preprint2025arXiv

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

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

preprint2022arXiv

Searching for cataclysmic variable stars in unidentified X-ray sources

We carry out a photometric search for new cataclysmic variable stars (CVs), with the goal of identification for candidates of AR~Scorpii-type binary systems. We select GAIA sources that are likely associated with unidentified X-ray sources, and analyze the light curves taken by the Zwicky Transient Facility, Transiting Exoplanet Survey Satellite, and Lulin One-meter Telescope in Taiwan. We investigate eight sources as candidates for CVs, among which six sources are new identifications. Another two sources have been recognized as CVs in previous studies, but no detailed investigations have been done. We identify two eclipsing systems that are associated with an unidentified XMM-Newton or Swift source, and one promising candidate for polar associated with an unidentified ASKA source. Two polar candidates may locate in the so-called period gap of a CV, and the other six candidates have an orbital period shorter than that of the period gap. Although we do not identify a promising candidate for AR~Scorpii-type binary systems, our study suggests that CV systems that have X-ray emission and do not show frequent outbursts may have been missed in previous surveys.

preprint2021arXiv

Modelling the correlated keV/TeV light curves of Be/gamma-ray binaries

Be/$γ$-ray binaries comprise a confirmed or presumptive pulsar orbiting a Be star and emit luminous $γ$-rays. Non-thermal emissions are thought to arise from synchrotron radiation and inverse-Compton (IC) scattering in the shock where the pulsar wind is terminated by the stellar outflow. We study wind interactions and shock radiations from such systems and show that the bimodal structures observed in keV/TeV light curves are caused by enhanced synchrotron radiation and IC scattering during disc passages. We use a simple radiation model to reproduce orbital modulations of keV X-ray and TeV $γ$-ray flux and compare with two confirmed pulsar/Be star binaries (i.e. PSR B1259-63/LS 2883 and PSR J2032+4127/MT91 213), and two candidates (i.e. HESS J0632+057 and LS I +61$^{\circ}$303). We find that the keV/TeV light curves of the former two binaries can be well explained by the inclined disc model, while modelling the modulated emissions of the latter two sources remains challenging with current orbital solutions. Therefore, we propose alternative orbital geometries for HESS J0632+057 and LS I +61$^{\circ}$303. We estimate the positions and inclination angles of Be discs by fitting correlated keV/TeV light curves. Our results could be beneficial for future measurements of orbital parameters and searches for radio pulsations from presumed pulsars.

preprint2020arXiv

A deep X-ray spectral imaging of the bow-shock pulsar wind nebula associated with PSR B1929+10

In this work, we study the X-ray bow-shock nebula powered by the mature pulsar PSR B1929+10 using data from XMM-Newton, with an effective exposure of $\sim$ 300 ks, offering the deepest investigation of this system thus far. We found the X-ray axial outflow extends as long as $\sim$ 8 arc minute behind the proper motion direction, which is a factor of two longer than the result reported in the previous study. Furthermore, we found evidence of two faint lateral outflows extending laterally with respect to the proper motion. We also found indications of spectral hardening along the axial outflow, suggesting that certain acceleration processes might occur along this feature.

preprint2020arXiv

An X-ray study for white dwarf binary AR Scorpii

We report a study of the X-ray emission from the white dwarf/M-type star binary system AR Scorpii using archival data taken in 2016-2020. It has been known that the X-ray emission is dominated by the optically thin thermal plasma emission, and its flux level varies significantly over the orbital phase. The X-ray emission also contains a component that modulates with the beat frequency between the white dwarf's spin frequency and orbital frequency. In this new analysis, the 2020 data taken by NICER shows that the X-ray emission is modulating with the spin frequency as well as the beat frequency, indicating that part of the X-ray emission is coming from the white dwarf's magnetosphere. It is found that the signal of the spin frequency appears only at a specific orbital phase, while the beat signal appears over the orbital phase. We interpret the X-ray emission modulating with the spin frequency and the beat frequency as a result of the synchrotron emission from electrons with a smaller and larger pitch angle, respectively. In a long-term evolution, the beat pulse profile averaged over the orbital phase changed from a single-peak structure in 2016/2018 to a double-peak structure in 2020. The observed X-ray flux levels measured in 2016/2017 are higher than those measured in 2018/2020. The plasma temperature and amplitude of the orbital waveform might vary with time too. These results indicate that the X-ray emission from AR Scorpii evolves on a timescale of years. This long-term evolution would be explained by a super-orbital modulation related to, for example, a precession of the white dwarf, or a fluctuation of the system related to activity of the companion star.

preprint2020arXiv

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

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.

preprint2020arXiv

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

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.

preprint2016arXiv

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

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.

preprint2016arXiv

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

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}$.

preprint2015arXiv

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

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.

preprint2015arXiv

Long-term X-ray Variability of Ultraluminous X-ray Sources

Long-term X-ray modulations on timescales from tens to hundreds of days have been widely studied for X-ray binaries located in the Milky Way and the Magellanic Clouds. For other nearby galaxies, only the most luminous X-ray sources can be monitored with dedicated observations. We here present the first systematic study of long-term X-ray variability of four ultraluminous X-ray sources (ESO 243-49 HLX-1, Holmberg IX X-1, M81 X-6, and NGC 5408 X-1) monitored with Swift. By using various dynamic techniques to analyse their light curves, we find several interesting low-frequency quasi-periodicities. Although the periodic signals may not represent any stable orbital modulations, these detections reveal that such long-term regular patterns may be related to superorbital periods and structure of the accretion discs. In particular, we show that the outburst recurrence time of ESO 243-49 HLX-1 varies over time and suggest that it may not be the orbital period. Instead, it may be due to some kinds of precession, and the true binary period is expected to be much shorter.

preprint2015arXiv

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

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.

preprint2015arXiv

Searches for millisecond pulsar candidates among the unidentified Fermi objects

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.

preprint2014arXiv

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

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.

preprint2014arXiv

Discovery of an ultracompact gamma-ray millisecond pulsar binary candidate

We report multi-wavelength observations of the unidentified Fermi object 2FGL J1653.6-0159. With the help of high-resolution X-ray observation, we have identified an X-ray and optical counterpart of 2FGL J1653.6-0159. The source exhibits a periodic modulation of 75 min in optical and possibly also in X-ray. We suggest that 2FGL J1653.6-0159 is a compact binary system with an orbital period of 75 min. Combining the gamma-ray and X-ray properties, 2FGL J1653.6-0159 is potentially a black widow/redback type gamma-ray millisecond pulsar (MSP). The optical and X-ray lightcurve profiles show that the companion is mildly heated by the high-energy emission and the X-rays are from intrabinary shock. Although no radio pulsation has been detected yet, we estimated that the spin period of the MSP is ~2ms based on a theoretical model. If pulsation can be confirmed in the future, 2FGL J1653.6-0159 will become the first ultracompact rotation-powered MSP.

preprint2014arXiv

Discovery of X-ray pulsations from "next Geminga" - PSR J1836+5925

We report the X-ray pulsation of ~173.3 ms for the "next Geminga", PSR J1836+5925, with recent XMM-Newton investigations. The X-ray periodicity is consistent wtih the gamma-ray ephemeris at the same epoch. The X-ray folded light curve has a sinusoidal structure which is different from the double-peaked gamma-ray pulse profile. We have also analysed the X-ray phase-averaged spectra which shows the X-ray emission from PSR J1836+5925 is thermal dominant. This suggests the X-ray pulsation mainly originates from the modulated hot spot on the stellar surface.

preprint2014arXiv

Emission mechanism of GeV-quiet soft gamma-ray pulsars; A case for peculiar geometry?

There is a growing new class of young spin-down powered pulsars called GeV-quiet soft gamma-ray pulsar; (1) spectral turnover appears around~10MeV, (2) the X-ray spectra of below 20 keV can be described by power law with photon index around 1.2 and (3) the light curve in X-ray/soft gamma-ray bands shows single broad pulse. Their emission properties are distinct from the normal gamma-ray pulsars, for which the spectral peak in $νF_ν$ appears in GeV energy bands and the X-ray/gamma-ray light curves show sharp and double (or more) peaks. In this paper, we discuss that X-ray/soft gamma-ray emissions of the GeV-quiet soft gamma-ray pulsars are caused bythe synchrotron radiation of the electron/positron pairs, which are created by the magnetic pair-creation process near the stellar surface. In our model, the viewing geometry is crucial factor to discriminate between the normal gamma-ray pulsars and soft gamma-ray pulsars. Our model suggests that the difference between the magnetic inclination angle ($α$) and the Earth viewing angle ($β$) of the soft gamma-ray pulsars is small, so that the synchrotron emissions from the high magnetic field region around the polar cap region dominates in the observed emissions. Furthermore, the inclination angle of the soft gamma-ray pulsar is relatively small, $α\leq 30$~degree, and our line of sight is out of the gamma-ray beam emitted via the curvature radiation process in the outer gap. We also analysis the six year $Fermi$ data for four soft gamma-ray pulsars to determine the upper limit of the GeV flux.

preprint2014arXiv

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

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.

preprint2014arXiv

High-Energy Emission of the First Millisecond Pulsar

We report on X-ray and gamma-ray observations of the millisecond pulsar B1937+21 taken with the Chandra X-ray Observatory, XMM-Newton, and the Fermi Large Area Telescope. The pulsar X-ray emission shows a purely non-thermal spectrum with a hard photon index of 0.9+/-0.1, and is nearly 100% pulsed. We found no evidence of varying pulse profile with energy as previously claimed. We also analyzed 5.5 yr of Fermi survey data and obtained much improved constraints on the pulsar's timing and spectral properties in gamma-rays. The pulsed spectrum is adequately fitted by a simple power-law with a photon index of 2.38+/-0.07. Both the gamma-ray and X-ray pulse profiles show similar two-peak structure and generally align with the radio peaks. We found that the aligned profiles and the hard spectrum in X-rays seem to be common properties among millisecond pulsars with high magnetic fields at the light cylinder. We discuss a possible physical scenario that could give rise to these features.

preprint2014arXiv

High-energy emissions from the gamma-ray binary LS 5039

We study mechanisms of multi-wavelength emissions (X-ray, GeV and TeV gamma-rays) from the gamma-ray binary LS~5039. This paper is composed of two parts. In the first part, we report on results of observational analysis using four year data of \fermi\ Large Area Telescope. Due to the improvement of instrumental response function and increase of the statistics, the observational uncertainties of the spectrum in $\sim$100-300 MeV bands and $>10$GeV bands are significantly improved. The present data analysis suggests that the 0.1-100GeV emissions from LS~5039 contain three different components; (i) the first component contributes to $<$1GeV emissions around superior conjunction, (ii) the second component dominates in 1-10GeV energy bands and (iii) the third component is compatible to lower energy tail of the TeV emissions. In the second part, we develop an emission model to explain the properties of the phase-resolved emissions in multi-wavelength observations. Assuming that LS~5039 includes a pulsar, we argue that both emissions from magnetospheric outer gap and inverse-Compton scattering process of cold-relativistic pulsar wind contribute to the observed GeV emissions. We assume that the pulsar is wrapped by two kinds of termination shock; Shock-I due to the interaction between the pulsar wind and the stellar wind and Shock-II due to the effect of the orbital motion. We propose that the X-rays are produced by the synchrotron radiation at Shock-I region and the TeV gamma-rays are produced by the inverse-Compton scattering process at Shock-II region.

preprint2014arXiv

High-energy observations of PSR B1259-63/LS 2883 through the 2014 periastron passage: connecting X-rays to the GeV flare

The binary system PSR B1259-63/LS 2883 is well sampled in radio, X-rays, and TeV gamma-rays, and shows orbital phase-dependent variability in these frequencies. The first detection of GeV gamma-rays from the system was made around the 2010 periastron passage. In this Letter, we present an analysis of X-ray and gamma-ray data obtained by the Swift/XRT, NuSTAR/FPM, and Fermi/LAT, through the recent periastron passage which occurred on 2014 May 4. While PSR B1259-63/LS 2883 was not detected by the LAT before and during this passage, we show that the GeV flares occurred at a similar orbital phase as in early 2011, thus establishing the repetitive nature of the post-periastron GeV flares. Multiple flares each lasting for a few days have been observed and short-term variability is seen as well. We also found X-ray flux variation contemporaneous with the GeV flare for the first time. A strong evidence of the keV-to-GeV connection came from the broadband high-energy spectra, which we interpret as synchrotron radiation from the shocked pulsar wind.

preprint2014arXiv

Multi-wavelength emissions from the millisecond pulsar binary PSR J1023+0038 during an accretion active state

Recent observations strongly suggest that the millisecond pulsar binary PSR J1023+0038 has developed an accretion disk since 2013 June. We present a multi-wavelength analysis of PSR J1023+0038, which reveals that 1) its gamma-rays suddenly brightened within a few days in June/July 2013 and has remained at a high gamma-ray state for several months; 2) both UV and X-ray fluxes have increased by roughly an order of magnitude, and 3) the spectral energy distribution has changed significantly after the gamma-ray sudden flux change. Time variabilities associated with UV and X-rays are on the order of 100-500 seconds and 50-100 seconds, respectively. Our model suggests that a newly formed accretion disk due to the sudden increase of the stellar wind could explain the changes of all these observed features. The increase of UV is emitted from the disk, and a new component in gamma-rays is produced by inverse Compton scattering between the new UV component and pulsar wind. The increase of X-rays results from the enhancement of injection pulsar wind energy into the intra-binary shock due to the increase of the stellar wind. We also predict that the radio pulses may be blocked by the evaporated winds from the disk and the pulsar is still powered by rotation.

preprint2014arXiv

NuSTAR observations and broadband spectral energy distribution modeling of the millisecond pulsar binary PSR J1023+0038

We report the first hard X-ray (3-79 keV) observations of the millisecond pulsar (MSP) binary PSR J1023+0038 using NuSTAR. This system has been shown transiting between a low-mass X-ray binary (LMXB) state and a rotation-powered MSP state. The NuSTAR observations were taken in both LMXB state and rotation-powered state. The source is clearly seen in both states up to ~79 keV. During the LMXB state, the 3-79 keV flux is about a factor of 10 higher that in the rotation-powered state. The hard X-rays show clear orbital modulation during the X-ray faint rotation-powered state but the X-ray orbital period is not detected in the X-ray bright LMXB state. In addition, the X-ray spectrum changes from a flat power-law spectrum during the rotation-powered state to a steeper power-law spectrum in the LMXB state. We suggest that the hard X-rays are due to the intra-binary shock from the interaction between the pulsar wind and the injected material from the low-mass companion star. During the rotation-powered MSP state, the X-ray orbital modulation is due to Doppler boosting of the shocked pulsar wind. At the LMXB state, the evaporating matter of the accretion disk due to the gamma-ray irradiation from the pulsar stops almost all the pulsar wind, resulting the disappearance of the X-ray orbital modulation.

preprint2014arXiv

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

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.

preprint2013arXiv

Discovery of X-ray Pulsation from the Geminga-like Pulsar PSR J2021+4026

We report the discovery of X-ray periodicity of ~265.3 ms from a deep XMM-Newton observation of the radio-quiet gamma-ray pulsar, PSR J2021+4026, located at the edge of the supernova remnant G78.2+2.1 (gamma-Cygni).The detected frequency is consistent with the gamma-ray pulsation determined by the observation of Fermi Gamma-ray Space Telescope at the same epoch. The X-ray pulse profile resembles the modulation of hot spot on the surface of the neutron star. The phase-averaged spectral analysis also suggests that the majority of the observed X-rays have a thermal origin. This is the third member in the class of radio-quiet pulsars with the significant pulsations detected from both X-rays and gamma-ray regimes.

preprint2013arXiv

Exploring the X-ray and gamma-ray properties of the redback millisecond pulsar PSR J1723-2837

We have investigated the X-ray and $γ$-ray properties of the redback millisecond pulsar PSR J1723-2837 with XMM-Newton, Chandra and Fermi. We have discovered the X-ray orbital modulation of this binary system with the minimum that coincides with the phases of radio eclipse. The X-ray emission is clearly non-thermal in nature which can be well described by a simple power-law with a photon index of $\sim1.2$. The phase-averaged luminosity is $\sim9\times10^{31}$ erg/s in 0.3-10 keV which consumes $\sim0.2\%$ of the spin-down power. We have detected the $γ-$ray emission in $0.1-300$ GeV from this system at a significance of $\sim6σ$ for the first time. The $γ-$rays in this energy range consumes $\sim2\%$ of the spin-down power and can be modeled by a power-law with a photon index of $\sim2.6$. We discuss the high energy properties of the new redback in the context of a intrabinary shock model.

preprint2013arXiv

Model on pulsed GeV radiation from magnetars

We discuss a possible scenario for radiation mechanism of pulsed GeV gamma-rays from magnetars. The magnetars have shown frequent X-ray bursts, which would be triggered by crust fractures and could release the energy of order of ~10^{41-42}erg. If the location of the crust cracking of the magnetic field is close to the magnetic pole, part of the released energy may excite the Alfevn wave that can propagate into outer magnetosphere. The oscillation of the magnetic field induces the available potential drop ~10^{15}Volts, which can accelerate the electrons and/or positrons to the Lorentz factor ~10^{7} in the outer magnetosphere. The curvature radiation process at outer magnetosphere can produce GeV gamma-rays. If the radiation process is occurred above r~5x 10^7cm from the stellar surface, the emitted GeV gamma-rays can escape from the pair-creation process with the X-rays and/or the magnetic field. The expected luminosity of the GeV emissions is order of L_γ< 10^{35} erg/s, and the radiation process will last for a temporal scale of years. The expected pulse profiles have a broad shape with sometimes sharp peaks. We apply the model to AXP 1E~2259+586.

preprint2013arXiv

Radiation Mechanism of the Soft Gamma-ray Pulsar PSR B1509-58

The outer gap model is used here to explain the spectrum and the energy dependent light curves of the X-ray and soft gamma-ray radiations of the spin-down powered pulsar PSR B1509-58.In the outer gap model, most pairs inside the gap are created around the null charge surface and the gap's electric field separates the two charges to move in opposite directions. Consequently, the region from the null charge surface to the light cylinder is dominated by the outflow of particles and that from the null charge surface to the star is dominated by the inflow of particles. The inflow and outflow of particles move along the magnetic field lines and emit curvature photons, and the incoming curvature photons are converted to pairs by the strong magnetic field of the star. These pairs emit synchrotron photons. We suggest that the X-rays and soft gamma-rays of PSR B1509-58 result from the synchrotron radiation of these pairs, and the viewing angle of PSR B1509-58 only receives the inflow radiation. The magnetic pair creation requires a large pitch angle, which makes the pulse profile of the synchrotron radiation distinct from that of the curvature radiation. We carefully trace the pulse profiles of the synchrotron radiation with different pitch angles. We find that the differences between the light curves of different energy bands are due to the different pitch angles of the secondary pairs, and the second peak appearing at E>10MeV comes from the region near the star, where the stronger magnetic field allows the pair creation to happen with a smaller pitch angle.

preprint2013arXiv

Search for pulsed gamma-ray emission from globular cluster M28

Using the data from the Large Area Telescope on board the Fermi Gamma-ray Space Telescope, we have searched for the gamma-ray pulsations from the direction of globular cluster M28 (NGC 6626). We report the discovery of a signal with the frequency consistent with that of the energetic millisecond pulsar (MSP) PSR B1821-24 in M28. A weighted H-test test statisic (TS) of 28.8 is attained which corresponds to a chance probability of ~1e-5 (4.3-sigma detection). With a phase-resolved analysis, the pulsed component is found to contribute ~25% of the total observed gamma-ray emission from the cluster. On the other hand, the unpulsed level provides a constraint for the underlying MSP population and the fundamental plane relations for the scenario of inverse Compton scattering. Follow-up timing observations in radio/X-ray are encouraged for further investigating this periodic signal candidate.

preprint2012arXiv

Discovery of an unidentified Fermi object as a black widow-like millisecond pulsar

The Fermi Gamma-ray Space Telescope has revolutionized our knowledge of the gamma-ray pulsar population, leading to the discovery of almost 100 gamma-ray pulsars and dozens of gamma-ray millisecond pulsars (MSPs). Although the outer-gap model predicts different sites of emission for the radio and gamma-ray pulsars, until now all of the known gamma-ray MSPs have been visible in the radio. Here we report the discovery of a "radio-quiet" gamma-ray emitting MSP candidate by using Fermi, Chandra, Swift, and optical observations. The X-ray and gamma-ray properties of the source are consistent with known gamma-ray pulsars. We also found a 4.63-hr orbital period in optical and X-ray data. We suggest that the source is a black widow-like MSP with a ~0.1 solar-mass late-type companion star. Based on the profile of the optical and X-ray light-curves, the companion star is believed to be heated by the pulsar while the X-ray emissions originate from pulsar magnetosphere and/or from intra-binary shock. No radio detection of the source has been reported yet and although no gamma-ray/radio pulsation has been found, we estimated that the spin period of the MSP is ~3-5 ms based on the inferred gamma-ray luminosity.

preprint2012arXiv

Discovery of γ-ray pulsation and X-ray emission from the black widow pulsar PSR J2051-0827

We report the discovery of pulsed γ-ray emission and X-ray emission from the black widow millisecond pulsar PSR J2051-0827 by using the data from the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope and the Advanced CCD Imaging Spectrometer array (ACIS-S) on the Chandra X-ray Observatory. Using 3 years of LAT data, PSR J2051-0827 is clearly detected in γ-ray with a signicance of \sim 8σ in the 0.2 - 20 GeV band. The 200 MeV - 20 GeV γ-ray spectrum of PSR J2051-0827 can be modeled by a simple power- law with a photon index of 2.46 \pm 0.15. Significant (\sim 5σ) γ-ray pulsations at the radio period were detected. PSR J2051-0827 was also detected in soft (0.3-7 keV) X-ray with Chandra. By comparing the observed γ-rays and X-rays with theoretical models, we suggest that the γ-ray emission is from the outer gap while the X-rays can be from intra-binary shock and pulsar magnetospheric synchrotron emissions.

preprint2012arXiv

Modeling high-energy light curves of the PSR B1259-63/LS 2883 binary based on 3-D SPH simulations

Temporal changes of X-ray to very-high-energy gamma-ray emissions from the pulsar-Be star binary PSR B1259-63/LS 2883 are studied based on 3-D SPH simulations of pulsar wind interaction with Be-disk and wind. We focus on the periastron passage of the binary and calculate the variation of the synchrotron and inverse-Compton emissions using the simulated shock geometry and pressure distribution of the pulsar wind. The characteristic double-peaked X-ray light curve from observations is reproduced by our simulation under a dense Be disk condition (base density ~10^{-9} g cm^{-3}). We interpret the pre- and post-periastron peaks as being due to a significant increase in the conversion efficiency from pulsar spin down power to the shock-accelerated particle energy at orbital phases when the pulsar crosses the disk before periastron passage, and when the pulsar wind creates a cavity in the disk gas after periastron passage, respectively. On the contrary, in the model TeV light curve, which also shows a double peak feature, the first peak appears around the periastron phase. The possible effects of cooling processes on the TeV light curve are briefly discussed.

preprint2012arXiv

Orbital Phase Dependent Gamma-ray emissions from the Black Widow Pulsar

We report on evidence for orbital phase-dependence of the gamma-ray emission from PSR B1957+20 black widow system by using the data of the Fermi Large Area Telescope. We divide an orbital cycle into two regions: a region containing the inferior conjunction, and the other region containing rest of the orbital cycle. We show that the observed spectra for the different orbital regions are fitted by different functional forms. The spectrum of the orbital region containing inferior conjunction can be described by a power-law with an exponential cutoff (PLE) model, which gives the best-fit model for the orbital phase that does not contain the inferior conjunction, plus an extra component above ~2.7 GeV. The emission above 3 GeV in this region is detected with a ~7-sigma confidence level. The gamma-ray data above ~2.7 GeV are observed to be modulated at the orbital period at the ~2.3-sigma level. We anticipate that the PLE component dominating below ~2.7 GeV originates from the pulsar magnetosphere. We also show that the inverse-Compton scattering of the thermal radiation of the companion star off a "cold" ultra-relativistic pulsar wind can explain the extra component above ~2.7 GeV. The black widow pulsar PSR B1957+20 may be the member of a new class of object, in the sense that the system is showing gamma-ray emission with both magnetospheric and pulsar wind origins.

preprint2012arXiv

Pulsed γ-ray emission from magnetar 1E 2259+586

Anomalous x-ray pulsars (AXPs) are thought to be magnetars which are young isolated neutron stars with extremely strong magnetic fields of > 10^14Gauss. Their tremendous magnetic fields inferred from the spin parameters provide a huge energy reservoir to power the observed x-ray emission. High-energy emission above 0.3 MeV has never been detected despite intensive search. Here, we present the possible Fermi Large Area Telescope (LAT) detection of γ-ray pulsations above 200 MeV from the AXP, 1E 2259+586, which puts the current theoretical models of γ-ray emission mechanisms of magnetars into challenge. We speculate that the high-energy γ-rays originate from the outer magnetosphere of the magnetar.

preprint2012arXiv

X-ray studies of the Black Widow Pulsar PSR B1957+20

We report on Chandra observations of the black widow pulsar, PSR B1957+20. Evidence for a binary-phase dependence of the X-ray emission from the pulsar is found with a deep observation. The binary-phase resolved spectral analysis reveals non-thermal X-ray emission of PSR B1957+20, confirming the results of previous studies. This suggests that the X-rays are mostly due to intra-binary shock emission which is strongest when the pulsar wind interacts with the ablated material from the companion star. The geometry of the peak emission is determined in our study. The marginal softening of the spectrum of the non-thermal X-ray tail may indicate that particles injected at the termination shock is dominated by synchrotron cooling.

preprint2011arXiv

Population study for $γ$-ray emitting Millisecond Pulsars and $Fermi$ unidentified sources

The $Fermi$-LAT has revealed that rotation powered millisecond pulsars (MSPs) are a major contributor to the Galactic $γ$-ray source population. We discuss the $γ$-ray emission process within the context of the outer gap accelerator model, and use a Monte-Calro method to simulate the Galactic population of the $γ$-ray emitting MSPs. We find that the outer gap accelerator controlled by the magnetic pair-creation process is preferable in explaining the possible correlation between the $γ$-ray luminosity and the spin down power. Our Monte-Calro simulation implies that most of the $γ$-ray emitting MSPs are radio quiet in the present sensitivity of the radio survey, indicating that most of the $γ$-ray MSPs have been unidentified. We argue that the Galactic $Fermi$ unidentified sources located at high latitudes should be dominated by MSPs, whereas the sources in the galactic plane are dominated by radio-quiet canonical pulsars.

preprint2011arXiv

Population study for gamma-ray pulsars: (III) Radiation characteristics and viewing geometry

We have performed a Monte-Calro simulation for Galactic population of pulsars and for the $γ$-ray observations. We apply outer gap model for the $γ$-ray emission process, and study the radiation characteristics as a function of the magnetic inclination angle $(α)$ and the Earth viewing angle ($ζ$). In our model, the $γ$-ray flux and the spectral cut-off energy tend to decrease as the inclination and viewing angles deviate from $90^{\circ}$. The emerging spectrum above 100 MeV becomes soft with a photon index $p\sim 1.8-2$ for $ζ\rightarrow 90^{\circ}$ and $p\sim 1.2-1.3$ for $ζ\ll 90^{\circ}$. Our simulation predicts that the pulsars with larger inclination angles ($α=70-90^{\circ}$) and larger viewing angles ($ζ=70-90^{\circ}$) have been preferentially detected by the $Fermi$ $γ$-ray telescope, and hence the observed pulse profiles of the $γ$-ray pulsars have the double peak structure rather than single peak. In the simulation, most $γ$-ray millisecond pulsars are categorized as the radio-quiet $γ$-ray pulsars, because its radio fluxes are under the sensitivities of the major radio surveys. Even we drastically increase the radio sensitivity by a factor of ten, the number of radio-selected millisecond pulsars detected by the $Fermi$ ten years observations is still much less than the expected $γ$-ray-selected millisecond pulsars, indicating the radio-quiet millisecond pulsars must contribute to the $Fermi$ unidentified sources and/or the $γ$-ray background radiations. The Unique radiation properties of the low-efficient $γ$-ray pulsar, PSR J0659+1414, can be explained by the present our gap model with a viewing geometry of $α\sim ζ=40^{\circ}-50^{\circ}$.

preprint2011arXiv

Three-dimensional Two-Layer Outer Gap Model: Fermi Energy Dependent Light Curves of the Vela Pulsar

We extend the two-dimensional two-layer outer gap model to a three-dimensional geometry and use it to study the high-energy emission of the Vela pulsar. In this model, the outer gap is divided into two parts, i.e. the main acceleration region on the top of last-open field lines and the screening region around the upper boundary of the gap. In the main acceleration region, the charge density is much lower than the Goldreich-Julian charge density and the charged particles are accelerated by the electric field along the magnetic field to emit multi-GeV photons. In the screening region, the charge density is larger than the Goldreich-Julian value to close the gap and particles in this region are responsible for multi-100MeV photon emission. We apply this three dimensional two-layer model to the Vela pulsar and compare the model light curves, the phase-averaged spectrum and the phase-resolved spectra with the recent Fermi observations, which also reveals the existence of the third peak between two main peaks. The phase position of the third peak moves with the photon energy, which cannot be explained by the geometry of magnetic field structure and the caustic effects of the photon propagation. We suggest that the existence of the third peak and its energy dependent movement results from the azimuthal structure of the outer gap.

preprint2011arXiv

X-ray and Gamma-ray Emissions from Different Evolutionary Stage of Rotation Powered Millisecond Pulsars

The $Fermi$-LAT has revealed that rotation powered millisecond pulsars (MSPs) are a major contributor to the Galactic $γ$-ray source population. Such pulsars may also be important in modeling the quiescent state of several low mass X-ray binaries (LMXBs), where optical observations of the companion star suggest the possible existence of rotation powered MSPs. To understand the observational properties of the different evolutionary stages of MSPs, the X-ray and $γ$-ray emission associated with the outer gap model is investigated. For rotation powered MSPs, the size of the outer gap and the properties of the high-energy emission are controlled by either the photon-photon pair-creation process or magnetic pair-creation process near the surface. For these pulsars, we find that the outer gap model controlled by the magnetic pair-creation process is preferable in explaining the possible correlations between the $γ$-ray luminosity or non-thermal X-ray luminosity versus the spin down power. For the accreting MSPs in quiescent LMXBs, the thermal X-ray emission at the neutron star surface resulting from deep crustal heating can control the conditions in the outer gap. We argue that the optical modulation observed in the quiescent state of several LMXBs originates from the irradiation of the donor star by $γ$-rays from the outer gap. In these systems, the irradiation luminosity required for the optical modulation of the source such as SAX J1808.4-3658 can be achieved for a neutron star of high mass. Finally, we discuss the high-energy emission associated with an intra-binary shock in black widow systems, e.g. PSR B1957+20.

preprint2011arXiv

XMM-Newton observation of PSR B2224+65 and its jet

We have investigated the pulsar PSR B2224+65 and its X-ray jet with XMM-Newton. Apart from the long X-ray jet which is almost perpendicular to the direction of proper motion, a putative extended feature at the pulsar position, which oriented in the opposite direction of the proper motion, is also suggested by this deep X-ray imaging. Non-detection of any coherent X-ray pulsation disfavors the magnetospheric origin of the X-rays observed from the position of PSR B2224+65 and hence suggest that the interpretation of pulsar wind nebula is more viable. We have also probed the origin of PSR B2224+65 and identified a runaway star, which possibly originated from the Cygnus OB9 association, as a candidate for the former binary companion of the neutron star's progenitor.

preprint2010arXiv

Alfven seismic vibrations of crustal solid-state plasma in quaking paramagnetic neutron star

Magneto-solid-mechanical model of two-component, core-crust, paramagnetic neutron star responding to quake-induced perturbation by differentially rotational, torsional, oscillations of crustal electron-nuclear solid-state plasma about axis of magnetic field frozen in the immobile paramagnetic core is developed. Particular attention is given to the node-free torsional crust-against-core vibrations under combined action of Lorentz magnetic and Hooke's elastic forces; the damping is attributed to Newtonian force of shear viscose stresses in crustal solid-state plasma. The spectral formulae for the frequency and lifetime of this toroidal mode are derived in analytic form and discussed in the context of quasi-periodic oscillations of the X-ray outburst flux from quaking magnetars. The application of obtained theoretical spectra to modal analysis of available data on frequencies of oscillating outburst emission suggests that detected variability is the manifestation of crustal Alfven's seismic vibrations restored by Lorentz force of magnetic field stresses.

preprint2010arXiv

Asteroseismic model of paramagnetic neutron star

We investigate an asteroseismic model of non-rotating paramagnetic neutron star with core-crust stratification of interior pervaded by homogeneous internal and dipolar external magnetic field, presuming that neutron degenerate Fermi-matter of the star core is in the state of Pauli's paramagnetic permanent magnetization caused by polarization of spin magnetic moments of neutrons along the axis of magnetic field of collapsed massive progenitor. The magnetic cohesion between metal-like crust and permanent-magnet-like core is considered as playing a main part in the dynamics of starquake. Focus is laid on the post-quake relaxation of the star by node-free torsional vibrations of highly conducting crustal solid-state plasma, composed of nuclei embedded in the degenerate Fermi-gas of relativistic electrons, about axis of magnetic field frozen in the immobile paramagnetic core. Two scenarios of these axisymmetric seismic vibrations are examined, in first of which these are considered as maintained by combined action of Lorentz magnetic and Hooke's elastic forces and in second one by solely Lorentz force. Based on the energy variational method of magneto-solid-mechanical theory of elastic continuous medium, the spectral formulae for the frequency and lifetime of this toroidal mode are obtained and discussed in the context of theoretical treatment of recently discovered quasi-periodic oscillations of the X-ray outburst flux from SGR 1806-20 and SGR 1900+14 as being produced by above seismic vibrations.

preprint2010arXiv

Multiwavelength properties of a new Geminga-like pulsar: PSR J2021+4026

In this paper, we report a detailed investigation of the multiwavelength properties of a newly detected gamma-ray pulsar, PSR J2021+4026, in both observational and theoretical aspects. We firstly identify an X-ray source in the XMM-Newton serendipitous source catalogue, 2XMM J202131.0+402645, located within the 95% confidence circle of PSR J2021+4026. With an archival Chandra observation, this identification provides an X-ray position with arcsecond accuracy which is helpful in facilitating further investigations. Searching for the pulsed radio emission at the position of 2XMM J202131.0+402645 with a 25-m telescope at Urumqi Astronomical Observatory resulted in null detection and places an upper-limit of 0.1~mJy for any pulsed signal at 18~cm. Together with the emission properties in X-ray and gamma-ray, the radio quietness suggests PSR J2021+4026 to be another member of Geminga-like pulsars. In the radio sky survey data, extended emission features have been identified in the gamma-ray error circle of PSR J2021+4026. We have also re-analyzed the gamma-ray data collected by FERMI's Large Area Telescope. We found that the X-ray position of 2XMM J202131.0+402645 is consistent with that of the optimal gamma-ray timing solution. We have further modeled the results in the context of outer gap model which provides us with constraints for the pulsar emission geometry such as magnetic inclination angle and the viewing angle. We have also discussed the possibility of whether PSR J2021+4026 has any physical association with the supernova remnant G78.2+2.1 (gamma-Cygni).

preprint2010arXiv

On the Transition from Accretion Powered to Rotation Powered Millisecond Pulsars

The heating associated with the deposition of $γ$-rays in an accretion disk is proposed as a mechanism to facilitate the transformation of a low mass X-ray binary to the radio millisecond pulsar phase. The $γ$-ray emission produced in the outer gap accelerator in the pulsar magnetosphere likely irradiates the surrounding disk, resulting in its heating and to the possible escape of matter from the system. We apply the model to PSR J1023+0038, which has recently been discovered as a newly born rotation powered millisecond pulsar. The predicted $γ$-ray luminosity $\sim 6 \times 10^{34}~\mathrm{erg~s^{-1}}$ can be sufficient to explain the disappearance of the truncated disk existing during the 8~month$\sim 2$~yr period prior to the 2002 observations of J1023+0038 and the energy input required for the anomalously bright optical emission of its companion star.

preprint2010arXiv

Population study for gamma-ray pulsars with the outer gap model

Inspired by increase of population of $γ$-ray emitting pulsars by the $Fermi$ telescope, we perform a population study for $γ$-ray emitting canonical pulsars. We use a Monte-Carlo technique to simulate the Galactic population of neutron stars and the radio pulsars. For each simulated neutron star, we consider the $γ$-ray emission from the outer gap accelerator in the magnetosphere. In our outer gap model, we apply the gap closure mechanism proposed by Takata et al., in which both photon-photon pair-creation and magnetic pair-creation processes are considered. Simulating the sensitivities of previous major radio surveys, our simulation predicts that there are $\sim 18-23$ radio loud and $\sim 26-34$ $γ$-ray-selected $γ$-ray pulsars, which can be detected with a $γ$-ray flux $F_γ\ge 10^{-10}~\mathrm{erg/cm^2 s}$. Applying the sensitivity of the six-month observation of the $Fermi$ telescope, 40-61 radio-selected and 36-75 $γ$-ray selected pulsars are detected within our simulation. We show that the distributions of various pulsar parameters for the simulated $γ$-ray pulsars can be consistent with the observed distribution of the $γ$-ray pulsars detected by the $Fermi$ telescope. We also predict that $\sim 64$ radio-loud and $\sim 340$ $γ$-ray-selected pulsars irradiate the Earth with a flux $F_γ\ge 10^{-11}~\mathrm{erg/cm^2 s}$, and most of those $γ$-ray pulsars are distributing with a distance more than 1~kpc and a flux $F_γ\sim 10^{-11}~\mathrm{erg/cm^2 s}$. The ration between the radio-selected and $γ$-ray-selected pulsars depend on the sensitivity of the radio surveys. We also discuss the Galactic distribution of the unidentified $Fermi$ sources and the canonical $γ$-ray pulsars.

preprint2007arXiv

A re-visit of the phase-resolved X-ray and γ-ray spectra of the Crab pulsar

We use a modified outer gap model to study the multi-frequency phase-resolved spectra of the Crab pulsar. The emissions from both poles contribute to the light curve and the phase-resolved spectra. Using the synchrotron self-Compton mechanism and by considering the incomplete conversion of curvature photons into secondary pairs, the observed phase-averaged spectrum from 100 eV - 10 GeV can be explained very well. The predicted phase-resolved spectra can match the observed data reasonably well, too. We find that the emission from the north pole mainly contributes to Leading Wing 1. The emissions in the remaining phases are mainly dominated by the south pole. The widening of the azimuthal extension of the outer gap explains Trailing Wing 2. The complicated phase-resolved spectra for the phases between the two peaks, namely Trailing Wing 1, Bridge and Leading Wing 2, strongly suggest that there are at least two well-separated emission regions with multiple emission mechanisms, i.e. synchrotron radiation, inverse Compton scattering and curvature radiation. Our best fit results indicate that there may exist some asymmetry between the south and the north poles. Our model predictions can be examined by GLAST.

preprint2006arXiv

Polarization of high-energy emissions from the Crab pulsar

We investigate polarization of high-energy emissions from the Crab pulsar in the frame work of the outer gap accelerator, following previous works of Cheng and coworkers. The recent version of the outer gap, which extends from inside the null charge surface to the light cylinder, is used for examining the synchrotron radiations from the secondary and the tertiary pairs, which are produced outside the gap. We calculate the light curve, the spectrum and the polarization characteristics, simultaneously, by taking into account gyration motion of the particles. The polarization position angle curve and the polarization degree are calculated to compare with the Crab optical data. We demonstrate that the radiations from inside the null charge surface make outer-wing and off-pulse emissions in the light curve, and the tertiary pairs contribute to bridge emissions. The emissions from the secondary pairs explain the main features of the observed light curve and spectrum. On the other hand, both emissions from inside the null charge surface and from the tertiary pairs are required to explain the optical polarization behavior of the Crab pulsar. The energy dependence of the polarization features is expected by the present model. For the Crab pulsar, the polarization position angle curve indicates that the viewing angle of the observer measured from the rotational axis is greater than $90^{\circ}$.

preprint2004arXiv

A pulsar outer gap model with trans-field structure

We investigate the electrodynamics of an outer gap in the meridional plane of the aligned-rotator. The charge depletion from the Goldreich-Julian charge density causes a large electric field along the magnetic field line. The electrons or the positrons are accelerated by the field-aligned electric field and radiate the $γ$-rays tangentially to the local magnetic field line. Some of such $γ$-rays collide with $X$-rays to materialize as the electron-positron pairs on different field lines from the field line on which they were emitted. As a result, the electric field structure is expected to change across the field lines. Including these trans-field effects, we solve the formation of the electric field self-consistently with the curvature radiation and the pair creation processes. The $γ$-ray emission and the pair creation are treated by use of Monte Carlo technique. We demonstrate that the distribution of the electric field along the field lines is affected by both the gap geometry and the external currents coming into the gap through the boundaries. In the electrodynamical model, it has been known that the solution disappears if the current density carried by the electron-positron pairs produced in the gap exceeds a critical value. We show that the critical current density is significantly increased when the trans-field structure is taken into account. We also find that the location of the inner boundary of the gap shifts toward the stellar surface from the conventional null surface as the current density increases. The reason for the shift is derived from the stability condition of the inner boundary. We also argue that the ideal-MHD condition holds outside of the gap only when the low energy particles coexist with the high energy particles migrating from the gap.