Source author record

K. S. Cheng

K. S. Cheng 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

102works
12topics
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

102 published item(s)

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

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

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

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.

preprint2019arXiv

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

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

preprint2016arXiv

Constraining warm dark matter mass with cosmic reionization and gravitational wave

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

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

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

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

preprint2016arXiv

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

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

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

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

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

preprint2015arXiv

Constraining the bulk Lorentz factor from the photosphere emission

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

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

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

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

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.

preprint2015arXiv

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

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

preprint2015arXiv

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

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

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

Chandra detection of a new diffuse X-ray component from the globular cluster 47 Tucanae

In re-analyzing the archival Chandra data of the globular cluster 47 Tucanae, we have detected a new diffuse X-ray emission feature within the half-mass radius of the cluster. The spectrum of the diffuse emission can be described by a power-law model plus a plasma component with photon index $Γ\sim1.0$ and plasma temperature $kT\sim0.2$ keV. While the thermal component is apparently uniform, the non-thermal contribution falls off exponentially from the core. The observed properties could possibly be explained in the context of multiple shocks resulted from the collisions among the stellar wind in the cluster and the inverse Compton scattering between the pulsar wind and the relic photons.

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

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 bubbles as a source of cosmic rays above 10^{15} eV

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

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

Multi-wavelength Emission from the Fermi Bubble I. Stochastic acceleration from Background Plasma

We analyse processes of electron acceleration in the Fermi Bubbles in order to define parameters and restrictions of the models, which are suggested for the origin of these giant radio and gamma-ray structures. In the case of leptonic origin of the nonthermal radiation from the Bubbles, these electrons should be produced somehow in-situ because of relatively short lifetime of high energy electrons, which lose their energy by synchrotron and inverse Compton processes. It has been suggested that electrons in Bubbles may be accelerated by shocks produced by tidal disruption of star accreting onto the central black hole or a process of re-acceleration of electrons ejected by supernova remnants. These processes will be investigated in subsequent papers. In this paper we focus to study in-situ stochastic (Fermi) acceleration by a hydromagnetic/supersonic turbulence, in which electrons can be directly accelerated from the background plasma. We showed that the acceleration from the background plasma is able to explain the observed fluxes of radio and gamma-ray emission from the Bubbles but the range of permitted parameters of the model is strongly restricted.

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

Origin of Hydrogen Ionization in the 1 pc Galactic Central Region

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

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

Gamma-ray Emission from Millisecond Pulsars - An Outergap Perspective

In this review paper we explain the following gamma-ray emission features from the millisecond pulsars. (1)Why is the dipolar field of millisecond pulsars so weak but the magnetic pair creation process may still be able to control the size of the outergap? (2)A sub-GeV pulse component could occur in the vicinity of the radio pulse of millisecond pulsars. (3)Orbital modulated gamma-rays should exist in the black widow systems for large viewing angle.

preprint2013arXiv

Long-term X-ray emission from Swift J1644+57

The X-ray emission from Swift J1644+57 is not steadily decreasing instead it shows multiple pulses with declining amplitudes. We model the pulses as reverse shocks from collisions between the late ejected shells and the externally shocked material, which is decelerated while sweeping the ambient medium. The peak of each pulse is taken as the maximum emission of each reverse shock. With a proper set of parameters, the envelope of peaks in the light curve as well as the spectrum can be modelled nicely.

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

Condensate dark matter stars

We investigate the structure and stability properties of compact astrophysical objects that may be formed from the Bose-Einstein condensation of dark matter. Once the critical temperature of a boson gas is less than the critical temperature, a Bose-Einstein Condensation process can always take place during the cosmic history of the universe. Therefore we model the dark matter inside the star as a Bose-Einstein condensate. In the condensate dark matter star model, the dark matter equation of state can be described by a polytropic equation of state, with polytropic index equal to one. We derive the basic general relativistic equations describing the equilibrium structure of the condensate dark matter star with spherically symmetric static geometry. The structure equations of the condensate dark matter stars are studied numerically. The critical mass and radius of the dark matter star are given by $M_{crit}\approx 2(l_a/1fm)^{1/2}(m_χ/1\;{\rm GeV})^{-3/2}M_{\odot}$ and $R_{crit}\approx 1.1 \times 10^6(l_a/1\;{\rm fm})^{1/2}(m_χ/1\;{\rm GeV})^{-3/2}$ cm respectively, where $l_a$ and $m_χ$ are the scattering length and the mass of dark matter particle, respectively.

preprint2012arXiv

Cosmic histories of star formation and reionization: An analysis with a power-law approximation

With a simple power-law approximation of high-redshift ($\gtrsim3.5$) star formation history, i.e., $\dotρ_*(z)\propto [(1+z)/4.5]^{-α}$, we investigate the reionization of intergalactic medium (IGM) and the consequent Thomson scattering optical depth for cosmic microwave background (CMB) photons. A constraint on the evolution index $α$ is derived from the CMB optical depth measured by the {\it Wilkinson Microwave Anisotropy Probe} (WMAP) experiment, which reads $α\approx2.18\lg{\mathscr{N}_γ}-3.89$, where the free parameter $\mathscr{N}_γ$ is the number of the escaped ionizing ultraviolet photons per baryon. Moreover, the redshift for full reionization, $z_f$, can also be expressed as a function of $α$ as well as $\mathscr{N}_γ$. By further taking into account the implication of the Gunn-Peterson trough observations to quasars for the full reionization redshift, i.e., $6\lesssim z_f \lesssim7$, we obtain $0.3\lesssimα\lesssim1.3$ and $80\lesssim\mathscr{N}_γ\lesssim230$. For a typical number of $\sim4000$ of ionizing photons released per baryon of normal stars, the fraction of these photons escaping from the stars, $f_{\rm esc}$, can be constrained to within the range of $(2.0-5.8)%$.

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

Gravitational effects of condensate dark matter on compact stellar objects

We study the gravitational effect of non-self-annihilating dark matter on compact stellar objects. The self-interaction of condensate dark matter can give high accretion rate of dark matter onto stars. Phase transition to condensation state takes place when the dark matter density exceeds the critical value. A compact degenerate dark matter core is developed and alter the structure and stability of the stellar objects. Condensate dark matter admixed neutron stars is studied through the two-fuid TOV equation. The existence of condensate dark matter deforms the mass-radius relation of neutron stars and lower their maximum baryonic masses and radii. The possible effects on the Gamma-ray Burst rate in high redshift are discussed.

preprint2012arXiv

Modelling the light curves of PSR B1259-63/LS 2883-II.The effects of anisotropic pulsar wind and Doppler-boosting

PSR B1259-63/LS 2883 is a binary system in which a 48-ms pulsar orbits around a Be star in a high eccentric orbit with a long orbital period of about 3.4 yr. It is special for having asymmetric two-peak profiles in both the X-ray and the TeV light curves. Recently, an unexpected GeV flare was detected by $Fermi$ gamma-ray observatory several weeks after the last periastron passage. In this paper, we show that this observed GeV flare could be produced by the Doppler-boosted synchrotron emission in the bow shock tail. An anisotropic pulsar wind model, which mainly affects the energy flux injection to the termination shock in different orbital phase, is also used in this paper, and we find that the anisotropy in the pulsar wind can play a significant role in producing the asymmetric two-peak profiles in both X-ray and TeV light curves. The X-ray and TeV photons before periastron are mainly produced by the shocked electrons around the shock apex and the light curves after periastron are contributed by the emission from the shock apex and the shock tail together, which result in the asymmetric two-peak light curves.

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

Probing Pre-galactic Metal Enrichment with High-Redshift Gamma-Ray Bursts

We explore high-redshift gamma-ray bursts (GRBs) as promising tools to probe pre-galactic metal enrichment. We utilize the bright afterglow of a Pop III GRB exploding in a primordial dwarf galaxy as a luminous background source, and calculate the strength of metal absorption lines that are imprinted by the first heavy elements in the intergalactic medium (IGM). To derive the GRB absorption line diagnostics, we use an existing highly-resolved simulation of the formation of a first galaxy which is characterized by the onset of atomic hydrogen cooling in a halo with virial temperature >10^4 K. We explore the unusual circumburst environment inside the systems that hosted Pop III stars, modeling the density evolution with the self-similar solution for a champagne flow. For minihalos close to the cooling threshold, the circumburst density is roughly proportional to (1+z) with values of about a few cm^{-3}. In more massive halos, corresponding to the first galaxies, the density may be larger, n>100 cm^{-3}. The resulting afterglow fluxes may be detectable with the JWST and VLA in the near-IR and radio wavebands, respectively, out to redshift z>20. We predict that the maximum of the afterglow emission shifts from near-IR to millimeter bands with peak fluxes from mJy to Jy at different observed times. GRBs are ideal tools for probing the metal enrichment in the early IGM, due to their high luminosities and featureless power-law spectra. The metals in the first galaxies produced by the first supernova (SN) explosions are likely to reside in low-ionization stages. We show that if the afterglow can be observed sufficiently early, analysis of the metal lines can distinguish whether the first heavy elements were produced in a PISN, or a core-collapse (Type II) SN, thus constraining the initial mass function of the first stars.

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

Weyl fluid dark matter model tested on the galactic scale by weak gravitational lensing

The higher dimensional Weyl curvature induces on the brane a new source of gravity. This Weyl fluid of geometrical origin (reducing in the spherically symmetric, static configuration to a dark radiation and dark pressure) modifies space-time geometry around galaxies and has been shown to explain the flatness of galactic rotation curves. Independent observations for discerning between the Weyl fluid and other dark matter models are necessary. Gravitational lensing could provide such a test. Therefore we study null geodesics and weak gravitational lensing in the dark radiation dominated region of galaxies in a class of spherically symmetric brane-world metrics. We find that the lensing profile in the brane-world scenario is distinguishable from dark matter lensing, despite both the brane-world scenario and dark matter models fitting the rotation curve data. In particular, in the asymptotic regions light deflection is 18% enhanced as compared to dark matter halo predictions. For a linear equation of state of the Weyl fluid we further find a critical radius, below which brane-world effects reduce, while above it they amplify light deflection. This is in contrast to any dark matter model, the addition of which always increases the deflection angle.

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

Discovery of gamma-ray emission from the supernova remnant Kes 17 with Fermi Large Area Telescope

We report the discovery of GeV emission at the position of supernova remnant Kes 17 by using the data from the Large Area Telescope on board the Fermi Gamma-ray Space Telescope. Kes 17 can be clearly detected with a significance of ~12 sigma in the 1 - 20 GeV range. Moreover, a number of gamma-ray sources were detected in its vicinity. The gamma-ray spectrum of Kes 17 can be well described by a simple power-law with a photon index of ~ 2.4. Together with the multi-wavelength evidence for its interactions with the nearby molecular cloud, the gamma-ray detection suggests that Kes 17 is a candidate acceleration site for cosmic-rays.

preprint2011arXiv

Discovery of GeV gamma-ray emission from PSR B1259-63/LS 2883

The binary system PSR B1259-63/LS 2883 consists of a 47.8 ms radio pulsar that orbits the companion Be star with a period of 3.4 years in a highly eccentric orbit. The system is well sampled in radio, X-rays, and TeV gamma-rays, and shows orbital-phase-dependent variability in all observed frequencies. Here we report on the discovery of >100 MeV gamma-rays from PSR B1259-63/LS 2883 through the 2010 pariastron passage. Using data collected with the Large Area Telescope aboard Fermi from 33 days before pariastron to 75 days after pariastron, PSR B1259-63/LS 2883 is detected at a significance of 13.6 standard deviations. The gamma-ray light curve is highly variable over the above period, with changing photon index that correlates with gamma-ray flux. In particular, two major flares that occur after the pariastron passage were observed. The onset of gamma-ray emission occurs close to, but not at the same orbital phases as, the two disk passages that occur ~1 month before and ~1 month after the pariastron passage. The fact that the GeV orbital light curve is different from that of the X-ray and TeV light curves strongly suggests that GeV gamma-ray emission originates from a different component. We speculate that the observed GeV flares may be resulting from Doppler boosting effects.

preprint2011arXiv

Exploring the dark accelerator HESS J1745-303 with Fermi Large Area Telescope

We present a detailed analysis of the gamma-ray emission from HESS J1745-303 with the data obtained by the Fermi Gamma-ray Space Telescope in the first ~29 months observation.The source can be clearly detected at the level of ~18-sigma and ~6-sigma in 1-20 GeV and 10-20 GeV respectively. Different from the results obtained by the Compton Gamma-ray Observatory, we do not find any evidence of variability. Most of emission in 10-20 GeV is found to coincide with the region C of HESS J1745-303. A simple power-law is sufficient to describe the GeV spectrum with a photon index of ~2.6. The power-law spectrum inferred in the GeV regime can be connected to that of a particular spatial component of HESS J1745-303 in 1-10 TeV without any spectral break. These properties impose independent constraints for understanding the nature of this "dark particle accelerator".

preprint2011arXiv

Gamma-ray emission from the globular clusters Liller 1, M80, NGC 6139, NGC 6541, NGC 6624, and NGC 6752

Globular clusters (GCs) are emerging as a new class of gamma-ray emitters, thanks to the data obtained from the Fermi Gamma-ray Space Telescope. By now, eight GCs are known to emit gamma-rays at energies >100~MeV. Based on the stellar encounter rate of the GCs, we identify potential gamma-ray emitting GCs out of all known GCs that have not been studied in details before. In this paper, we report the discovery of a number of new gamma-ray GCs: Liller 1, NGC 6624, and NGC 6752, and evidence for gamma-ray emission from M80, NGC 6139, and NGC 6541, in which gamma-rays were found within the GC tidal radius. With one of the highest metallicity among all GCs in the Milky Way, the gamma-ray luminosity of Liller 1 is found to be the highest of all known gamma-ray GCs. In addition, we confirm a previous report of significant gamma-ray emitting region next to NGC 6441. We briefly discuss the observed offset of gamma-rays from some GC cores. The increasing number of known gamma-ray GCs at distances out to ~10 kpc is important for us to understand the gamma-ray emitting mechanism and provides an alternative probe to the underlying millisecond pulsar populations of the GCs.

preprint2011arXiv

Gamma-ray properties of globular clusters and the "fundamental planes"

We report on the discovery of gamma-ray emission from several globular clusters (GCs), including Terzan 5, the second known gamma-ray GCs. By now, more than a dozen GCs are known to emit gamma-rays of energies above 100 MeV, thus enabling us to carry out the first detailed correlation study with several cluster properties. We found strong correlations between the observed gamma-ray luminosities and four cluster parameters: stellar encounter rate, metallicity [Fe/H], and energy densities of the soft photons at the cluster locations. These "fundamental planes" of gamma-ray GCs put an intimate relation of the observed gamma-rays to the underlying millisecond pulsar population and have important implications on the origin of the gamma-ray emission of GCs.

preprint2011arXiv

Inflating wormholes in the braneworld models

The braneworld model, in which our Universe is a three-brane embedded in a five-dimensional bulk, allows the existence of wormholes, without any violation of the energy conditions. A fundamental ingredient of traversable wormholes is the violation of the null energy condition (NEC). However, in the brane world models, the stress energy tensor confined on the brane, threading the wormhole, satisfies the NEC. In conventional general relativity, wormholes existing before inflation can be significantly enlarged by the expanding spacetime. We investigate the evolution of an inflating wormhole in the brane world scenario, in which the wormhole is supported by the nonlocal brane world effects. As a first step in our study we consider the possibility of embedding a four-dimensional brane world wormhole into a five dimensional bulk. The conditions for the embedding are obtained by studying the junction conditions for the wormhole geometry, as well as the full set of the five dimensional bulk field equations. For the description of the inflation we adopt the chaotic inflation model. We study the dynamics of the brane world wormholes during the exponential inflation stage, and in the stage of the oscillating scalar field. A particular exact solution corresponding to a zero redshift wormhole is also obtained. The resulting evolution shows that while the physical and geometrical parameters of a zero redshift wormhole decay naturally, a wormhole satisfying some very general initial conditions could turn into a black hole, and exist forever.

preprint2011arXiv

Internal shock model for the X-ray flares of Swift J1644+57

Swift J1644+57 is an unusual transient event, likely powered by the tidal disruption of a star by a massive black hole. There are multiple short timescales X-ray flares were seen over a span of several days. We propose that these flares could be produced by internal shocks. In the internal shock model, the forward and reverse shocks are produced by collisions between relativistic shells ejected from central engine. The synchrotron emission from the forward and reverse shocks could dominate at two quite different energy bands under some conditions, the relativistic reverse shock dominates the X-ray emission and the Newtonian forward shock dominates the infrared and optical emission. We show that the spectral energy distribution of Swift J1644+57 could be explained by internal shock model.

preprint2011arXiv

Modeling the multi-wavelength light curves of PSR B1259-63/SS 2883

PSR B1259-63/SS 2883 is a binary system in which a 48-ms pulsar orbits around a Be star in a high eccentric orbit with a long orbital period of about 3.4 yr. Extensive broadband observational data are available for this system from radio band to very high energy (VHE) range. The multi-frequency emission is unpulsed and nonthermal, and is generally thought to be related to the relativistic electrons accelerated from the interaction between the pulsar wind and the stellar wind, where X-ray emission is from the synchrotron process and the VHE emission is from the inverse Compton (IC) scattering process. Here a shocked wind model with variation of the magnetic parameter $σ$ is developed for explaining the observations. By choosing proper param- eters, our model could reproduce two-peak profile in X-ray and TeV light curves. The effect of the disk exhibits an emission and an absorption components in the X-ray and TeV bands respectively. We suggest that some GeV flares will be produced by Doppler boosting the synchrotron spectrum. This model can possibly be used and be checked in other similar systems such as LS I+61o303 and LS 5039.

preprint2011arXiv

Origin of the Fermi Bubble

Fermi has discovered two giant gamma-ray-emitting bubbles that extend nearly 10kpc in diameter north and south of the galactic center (GC). The existence of the bubbles was first evidenced in X-rays detected by ROSAT and later WMAP detected an excess of radio signals at the location of the gamma-ray bubbles. We propose that periodic star capture processes by the galactic supermassive black hole, Sgr A$^*$, with a capture rate $3\times 10^{-5}$yr$^{-1}$ and energy release $\sim 3\times 10^{52}$erg per capture can produce very hot plasma $\sim 10$keV with a wind velocity $\sim 10^8$cm/s injected into the halo and heat up the halo gas to $\sim 1$keV, which produces thermal X-rays. The periodic injection of hot plasma can produce shocks in the halo and accelerate electrons to $\sim$TeV, which produce radio emission via synchrotron radiation, and gamma-rays via inverse Compton scattering with the relic and the galactic soft photons.

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

Probing the tidal disruption flares of massive black holes with high-energy neutrinos

The recently discovered high-energy transient Swift J164449.3+573451 (Sw J1644+57) is thought to arise from the tidal disruption of a passing star by a dormant massive black hole. Modeling of the broadband emission suggests the presence of a powerful relativistic jet, which contributes dominantly to the observed X-ray emission. Here we suggest that protons can be accelerated to ultra-high energies by internal shocks occurring in the jets, but their flux is insufficient to account for the observed flux of ultra-high energy cosmic rays. High energy protons can produce ~0.1-10 PeV neutrinos through photomeson interactions with X-ray photons. The large X-ray fluence (7x10^-4 erg cm^-2) and high photopion efficiency, together with the insignificant cooling of secondary mesons, result in bright neutrino emission expected from Sw J1644+57 if the jet composition is matter-dominated. One to several neutrinos may be detected by a Km^3-scale detector from one tidal disruption event similar to Sw J1644+57, thereby providing a powerful probe of the composition of the jets.

preprint2011arXiv

The Fermi Bubble as a Source of Cosmic Rays in the Energy Range > 10E15 eV

The {\it Fermi} Large Area Telescope has recently discovered two giant gamma-ray bubbles which extend north and south of the Galactic center with diameters and heights of the order of $H\sim 10$ kpc. We suggest that the periodic star capture processes by the Galactic supermassive black hole Sgr A$^*$, with a capture rate of $τ_{\rm cap}^{-1}\sim 3\times 10^{-5}$ yr$^{-1}$ and an energy release of $W\sim 3\times 10^{52}$ erg per capture, can result in hot plasma injecting into the Galactic halo at a wind velocity of $u\sim 10^8$ cm s$^{-1}$. The periodic injection of hot plasma can produce a series of shocks. Energetic protons in the bubble are re-accelerated when they interact with these shocks. We show that for energy larger than $E> 10^{15}$ eV, the acceleration process can be better described by the stochastic second-order Fermi acceleration. We propose that hadronic cosmic rays (CRs) within the ``knee'' of the observed CR spectrum are produced by Galactic supernova remnants distributed in the Galactic disk. Re-acceleration of these particles in the Fermi Bubble produces CRs beyond the knee. With a mean CR diffusion coefficient in this energy range in the bubble $D_B\sim 3\times 10^{30}$ cm$^2$ s$^{-1}$, we can reproduce the spectral index of the spectrum beyond the knee and within. The conversion efficiency from shock energy of the bubble into CR energy is about 10\%. This model provides a natural explanation of the observed CR flux, spectral indices, and matching of spectra at the knee.

preprint2011arXiv

The Fundamental Plane of Gamma-ray Globular Clusters

We have investigated the properties of a group of $γ$-ray emitting globular clusters (GCs) which have recently been uncovered in our Galaxy. By correlating the observed $γ$-ray luminosities $L_γ$ with various cluster properties, we probe the origin of the high energy photons from these GCs. We report $L_γ$ is positively correlated with the encounter rate $Γ_{c}$ and the metalicity $\left[{\rm Fe/H}\right]$ which place an intimate link between the gamma-ray emission and the millisecond pulsar population. We also find a tendency that $L_γ$ increase with the energy densities of the soft photon at the cluster location. Furthermore, the two-dimensional regression analysis suggests that $L_γ$, soft photon densities, and $Γ_{c}$/$\left[{\rm Fe/H}\right]$ possibly span fundamental planes which potentially provide better predictions for the $γ$-ray properties of GCs.

preprint2011arXiv

The luminosity function of Swift long gamma-ray bursts

The accumulation of {\it Swift} observed gamma-ray bursts (GRBs) gradually makes it possible to directly derive a GRB luminosity function (LF) from observational luminosity distribution, where however two complexities must be involved as (i) the evolving connection between GRB rate and cosmic star formation rate and (ii) observational selection effects due to telescope thresholds and redshift measurements. With a phenomenological investigation on these two complexities, we constrain and discriminate two popular competitive LF models (i.e., broke-power-law LF and single-power-law LF with an exponential cutoff at low luminosities). As a result, we find that the broken-power-law LF could be more favored by the observation, with a break luminosity $L_b=2.5\times10^{52}\rm erg s^{-1}$ and prior- and post-break indices $ν_1=1.72$ and $ν_2=1.98$. For an extra evolution effect expressed by a factor $(1+z)^δ$, if the matallicity of GRB progenitors is lower than $\sim0.1Z_{\odot}$ as expected by some collapsar models, then there may be no extra evolution effect other than the metallicity evolution (i.e., $δ$ approaches to be zero). Alternatively, if we remove the theoretical metallicity requirement, then a relationship between the degenerate parameters $δ$ and $Z_{\max}$ can be found, very roughly, $δ\sim2.4(Z_{\max}/Z_{\odot}-0.06)$. This indicates that an extra evolution could become necessary for relatively high metallicities.

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

Three-dimensional Two-Layer Outer Gap Model: the Third Peak of 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. 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 effect 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

A past capture event at Sagittarius A* inferred from the fluorescent X-ray emission of Sagittarius B clouds

The fluorescent X-ray emission from neutral iron in the molecular clouds (Sgr B) indicates that the clouds are being irradiated by an external X-ray source. The source is probably associated with the Galactic central black hole (Sgr A*), which triggered a bright outburst one hundred years ago. We suggest that such an outburst could be due to a partial capture of a star by Sgr A*, during which a jet was generated. By constraining the observed flux and the time variability ($\sim$ 10 years) of the Sgr B's fluorescent emission, we find that the shock produced by the interaction of the jet with the dense interstellar medium represents a plausible candidate for the X-ray source emission.

preprint2010arXiv

Dynamical Formation of Millisecond Pulsars in Globular Clusters

The cumulative luminosity distribution functions (CLFs) of radio millisecond pulsars (MSPs) in globular clusters (GCs) and in the Galactic field at a frequency of 1.4 GHz have been examined. Assuming a functional form, $N \propto L^q$ where $N$ is the number of MSPs and $L$ is the luminosity at 1.4 GHz, it is found that the CLFs significantly differ with a steeper slope, $q=-0.83 \pm 0.05$, in GCs than in the Galactic field ($q=-0.48 \pm 0.04$), suggesting a different formation or evolutionary history of MSPs in these two regions of the Galaxy. To probe the production mechanism of MSPs in clusters, a search of the possible relationships between the MSP population and cluster properties was carried out. The results of an investigation of 9 GCs indicate positive correlations between the MSP population and the stellar encounter rate and metallicity. This provides additional evidence suggesting that stellar dynamical interactions are important in the formation of the MSP population in GCs.

preprint2010arXiv

Electron-positron energy deposition rate from neutrino pair annihilation on the rotation axis of neutron and quark stars

We investigate the deposition of energy due to the annihilations of neutrinos and antineutrinos on the rotation axis of rotating neutron and quark stars, respectively. The source of the neutrinos is assumed to be a neutrino-cooled accretion disk around the compact object. Under the assumption of the separability of the neutrino null geodesic equation of motion we obtain the general relativistic expression of the energy deposition rate for arbitrary stationary and axisymmetric space-times. The neutrino trajectories are obtained by using a ray tracing algorithm, based on numerically solving the Hamilton-Jacobi equation for neutrinos by reversing the proper time evolution. We obtain the energy deposition rates for several classes of rotating neutron stars, described by different equations of state of the neutron matter, and for quark stars, described by the MIT bag model equation of state and in the CFL (Color-Flavor-Locked) phase, respectively. The electron-positron energy deposition rate on the rotation axis of rotating neutron and quark stars is studied for two accretion disk models (isothermal disk and accretion disk in thermodynamical equilibrium). Rotation and general relativistic effects modify the total annihilation rate of the neutrino-antineutrino pairs on the rotation axis of compact stellar, as measured by an observer at infinity. The differences in the equations of state for neutron and quark matter also have important effects on the spatial distribution of the energy deposition rate by neutrino-antineutrino annihilation.

preprint2010arXiv

Evidence for gamma-ray emission from the low-mass X-ray binary system FIRST J102347.6+003841

The low-mass X-ray binary (LMXB) system FIRST J102347.6+003841 hosts a newly born millisecond pulsar (MSP) PSR J1023+0038 that was revealed as the first and only known rotation-powered MSP in a quiescent LMXB. While the system is shown to have an accretion disk before 2002, it remains unclear how the accretion disk has been removed in order to reveal the radio pulsation in 2007. In this Letter, we report the discovery of gamma-rays spatially consistent with FIRST J102347.6+003841, at a significance of 7 standard deviations, using data obtained by the Fermi Gamma-ray Space Telescope. The gamma-ray spectrum can be described by a power law (PL) with a photon index of 2.9+-0.2, resulting in an energy flux above 200 MeV of (5.5+-0.9)x10^{-12} erg cm^{-2}s^{-1}. The gamma-rays likely originate from the MSP PSR J1023+0038, but also possibly from an intrabinary shock between the pulsar and its companion star. To complement the gamma-ray study, we also re-investigate the XMM-Newton data taken in 2004 and 2008. Our X-ray spectral analysis suggests that a broken PL with two distinct photon indices describes the X-ray data significantly better than a single PL. This indicates that there exists two components and that both components appear to vary with the orbital phase. The evidence for gamma-ray emission conforms with a recent suggestion that gamma-rays from PSR J1023+0038 may be responsible for ejecting the disk material out of the system.

preprint2010arXiv

Fermi Discovery of Gamma-ray Emission from the Globular Cluster Terzan 5

We report the discovery of gamma-ray emission from the Galactic globular cluster Terzan 5 using data taken with the Fermi Gamma-ray Space Telescope, from 2008 August 8 to 2010 January 1. Terzan 5 is clearly detected in the 0.5-20 GeV band by Fermi at ~27σlevel. This makes Terzan 5 as the second gamma-ray emitting globular cluster seen by Fermi after 47 Tuc. The energy spectrum of Terzan 5 is best represented by an exponential cutoff power-law model, with a photon index of ~1.9 and a cutoff energy at ~3.8 GeV. By comparing to 47 Tuc, we suggest that the observed gamma-ray emission is associated with millisecond pulsars, and is either from the magnetospheres or inverse Compton scattering between the relativistic electrons/positrons in the pulsar winds and the background soft photons from the Galactic plane. Furthermore, it is suggestive that the distance to Terzan 5 is less than 10 kpc and > 10 GeV photons can be seen in the future.

preprint2010arXiv

High Redshift Gamma-Ray Bursts: Observational Signatures of Superconducting Cosmic Strings?

The high-redshift gamma-ray bursts (GRBs), GRBs 080913 and 090423, challenge the conventional GRB progenitor models by their short durations, typical for short GRBs, and their high energy releases, typical for long GRBs. Meanwhile, the GRB rate inferred from high-redshift GRBs also remarkably exceeds the prediction of the collapsar model, with an ordinary star formation history. We show that all these contradictions could be eliminated naturally, if we ascribe some high-redshift GRBs to electromagnetic bursts of superconducting cosmic strings. High-redshift GRBs could become a reasonable way to test the superconducting cosmic string model, because the event rate of cosmic string bursts increases rapidly with increasing redshifts, whereas the collapsar rate decreases.

preprint2010arXiv

Inflation and late time acceleration in braneworld cosmological models with varying brane tension

Braneworld models with variable brane tension $λ$ introduce a new degree of freedom that allows for evolving gravitational and cosmological constants, the latter being a natural candidate for dark energy. We consider a thermodynamic interpretation of the varying brane tension models, by showing that the field equations with variable $λ$ can be interpreted as describing matter creation in a cosmological framework. The particle creation rate is determined by the variation rate of the brane tension, as well as by the brane-bulk energy-matter transfer rate. We investigate the effect of a variable brane tension on the cosmological evolution of the Universe, in the framework of a particular model in which the brane tension is an exponentially dependent function of the scale factor. The resulting cosmology shows the presence of an initial inflationary expansion, followed by a decelerating phase, and by a smooth transition towards a late accelerated de Sitter type expansion. The varying brane tension is also responsible for the generation of the matter in the Universe (reheating period). The physical constraints on the model parameters, resulted from the observational cosmological data, are also investigated.

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

Neutrino and Electron-positron Pair Emission from Phase-induced Collapse of Neutron Stars to Quark Stars

We study the energy released from phase-transition induced collapse of neutron stars, which results in large amplitude stellar oscillations. To model this process we use a Newtonian hydrodynamic code, with a high resolution shock-capturing scheme. The physical process considered is a sudden phase transition from normal nuclear matter to a mixed phase of quark and nuclear matter. We show that both the temperature and the density at the neutrinosphere oscillate with time. However, they are nearly 180 degree out of phase. Consequently, extremely intense, pulsating neutrino/antineutrino and leptonic pair fluxes will be emitted. During this stage several mass ejecta can be ejected from the stellar surface by the neutrinos and antineutrinos. These ejecta can be further accelerated to relativistic speeds by the electron/positron pairs, created by the neutrino and antineutrino annihilation outside the stellar surface. We suggest that this process may be a possible mechanism for short Gamma-Ray Bursts.

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.

preprint2010arXiv

The Origin of Gamma-Rays from Globular Clusters

Fermi has detected gamma-ray emission from eight globular clusters. We suggest that the gamma-ray emission from globular clusters may result from the inverse Compton scattering between relativistic electrons/positrons in the pulsar wind of MSPs in the globular clusters and background soft photons including cosmic microwave/relic photons, background star lights in the clusters, the galactic infrared photons and the galactic star lights. We show that the gamma-ray spectrum from 47 Tuc can be explained equally well by upward scattering of either the relic photons, the galactic infrared photons or the galactic star lights whereas the gamma-ray spectra from other seven globular clusters are best fitted by the upward scattering of either the galactic infrared photons or the galactic star lights. We also find that the observed gamma-ray luminosity is correlated better with the combined factor of the encounter rate and the background soft photon energy density. Therefore the inverse Compton scattering may also contribute to the observed gamma-ray emission from globular clusters detected by Fermi in addition to the standard curvature radiation process. Furthermore, we find that the emission region of high energy photons from globular cluster produced by inverse Compton scattering is substantially larger than the core of globular cluster with a radius >10pc. The diffuse radio and X-rays emitted from globular clusters can also be produced by synchrotron radiation and inverse Compton scattering respectively. We suggest that future observations including radio, X-rays, and gamma-rays with energy higher than 10 GeV and better angular resolution can provide better constraints for the models.

preprint2010arXiv

The role of newly born magnetars in gamma-ray burst X-ray afterglow emission: Energy injection and internal emission

Swift observations suggest that the central compact objects of some gamma-ray bursts (GRBs) could be newly born millisecond magnetars. Therefore, by considering the spin evolution of the magnetars against r-mode instability, we investigate the role of the magnetars in GRB X-ray afterglow emission. Besides modifying the conventional energy injection model, we pay particular attention to the internal X-ray afterglow emission, whose luminosity is assumed to track the magnetic dipole luminosity of the magentars with a certain fraction. Following a comparison between the model and some selected observational samples, we suggest that some so-called "canonical" X-ray afterglows including the shallow decay, normal decay, and steeper-than-normal decay phases could be internally produced by the magnetars (possibly through some internal dissipations of the magnetar winds), while the (energized) external shocks are associated with another type of X-ray afterglows. If this is true, from those internal X-ray afterglows, we can further determine the magnetic field strengths and the initial spin periods of the corresponding magnetars.

preprint2009arXiv

Can stellar mass black holes be quark stars?

We investigate the possibility that stellar mass black holes, with masses in the range of $3.8M_{\odot}$ and $6M_{\odot}$, respectively, could be in fact quark stars in the Color-Flavor-Locked (CFL) phase. Depending on the value of the gap parameter, rapidly rotating CFL quark stars can achieve much higher masses than standard neutron stars, thus making them possible stellar mass black hole candidates. Moreover, quark stars have a very low luminosity and a completely absorbing surface - the infalling matter on the surface of the quark star is converted into quark matter. A possibility of distinguishing CFL quark stars from stellar mass black holes could be through the study of thin accretion disks around rapidly rotating quark stars and Kerr type black holes, respectively. Furthermore, we show that the radiation properties of accretion disks around black holes and CFL quark stars are also very similar. However, strange stars exhibit a low luminosity, but high temperature bremsstrahlung spectrum, which, in combination with the emission properties of the accretion disk, may be the key signature to differentiate massive strange stars from black hole.

preprint2009arXiv

Could the compact remnant of SN 1987A be a quark star?

The standard model for Type II supernovae explosion, confirmed by the detection of the neutrinos emitted during the supernova explosion, predicts the formation of a compact object, usually assumed to be a neutron star. However, the lack of the detection of a neutron star or pulsar formed in the SN 1987A still remains an unsolved mystery. In this paper we suggest that the newly formed neutron star at the center of SN1987A may undergo a phase transition after the neutrino trapping time scale (~10 s). Consequently the compact remnant of SN 1987A may be a strange quark star, which has a softer equation of state than that of neutron star matter. Such a phase transition can induce the stellar collapse and result in a large amplitude stellar oscillations. We use a three dimensional Newtonian hydrodynamic code to study the time evolution of the temperature and density at the neutrinosphere. Extremely intense pulsating neutrino fluxes, with submillisecond period and with neutrino energy (> 30 MeV) can be emitted because the oscillations of the temperature and density are out of phase almost 180 degree. If this is true we predict that the current X-ray emission from the compact remnant of SN 1987A will be lower than 10^34 erg s-1, and it should be a thermal bremsstrahlung spectrum for a bare strange star with surface temperature of around ~10^7 K.

preprint2009arXiv

Electron-positron energy deposition rate from neutrino pair annihilation in the equatorial plane of rapidly rotating neutron and quark stars

The neutrino-antineutrino annihilation into electron-positron pairs near the surface of compact general relativistic stars could play an important role in supernova explosions, neutron star collapse, or for close neutron star binaries near their last stable orbit. General relativistic effects increase the energy deposition rates due to the annihilation process. We investigate the deposition of energy and momentum due to the annihilations of neutrinos and antineutrinos in the equatorial plane of the rapidly rotating neutron and quark stars, respectively. We analyze the influence of general relativistic effects, and we obtain the general relativistic corrections to the energy and momentum deposition rates for arbitrary stationary and axisymmetric space-times. We obtain the energy and momentum deposition rates for several classes of rapidly rotating neutron stars, described by different equations of state of the neutron matter, and for quark stars, described by the MIT bag model equation of state and in the CFL (Color-Flavor-Locked) phase, respectively. Compared to the Newtonian calculations, rotation and general relativistic effects increase the total annihilation rate measured by an observer at infinity. The differences in the equations of state for neutron and quark matter also have important effects on the spatial distribution of the energy deposition rate by neutrino-antineutrino annihilation.

preprint2009arXiv

Modeling the radio and optical/NIR afterglows of GRB 980703: a numerical study

Extensive multi-band afterglow data are available for GRB 980703. Especially, its radio afterglow was very bright and was monitored until more than 1000 days after the trigger time. Additionally, there is no obvious special feature, i.e., no rebrightenings, no plateau, and no special steep decay or slow decay in the multi-band afterglow light curves. All these conditions make GRB 980703 a precious sample in gamma-ray burst research. Here we use the observational data of GRB 980703 to test the standard fireball model in depth. It is found that the model can give a satisfactory explanation to the multi-band and overall afterglow light curves. The beaming angle of GRB 980703 is derived as ~0.23 radian, and the circum-burst medium density is ~ 27 cm-3. The total isotropic equivalent kinetic energy of the ejecta is ~ 3.8E52 ergs. A rest-frame extinction of Av ~ 2.5 mag in the host galaxy is also derived.

preprint2009arXiv

Nuclear Interaction Gamma-Ray Lines from the Galactic Center Region

Aims. The accretion of stars onto the central supermassive black hole at the center of the Milky Way is predicted to generate large fluxes of subrelativistic ions in the Galactic center region. We analyze the intensity, shape and spatial distribution of de-excitation gamma-ray lines produced by nuclear interactions of these energetic particles with the ambient medium. Methods. We first estimate the amount and mean kinetic energy of particles released from the central black hole during star disruption. We then calculate from a kinetic equation the energy and spatial distributions of these particles in the Galactic center region. These particle distributions are then used to derive the characteristics of the main nuclear interaction gamma-ray lines. Results. Because the time period of star capture by the supermassive black hole is expected to be shorter than the lifetime of the ejected fast particles against Coulomb losses, the gamma-ray emission is predicted to be stationary. We find that the nuclear de-excitation lines should be emitted from a region of maximum 5$^\circ$ angular radius. The total gamma-ray line flux below 8 MeV is calculated to be $\approx10^{-4}$ photons cm$^{-2}$ s$^{-1}$. The most promising lines for detection are those at 4.44 and $\sim$6.2 MeV, with a predicted flux in each line of $\approx$$10^{-5}$ photons cm$^{-2}$ s$^{-1}$. Unfortunately, it is unlikely that this emission can be detected with the INTEGRAL observatory. But the predicted line intensities appear to be within reach of future gamma-ray space instruments. A future detection of de-excitation gamma-ray lines from the Galactic center region would provide unique information on the high-energy processes induced by the central supermassive black hole and the physical conditions of the emitting region.

preprint2009arXiv

Relativistic Mass Ejecta from Phase-transition-induced Collapse of Neutron Stars

We study the dynamical evolution of a phase-transition-induced collapse neutron star to a hybrid star, which consists of a mixture of hadronic matter and strange quark matter. The collapse is triggered by a sudden change of equation of state, which result in a large amplitude stellar oscillation. The evolution of the system is simulated by using a 3D Newtonian hydrodynamic code with a high resolution shock capture scheme. We find that both the temperature and the density at the neutrinosphere are oscillating with acoustic frequency. However, they are nearly 180$^{\circ}$ out of phase. Consequently, extremely intense, pulsating neutrino/antineutrino fluxes will be emitted periodically. Since the energy and density of neutrinos at the peaks of the pulsating fluxes are much higher than the non-oscillating case, the electron/positron pair creation rate can be enhanced dramatically. Some mass layers on the stellar surface can be ejected by absorbing energy of neutrinos and pairs. These mass ejecta can be further accelerated to relativistic speeds by absorbing electron/positron pairs, created by the neutrino and antineutrino annihilation outside the stellar surface. The possible connection between this process and the cosmological Gamma-ray Bursts is discussed.

preprint2009arXiv

Thin accretion disks around neutron and quark stars

The possibility of observationally discriminating between various types of neutron stars, described by different equations of state of the nuclear matter, as well as differentiating neutron stars from other types of exotic objects, like, for example, quark stars, is one of the fundamental problems in contemporary astrophysics. We consider and investigate carefully the possibility that different types of rapidly rotating neutron stars, as well as other type of compact general relativistic objects, can be differentiated from the study of the emission properties of the accretion disks around them. We obtain the energy flux, the temperature distribution and the emission spectrum from the accretion disks around several classes of rapidly rotating neutron stars, described by different equations of state of the neutron matter, and for quark stars, described by the MIT bag model equation of state and in the CFL (Color-Flavor-Locked) phase, respectively. Particular signatures appear in the electromagnetic spectrum, thus leading to the possibility of directly testing the equation of state of the dense matter by using astrophysical observations of the emission spectra from accretion disks.

preprint2009arXiv

Variation of microphysics in wind bubbles: an alternative mechanism for explaining the rebrightenings in GRB afterglows

Conventionally, long GRBs are thought to be caused by the core collapses of massive stars. During the lifetime of a massive star, a stellar wind bubble environment should be produced. Furthermore, the microphysics shock parameters may vary along with the evolution of the fireball. Here we investigate the variation of the microphysics shock parameters under the condition of wind bubble environment, and allow the microphysics shock parameters to be discontinuous at shocks in the ambient medium. It is found that our model can acceptably reproduce the rebrightenings observed in GRB afterglows, at least in some cases. The effects of various model parameters on rebrightenings are investigated. The rebrightenings observed in both the R-band and X-ray afterglow light curves of GRB 060206, GRB 070311 and GRB 071010A are reproduced in this model.

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

Gravitational Waves from Phase-Transition Induced Collapse of Neutron Stars

(Abridged) We study the gravitational radiation from gravitational collapses of rapidly rotating neutron stars induced by a phase-transition from normal nuclear matter to a mixed phase of quark and nuclear matter in the core of the stars. The study is based on self-consistent three dimensional hydrodynamic simulations with Newtonian gravity and a high resolution shock capturing scheme, and the quadrupole formula of gravitational radiation. The quark matter of the mixed phase is described by the MIT bag model and the normal nuclear matter is described by an ideal fluid equation of state (EOS). 1. We determined the magnitudes of the emitted gravitational waves for several collapse scenarios. 2. We determined the types and frequencies of the fluid oscillation modes excited by the process. In particular, we find that the gravitational wave signals produced by the collapses are dominated by the fundamental quadrupole and quasi-radial modes of the final equilibrium configurations. In some collapse scenarios, we find that the oscillations are damped out within a few dynamical timescales due to the growth of differential rotations and the formation of strong shock waves. 3. We showed that the spectrum of the gravitational wave signals is sensitive to the EOS, implying that the detection of such gravitational waves could provide useful constraints on the EOS of newly born quark stars. 4. For the range of rotation periods we have studied, we found no sign of the development of nonaxisymmetric dynamical instabilities in the collapse process.

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

preprint2005arXiv

Galactic metric, dark radiation, dark pressure and gravitational lensing in brane world models

In the braneworld scenario, the four dimensional effective Einstein equation has extra terms which arise from the embedding of the 3-brane in the bulk. These non-local effects, generated by the free gravitational field of the bulk, may provide an explanation for the dynamics of the neutral hydrogen clouds at large distances from the galactic center, which is usually explained by postulating the existence of the dark matter. We obtain the exact galactic metric, the dark radiation and the dark pressure in the flat rotation curves region in the brane world scenario. Due to the presence of the bulk effects, the flat rotation curves could extend several hundred kpc. The limiting radius for which bulk effects are important is estimated and compared with the numerical values of the truncation parameter of the dark matter halos, obtained from weak lensing observations. There is a relatively good agreement between the predictions of the model and observations. The deflection of photons is also considered and the bending angle of light is computed. The bending angle predicted by the brane world models is much larger than that predicted by standard general relativistic and dark matter models. The angular radii of the Einstein rings are obtained in the small angles approximation. The predictions of the brane world model for the tangential shear are compared with the observational data obtained in the weak lensing of galaxies in the Red-Sequence Cluster Survey. Therefore the study of the light deflection by galaxies and the gravitational lensing could discriminate between the different dynamical laws proposed to model the motion of particles at the galactic level and the standard dark matter models.

preprint2005arXiv

GRB afterglows: deep Newtonian phase and its application

Gamma-ray burst afterglows have been observed for months or even years in a few cases. It deserves noting that at such late stages, the remnants should have entered the deep Newtonian phase, during which the majority of shock-accelerated electrons will no longer be highly relativistic. To calculate the afterglows, we must assume that the electrons obey a power-law distribution according to their kinetic energy, not simply the Lorentz factor.

preprint2005arXiv

Modeling the Optical Afterglow of GRB 030329

The best-sampled afterglow light curves are available for GRB 030329. A distinguishing feature of this event is the obvious rebrightening at around 1.6 days after the burst. Proposed explanations for the rebrightening mainly include the two-component jet model and the refreshed shock model, although a sudden density-jump in the circumburst environment is also a potential choice. Here we re-examine the optical afterglow of GRB 030329 numerically in light of the three models. In the density-jump model, no obvious rebrightening can be produced at the jump moment. Additionally, after the density jump, the predicted flux density decreases rapidly to a level that is significantly below observations. A simple density-jump model thus can be excluded. In the two-component jet model, although the observed late afterglow (after 1.6 days) can potentially be explained as emission from the wide-component, the emergence of this emission actually is too slow and it does not manifest as a rebrightening as previously expected. The energy-injection model seems to be the most preferred choice. By engaging a sequence of energy-injection events, it provides an acceptable fit to the rebrightening at $\sim 1.6$ d, as well as the whole observed light curve that extends to $\sim 80$ d. Further studies on these multiple energy-injection processes may provide a valuable insight into the nature of the central engines of gamma-ray bursts.

preprint2004arXiv

Constraints on extra-dimensions and variable constants from cosmological gamma ray bursts

The observation of the time delay between the soft emission and the high-energy radiation from cosmological gamma ray bursts can be used as an important observational test of multi-dimensional physical theories. The main source of the time delay is the variation of the electromagnetic coupling, due to dimensional reduction, which induces an energy dependence of the speed of light. For photons with energies around 1 TeV, the time delay could range from a few seconds in the case of Kaluza-Klein models to a few days for models with large extra-dimensions. Based on these results we suggest that the detection of the 18-GeV photon $ \sim $4500 s after the keV/MeV burst in GRB 940217 provides a strong evidence for the existence of extra-dimensions. The time delay of photons, if observed by the next generation of high energy detectors, like, for example, the SWIFT and GLAST satellite based detectors, or the VERITAS ground-based TeV gamma-ray instrument, could differentiate between the different models with extra-dimensions.

preprint2004arXiv

Nucleation of quark matter in neutron stars cores

We consider the general conditions of quark droplets formation in high density neutron matter. The growth of the quark bubble (assumed to contain a sufficiently large number of particles) can be described by means of a Fokker-Planck equation. The dynamics of the nucleation essentially depends on the physical properties of the medium it takes place. The conditions for quark bubble formation are analyzed within the frameworks of both dissipative and non-dissipative (with zero bulk and shear viscosity coefficients) approaches. The conversion time of the neutron star to a quark star is obtained as a function of the equation of state of the neutron matter and of the microscopic parameters of the quark nuclei. As an application of the obtained formalism we analyze the first order phase transition from neutron matter to quark matter in rapidly rotating neutron stars cores, triggered by the gravitational energy released during the spinning down of the neutron star. The endothermic conversion process, via gravitational energy absorption, could take place, in a very short time interval, of the order of few tens seconds, in a class of dense compact objects, with very high magnetic fields, called magnetars.

preprint2004arXiv

Pulsar Wind Nebulae and the X-Ray Emission of Non-Accreting Neutron Stars

The general properties of the non-thermal non-pulsed X-ray emission of rotation powered pulsars are investigated in the context of a pulsar wind nebula model. An examination of the observed X-ray emission from a sample of 23 pulsars from ASCA in the energy range between 2-10 keV reveals that the relation of X-ray luminosity, $L_x$, to the pulsar spin down power, $\dot{E}$, is steeper for the non-pulsed component than for the pulsed component. Specifically, $L_{x}^{npul} \propto \dot{E}^{1.4 \pm 0.1}$ for the non-pulsed component, whereas $L_{x}^{pul} \propto \dot{E}^{1.2 \pm 0.08}$ for the pulsed component. The former relation is consistent with emission from a pulsar wind nebula model in which $L_{x}^{npul} \propto \dot{E}^{p/2}$ where $p$ is the power law index of the electron energy distribution. The relation for the pulsed component, on the other hand, is consistent with a magnetospheric emission model.

preprint2003arXiv

The radio afterglow from the giant flare of SGR 1900+14: the same mechanism as afterglows from classic gamma-ray bursts?

A radio afterglow was detected following the 1998 August 27 giant flare from the soft gamma repeater (SGR) 1900+14. Its short-lived behavior is quite different from the radio nebula of SGR 1806-20, but very similar to radio afterglows from classic gamma-ray bursts (GRBs). Motivated by this, we attempt to explain it with the external shock model as invoked in the standard theory of GRB afterglows. We find that the light curve of this radio afterglow is not consistent with the forward shock emission of an ultra-relativistic outflow, which is suggested to be responsible for the initial hard spike of the giant flare. Nevertheless, shock emission from a mildly or sub-relativistic outflow expanding into the interstellar medium could fit the observations. The possible origin for this kind of outflow is discussed, based on the magnetar model for SGRs. Furthermore, we suggest that the presence of an ultra-relativistic fireball from SGR giant flares could be tested by rapid radio to optical follow-up observations in the future.

preprint2002arXiv

Beaming Effects in Gamma-Ray Bursts

Based on a refined generic dynamical model, we investigate afterglows from jetted gamma-ray burst (GRB) remnants numerically. In the relativistic phase, the light curve break could marginally be seen. However, an obvious break does exist at the transition from the relativistic phase to the non-relativistic phase, which typically occurs at time 10 to 30 days. It is very interesting that the break is affected by many parameters, especially by the electron energy fraction (xi_e), and the magnetic energy fraction (xi_B^2). Implication of orphan afterglow surveys on GRB beaming is investigated. The possible existence of a kind of cylindrical jets is also discussed.

preprint2001arXiv

The Radio and Gamma-Ray Luminosities of Blazars

Based on the $γ$-ray data of blazars in the third EGRET catalog and radio data at 5 GHz, we studied the correlation between the radio and $γ$-ray luminosities using two statistical methods. The first method was the partial correlation analysis method, which indicates that there exist correlations between the radio and $γ$-ray luminosities in both high and low states as well as in the average case. The second method involved a comparison of expected $γ$-ray luminosity distribution with the observed data using the Kolmogorov-- Smirnov (KS) test. In the second method, we assumed that there is a correlation between the radio and $γ$-ray luminosities and that the $γ$-ray luminosity function is proportional to the radio luminosity function. The KS test indicates that the expected gamma-ray luminosity distributions are consistent with the observed data in a reasonable parameter range. Finally, we used different $γ$-ray luminosity functions to estimate the possible 'observed' $γ$-ray luminosity distributions by GLAST.

preprint1995arXiv

Statistics Of The Burst Model At Super-critical Phase

We investigate the statistics of a model of type-I X-ray burst [Phys. Rev. E, {\bf 51}, 3045 (1995)] in its super-critical phase. The time evolution of the burnable clusters, places where fire can pass through, is studied using simple statistical arguments. We offer a simple picture for the time evolution of the percentage of space covered by burnable clusters. A relation between the time-average and the peak percentage of space covered by burnable clusters is also derived.