NOVEL CELESTIAL NAVIGATION FOR SATELLITE USING X-RAY PULSARS

The X-ray pulsar-based navigation is a novel technology for the satellite autonomous navigation. The position and the velocity of the satellite are deterimined by using the pulse phases detected at the satellite and predicted by the pulse timing models. With the detected pulse phase, the satellite position with respect to the Earth center can be calculated along the line-of-sight to the pulsar. Using three pulsars, the satellite position in the in- ertial frame can be resolved. The extended Kalman filter （EKF） algorithm is designed to incorporate the range measurements with the satellite dynamics. Simulation verification shows that the proposed algorithm can accu- rately determine the satellite orbit, with the position error less than 100 m. Furthermore, the factors influencing the navigation performance are also discussed.

On time transfer in X-ray pulsar navigation

X-ray pulsar navigation(XPNAV) is a new approach for spacecraft autonomous navigation.The system gets position information utilizing accurate timing methods.Among the timing models,the high-order relativistic effects on the propagated signal must be incorporated to attain precise timing.The time transfer model is provided in detail here in two parts:the time frame transformation and the relativistic effects.

Pulsar/CNS integrated navigation based on federated UKF

In order to improve the autonomous navigation capability of satellite,a pulsar/CNS（celestial navigation system） integrated navigation method based on federated unscented Kalman filter（UKF） is proposed.The celestial navigation is a mature and stable navigation method.However,its position determination performance is not satisfied due to the low accuracy of horizon sensor.Single pulsar navigation is a new navigation method,which can provide highly accurate range measurements.The major drawback of single pulsar navigation is that the system is completely unobservable.As two methods are complementary to each other,the federated UKF is used here for fusing the navigation data from single pulsar navigation and CNS.Compared to the traditional celestial navigation method and single pulsar navigation,the integrated navigation method can provide better navigation performance.The simulation results demonstrate the feasibility and effectiveness of the navigation method.

In-Orbit Calibration and Data Analysis of China＇s First Grazing Incidence Focusing X-Ray Pulsar Telescope

The grazing incidence focusing X-ray pulsar telescope(iFXPT), as the main payload of the X-ray Pulsar Navigation Test Satellite(XPNAV-1), will have great significance on China's space scientific exploration and X-ray pulsar navigation. With PSR B0531+21(Crab Pulsar) as the observation target, the pulsar profile has been recovered based on the data obtained by iFXPT, realizing the main objective of "observing" PSR B0531+21 for the first time in China. This payload mainly consists of the Wolter-I X-ray optics, silicon drift detector, magnetic deflector, electronics, high-energy particle shield and high-stability structures. Currently, the iFXPT, with its good in-orbit performance, has obtained a considerable observation data. The effective area, sensitivity and energy response have been calibrated both on ground and in-orbit, demonstrating a high degree of consistency. Meanwhile, the in-orbit observation data and information for pulsar navigation has also been analyzed simultaneously. As a result, the feasibility of the exploration scheme and the performance of the telescope have been fully validated.

Effects of dynamo magnetic fields on observational properties of accreting millisecond X-ray pulsars

In this paper,we have investigated accreting millisecond X-ray pulsars,which are rapidly rotating neutron stars in low-mass X-ray binaries.These systems exhibit coherent X-ray pulsations that arise when the accretion flow is magnetically channeled to the stellar surface.Here,we have developed the fundamental equations for an accretion disk around accreting millisecond X-ray pulsars in the presence of a dynamo generated magnetic field in the inner part of the disk.We have also formulated the numerical method for the structure equations in the inner region of the disk and the highest accretion rate is enough to form the inner region of the disk,which is overpowered by radiation pressure and electron scattering.Finally,we have examined our results with the effects of dynamo magnetic fields on accreting millisecond X-ray pulsars.

Denoising of X-ray pulsar observed profile using biorthogonal lifting wavelet transform

In X-ray pulsar-based navigation, strong X-ray background noise leads to a low signal-to-noise ratio（SNR） of the observed profile, which consequently makes it very difficult to obtain an accurate pulse phase that directly determines the navigation precision. This signifies the necessity of denoising of the observed profile. Considering that the ultimate goal of denoising is to enhance the pulse phase estimation, a profile denoising algorithm is proposed by fusing the biorthogonal lifting wavelet transform of the linear phase characteristic with the thresholding technique. The statistical properties of X-ray background noise after epoch folding are studied. Then a wavelet-scale dependent threshold is introduced to overcome correlations between wavelet coefficients. Moreover, a modified hyperbola shrinking function is presented to remove the impulsive oscillations of the observed profile. The results of numerical simulations and real data experiments indicate that the proposed method can effectively improve SNR of the observed profile and pulse phase estimation accuracy, especially in short observation durations. And it also outperforms the Donoho thresholding strategy normally used in combination with the orthogonal discrete wavelet transform.

Suzaku observation of Be / X-ray binary pulsar EXO 2030+375

We study the timing and spectral properties of Be/X-ray binary pulsar EXO 2030＋375 using a Suzaku observation taken on 2012 May 23, during a less intense Type I outburst. Pulsations were clearly detected in the X-ray light curves at a barycentric period of 41.2852 s, which suggest that the pulsar is spinning-up. The pulse profiles were found to be peculiar, e.g. unlike those obtained from the earlier Suzaku observation acquired on 2007 May 14. A single-peaked narrow profile at soft X-rays （0.5-10 keV range） changed to a double-peaked broad profile in the 12-55 keV energy range and again reverted back to a smooth single-peaked profile at hard X-rays （55-70 keV range）. The 1.0-100.0 keV broadband spectrum of the pulsar was found to be well described by three continuum models described as （i） a partial covering high energy cut-off power-law model, （ii） a partially absorbed power-law with high- energy exponential rolloff and （iii） a partial covering Negative and Positive power law with EXponential （NPEX） continuum model. Unlike the earlier Suzaku observation during which several low energy emission lines were detected, a weak and narrow Iron Ka emission line at 6.4 keV was only present in the pulsar spectrum during the 2012 May outburst. Non-detection of any absorption like feature in the 1-100 keV energy range supports the claim of the absence of the cyclotron resonance scattering feature in EXO 2030＋375 from the earlier Suzaku observation. Pulse-phase resolved spectroscopy revealed the presence of additional dense matter causing the absence of a second peak from the soft X-ray pulse profiles. The details of the results are described in the paper.

Modification of gravitational redshift of x-ray burst produced by pulsar surface magnetoplasma

In this paper, the propagation of x-ray bursts in the magnetoplasma of pulsar magnetosphere is discussed. The electromagnetic interaction between x-ray bursts and magnetoplasma is described as some geometry. The electromagnetic effects of surface superstrong magnetic field and dynamic effects of outflowing magnetoplasma of pulsars are treated as an optical metric. The Gordon metric is introduced to represent the gravitational metric and optical metric. So the propagation of x-ray bursts in magnetoplasma of pulsars can be described as x-ray bursts transmitting in an effective space characterized by Gordon metric. The modification of gravitational redshift, attributed to the flowing magnetoplasma of pulsars, is obtained and it is shown that the modification is of redshift and can reach the same magnitude as the gravitational redshift for ordinary pulsars.

Detectability of rotation-powered pulsars in future hard X-ray surveys

Recent INTEGRAL/IBIS hard X-ray surveys have detected about 10 young pulsars. We show hard X-ray properties of these 10 young pulsars, which have a luminosity of 10^33 - 10^37 erg s^-1 and a photon index of 1.6-2.1 in the energy range of 20-100 keV. The correlation between X-ray luminosity and spin-down power of Lx ∝ Lsd^1.31 suggests that the hard X-ray emission in rotation-powered pulsars is dominated by the pulsar wind nebula （PWN） component. Assuming spectral properties are similar in 20-100keV and 2-10 keV for both the pulsar and PWN components, the hard X-ray luminosity and flux of 39 known young X-ray pulsars and 8 millisecond pulsars are obtained, and a correlation of Lx ∝ Lsd^1.5 is derived. About 20 known young X-ray pulsars and 1 millisecond pulsars could be detected with future INTEGRAL and HXMT surveys. We also carry out Monte Carlo simulations of hard X-ray pulsars in the Galaxy and the Gould Belt, assuming values for the pulsar birth rate, initial position, proper motion velocity, period, and magnetic field distribution and evolution based on observational statistics and the Lx - Lsd reltions：Lx ∝Lsd^1.31 and Lx∝Lsd^1.5 More than 40 young pulsars （mostly in the Galactic plane） could be detected after ten years of INTEGRAL surveys and the launch of HXMT. So, the young pulsars would be a significant part of the hard X-ray source population in the sky, and will contribute to unidentified hard X-ray sources in present and future hard X-ray surveys by INTEGRAL and HXMT.

Spectral Properties of Anomalous X-ray Pulsars

We examine the spectra of the persistent emission from anomalous X-ray pulsars (AXPs) and their variation with the spin-down rate . Based on an accretion-powered model, the influences of both the magnetic field and the mass accretion rate on the spectral properties of AXPs are addressed. We then investigate the relation between the spectral property of AXPs and mass accretion rate . The result shows that there exists a linear correlation between the photon index and the mass accretion rate: the spectral hardness increases with increasing . A possible emission mechanism for the explanation of the spectral properties of AXPs is also discussed.

Understanding the X-ray spectrum of anomalous X-ray pulsars and soft gamma-ray repeaters

Hard X-rays above 10 keV are detected from several anomalous X-ray pul- sars （AXPs） and soft gamma-ray repeaters （SGRs）, and different models have been proposed to explain the physical origin within the frame of either a magnetar model or a fallback disk system. Using data from Suzaku and INTEGRAL, we study the soft and hard X-ray spectra of four AXPs/SGRs： 1RXS J 170849-400910, 1E 1547.0- 5408, SGR 1806-20 and SGR 0501＋4516. It is found that the spectra could be well reproduced by the bulk-motion Comptonization （BMC） process as was first suggested by Triimper et al., showing that the accretion scenario could be compatible with X- ray emission from AXPs/SGRs. Simulated results from the Hard X-ray Modulation Telescope using the BMC model show that the spectra would have discrepancies from the power-law, especially the cutoff at -200 keV. Thus future observations will allow researchers to distinguish different models of the hard X-ray emission and will help us understand the nature of AXPs/SGRs.

A Two-Temperature Supernova Fallback Disk Model for Anomalous X-ray Pulsars

We present a case study of the relevance of the radially pulsational instability of a two-temperature accretion disk around a neutron star to anomalous X-ray pulsars (AXPs). Our estimates are based on the approximation that such a neutron star disk with mass in the range of 10-6-10-5 M⊙ is formed by supernova fallback. We derive several peculiar properties of the accretion disk instability: a narrow interval of X-ray pulse periods; lower X-ray luminosities; a period derivative and an evolution time scale. All these results are in good agreement with the observations of the AXPs.

Broad-band spectroscopy of the transient X-ray binary pulsar KS 1947+300 during 2013 giant outburst: Detection of pulsating soft X-ray excess component

We present the results obtained from detailed timing and spectral studies of the Be/X-ray binary pulsar KS 1947＋300 during its 2013 giant outburst. We used data from Suzaku observations of the pulsar at two epochs, i.e. on 2013 October 22 （close to the peak of the outburst） and 2013 November 22. X- ray pulsations at - 18.81 s were clearly detected in the light curves obtained from both observations. Pulse periods estimated during the outburst showed that the pulsar was spinning up. The pulse profile was found to be single-peaked up to -10 keV beyond which a sharp peak followed by a dip-like feature appeared at hard X-rays. The dip-like feature has been observed up to -70 keV. The 1-110 keV broad-band spectroscopy of both observations revealed that the best-fit model was comprised of a partially absorbed Negative and Positive power law with EXponential cutoff （NPEX） continuum model along with a blackbody component for the soft X-ray excess and two Gaussian functions at 6.4 and 6.7 keV for emission lines. Both the lines were identified as emission from neutral and He-like iron atoms. To fit the spectra, we included the previously reported cyclotron absorption line at 12.2 keV. From the spin-up rate, the magnetic field of the pulsar was estimated to be -1.2 x 10^12 G and found to be comparable to that obtained from the detection of the cyclotron absorption feature. Pulse-phase resolved spectroscopy revealed the pulsating nature of the soft X-ray excess component in phase with the continuum flux. This confirms that the accretion column and/or accretion stream are the most probable regions of the soft X-ray excess emission in KS1947＋300. The presence of the pulsating soft X-ray excess in phase with continuum emission may be the possible reason for not observing the dip at soft X-rays.

On the magnetic fields of ultraluminous X-ray pulsars

So far quite a few ultraluminous X-ray(ULX) pulsars have been discovered.In this work,we construct a super-Eddington,magnetic accretion disk model to estimate the dipole magnetic field of eight ULX pulsars based on their observed spin-up variations and luminosities.We obtain two branches of dipole magnetic field solutions.They are distributed in the range of B-(0.156-64.5) × 10^(10) G and-(0.275-79.0) × 10^(13) G corresponding to the low-and high-B solutions respectively.The low magnetic field solutions correspond to the state that the neutron stars are far away from the spin equilibrium,and the high magnetic field solutions are close to the spin equilibrium.The ultra-strong magnetic fields derived in Be-type ULX pulsars imply that the accretion mode in Be-type ULX pulsars could be more complicated than in the persistent ULX pulsars and may not be accounted for by the magnetized accretion disk model.We suggest that the transition between the accretor and the propeller regimes may be used to distinguish between the low-and high-B magnetic field solutions in addition to the detection of the cyclotron resonance scattering features.

Pulse phase and doppler frequency estimation of X-ray pulsars under conditions of spacecraft and binary motion and its application in navigation

The pulse phase and doppler frequency estimation of X-ray pulsars in dynamic situations and its application in navigation is a problem that has not been fully investigated. In this paper, solutions are proposed to solve this problem under conditions of spacecraft and binary motion. A high-precision doppler frequency (velocity) measurement model as well as a phase (range) measurement model is established. The averaged maximum-likelihood estimator is developed for the dynamic pulse phase estimation. The pulse phase tracking technique is used in the doppler frequency determination. The tracking filter is redesigned and compared with the existing algorithms. The comparison verifies the advantage of the filter algorithm presented in this pa- per. Unlike traditional views, it is found that in dynamic situations, shorter observation interval lengths will result in higher-accuracy phase and frequency estimates as the tracking filter outputs. A photon-level integrated numerical simulation is performed. Simulation results testify to the validity of the proposed phase and doppler frequency estimation scheme, and show that incorporation of velocity measurements as well as the range ones into the navigation estimator will improve the navigation steady-state performance.

Review of X-ray pulsar spacecraft autonomous navigation

This article provides a review on X-ray pulsar-based navigation(XNAV).The review starts with the basic concept of XNAV,and briefly introduces the past,present and future projects concerning XNAV.This paper focuses on the advances of the key techniques supporting XNAV,including the navigation pulsar database,the X-ray detection system,and the pulse time of arrival estimation.Moreover,the methods to improve the estimation performance of XNAV are reviewed.Finally,some remarks on the future development of XNAV are provided.

X-ray pulsar-based navigation using pulse phase and Doppler frequency measurements

In order to eliminate the impact of the Doppler effects caused by the motion of the spacecraft on the X-ray pulsar-based navigation, an innovative navigation method using the pulse phase and Doppler frequency measurements of the X-ray pulsars is proposed. Given the initial estimate of the spacecraft's state,the real-time photon arrival model is established at the spacecraft with respect to the spacecraft's position and velocity predicted by the orbit dynamic model and their estimation errors. On this basis, a maximum likelihood estimation algorithm directly using the observed photon event timestamps is developed to extract a single pair of pulse phase and Doppler frequency measurements caused by the spacecraft's state estimation error. Since the phase estimation error increases as the observation time increases, we propose a new measurement updating scheme of referring the measurements to the middle time of an observation interval. By using the ground-based simulation system of X-ray pulsar signals, a series of photon-level simulations are performed. The results testify to the feasibility and real-timeliness of the proposed navigation method, and show that the incorporation of the Doppler measurement as well as the pulse phase into the navigation filter can improve the navigation accuracy.

Grouping bi-chi-squared method for pulsar navigation experiment using observations of Rossi X-ray timing explorer

X-ray pulsar-based navigation is a revolutionary technology which is capable of providing the required navigation information in the solar system. Performing as an important pulsar-based navigation technique, the Significance Enhancement of Pulse-profile with Orbit-dynamics(SEPO)can estimate orbital elements by using the distortion of pulse profile. Based on the SEPO technique,we propose a grouping bi-chi-squared technique and adopt the observations of Rossi X-ray Timing Explorer(RXTE) to carry out experimental verification. The pulsar timing is refined to determine spin parameters of the Crab pulsar(PSR B0531+21) during the observation periods, and a short-term dynamic model for RXTE satellite is established by considering the effects of diverse perturbations.Experimental results suggest that the position and velocity errors are 16.3 km(3σ) and 13.3 m/s(3σ) with two adjacent observations split by one day(exposure time of 1.5 ks), outperforming those of the POLAR experiment whose results are 19.2 km(3σ) and 432 m/s(3σ). The approach provided is particularly applicable to the estimation of orbital elements via using high-flux observations.

A report on TypeⅡX-ray bursts from SMC X-1

We study RXTE PCA data for the high mass X-ray binary source SMC X-1 between 2003–10 and 2003–12 when the source was in its high states.The source is found to be frequently bursting which can be seen as flares in lightcurves that occur at a rate of one every 800 s, with an average of 4–5 Type Ⅱ X-ray bursts per hour.We note that typically a burst was short, lasting for a few tens of seconds in addition to a few long bursts spanning more than a hundred seconds that were also observed.The flares apparently occupied 2.5% of the total observing time of 225.5 ks.We note a total of 272 flares with mean FWHM of the flare ～21 s.The rms variability and aperiodic variability are independent of flares.As observed, the pulse profiles of the lightcurves do not change their shape, implying that there is no change in the geometry of an accretion disk due to a burst.The hardness ratio and rms variability of lightcurves exhibit no correlation with the flares.The flare fraction shows a positive correlation with the peak-to-peak ratio of the primary and secondary peaks of the pulse profile.The observed hardening or softening of the spectrum cannot be correlated with the flaring rate but may be due to the interstellar absorption of X-rays as evident from the change in hydrogen column density(n_H).It is found that the luminosity of the source increases with the flaring rate.Considering that the viscous timescale is equal to the mean recurrence time of flares, we fixed the viscosity parameter α ～ 0.16.

Variations in the X-ray eclipse transitions of Cen X-3

We report here an investigation of the X-ray eclipse transitions of the high mass X-ray binary pulsar Cen X-3 in different intensity states. The long term light curve of Cen X-3 obtained with RXTE-ASM spanning more than 5000 d shows strong aperiodic flux variations with low and high states. We have investigated the eclipse transitions of Cen X-3 in different intensity states with data obtained from pointed observations with the more sensitive instruments onboard ASCA, BeppoSAX, XMM- Newton, Chandra and RXTE. We found a very clear trend of sharp eclipse transitions in the high state and longer transitions in the low state. This is a confirmation of this feature first observed with the RXTE-ASM but now observed with much better clarity. From the light curves obtained from several missions, it is seen that the eclipse egress in the low state starts earlier by an orbital phase of 0.02 compared to the eclipse egress in the high state indicating that the observed X-rays originate from a much larger region. We have also performed spectral analysis of the post-eclipse part of each observation. From BeppoSAX observations, the out-of-eclipse X-ray flux is found to differ by a factor of -26 during the high and low intensity states while the eclipse count rates differ by a factor of only ~ 4.7. This indicates that in the low state, there is an additional scattering medium which scatters some of the source photons towards the observer even when the neutron star is completely eclipsed. We could also resolve the three iron line components using XMM-Newton observation in the low state. By comparing the iron line equivalent width during the high and low states, it is seen that the width of the iron line is relatively large during the low state which supports the fact that significant reprocessing and scattering of X-rays takes place in the low state.