Source author record

R. X. Xu

R. X. Xu 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

77works
10topics
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

77 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

An Ultrahigh-energy $γ$-ray Bubble Powered by a Super PeVatron

We report the detection of a $γ$-ray bubble spanning at least 100$\rm deg^2$ in ultra high energy (UHE) up to a few PeV in the direction of the star-forming region Cygnus X, implying the presence Super PeVatron(s) accelerating protons to at least 10 PeV. A log-parabola form with the photon index $Γ(E) = (2.71 \pm 0.02) + (0.11 \pm 0.02) \times \log_{10} (E/10 \ {\rm TeV})$ is found fitting the gamma-ray energy spectrum of the bubble well. UHE sources, `hot spots' correlated with very massive molecular clouds, and a quasi-spherical amorphous $γ$-ray emitter with a sharp central brightening are observed in the bubble. In the core of $\sim 0.5^{\circ}$, spatially associating with a region containing massive OB association (Cygnus OB2) and a microquasar (Cygnus X-3), as well as previously reported multi-TeV sources, an enhanced concentration of UHE $γ$-rays are observed with 2 photons at energies above 1 PeV. The general feature of the bubble, the morphology and the energy spectrum, are reasonably reproduced by the assumption of a particle accelerator in the core, continuously injecting protons into the ambient medium.

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

Constraints on heavy decaying dark matter from 570 days of LHAASO observations

The Kilometer Square Array~(KM2A) of the Large High Altitude Air Shower Observatory (LHAASO) aims at surveying the northern gamma-ray sky at energies above 10 TeV with unprecedented sensitivity. Gamma-ray observations have long been one of the most powerful tools for dark matter searches, as e.g., high-energy gamma-rays could be produced by the decays of heavy dark matter particles. In this letter, we present the first dark matter analysis with LHAASO-KM2A, using the first 340~days of data from 1/2-KM2A and 230~days of data from 3/4-KM2A. Several regions of interest are used to search for a signal and account for the residual cosmic-ray background after gamma/hadron separation. We find no excess of dark matter signals, and thus place some of the strongest gamma-ray constraints on the lifetime of heavy dark matter particles with mass between 10^5 and 10^9~GeV. Our results with LHAASO are robust, and have important implications for dark matter interpretations of the diffuse astrophysical high-energy neutrino emission.

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

Discovery of a new $γ$-ray source LHAASO J0341+5258 with emission up to 200TeV

We report the discovery of a new unidentified extended $γ$-ray source in the Galactic plane named LHAASO J0341+5258 with a pre-trial significance of 8.2 standard deviations above 25 TeV. The best fit position is R.A.$=55.34^{\circ}\pm0.11^{\circ}$ and Dec$=52.97^{\circ}\pm0.07^{\circ}$. The angular size of LHAASO J0341+5258 is $0.29^\circ \pm 0.06^\circ_{stat} \pm0.02^\circ_{sys}$. The flux above 25 TeV is about $20\%$ of the flux of Crab Nebula. Although a power-law fit of the spectrum from 10 TeV to 200 TeV with the photon index $α=2.98 \pm 0.19_{stat} \pm 0.02_{sys}$ is not excluded, the LHAASO data together with the flux upper limit at 10 GeV set by the Fermi LAT observation, indicate a noticeable steepening of an initially hard power-law spectrum %($α\leq 1.75$) spectrum with a cutoff at $\approx 50$ TeV. We briefly discuss the origin of UHE gamma-rays. The lack of an energetic pulsar and a young SNR inside or in the vicinity of LHAASO J0341+5258 challenge, but do not exclude both the leptonic and hadronic scenarios of gamma-ray production.

preprint2026arXiv

Discovery of the Ultra-high energy gamma-ray source LHAASO J2108+5157

We report the discovery of a UHE gamma-ray source, LHAASO J2108+5157, by analyzing the LHAASO-KM2A data of 308.33 live days. Significant excess of gamma-ray induced showers is observed in both energy bands of 25-100 TeV and $\gt$100 TeV with 9.5 sigma and 8.5 sigma, respectively. This source is not significantly favored as an extensive source with the angular extension smaller than the point-spread function of KM2A. The measured energy spectrum from 20 to 200 TeV can be approximately described by a power-law function with an index of -2.83$\pm$ 0.18stat. A harder spectrum is demanded at lower energies considering the flux upper limit set by Fermi-LAT observations. The position of the gamma-ray emission is correlated with a giant molecular cloud, which favors a hadronic origin. No obvious counterparts have been found, deeper multiwavelength observations will help to shed new light on this intriguing UHE source.

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

Exploring Lorentz Invariance Violation from Ultra-high-energy Gamma Rays Observed by LHAASO

Recently the LHAASO Collaboration published the detection of 12 ultra-high-energy gamma-ray sources above 100 TeV, with the highest energy photon reaching 1.4 PeV. The first detection of PeV gamma rays from astrophysical sources may provide a very sensitive probe of the effect of the Lorentz invariance violation (LIV), which results in decay of high-energy gamma rays in the superluminal scenario and hence a sharp cutoff of the energy spectrum. Two highest energy sources are studied in this work. No signature of the existence of LIV is found in their energy spectra, and the lower limits on the LIV energy scale are derived. Our results show that the first-order LIV energy scale should be higher than about 10^5 times the Planck scale M_{pl} and that the second-order LIV scale is >10^{-3}M_{pl}. Both limits improve by at least one order of magnitude the previous results.

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.

preprint2024arXiv

DeepSeek LLM: Scaling Open-Source Language Models with Longtermism

The rapid development of open-source large language models (LLMs) has been truly remarkable. However, the scaling law described in previous literature presents varying conclusions, which casts a dark cloud over scaling LLMs. We delve into the study of scaling laws and present our distinctive findings that facilitate scaling of large scale models in two commonly used open-source configurations, 7B and 67B. Guided by the scaling laws, we introduce DeepSeek LLM, a project dedicated to advancing open-source language models with a long-term perspective. To support the pre-training phase, we have developed a dataset that currently consists of 2 trillion tokens and is continuously expanding. We further conduct supervised fine-tuning (SFT) and Direct Preference Optimization (DPO) on DeepSeek LLM Base models, resulting in the creation of DeepSeek Chat models. Our evaluation results demonstrate that DeepSeek LLM 67B surpasses LLaMA-2 70B on various benchmarks, particularly in the domains of code, mathematics, and reasoning. Furthermore, open-ended evaluations reveal that DeepSeek LLM 67B Chat exhibits superior performance compared to GPT-3.5.

preprint2020arXiv

An in-depth investigation of 11 pulsars discovered by FAST

We present timing solutions and analyses of 11 pulsars discovered by the Five-hundred-meter Aperture Spherical radio Telescope (FAST). These pulsars were discovered using an ultra-wide bandwidth receiver in drift-scan observations made during the commissioning phase of FAST, and were then confirmed and timed using the 64-m Parkes Radio Telescope. Each pulsar has been observed over a span of at least one year. Highlighted discoveries include PSR J0344-0901, which displays mode-changing behaviour and may belong to the class of so-called `swooshing' pulsars (alongside PSRs B0919+06 and B1859+07); PSR J0803-0942, whose emission is almost completely linearly polarised; and PSRs J1900-0134 and J1945+1211, whose well defined polarisation angle curves place stringent constraints on their emission geometry. We further discuss the detectability of these pulsars by earlier surveys, and highlight lessons learned from our work in carrying out confirmation and monitoring observations of pulsars discovered by a highly sensitive telescope, many of which may be applicable to next-generation pulsar surveys. This paper marks one of the first major releases of FAST-discovered pulsars, and paves the way for future discoveries anticipated from the Commensal Radio Astronomy FAST Survey (CRAFTS).

preprint2020arXiv

The observation of the Crab Nebula with LHAASO-KM2A for the performance study

As a sub-array of the Large High Altitude Air Shower Observatory (LHAASO), KM2A is mainly designed to cover a large fraction of the northern sky to hunt for gamma-ray sources at energies above 10 TeV. Even though the detector construction is still underway, a half of the KM2A array has been operating stably since the end of 2019. In this paper, we present the pipeline of KM2A data analysis and the first observation on the Crab Nebula, a standard candle in very high energy gamma-ray astronomy. We detect gamma-ray signals from the Crab Nebula in both energy ranges of 10$-$100 TeV and $>$100 TeV with high significance, by analyzing the KM2A data of 136 live days between December 2019 and May 2020. With the observations, we test the detector performance including angular resolution, pointing accuracy and cosmic ray background rejection power. The energy spectrum of the Crab Nebula in the energy range 10-250 TeV fits well with a single power-law function dN/dE =(1.13$\pm$0.05$_{stat}$$\pm$0.08$_{sys}$)$\times$10$^{-14}$$\cdot$(E/20TeV)$^{-3.09\pm0.06_{stat}\pm0.02_{sys}}$ cm$^{-2}$ s$^{-1}$ TeV$^{-1}$. It is consistent with previous measurements by other experiments. This opens a new window of gamma-ray astronomy above 0.1 PeV through which ultrahigh-energy gamma-ray new phenomena, such as cosmic PeVatrons, might be discovered.

preprint2016arXiv

eXTP -- enhanced X-ray Timing and Polarimetry Mission

eXTP is a science mission designed to study the state of matter under extreme conditions of density, gravity and magnetism. Primary targets include isolated and binary neutron stars, strong magnetic field systems like magnetars, and stellar-mass and supermassive black holes. The mission carries a unique and unprecedented suite of state-of-the-art scientific instruments enabling for the first time ever the simultaneous spectral-timing-polarimetry studies of cosmic sources in the energy range from 0.5-30 keV (and beyond). Key elements of the payload are: the Spectroscopic Focusing Array (SFA) - a set of 11 X-ray optics for a total effective area of about 0.9 m^2 and 0.6 m^2 at 2 keV and 6 keV respectively, equipped with Silicon Drift Detectors offering <180 eV spectral resolution; the Large Area Detector (LAD) - a deployable set of 640 Silicon Drift Detectors, for a total effective area of about 3.4 m^2, between 6 and 10 keV, and spectral resolution <250 eV; the Polarimetry Focusing Array (PFA) - a set of 2 X-ray telescope, for a total effective area of 250 cm^2 at 2 keV, equipped with imaging gas pixel photoelectric polarimeters; the Wide Field Monitor (WFM) - a set of 3 coded mask wide field units, equipped with position-sensitive Silicon Drift Detectors, each covering a 90 degrees x 90 degrees FoV. The eXTP international consortium includes mostly major institutions of the Chinese Academy of Sciences and Universities in China, as well as major institutions in several European countries and the United States. The predecessor of eXTP, the XTP mission concept, has been selected and funded as one of the so-called background missions in the Strategic Priority Space Science Program of the Chinese Academy of Sciences since 2011. The strong European participation has significantly enhanced the scientific capabilities of eXTP. The planned launch date of the mission is earlier than 2025.

preprint2016arXiv

Rotational evolution of magnetars in the presence of a fallback disk

Magnetars may have strong surface dipole field. Observationally, two magnetars may have passive fallback disks. In the presence of a fallback disk, the rotational evolution of magnetars may be changed. In the self-similar fallback disk model, it is found that: (1) When the disk mass is significantly smaller than $10^{-6} \,\rm M_{\odot}$, the magnetar is unaffected by the fallback disk and it will be a normal magnetar. (2) When the disk mass is large, but the magnetar's surface dipole field is about or below $10^{14} \,\rm G$, the magnetar will also be a normal magnetar. A magnetar plus a passive fallback disk system is expected. This may correspond to the observations of magnetars 4U 0142$+$61, and 1E 2259$+$586. (3) When the disk mass is large, and the magnetar's surface dipole field is as high as $4\times 10^{15} \,\rm G$, the magnetar will evolve from the ejector phase to the propeller phase, and then enter rotational equilibrium. The magnetar will be slowed down quickly in the propeller phase. The final rotational period can be as high $2\times 10^4 \,\rm s$. This may correspond to the super-slow magnetar in the supernova remnant RCW 103. Therefore, the three kinds of magnetars can be understood in a unified way.

preprint2016arXiv

Spontaneous Magnetization of Solid Quark-cluster Stars

Pulsar-like compact stars usually have strong magnetic fields, with the strength from $\sim 10^8$ to $\sim 10^{12}$ Gauss on surface. How such strong magnetic fields can be generated and maintained is still an unsolved problem, which is, in principle, related to the interior structure of compact stars, i.e., the equation of state of cold matter at supra-nuclear density. In this paper we are trying to solve the problem in the regime of solid quark-cluster stars. Inside quark-cluster stars, the extremely low ratio of number density of electrons to that of baryons $n_e/n_b$ and the screening effect from quark-clusters could reduce the long-range Coulomb interaction between electrons to short-range interaction. In this case, the Stoner's model could apply, and we find that the condition for ferromagnetism is consistent with that for validity of Stoner's model. Under the screened Coulomb repulsion, the electrons inside the stars could spontaneously magnetized and become ferromagnetic, and hence would contribute non-zero net magnetic momentum to the whole star. We conclude that, for most cases in solid quark-cluster stars, the amount of net magnetic momentum, which is proportional to the amount of unbalanced spins $ξ=(n_+-n_-)/n_e$ and depends on the number density of electrons $n_e=n_++n_-$, could be significant with non-zero $ξ$. The net magnetic moments of electron system in solid quark-cluster stars could be large enough to induce the observed magnetic fields for pulsars with $B\sim 10^{11}$ to $\sim 10^{13}$ Gauss.

preprint2015arXiv

A Study of Multi-frequency Polarization Pulse Profiles of Millisecond Pulsars

We present high signal-to-noise ratio, multi-frequency polarization pulse profiles for 24 millisecond pulsars that are being observed as part of the Parkes Pulsar Timing Array (PPTA) project. The pulsars are observed in three bands, centred close to 730, 1400 and 3100 MHz, using a dual-band 10 cm/50 cm receiver and the central beam of the 20 cm multibeam receiver. Observations spanning approximately six years have been carefully calibrated and summed to produce high S/N profiles. This allows us to study the individual profile components and in particular how they evolve with frequency. We also identify previously undetected profile features. For many pulsars we show that pulsed emission extends across almost the entire pulse profile. The pulse component widths and component separations follow a complex evolution with frequency; in some cases these parameters increase and in other cases they decrease with increasing frequency. The evolution with frequency of the polarization properties of the profile is also non-trivial. We provide evidence that the pre- and post-cursors generally have higher fractional linear polarization than the main pulse. We have obtained the spectral index and rotation measure for each pulsar by fitting across all three observing bands. For the majority of pulsars, the spectra follow a single power-law and the position angles follow a $λ^2$ relation, as expected. However, clear deviations are seen for some pulsars. We also present phase-resolved measurements of the spectral index, fractional linear polarization and rotation measure. All these properties are shown to vary systematically over the pulse profile.

preprint2015arXiv

Pulsar wind model for the spin-down behavior of intermittent pulsars

Intermittent pulsars are part-time radio pulsars. They have higher slow down rate in the on state (radio-loud) than in the off state (radio-quiet). This gives the evidence that particle wind may play an important role in pulsar spindown. The effect of particle acceleration is included in modeling the rotational energy loss rate of the neutron star. Applying the pulsar wind model to the three intermittent pulsars (PSR B1931+24, PSR J1841-0500, and PSR J1832+0029), their magnetic field and inclination angle are calculated simultaneously. The theoretical braking indices of intermittent pulsars are also given. In the pulsar wind model, the density of the particle wind can always be the Goldreich-Julian density. This may ensure that different on states of intermittent pulsars are stable. The duty cycle of particle wind can be determined from timing observations. It is consistent with the duty cycle of the on state. Inclination angle and braking index observations of intermittent pulsars may help to test different models of particle acceleration. At present, the inverse Compton scattering induced space charge limited flow with field saturation model can be ruled out.

preprint2015arXiv

Small glitches: the role of strange nuggets?

Pulsar glitches, i.e. the sudden spin-ups of pulsars, have been detected for most pulsars that we known. The mechanism giving rise to this kind of phenomenon is uncertain, although a large data set has been built. In the framework of star-quake model, based on~\cite{Baym1971}, the glitch-sizes (the relative increases of spin-frequencies during glitches) $ΔΩ/Ω$ depend on the released energies during glitches, with less released energies corresponding to smaller glitch sizes. On the other hand, as one of dark matter candidates, our Galaxy might be filled with the so called strange nuggets (SNs) which are the relics from the early Universe. In this case the collisions between pulsars and SNs is inevitable, and these collisions would lead to glitches when enough elastic energy has been accumulated during the spin-down process. The SNs-triggered glitches could release less energy, because the accumulated elastic energy would be less than that in the scenario of glitches without SNs. Therefore, if a pulsar is hit frequently by SNs, it would tend to have more small size glitches, whose values of $ΔΩ/Ω$ are smaller than that in standard star-quake model (with larger amounts of released energy). Based on the assumption that in our Galaxy the distribution of SNs is similar to that of dark matter, as well as on the glitch data in ATNF Pulsar Catalogue and Jodrell Bank glitch table, we find that in our Galaxy the incidences of small size glitches exhibit tendencies consistent with the collision rates of pulsars and strange nuggets. Further test of this scenario is expected by detecting more small glitches (e.g., by the Square Kilometre Array).

preprint2015arXiv

Structures of the Vela pulsar and the glitch crisis from the Brueckner theory

Detailed structures of the Vela pulsar (PSR B0833-45, with a period of $89.33$ milliseconds) are predicted by adopting a recently-constructed unified treatment of all parts of neutron stars: the outer crust, the inner crust and the core based on modern microscopic Brueckner-Hartree-Fock calculations. To take the pulsar mass ranging from $1.0M_{\odot}$ to $2.0M_{\odot}$, we calculate the central density, the core/crust radii, the core/crust mass, the core/crustal thickness, the moment of inertia, and the crustal moment of inertia. Among them, the crustal moment of inertia could be effectively constrained from the accumulated glitch observations, which has been a great debate recently, known as "glitch crisis". Namely, superfluid neutrons contained in the inner crust, which are regarded as the origin of the glitch in the standard two-component model, could be largely entrained in the nuclei lattices, then there may not be enough superfluid neutrons ($\sim 4/5$ less than the previous value) to trigger the large glitches ($Δν/ν_0 \sim 10^{-6}$) in the Vela pulsar. We then provide the first analysis of the crisis based on the microscopic basis, by confronting the glitch observations with the theoretical calculations for the crustal moment of inertia. We find that despite some recent opposition to the crisis argument, the glitch crisis is still present, which means that besides the crust superfluid neutrons, core neutrons might be necessary for explaining the large glitches of the Vela pulsar.

preprint2015arXiv

The extremely long period X-ray source in a young supernova remnant: a Thorne-Zytkow Object descendant?

The origin of the 6.67 hr period X-ray source, 1E161348-5055, in the young supernova remnant RCW 103 is puzzling. We propose that it may be the descendant of a Thorne-Zytkow Object (TZO). A TZO may at its formation have a rapidly spinning neutron star as a core, and a slowly rotating envelope. We found that the core could be braked quickly to an extremely long spin period by the coupling between its magnetic field and the envelope, and that the envelope could be disrupted by some powerful bursts or exhausted via stellar wind. If the envelope is disrupted after the core has spun down, the core will become an extremely long-period compact object, with a slow proper motion speed, surrounded by a supernova-remnant-like shell. These features all agree with the observations of 1E161348-5055. TZOs are expected to have produced extraordinary high abundances of lithium and rapid proton process elements that would remain in the remnants and could be used to test this scenario.

preprint2014arXiv

Wind braking of magnetars: to understand magnetar's multiwave radiation properties

Magnetars are proposed to be peculiar neutron stars powered by their super strong magnetic field. Observationally, anomalous X-ray pulsars and soft gamma-ray repeaters are believed to be magnetar candidates. While more and more multiwave observations of magnetars are available, unfortunately, we see accumulating failed predictions of the traditional magnetar model. These challenges urge rethinking of magnetar. Wind braking of magnetars is one of the alternative modelings. The release of magnetic energy may generate a particle outflow (i.e., particle wind), that results in both an anomalous X-ray luminosity and significantly high spindown rate. In this wind braking scenario, only strong multipole field is necessary for a magnetar (a strong dipole field is no longer needed). Wind braking of magnetars may help us to understand their multiwave radiation properties, including (1) Non-detection of magnetars in Fermi-LAT observations, (2) The timing behaviors of low magnetic field magnetars, (3) The nature of anti-glitches, (4) The criterion for magnetar's radio emission, etc. In the wind braking model of magentars, timing events of magnetars should always be accompanied by radiative events. It is worth noting that the wind engine should be the central point in the research since other efforts with any reasonable energy mechanism may also reproduce the results.

preprint2013arXiv

A corresponding-state approach to quark-cluster matter

The state of super-dense matter is essential for us to understand the nature of pulsars, but the non- perturbative quantum chromodynamics (QCD) makes it very difficult for direct calculations of the state of cold matter at realistic baryon number densities inside compact stars. Nevertheless, from an observational point of view, it is conjectured that pulsars could be made up of quark clusters since the strong coupling between quarks might render quarks grouped in clusters. We are trying an effort to find an equation of state of condensed quark-cluster matter in a phenomenological way. Supposing that the quark-clusters could be analogized to inert gases, we apply here the corresponding-state approach to derive the equation of state of quark-cluster matter, as was similarly demonstrated for nuclear and neutron-star matter in 1970s. According to the calculations presented, the quark-cluster stars, which are composed of quark-cluster matter, could then have high maximum mass that is consistent with observations and, in turn, further observations of pulsar mass would also put constraints to the properties of quark-cluster matter. Moreover, the melting heat during solid-liquid phase conversion and the related astrophysical consequences are also briefly discussed.

preprint2013arXiv

H-cluster stars

The study of dense matter at ultra-high density has a very long history, which is meaningful for us to understand not only cosmic events in extreme circumstances but also fundamental laws of physics. It is well known that the state of cold matter at supra-nuclear density depends on the non-perturbative nature of quantum chromo-dynamics (QCD) and is essential for modeling pulsars. A so-called H-cluster matter is proposed in this paper as the nature of dense matter in reality. In compact stars at only a few nuclear densities but low temperature, quarks could be interacting strongly with each other there. That might render quarks grouped in clusters, although the hypothetical quark-clusters in cold dense matter has not been confirmed due to the lack of both theoretical and experimental evidence. Motivated by recent lattice QCD simulations of the H-dibaryons (with structure uuddss), we are therefore considering here a possible kind of quark-clusters, H-clusters, that could emerge inside compact stars during their initial cooling, as the dominant components inside (the degree of freedom could then be H-clusters there).Taking into account the in-medium stiffening effect, we find that at baryon densities of compact stars $H$-cluster matter could be more stable than nuclear matter. We also find that for the H-cluster matter with lattice structure, the equation of state could be so stiff that it would seem to be "superluminal" in most dense region. However, the real sound speed for H-cluster matter is in fact hard to calculate, so at this stage we do not put constraints on our model from the usual requirment of causality.

preprint2013arXiv

The timing behavior of magnetar Swift J1822.3-1606: timing noise or a decreasing period derivative?

The different timing results of the magnetar Swift J1822.3-1606 is analyzed and understood theoretically. It is pointed that different timing solutions are caused not only by timing noise, but also that the period derivative is decreasing after outburst. Both the decreasing period derivative and the large timing noise may be originated from wind braking of the magnetar. Future timing of Swift J1822.3-1606 will help us make clear whether its period derivative is decreasing with time or not.

preprint2013arXiv

Wind braking of magnetars

(adapted)Considering recent observations challenging the traditional magnetar model, we explore the wind braking of magnetars. There is evidence for strong multipole magnetic fields in active magnetars, but the dipole field inferred from spin down measurements may be strongly biased by a particle wind. Recent challenging observations of magnetars may be explained naturally in the wind braking scenario: (1) The supernova energies of magnetars are of normal value; (2) The non-detection in Fermi observations of magnetars; (3) The problem posed by the low-magnetic field soft gamma-ray repeaters; (4) The relation between magnetars and high magnetic field pulsars; (5) A decreasing period derivative during magnetar outbursts. Transient magnetars may still be magnetic dipole braking. This may explain why low luminosity magnetars are more likely to have radio emissions. In the wind braking scenario, magnetars are neutron stars with strong multipole field. For some sources, a strong dipole field may be no longer needed. A magnetism-powered pulsar wind nebula and a braking index smaller than three are the two predictions of the wind braking model.

preprint2013arXiv

X-ray and Optical Plateau Following the Main Bursts in Gamma-Ray Bursts and SNe II-P: A hint to the similar late injection behavior?

We analyze the emission plateaus in the X-ray afterglow lightcurves of gamma-ray bursts (GRBs) and in the optical lightcurves of Type II superpernovae (SNe IIP) in order to study whether they have similar late energy injection behaviors. We show that correlations of bolometric energies (or luminosities) between the prompt explosions and the plateaus for the two phenomena are similar. The Type II SNe are in the low energy end of the GRBs. The bolometric energies (or luminosities) in prompt phase E_{\rm expl} (or L_{\rm expl}) and in plateau phase E_{\rm plateau} (or L_{\rm plateau}) share relations of E_{\rm expl} \propto E_{\rm plateau}^{0.73\pm 0.14} and L_{\rm expl} \propto L_{\rm plateau}^{\sim 0.70}. These results may indicate a similar late energy injection behavior to reproduce the observed plateaus of the plateaus in the two phenomena.

preprint2012arXiv

Absorption features caused by oscillations of electrons on the surface of a quark star

If quark stars exist, they may be enveloped in thin electron layers (electron seas), which uniformly surround the entire star. These layers will be affected by the magnetic fields of quark stars in such a way that the electron seas would transmit hydromagnetic cyclotron waves, as studied in this paper. Particular attention is devoted to vortex hydrodynamical oscillations of the electron sea. The frequency spectrum of these oscillations is derived in analytic form. If the thermal X-ray spectra of quark stars are modulated by vortex hydrodynamical vibrations, the thermal spectra of compact stars, foremost central compact objects (CCOs) and X-ray dim isolated neutron stars (XDINSs), could be used to verify the existence of these vibrational modes observationally. The central compact object 1E 1207.4-5209 appears particularly interesting in this context, since its absorption features at 0.7 keV and 1.4 keV can be comfortably explained in the framework of the hydro-cyclotron oscillation model.

preprint2012arXiv

Detection of 107 glitches in 36 southern pulsars

Timing observations from the Parkes 64-m radio telescope for 165 pulsars between 1990 and 2011 have been searched for period glitches. A total of 107 glitches were identified in 36 pulsars, where 61 have previously been reported and 46 are new discoveries. Glitch parameters were measured by fitting the timing residual data. Observed relative glitch sizes Δν_g/νrange between 10^-10 and 10^-5, where ν= 1/P is the pulse frequency. We confirm that the distribution of Δν_g/νis bimodal with peaks at approximately 10^-9 and 10^-6. Glitches are mostly observed in pulsars with characteristic ages between 10^3 and 10^5 years, with large glitches mostly occurring in the younger pulsars. Exponential post-glitch recoveries were observed for 27 large glitches in 18 pulsars. The fraction Q of the glitch that recovers exponentially also has a bimodal distribution. Large glitches generally have low Q, typically a few per cent, but large Q values are observed in both large and small glitches. Observed time constants for exponential recoveries ranged between 10 and 300 days with some tendency for longer timescales in older pulsars. Shorter timescale recoveries may exist but were not revealed by our data which typically have observation intervals of 2 - 4 weeks. For most of the 36 pulsars with observed glitches, there is a persistent linear increase in \dotνin the inter-glitch interval. Where an exponential recovery is also observed, the effects of this are superimposed on the linear increase in \dotν. In some cases, the slope of the linear recovery changes at the time of a glitch. The \ddotνvalues characterising the linear changes in \dotνare almost always positive and, after subtracting the magnetospheric component of the braking, are approximately proportional to the ratio of |\dotν| and the inter-glitch interval, as predicted by vortex-creep models.

preprint2012arXiv

Is magnetar a fact or fiction to us?

The key point of studying AXPs/SGRs (anomalous X-ray pulsars/soft gamma-ray repeaters) is relevant to the energy budget. Historically, rotation was thought to be the only free energy of pulsar until the discovery of accretion power in X-ray binaries. AXPs/SGRs could be magnetars if they are magnetism-powered, but would alternatively be quark-star/fallback-disk systems if more and more observations would hardly be understood in the magnetar scenario.

preprint2012arXiv

SGR 0418+5729: a small inclination angle resulting in a not so low dipole magnetic field?

The spin down behaviors of SGR 0418+5729 are investigated. The pulsar spin down model of Contopoulos & Spitkovsky (2006) is applied to SGR 0418+5729. It is shown that SGR 0418+5729 lies below the pulsar death line and its rotation-powered magnetospheric activities may therefore have stopped. The compact star is now spun down by the magnetic dipole moment perpendicular to its rotation axis. Our calculations show that under these assumption there is the possibility of SGR 0418+5729 having a strong dipole magnetic field, if there is a small magnetic inclination angle. Its dipole magnetic field may be much higher than the characteristic magnetic field. Therefore, SGR 0418+5729 may be a normal magnetar instead of a low magnetic field magnetar.

preprint2012arXiv

The missing compact star of SN1987A: a solid quark star?

To investigate the missing compact star of Supernova 1987A, we analyzed both the cooling and the heating processes of a possible compact star based on the upper limit of observational X-ray luminosity. From the cooling process we found that a solid quark-cluster star, which has a stiffer equation of state than that of conventional liquid quark star, has a heat capacity much smaller than a neutron star. It can cool down quickly, which can naturally explain the non-detection of a point source (neutron star or quark star) in X-ray band. On the other hand, we consider the heating process from magnetospheric activity and possible accretion, and obtain some constraints to the parameters of a possible pulsar. We conclude that a solid quark-cluster star can be fine with the observational limit in a large and acceptable parameter space. A pulsar with a short period and a strong magnetic field (or with a long period and a weak field) would has luminosity higher than the luminosity limit if the optical depth is not large enough to hide the compact star. The constraints of the pulsar parameters can be tested if the central compact object in 1987A is discovered by advanced facilities in the future.

preprint2012arXiv

What can Fermi tell us about magnetars?

We have analyzed the physical implications of Fermi observations of magnetars. Observationally, no significant detection is reported in Fermi observations of all magnetars. Then there are conflicts between outer gap model in the case of magnetars and Fermi observations. One possible explanation is that magnetars are wind braking instead of magnetic dipole braking. In the wind braking scenario, magnetars are neutron stars with strong multipole field. A strong dipole field is no longer required. A magnetism-powered pulsar wind nebula and a braking index smaller than three are the two predictions of wind braking of magnetars. Future deeper Fermi observations will help us make clear whether they are wind braking or magnetic dipole braking. It will also help us to distinguish between the magnetar model and the accretion model for AXPs and SGRs.

preprint2011arXiv

AXPs and SGRs in the outer gap model: confronting Fermi observations

Anomalous X-ray pulsars (AXPs) and soft gamma-ray repeaters (SGRs) are magnetar candidates, i.e., neutron stars powered by strong magnetic field. If they are indeed magnetars, they will emit high-energy gamma-rays which are detectable by Fermi-LAT according to the outer gap model. However, no significant detection is reported in recent Fermi-LAT observations of all known AXPs and SGRs. Considering the discrepancy between theory and observations, we calculate the theoretical spectra for all AXPs and SGRs with sufficient observational parameters. Our results show that most AXPs and SGRs are high-energy gamma-ray emitters if they are really magnetars. The four AXPs 1E 1547.0-5408, XTE J1810-197, 1E 1048.1-5937, and 4U 0142+61 should have been detected by Fermi-LAT. Then there is conflict between out gap model in the case of magnetars and Fermi observations. Possible explanations in the magnetar model are discussed. On the other hand, if AXPs and SGRs are fallback disk systems, i.e., accretion-powered for the persistent emissions, most of them are not high-energy gamma-ray emitters. Future deep Fermi-LAT observations of AXPs and SGRs will help us make clear whether they are magnetars or fallback disk systems.

preprint2011arXiv

Constraint on parameters of the Inverse Compton Scattering model for radio pulsars

The inverse Compton scattering (ICS) model can explain various pulse profile shapes and diversity of pulse profile evolution based on the mechanism that the radio emission is generated through inverse Compton scattering between secondary relativistic particles and radio waves from polar gap avalanches. In this paper, we study the parameter space of ICS model for 15 pulsars, which share the common pulse profile evolution phenomena that the pulse profiles are narrower at higher observing frequencies. Two key parameters, the initial Lorentz factor and the energy loss factor of secondary particles are constrained using the least square fitting method, where we fit the theoretical curve of the "beam-frequency mapping" of the ICS model to the observed pulse widths at multiple frequencies. The uncertainty of the inclination and viewing angles are taken into account in the fitting process. It is found that the initial Lorentz factor is larger than 4000, and the energy loss factor is between 20 and 560. The Lorentz factor is consistent with the prediction of the inner vacuum gap model. Such high energy loss factors suggest significant energy loss for secondary particles at altitudes of a few tens to hundreds of kilometers.

preprint2011arXiv

Magnetars: fact or fiction?

Anomalous X-ray pulsars (AXPs) and soft gamma-ray repeaters (SGRs) are enigmatic pulsar-like objects. The energy budget is the fundamental problem in their studies. In the magnetar model, they are supposed to be powered by the extremely strong magnetic fields (>~ 10^14 G) of neutron stars. Observations for and against the magnetar model are both summarized. Considering the difficulties encountered by the magnetar model to comfortably understand more and more observations, one may doubt that AXPs and SGRs are really magnetars. If they are not magnetar candidates (including magnetar-based models), then they must be "quark star/fallback disk" systems.

preprint2011arXiv

On Hoyle-Narlikar-Wheeler mechanism of vibration energy powered magneto-dipole emission of neutron stars

We revisit the well-known Hoyle-Narlikar-Wheeler proposition that neutron star emerging in the magnetic-flux-conserving process of core-collapse supernova can convert the stored energy of Alfven vibrations into power of magneto-dipole radiation. We show that the necessary requirement for the energy conversion is the decay of internal magnetic field. In this case the loss of vibration energy of the star causes its vibration period, equal to period of pulsating emission, to lengthen at a rate proportional to the rate of magnetic field decay. These prediction of the model of vibration powered neutron star are discussed in juxtaposition with data on pulsating emission of magnetars whose radiative activity is generally associated with the decay of ultra strong magnetic field.

preprint2011arXiv

Pulsars and Gravitational Waves

The relationship between pulsar-like compact stars and gravitational waves is briefly reviewed. Due to regular spins, pulsars could be useful tools for us to detect ~nano-Hz low-frequency gravitational waves by pulsar-timing array technique; besides, they would also be ~kilo-Hz high-frequency gravitational wave radiators because of their compactness. The wave strain of an isolate pulsar depends on the equation state of cold matter at supra-nuclear densities. Therefore, a real detection of gravitational wave should be very meaningful in gravity physics, micro-theory of elementary strong interaction, and astronomy.

preprint2011arXiv

Pulsating magneto-dipole radiation of a quaking neutron star powered by energy of Alfven seismic vibrations

We compute the characteristic parameters of magnetic dipole radiation of a neutron star undergoing torsional seismic vibrations under the action of Lorentz restoring force about axis of a dipolar magnetic field experiencing decay. After brief outline of general theoretical background of the model of vibration powered neutron star, we present numerical estimates of basic vibration and radiation characteristics, such as the oscillation frequency, lifetime, luminosity of radiation, and investigate their time dependence upon magnetic field decay. The presented analysis suggests that gradual decrease in frequencies of pulsating high-energy emission detected from a handful of currently monitored AXP/SGR-like X-ray sources can be explained as being produced by vibration powered magneto-dipole radiation of quaking magnetars. \keywords{neutron stars, torsion Alfvén vibrations, vibration powered radiation, magnetic field decay, magnetars}

preprint2011arXiv

Radiative activity of magnetic white dwarf undergoing Lorentz-force-driven torsional vibrations

We study radiative activity of magnetic white dwarf undergoing torsional vibrations about axis of its own dipole magnetic moment under the action of Lorentz restoring force. It is shown that pulsating white dwarf can convert its vibration energy into the energy of magneto-dipole emission, oscillating with the frequency equal to the frequency of Alfvén torsional vibrations, provided that internal magnetic field is decayed. The most conspicuous feature of the vibration energy powered radiation in question is the lengthening of periods of oscillating emission; the rate of period elongation is determined by the rate magnetic field decay.

preprint2011arXiv

Strong pulses detected from a rotating radio transient J1819$-$1458

We analyze individual pulses detected from RRAT J1819$-$1458. From April 2007 to April 2010, we carried out observations using the Nanshan 25-m radio telescope of Urumqi Observatory at a central frequency of 1541.25 MHz. We obtain a dispersion measure $DM=195.7\pm0.3$ pc cm^{-3} by analyzing all the 423 detected bursts. The tri-band pattern of arrival time residuals is confirmed by a single pulse timing analysis. Twenty-seven bimodal bursts located in the middle residual band are detected, and, profiles of two typical bimodal bursts and two individual single-peak pulses are presented. We determine the statistical properties of SNR and W$_{50}$ of bursts in different residual bands. The W$_{50}$ variation with SNR shows that the shapes of bursts are quite different from each other. The cumulative probability distribution of intensity for a possible power law with index $α=1.6\pm0.2$ is inferred from the number of those bursts with $SNR\ge6$ and high intensities.

preprint2011arXiv

The optical/UV excess of isolated neutron stars in the RCS model

The X-ray dim isolated neutron stars (XDINSs) are peculiar pulsar-like objects, characterized by their very well Planck-like spectrum. In studying their spectral energy distributions, the optical/UV excess is a long standing problem. Recently, Kaplan et al. (2011) have measured the optical/UV excess for all seven sources, which is understandable in the resonant cyclotron scattering (RCS) model previously addressed. The RCS model calculations show that the RCS process can account for the observed optical/UV excess for most sources . The flat spectrum of RX J2143.0+0654 may due to contribution from bremsstrahlung emission of the electron system in addition to the RCS process.

preprint2011arXiv

Vibration Powered Radiation of Quaking Magnetar

In juxtaposition with the standard model of rotation powered pulsar, the model of vibration powered magnetar undergoing quake-induced torsional Alfven vibrations in its own ultra strong magnetic field experiencing decay is considered. The presented line of argument suggests that gradual decrease of frequencies (lengthening of periods) of long-periodic pulsed radiation detected from a set of X-ray sources can be attributed to magnetic-field-decay induced energy conversion from seismic vibrations to magneto-dipole radiation of quaking magnetar.

preprint2010arXiv

A new parametric equation of state and quark stars

It is still a matter of debate to understand the equation of state of cold supra-nuclear matter in compact stars because of unknown on-perturbative strong interaction between quarks. Nevertheless, it is speculated from an astrophysical view point that quark clusters could form in cold quark matter due to strong coupling at realistic baryon densities. Although it is hard to calculate this conjectured matter from first principles, one can expect the inter-cluster interaction to share some general features to nucleon-nucleon interaction. We adopt a two-Gaussian component soft-core potential with these general features and show that quark clusters can form stable simple cubic crystal structure if we assume Gaussian form wave function. With this parameterizing, Tolman-Oppenheimer-Volkoff equation is solved with reasonable constrained parameter space to give mass-radius relation of crystalline solid quark star. With baryon densities truncated at 2 times nuclear density at surface and range of interaction fixed at 2fm we can reproduce similar mass-radius relation to that obtained with bag model equations of state. The maximum mass ranges from about 0.5 to 3 solar mass. Observed maximum pulsar mass (about 2 solar mass) is then used to constrain parameters of this simple interaction potential.

preprint2010arXiv

A note on the discovery of a 2M_\odot pulsar

It is conventionally thought that the state equation of dense matter softens and thus cannot result in high maximum mass if pulsars are quark stars, and that a recently discovered $2M_\odot$ pulsar (PSR J1614-2230) may make quark stars to be unlikely. However, this standard point of view would be revisited and updated if quark clustering could occur in cold quark matter because of the strong coupling between quarks at realistic baryon densities of compact stars, and it is addressed that the state equation of clustering quark matter stiffs to support compact stars with maximum mass $M_{\rm max}>2M_\odot$. In this brief note, it is demonstrated that large parameter spaces are allowed for $M_{\rm max}>2M_\odot$ in a Lennard-Jones model of clustered quark matter, and the newly measured highest mass of PSR J1614-2230 would be meaningful to constrain the number of quarks inside single quark-cluster, to be $N_q<\sim 10^3$.

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

Are AXPs/SGRs magnetars?

Anomalous X-ray Pulsars and Soft Gamma-Ray Repeaters have been generally recognized as neutron stars with super strong magnetic fields, namely "magnetars". The "magnetars" manifest that the luminosity in X-ray band are larger than the rotational energy loss rate, i.e. $L_{X}>\dot {E}_{\rm rot}$, and then the radiation energy is coming from the energy of magnetic field. Here it is argued that magnetars may not really exist. Some X-ray and radio observational results are contradicted with the magnetar model. (1) The X-ray luminosity of PSR J1852+0040 is much larger than the rotational energy loss rate ($L_{X}/\dot {E}_{\rm rot}\simeq 18)$, but the magnetic field is just $3.1\times 10^{11}$ G. Does this X-ray radiation energy come from the magnetic field? (2) In contrast to the above, the magnetic fields of radio pulsars J1847-0130 and PSR J1718-3718 are higher than that of AXP 1E 2259+586, why is the radiation energy of those two radio pulsars still coming from rotational energy? Furthermore, the magnetic field of the newly discovered SGR 0418+5729 with the lowest magnetic field is 3.0 e13 G, lower than the critical magnetic field $B_{\rm C}=4.414$ e13 G) (Esposito et al. 2010). (3) Some "magnetars" also emit normal transient radio pulses, what is the essential difference between radio pulsars and the "magnetars"? The observational fact arguments will be presented at first, then we discuss in what situation the conventional method to obtain magnetic field could not be correct.

preprint2010arXiv

Lower bounds of altitudes for pulsar $γ$-ray radiation

Determining radiation location observationally plays a very important role in testing the pulsar radiation models. One-photon pair production in the strong magnetic field, $γ-e^{+}e^{1}$, is one of the important physical processes in pulsar radiation mechanisms. Photons near pulsar surface with sufficient energy will be absorbed in the magnetosphere and the absorption optical depth for these GeV $γ$-ray photons is usually large. In this paper, we include the aberrational, rotational and general relativistic effects and calculate the $γ$-B optical depth for $γ$-ray photons. Then we use the derived optical depth to determine the radiation altitude lower bounds for photons with given energies. As a case study, we calculate the lower bounds of radiation altitudes of Crab pulsar for photons with energy from 5 GeV to 1 TeV.

preprint2010arXiv

Microlensing pulsars

We investigate the possibilities that pulsars act as the lens in gravitational microlensing events towards the galactic bulge or a spiral arm. Our estimation is based on expectant survey and observations of FAST (Five hundred meter Aperture Spherical Telescope) and SKA (Square Kilometer Array), and two different models of pulsar distribution are used. We find that the lensing rate is > 1 event/decade, being high enough to search the real events. Therefore, the microlensing observations focusing on pulsars identified by FAST or SKA in the future are meaningful. As an independent determination of pulsar mass, a future detection of microlensing pulsars should be significant in the history of studying pulsars, especially in constraining the state of matter (either hadronic or quark matter) at supra-nuclear densities. The observations of such events by using advanced optical facilities (e.g., the James Webb Space Telescope and the Thirty Meter Telescope) in future are highly suggested.

preprint2010arXiv

Non-detection in a Fermi/LAT observation of AXP 4U 0142+61: magnetars?

Significant research in compact stars is currently focused on two kinds of enigmatic sources: anomalous X-ray pulsars (AXPs) and soft gamma-ray repeaters (SGRs). Although AXPs and SGRs are popularly thought to be magnetars, other models (e.g. the accretion model) to understand the observations can still not be ruled out. It is worth noting that a non-detection in a Fermi/LAT observation of AXP 4U 0142+61 has been reported recently by Sasmaz Mus & Gogus. We propose here that Fermi/LAT observations may distinguish between the magnetar model and the accretion model for AXPs and SGRs. We explain how this null observation of AXP 4U 0142+61 favors the accretion model. Future Fermi/LAT observations of AXP 1E 1547.0-5408 and AXP 1E 1048.1-5937 are highly recommended.

preprint2010arXiv

Quaking neutron star deriving radiative power of oscillating magneto-dipole emission from energy of Alfvén seismic vibrations

It is shown that depletion of the magnetic field pressure in a quaking neutron star undergoing Lorentz-force-driven torsional seismic vibrations about axis of its dipole magnetic moment is accompanied by the loss of vibration energy of the star that causes its vibration period to lengthen at a rate proportional to the rate of magnetic field decay. Highlighted is the magnetic-field-decay induced conversion of the energy of differentially rotational Alfvén vibrations into the energy of oscillating magneto-dipole radiation. A set of representative examples of magnetic field decay illustrating the vibration energy powered emission with elongating periods produced by quaking neutron star are considered and discussed in the context of theory of magnetars.

preprint2010arXiv

The annular gap model for gamma-ray emission from young and millisecond pulsars

Pulsed high energy radiation from pulsars is not yet completely understood. In this paper, we use the 3D self-consistent annular gap model to study light curves for both young and millisecond pulsars observed by the Fermi Gamma-ray Space Telescope. The annular gap can generate high energy emission for short-period pulsars. The annular gap regions are so large that they have enough electric potential drop to accelerate charged particles to produce gamma-ray photons. For young pulsars, the emission region is from the neutron star surface to about half of the light cylinder radius, and the peak emissivity is in the vicinity of the null charge surface. The emission region for the millisecond pulsars is located much lower than that of the young pulsars. The higher energy gamma-ray emission comes from higher altitudes in the magnetosphere. We present the simulated light curves for three young pulsars (the Crab, the Vela, the Geminga) and three millisecond pulsars (PSR J0030+0451, PSR J0218+4232, PSR J0437-3715) using the annular gap model. Our simulations can reproduce the main properties of observed light curves.

preprint2010arXiv

Toward an understanding of thermal X-ray emission of pulsars

We present a theoretical model for the thermal X-ray emission and cooling of isolated pulsars, assuming that pulsars are solid quark stars. We calculate the heat capacity for such a quark star, and the results show that the residual thermal energy cannot sustain the observed thermal X-ray luminosities seen in typical isolated X-ray pulsars. We conclude that other heating mechanisms must be in operation if the pulsars are in fact solid quark stars. Two possible heating mechanisms are explored. Firstly, for pulsars with little magnetospheric activities, accretion from the interstellar medium or from the material in the associated supernova remnants may power the observed thermal emission. In the propeller regime, a disk-accretion rate ${\dot M}\sim$1% of the Eddington rate with an accretion onto the stellar surface at a rate of $\sim 0.1% {\dot M}$ could explain the observed emission luminosities of the dim isolated neutron stars and the central compact objects. Secondly, for pulsars with significant magnetospheric activities, the pulsar spindown luminosities may have been as the sources of the thermal energy via reversing plasma current flows. A phenomenological study between pulsar bolometric X-ray luminosities and the spin energy loss rates presents the probable existence of a 1/2-law or a linear law, i.e. $L_{\rm bol}^{\infty}\propto\dot{E}^{1/2}$ or $L_{\rm bol}^{\infty}\propto\dot{E}$. This result together with the thermal properties of solid quark stars allow us to calculate the thermal evolution of such stars. Thermal evolution curves, or cooling curves, are calculated and compared with the `temperature-age' data obtained from 17 active X-ray pulsars. It is shown that the bolometric X-ray observations of these sources are consistent with the solid quark star pulsar model.

preprint2009arXiv

Lennard-Jones quark matter and massive quark stars

Quark clustering could occur in cold quark matter because of the strong coupling between quarks at realistic baryon densities of compact stars. Although one may still not be able to calculate this conjectured matter from first principles, the inter-cluster interaction might be analogized to the interaction between inert molecules. Cold quark matter would then crystallize in a solid state if the inter-cluster potential is deep enough to trap the clusters in the wells. We apply the Lennard-Jones potential to describe the inter-cluster potential, and derive the equations of state, which are stiffer than that derived in conventional models (e.g., MIT bag model). If quark stars are composed of Lennard-Jones matter, they could have high maximum masses ($>2M_{\odot}$) as well as very low masses ($<10^{-3}M_{\odot}$). These features could be tested by observations.

preprint2007arXiv

An annular gap acceleration model for $γ$-ray emission of pulsars

If the binding energy of the pulsar's surface is not so high (the case of a neutron star), both the negative and positive charges will flow out freely from the surface of the star. The annular free flow model for $γ$-ray emission of pulsars is suggested in this paper. It is emphasized that: (1). Two kinds of acceleration regions (annular and core) need to be taken into account. The annular acceleration region is defined by the magnetic field lines that cross the null charge surface within the light cylinder. (2). If the potential drop in the annular region of a pulsar is high enough (normally the cases of young pulsars), charges in both the annular and the core regions could be accelerated and produce primary gamma-rays. Secondary pairs are generated in both regions and stream outwards to power the broadband radiations. (3). The potential drop in the annular region grows more rapidly than that in the core region. The annular acceleration process is a key point to produce wide emission beams as observed. (4). The advantages of both the polar cap and outer gap models are retained in this model. The geometric properties of the $γ$-ray emission from the annular flow is analogous to that presented in a previous work by Qiao et al., which match the observations well. (5). Since charges with different signs leave the pulsar through the annular and the core regions, respectively, the current closure problem can be partially solved.

preprint2006arXiv

Is PSR B0943+10 a low-mass quark star?

A recent X-ray observation has shown that the radio pulsar PSR B0943+10, with clear drifting subpulses, has a much smaller polar cap area than that of conventional pulsars with mass of $\sim\msun$ and radius of $\sim10$ km. Zhang et al. (2005) addressed then that this new result conflicts with the standard vacuum gap model. Nonetheless, the discrepancy could be explained if PSR B0943+10 is actually a low-mass quark star. It is found that the potential drop in the open-field-line region of oblique pulsars (i.e., inclination angle $α\neq 0$) might be $\sim 10^2$ times that of aligned pulsars, and that PSR B0943+10 with $α= 12.4^{\rm o}$ could be well above the deathline. We thus conclude that the Ruderman-Sutherland-type vacuum gap model still works well for this pulsar if it is a bare quark star with a mass of $\sim 0.02M_\odot$ and a radius of $\sim 2.6$ km.

preprint2006arXiv

Toward an understanding of the periastron puzzle of PSR B1259-63

Efforts are made to understand the timing behaviors (e.g., the jumps of the projected pulsar semimajor axis in the periastron passages) observed via 13-year timing for PSR B1259-63. In the first, planet-like objects are suggested to orbit around the Be star, which may gravitationally perturb the (probably low mass) pulsar when it passes through periastron. Second, an accretion disk should exist outside the pulsar's light cylinder, which creates a spindown torque on the pulsar due to the propeller effect. The observed negative braking index and the discrepant timing residuals close to periastron could be relevant to the existence of the disk with varying accretion rates. Third, a speculation is presented that the accretion rate may increase on a long timescale in order to explain the negative braking index.

preprint2004arXiv

Quakes in Solid Quark Stars

A starquake mechanism for pulsar glitches is developed in the solid quark star model. It is found that the general glitch natures (i.e., the glitch amplitudes and the time intervals) could be reproduced if solid quark matter, with high baryon density but low temperature, has properties of shear modulus μ= 10^{30~34} erg/cm^3 and critical stress σ_c = 10^{18~24} erg/cm^3. The post-glitch behavior may represent a kind of damped oscillations.

preprint2002arXiv

Bare strange quark stars: formation and emission

Recent achievements of bare strange stars are briefly reviewed. A nascent protostrange star should be bare because of strong mass ejection and high temperature after the supernova detonation flame, and a crust can also hardly form except for a super-Eddington accretion. The magnetosphere of a bare strange star is composed mainly of electron-positron pair plasma, where both inner and outer vacuum gaps work for radio as well as high energy nonthermal emission. A featureless thermal spectrum is expected since no ion is above the quark surface, whilst electron cyclotron lines could appear in some bare strange stars with suitable magnetic fields. Various astrophysical implications of bare strange stars are discussed.

preprint2002arXiv

What can the redshift observed in EXO 0748-676 tell us?

The mass-radius relations for bare and crusted strange stars are calculated with the bag model. Comparing these relations with the observed one derived from redshift, we address that the conclusion, that EXO 0748-676 can not be a strange star, is incorrect. Various strange star models can show that EXO 0748-676 could have a mass of (1.3 ~ 1.7) M_sun and a radius of (8.4 ~ 11.4) km. It is proposed that some nascent strange stars could be bare and have mass \~ 0.1 M_sun, whereas their masses increase during a long accretion history.

preprint2001arXiv

The birth of strange stars and their dynamo-originated magnetic fields

It is shown that protostrange stars (PSSs) can be convective and that there are two possible scenarios describing their turbulence. Besides the local turbulence on the scale which is less than the mean free path of neutrinos, large-scale (~1 km) convection also may occur with properties that are similar to those of convection in protoneutron stars (PNSs). We thus suggest that strange stars can also create dynamo-originated magnetic fields during the deleptonization episode soon after a supernova explosion. Further detailed investigations are needed to see whether or not strange stars and neutron stars can be distinguished according to the differences in dynamo actions in strange quark matter and in neutron matter. The magnetic fields of strange stars and neutron stars may also behave very differently during the accretion-phase when the fields decay.

preprint2001arXiv

The fastest rotating pulsar: A strange star?

According to the observational limits on the radius and mass, the fastest rotating pulsar (PSR 1937+21) is probably a strange star, or at least some neutron star equations of state should be ruled out, if we suggest that a dipole magnetic field is relevant to its radio emission. We presume that the millisecond pulsar is a strange star with {\em much} low mass, small radius, and weak magnetic moment.

preprint2000arXiv

Are there real orthogonal polarization modes in pulsar radio emission?

The Orthogonal Polarization Modes (OPM) have been reported observationally and accepted widely by pulsar researchers. However, no acceptable theory can show the origin of the OPM, which becomes a mystery in pulsar research field. Here a possible way to solve this mystery is presented. We ask a question: Does there exist any real so-called OPM in pulsar radiation? It is proposed in this paper that the `observed OPM' in individual pulses could be the results of depolarization of pulsar radiation and the observational uncertainties originated from polarimeter in observation. A possible way to check this idea is suggested. If the idea is verified, the pulsar research would be influenced significantly in theory and in observation.

preprint2000arXiv

On the inverse Compton scattering model of radio pulsars

Some characteristics of the inverse Compton scattering (ICS) model are reviewed. At least the following properties of radio pulsars can be reproduced in the model: core or central emission beam, one or two hollow emission cones, different emission heights of these components, diverse pulse profiles at various frequencies, linear and circular polarization features of core and cones.

preprint2000arXiv

What if pulsars are born as strange stars?

The possibility and the implications of the idea, that pulsars are born as strange stars, are explored. Strange stars are very likely to have atmospheres with typical mass of $\sim 5\times 10^{-15}M_\odot$ but bare polar caps almost throughout their lifetimes, if they are produced during supernova explosions. A direct consequence of the bare polar cap is that the binding energies of both positively and negatively charged particles at the bare quark surface are nearly infinity, so that the vacuum polar gap sparking scenario as proposed by Ruderman & Sutherland should operate above the cap, regardless of the sense of the magnetic pole with respect to the rotational pole. Heat can not accumulate on the polar cap region due to the large thermal conductivity on the bare quark surface. We test this ``bare polar cap strange star'' (BPCSS) idea with the present broad band emission data of pulsars, and propose several possible criteria to distinguish BPCSSs from neutron stars.