Effect of Ephemeris on Pulsar Timing and Navigation Accuracy Based on X-ray Pulsar Navigation-I Data

: Solar system ephemeris is very important for pulsar timing and navigation. In order to explore the effect of different precision ephemerides on X-ray pulsar timing and navigation, the differences between timing and navigation results with four JPL Development Ephemerides based on the data of X-ray pulsar navigation-I (XPNAV-I) were compared and analyzed in this paper. For pulsar timing, the ephemeris has a systematic effect on time scale conversion (nanosecond difference), light-travel delay (millisecond difference) and timing residuals (microsecond difference)


Introduction
Pulsars are fast-rotating neutron stars in the universe with stable rotation periods, and the long-term stability is higher than that of the most precise atomic clock at present [1,2].Since the first pulsar was discoved in 1967, the study of its physical properties has been a hot spot in the fields of radio astronomy and space high-energy radiation.Many radioobservation arrays and space high-energy radiation detectors around the world have long-term observation records of pulsars.The study of pulsars can be carried out to verify the relativistic effects, detect low-frequency gravitational waves and measure the mass of planets in the solar system [2][3][4]; At the same time, the pulsar time scale based on the ground radio pulsar timing array and spatial X-ray timing data is an important supplement to the integrated positioning, navigation and timing system (PNT), and the X-ray pulsar based navigation (XPNAV) with pulsar X-ray signal as the source is the most promising technology to achieve the long-time and high-precision Autonomous Navigation of deep space probe [5][6][7][8][9][10][11].The research on pulsar timing and XPNAV is of great significance for space exploration by attracting the attention of the international community.
Since the concepts of pulsar timing and navigation were proposed in the 1970s, many researchers have demonstrated their specific implementation methods and feasibility by developing a set of relatively mature theoretical methods [9][10][11][12][13][14][15].At the same time, related satellites in orbit experiments were also gradually carried out [7,11].At present, the most effective is the Neutron-star Interior Composition ExploreR/Station Explorer for X-ray Timing and Nav-igation Technology (NICER/SEXTANT) project carried out by NASA on the International Space Station in 2018, with a positioning error of better than 5 km [11]; Chinese institutions and researchers also carried out timing or satellite orbit determination experiments based on the observation data of the Insight-Hard X-ray Modulation Telescope (Insight-HXMT), the POLAR Gamma Ray Burst Polarimeter on the TianGong-2 Spacelab and the XPNAV-I [7,[16][17][18][19].The results verified the feasibility of X-ray pulsar navigation.
For a spacecraft flying in the solar system, to achieve pulsar navigation or improve its positioning and velocity measurement accuracy, precise solar system ephemeris is extremely necessary to correct the error terms caused by the relativistic effect and Doppler frequency shift [14].At present, many institutions around the world have released the solar system ephemerides for different purposes, among which the Development Ephemerides (DE) released by NASA's Jet Propulsion Laboratory (JPL) has the highest comprehensive accuracy and is the most widely used [20].Since JPL released the first ephemeris DE102 in 1981, 18 different versions of ephemerides have been released successively.Each ephemeris has certain differences in the observation data, inertial reference frame, time span, usage, etc., which leads to differences in the use of different ephemerides to calculate the position or speed of celestial bodies, and the position differences may reach tens of kilometers [21].
For pulsar navigation or timing, the effect of ephemeris on the accuracy of results is complex, and the noise components caused by different celestial errors are not independent but red noise coupled with each other.In recent years, there have still not been many studies in this field.Some scholars regard the ephemeris error as a geometric error, study its impact on the near-Earth satellite navigation and orbit determination system based on X-ray pulsar, and regard it as a part of the state vector processed by the filtering method [22,23]; Ning et al. [24] proposed a celestial/pulsar navigation method using time and ephemeris estimation.In order to eliminate the influence of the ephemeris errors of a planet's moon, the ephemeris of the planet's moons is augmented to the state and estimated in real time.Some of the latest studies also consider the influence of ephemeris on timing and positioning and propose improved methods [16,25].In the research and application of pulsar navigation or timing, one of the important basic works is to determine the relevant pulsar rotation parameters and establish and update the navigation database in real time (especially for relatively young pulsars such as Crab pulsar, undetectable micro-jumps in pulsar rotation may cause serious "timing noise", so it is very necessary to obtain the latest timing parameters [26]).However, the parameters released in different periods or institutions may be fitted based on different ephemerides.When using these parameters, the internal self-consistency of the system should also be considered.
This paper aims at the problems and effects of ephemeris on pulsar navigation and timing.Firstly, the commonly used X-ray pulsar navigation and timing methods and data processing processes were summarized and analyzed, including time scale conversion, light-travel delay calculation, etc.On this basis, based on the observation data of XPNAV-I, the differences between the four commonly used JPL ephemerides DE200, DE405, DE421 and DE430 used in time scale conversion and light-travel delay calculation were calculated and compared, and the effect of ephemeris difference on pulse profiles and timing residuals were analyzed (For the description of XPNAV-I data and the four ephemerides, see reference [27] and Appendix A).Because there is no observation data that can realize navigation measurement at present, based on the orbit determination method of single pulsar observation data, the influence of each ephemeris on orbit determination was calculated and analyzed, and some suggestions on the use of ephemeris in pulsar navigation and timing are given according to the experimental results.

Effect of Ephemeris on Pulsar Timing
Using the high long-term stability of pulsar rotation to establish pulsar time to assist ground atomic time thereby improves the stability and reliability of ground time reference, which has a high practical value [8].At the same time, pulsar timing is also the basis of pulsar navigation, and quality of timing results directly affects the accuracy of spacecraft positioning and velocity measurement [7].At present, the processing flow of observation data is generally the first to convert the detector's local recording time (Universal Time Coordinated(UTC) or Terrestrial Time(TT)) to the Barycentric Coordinate Time (TCB) or Barycentric Dynamical Time (TDB).The next step is to calculate the time delay from the detector to the solar system barycentric (SSB) at the photon arrival time according to the space-time coordinate information of the detector.Then, the phase of the corresponding time according to the pulsar's timing model is calculated, and the pulse profile is obtained by folding the discrete phase.Finally, the phase delay is obtained by cross-correlation operation with the standard pulse profile, and then the time of arrival (TOA) of the pulse is calculated [5,17].The time scale conversion and time delay calculation can be expressed by the following formula: where t SSB is the time when the photon or pulse reaches the SSB; t SC is the local recording time of the detector; ∆E represents the Einstein delay caused by time scale transformation; ∆R represents Romer delay, i.e., the geometric delay between detector and SSB; ∆S represents the Shapiro delay, which is the relativistic effect delay caused by the gravity of all celestial bodies in the solar system; ∆P is the parallax delay caused by the change in spacecraft position.Since the subsequent experimental analysis is mainly based on the Crab observation data of XPNAV-I, for the observation data of a single pulsar, only the time delay between the detector and SSB needs to be considered, rather than the extrasolar delay such as a double pulsar system delay and interstellar delay.In addition, the standard delay formula should also include the dispersion delay of the solar system, but it has little impact on high-energy X-ray observation, and the delay generated by different ephemerides can be considered equal, so it is ignored.

Effect on Time Scale Conversion
Fields such as modern astronomy and geodesy use different time scales for different purposes.For pulsar observations, TT or UTC with the same scale as International Atomic Time(TAI) and only linear transformation is generally adopted; SSB is the inertial reference origin for the study of the motion of various celestial bodies in the solar system, and the corresponding time scale is TCB or TDB.For XPNAV-I satellite and JPL ephemeris, it is the conversion between UTC time and TDB time, and the standard conversion equation is as follows [5,28,29], i. UTC→TAI: where LS represents the leap seconds between UTC and TAI, which can be obtained from the IERS Bulletin.
ii. TAI→TT: TT = TAI + 32.184 s (3) 32.184 s is the origin difference between TAI and TT. iii.TT→TCG: where L G represents the change rate of TCG relative to TT; JD represents Julian day; T 0 = 2, 443, 144.5003725, which is the origin of time (i.e., TT at 1977 January 1.0 TAI).iv.TCG→TCB: where v e is the velocity vector of the Earth centroid relative to SSB, r SC is the position vector of the satellite relative to the Earth centroid, and U ext (r e ) is the scalar potential of other celestial bodies in the solar system except the Earth at the Earth centroid; O c −4 represents the high-order delay term, and its calculation formula is where U ext (r e ) represents the vector potential corresponding to U ext (r e ); O c −4 is less than the order of 10 −16 , which can be ignored for calculations that do not require very high precision.The above transformations are in the first-order post-Newtonian system, and the asteroid correction term is considered.The integral operation is very time-consuming, and the following linear formula can be used for approximate calculation, where L C represents the average change rate of TCB relative to TCG, and P(TT)−P(T 0 ) can be calculated according to the corresponding time ephemeris.v. TCB→TDB where L B represents the average rate of change in TDB relative to TCB.The final value obtained by TDB−TT is the delay generated by the time scale transformation, which is the Einstein delay term in Formula (1).It can be found that there are differences in the conversion process between TCG and TCB (when two quasi-inertial system coordinates) due to the use of different ephemerides; in general, the higher the ephemeris accuracy, the more accurate the calculated delay.For the values of the above constants, there are slight differences in different ephemerides; for example, L TE200 C = 1.48082685703 × 10 −8 , L TE405 C = 1.48082685594 × 10 −8 (the difference value is 1.09 × 10 −17 , superscript TE200 and TE405, respectively, represent the time ephemeris calculated based on DE200 and DE405).
The International Earth Rotation Service (IERS) released the time ephemeris TE200 and TE405, but the delay of DE421 and DE430 needs to be obtained by integration according to Formula (5).Table 1 shows the values at different time scales and the Einstein delay of the start time of the observation data (using the coordinate value corresponding to the time of XPNAV-I) calculated according to the four ephemerides.Figure 1 shows the difference between the Einstein delay calculated by DE200 at the initial moment of the 35 groups of photon recording files and the corresponding values of the other three ephemerides.
It was found that the ephemeris has little influence on the calculation of the Einstein delay.During the data observation period adopted in this paper, the difference between the Einstein delays calculated by each ephemeris is less than 1 ns, among which the difference between DE405, DE421 and DE430 is smaller, and the difference between the latest two ephemeris DE421 and 430 is less than 10 −10 s.In reference [30], the long-term difference between TE200 and TE405 time ephemeris was analyzed, and the maximum difference was less than 5 ns in this century.The celestial mass parameters used in the above results are the actual parameters of each ephemeris.If unified mass parameters are adopted, the difference value will be smaller.It was found that the ephemeris has little influence on the calculation of the Einstein delay.During the data observation period adopted in this paper, the difference between the Einstein delays calculated by each ephemeris is less than 1 ns, among which the difference between DE405, DE421 and DE430 is smaller, and the difference between the latest two ephemeris DE421 and 430 is less than 10 −10 s.In reference [30], the long-term difference between TE200 and TE405 time ephemeris was analyzed, and the maximum difference was less than 5 ns in this century.The celestial mass parameters used in the above results are the actual parameters of each ephemeris.If unified mass parameters are adopted, the difference value will be smaller.

Effect on Light-Travel Delay
After the time scales at the spacecraft and SSB are unified, the spatial light-travel delay can be calculated according to the geometric relationship between them, that is, ∆R, ∆S and ∆P in Formula (1) (observation formula in pulsar navigation or timing).Formula ( 9) is the solar system light-travel delay for X-ray pulsar observation [5,10,12] where n is the unit vector of the pulsar in the Barycentric Celestial Reference System (BCRS); r is the position vector of the satellite relative to the SSB at the time of observation; v is the velocity vector of the pulsar relative to SSB; ∆t = t-t 0 is the difference between the observation time and the reference epoch; R 0 is the distance between the pulsar and the SSB at the reference epoch; G is Newton's gravitational constant and M k is the mass of the kth solar system celestial body; r k is the position vector of the satellite relative to the kth solar system celestial body at the observation time, and r k is the modulus of r k ; N is the number of solar system celestial bodies for which Shapiro delay is to be calculated; r is the distance from the satellite to the Sun at the observation time; M S is the solar mass; ψ is the angle of the Sun and pulsar relative to the satellite at the observation time.The descriptions of the items on the right side of Formula ( 9) are as follows: ⅰ.Solar system Roemer delay

Effect on Light-Travel Delay
After the time scales at the spacecraft and SSB are unified, the spatial light-travel delay can be calculated according to the geometric relationship between them, that is, ∆R, ∆S and ∆P in Formula (1) (observation formula in pulsar navigation or timing).Formula ( 9) is the solar system light-travel delay for X-ray pulsar observation [5,10,12] where n is the unit vector of the pulsar in the Barycentric Celestial Reference System (BCRS); r is the position vector of the satellite relative to the SSB at the time of observation; v is the velocity vector of the pulsar relative to SSB; ∆t = t − t 0 is the difference between the observation time and the reference epoch; R 0 is the distance between the pulsar and the SSB at the reference epoch; G is Newton's gravitational constant and M k is the mass of the kth solar system celestial body; r k is the position vector of the satellite relative to the kth solar system celestial body at the observation time, and r k is the modulus of r k ; N is the number of solar system celestial bodies for which Shapiro delay is to be calculated; r s is the distance from the satellite to the Sun at the observation time; M S is the solar mass; ψ is the angle of the Sun and pulsar relative to the satellite at the observation time.The descriptions of the items on the right side of Formula (9) are as follows: i. Solar system Roemer delay The solar system Roemer delay includes the geometric delay term caused by the spacecraft position and the pulsar proper motion.Now define r P to represent the projection of the vector r in the n direction, that is, r P = n•r; r V represents the vector perpendicular to n, r V = r − r p n.If the pulsar 3D proper motion is assumed to be constant, i.e., v∆t ≡ l, and the items with delay less than 1 ns are ignored, the Roemer delay can be simplified as where the first item of the above formula corresponds to the first item on the right of Formula ( 9), which represents the zero-order Roemer delay caused by the radial position motion of the spacecraft.The second item corresponds to the second item of Formula (9), which is the time delay of the solar system caused by the pulsar radial proper motion.The third item corresponds to the third item of Formula ( 9), which is the second-order Roemer delay of the solar system caused by the transverse direction induced by radial motion and can be regarded as the correction of the second item.
ii. Solar system Parallax delay Parallax delay corresponds to the fourth item on the right side of Formula (9), which is the curvature delay term caused by the inconsistency of radial components caused by the position change of the spacecraft during continuous observation, which can be simplified as iii.Solar system Shapiro delay The fifth and sixth items on the right side of Formula ( 9) represent the first-order and second-order Shapiro delay of the solar system, respectively, which are the relativistic effect delay caused by the space-time curvature caused by the gravitational field of all celestial bodies in the solar system.If only the gravitational delay caused by the Sun is considered, the relevant items can be simplified as follows It was found that each item of calculating light-travel delay involves the use of ephemeris, so the influence of ephemeris is particularly significant.For near-Earth satellite or ground-based radio observations, the spatial delay error caused by the Earth ephemeris error and the relativistic effect caused by massive celestial bodies such as the Sun and Jupiter is the error items that must be considered.A large error leads to the inability to extract the pulse profile that conforms to the periodic sequence.Table 2 shows the delay items of the XPNAV-I satellite at the start time of observation calculated based on four ephemerides.Figure 2 shows the difference between the light-travel delay calculated by the DE200 at the initial moment of the 35 groups of photons recording files and the corresponding values of the other three ephemerides.By comparing the results, it was found that Roemer delay is the main component of light-travel delay and the delay with the largest difference in calculation results between different ephemerides.The difference between the Roemer delay calculated based on DE200 and the values of the other three ephemerides is greater than 1 ms.The difference between Roemer delay calculated by DE405, DE421 and DE430 is less than 10 µs.The difference between the two latest ephemeris Roemer delay is 0.1 microsecond.The ephemeris has little effect on parallax delay and Shapiro delay, and the difference between the results is in the order of picoseconds.
By comprehensively comparing Figures 1 and 2, it can be found that the time difference shows a linear trend as a whole during the observation period, but the difference between Einstein delay in Figure 1 shows a small jump at a few points, which may be related to the change in satellite position (or there is some error in the satellite position files).It is a pity that this paper does not compare the long-term change trend of the difference, but the simulation analysis of similar satellites shows that Figures 1 and 2 conform to the long-term change trend.

Comprehensive Impact Analysis on Timing Results
After time scale conversion and light-travel delay calculation of the original photons sequence observed by the spacecraft, the phase of the corresponding time can be calculated according to the periodic information of the observed pulsar.The phase calculation formula is the following Taylor expansion [31]: where t 0 is the reference epoch; ϕ 0 is the rotation phase of the reference epoch; f 0 , .f 0 and .. f 0 are the rotation frequency and its first and second derivatives.After calculating the phase, the pulse profile in the observation period can be obtained according to the epoch folding method.By comparing the folded pulse profile p with the standard pulse profile p, the phase delay of the pulse can be obtained, and then the initial phase can be obtained.The cross-correlation method is the most commonly used phase comparison method at present.For discrete pulse, its cross-correlation function is: where ϕ is the phase difference between two profiles; N b is the number of pulse bins.The phase delay shall maximize the cross-correlation function.If the starting point of the folding period is the starting point of the observation period, the estimation of the initial phase is: Based on the above analysis, it was found that the ephemeris is an important variable to be considered in time scale conversion and spatial light-travel delay calculation, and its influence does not need to be considered too much in phase calculation and phase comparison.In the actual timing applications, if the actual period of pulsars is not searched according to the photon sequence, but the timing parameters published by some institutions (including period, phase and its derivatives), it is necessary to pay attention to the ephemeris used in the data.For example, the Crab timing parameters published by the British Jodrell Bank Observatory use DE200, and the ephemeris used for similar data of the Australian Telescope National Facility (ATNF) is DE405 [32,33].The mixing of several ephemerides in the timing process may produce a certain systematic deviation.

Effect of Ephemeris on Satellite Navigation and Orbit Determination
According to the traditional navigation and positioning theory, in Formula (9), if the timing result on the right side is known, r can be calculated iteratively.Without considering the clock difference, the absolute spatial positioning can be achieved by observing three pulsars simultaneously [10,12].However, considering the limited relevant research at present, there is no space experiment using this ideal model for navigation verification.In general, a single detector can be used to observe a pulsar for a long time, and the satellite orbit can be determined based on the orbital-dynamics theory, or a single detector sequentially observes multiple pulsars and then synchronizes the observations to the same epoch, which can also achieve positioning in theory.A similar approach was used in a successful verification experiment conducted by NASA in 2018 [11].
Considering the influence of ephemeris, if the first model is applied, its influence is similar to pulsar timing, and the error value it generates in the positioning equation can be calculated according to the law of error propagation.Wang and Xu et al. conducted some research and analysis on this.For the single satellite observation model, considering that the observation period is generally long, it is difficult to model the error caused by the ephemeris, but it is also possible to conduct experiments to compare the effects of different ephemerides on the final orbit determination results to make an effective evaluation of the ephemeris.

Positioning Error Modeling
Combined with the analysis of reference [22,23], considering the systematic deviation caused by ephemeris error, the measurement model of the X-ray pulsar navigation system can be modified as: where ∆τ SC is the difference between the measured TOA and the predicted TOA, and ∆τ is the time delay caused by the deviation of the real position of the satellite; B k is the systematic deviation caused by the ephemeris error of the kth planet or celestial body (the serial number similar to the JPL ephemeris), and N is the number of celestial bodies considered; ∆ε is the time delay caused by other measurement errors.The specific expression of B k is related to the actual navigation situation and the influence quantity of the celestial body, which can be derived from the timing formula.
The basis of the analysis in this paper is to know the position error of celestial bodies in three-dimensional space coordinates s and analyze them one by one.However, to achieve high-precision timing or positioning, the disturbance of large mass celestial bodies and even some asteroids in the solar system needs to be considered.Moreover, the comprehensive influence of various errors is not a simple linear superposition but a nonlinear system problem.In order to illustrate its impact, more in-depth research is needed.

Effect of Ephemeris on Single Pulsar Navigation or Orbit Determination
Considering the limited observation data of X-ray pulsar and the immature navigation technology, the orbit determination of detector based on single X-ray pulsar observation data combined with spacecraft orbit dynamics model has certain application value in the verification of X-ray pulsar navigation theory and the performance test of the detector.For near-Earth satellites, the kinematic equation of spacecraft in an inertial coordinate system is [34]: ..
where r SC/E , .
r SC/C and ..
r SC/E are the position, velocity and acceleration vectors of the spacecraft relative to the Geocentric Celestial Reference System (GCRS) at time t, respectively; p is the kinetic model parameter adopted.The position vector of spacecraft in BCRS can be deduced from r SC/E and ephemeris, that is, r in Formula (9).The two formulas can be combined to solve the position vector according to the pulsar timing observations.The principal term of the formula of motion is the two-body motion of the first item on the right side; the second item represents the sum of the perturbation accelerations acting on the Universe 2022, 8, 360 9 of 18 spacecraft, in which the gravitational perturbation produced by some celestial bodies is a non-negligible amount.
In most dynamic orbit determination analyses and experimental simulations, pulse TOA is still taken as the observation measurement.However, for most detectors, including XPNAV-I, the observation pulse with a high signal-to-noise ratio (SNR) cannot be obtained in a short time, so a stable and continuous TOA observation measurement cannot be obtained.In order to solve this problem, an autonomous orbit determination algorithm based on normal plane constraint was proposed in reference [17], and the navigation test was carried out by using the three-month Crab observation data of XPNAV-I; the accuracy of orbit determination at the control point was 38.4 km.Since most of the published data were distributed in December 2016, the line-of-sight direction of the pulsar is almost perpendicular to the satellite orbit plane, which is not conducive to the orbit correction of this method.
Considering the observational characteristics of the experimental data, the pulsar navigation method "Significance Enhancement of Pulse-profile with Orbit-dynamics" (SEPO) in reference [18] was adopted, instead of selecting the pulse TOA as the measurement for pulsar navigation, the observed pulse profile was directly analyzed.Since the pulse profile and the orbit are closely coupled, and the simulated profile is highly similar to the actual profile and is also modulated by the orbit, the deformation of the profile can reflect the significant change in the orbit.The significance χ 2 of the pulse profile is defined as follows: among them, p i is the pulse profile generated by simulation, and p i is the standard pulse profile of the original observation data.When the predicted orbit is closer to the real orbit, the closer the pulse profile is to the real profile, and the greater the significance χ 2 is.Zheng et al. applied this method to the orbit determination analysis of high-energy space detectors such as POLAR of TG-2 space station, Insight-HXMT Satellite, etc., and verified the effectiveness of this method [18,19].

Experimental Results and Discussion
The above mainly expounds on the commonly used X-ray pulsar navigation and timing methods.When analyzing the influence of ephemeris on timing, the differences in different ephemerides used in time scale and optical delay calculation are analyzed based on XPNAV-I observation data.By using the same data, this section mainly calculates and compares the effects of different ephemerides on pulse profile, timing residuals and single pulsar orbit determination and gives some suggestions on the use of ephemeris based on the results.

Pulse Profile
The pulse profile is obtained by the epoch folding method, and the timing parameters are the crab timing parameters regularly released by Jodrell Bank Observatory.(In order to minimize timing noise, it is necessary to use parameters that are as consistent as possible with the date of the observed data.During the observation period, 2016/11/17~12/19, Jodrell has the data of reference epoch 11/23 and 12/16, while the reference epoch of ATNF similar data is very old.See Appendix A Table A3 for details) The standard pulse profile obtained from all 35 groups of observation data is shown in Figure 3.
It can be found that there is a certain systematic deviation between the standard profile obtained by DE200 and the profiles obtained by the other three ephemerides.The calculated deviation value is about 1.3 ms.The deviation value is approximately equal to the timing deviation of each ephemeris; the difference between the profiles obtained by DE405, 421 and 430 is small, and the photon number difference of discrete phase points is generally less than 100.Due to the small number of profile bins, the systematic deviation between the three cannot be displayed through the profiles.For the systematic deviation between the Universe 2022, 8, 360 10 of 18 DE200 profile and the other three, the possible reason for the deviation is that the DE200 is transformed from the B1950 coordinate system to the J2000 coordinate system through coordinate transformation according to the older ephemeris, while the newer ephemeris is compiled according to the latest observation data at that time and based on ICRF.The two methods are likely to produce some spatial deviations.
compares the effects of different ephemerides on pulse profile, timing residuals and single pulsar orbit determination and gives some suggestions on the use of ephemeris based on the results.

Pulse Profile
The pulse profile is obtained by the epoch folding method, and the timing parameters are the crab timing parameters regularly released by Jodrell Bank Observatory.(In order to minimize timing noise, it is necessary to use parameters that are as consistent as possible with the date of the observed data.During the observation period, 2016/11/17~12/19, Jodrell has the data of reference epoch 11/23 and 12/16, while the reference epoch of ATNF similar data is very old.See Appendix A Table A3 for details) The standard pulse profile obtained from all 35 groups of observation data is shown in Figure 3.It can be found that there is a certain systematic deviation between the standard profile obtained by DE200 and the profiles obtained by the other three ephemerides.The calculated deviation value is about 1.3 ms.The deviation value is approximately equal to the timing deviation of each ephemeris; the difference between the profiles obtained by DE405, 421 and 430 is small, and the photon number difference of discrete phase points is generally less than 100.Due to the small number of profile bins, the systematic deviation between the three cannot be displayed through the profiles.For the systematic deviation between the DE200 profile and the other three, the possible reason for the deviation is that the DE200 is transformed from the B1950 coordinate system to the J2000 coordinate system through coordinate transformation according to the older ephemeris, while the newer ephemeris is compiled according to the latest observation data at that time and based on ICRF.The two methods are likely to produce some spatial deviations.
If the systematic deviation of the DE200 profile relative to the other three is corrected, the resulting standard pulse profile is shown in Figure 4.The Pearson correlation coefficient between each profile is calculated according to the profiles after error correction, as shown in Table 3.If the systematic deviation of the DE200 profile relative to the other three is corrected, the resulting standard pulse profile is shown in Figure 4.The Pearson correlation coefficient between each profile is calculated according to the profiles after error correction, as shown in Table 3.It can be found that the difference between the profiles is small after correcting the system deviation, and the Pearson correlation coefficients between the profiles are higher than 0.99.The profiles of DE200 can still see obvious photon number differences at some discrete phase points.The profiles of the two latest ephemerides almost coincide, and the photon number difference at most discrete phase points is less than 10.It can be found that the difference between the profiles is small after correcting the system deviation, and the Pearson correlation coefficients between the profiles are higher than 0.99.The profiles of DE200 can still see obvious photon number differences at some discrete phase points.The profiles of the two latest ephemerides almost coincide, and the photon number difference at most discrete phase points is less than 10.

Timing Residuals
Each group of observation data was folded to obtain a pulse profile, and its phase starting point is the phase corresponding to the arrival time of the first photon of the data, and then the phase difference was obtained by cross-correlation operation with the standard pulse profile.Without considering the integer ambiguity, the obtained phase difference is the normalized phase delay, and the difference between the phase calculated by the standard timing parameters and multiplied by the period is the timing residual.
Here, the fast Fourier transform (FFT) was used to calculate the cross-correlation, and the timing residuals calculated based on each ephemeris are shown in Figure 5. Table 4 shows the corresponding RMSE and MAE values.It can be found that the difference between the profiles is small after correcting the system deviation, and the Pearson correlation coefficients between the profiles are higher than 0.99.The profiles of DE200 can still see obvious photon number differences at some discrete phase points.The profiles of the two latest ephemerides almost coincide, and the photon number difference at most discrete phase points is less than 10.

Timing Residuals
Each group of observation data was folded to obtain a pulse profile, and its phase starting point is the phase corresponding to the arrival time of the first photon of the data, and then the phase difference was obtained by cross-correlation operation with the standard pulse profile.Without considering the integer ambiguity, the obtained phase difference is the normalized phase delay, and the difference between the phase calculated by the standard timing parameters and multiplied by the period is the timing residual.Here, the fast Fourier transform (FFT) was used to calculate the cross-correlation, and the timing residuals calculated based on each ephemeris are shown in Figure 5. Table 4 shows the corresponding RMSE and MAE values.It was found that the timing residuals of each group calculated based on the four ephemerides are less than 150 microseconds, and the RMSE and MAE values are also very close; the difference is less than 2 µs.Similarly, the timing residuals obtained by DE200 are quite different from those obtained by the other three ephemerides (mainly expressed in the timing residuals at a single point); the difference between the timing residuals of each group of the results obtained from DE405, DE421 and DE430 is very small, and their MAE values of is slightly better than that of DE200, indicating that their internal accuracy is higher.The systematic trend of timing residuals due to different ephemerides is the same as that of time conversion and light-travel delay calculation.
In order to verify the timing performance when different ephemerides are mixed, the timing residuals were obtained by cross-correlation operation between the DE405 standard pulse profile and other groups of pulse profiles obtained from other ephemerides, and  It was found that the timing residuals of each group calculated based on the four ephemerides are less than 150 microseconds, and the RMSE and MAE values are also very close; the difference is less than 2 µs.Similarly, the timing residuals obtained by DE200 are quite different from those obtained by the other three ephemerides (mainly expressed in the timing residuals at a single point); the difference between the timing residuals of each group of the results obtained from DE405, DE421 and DE430 is very small, and their MAE values of is slightly better than that of DE200, indicating that their internal accuracy is higher.The systematic trend of timing residuals due to different ephemerides is the same as that of time conversion and light-travel delay calculation.
In order to verify the timing performance when different ephemerides are mixed, the timing residuals were obtained by cross-correlation operation between the DE405 standard pulse profile and other groups of pulse profiles obtained from other ephemerides, and the results are shown in Figure 6; Table 5 shows the corresponding RMSE and MAE values.
It was found that the timing residuals of each group calculated based on the four ephemerides are less than 150 microseconds, and the RMSE and MAE values are also very close; the difference is less than 2 µs.Similarly, the timing residuals obtained by DE200 are quite different from those obtained by the other three ephemerides (mainly expressed in the timing residuals at a single point); the difference between the timing residuals of each group of the results obtained from DE405, DE421 and DE430 is very small, and their MAE values of is slightly better than that of DE200, indicating that their internal accuracy is higher.The systematic trend of timing residuals due to different ephemerides is the same as that of time conversion and light-travel delay calculation.
In order to verify the timing performance when different ephemerides are mixed, the timing residuals were obtained by cross-correlation operation between the DE405 standard pulse profile and other groups of pulse profiles obtained from other ephemerides, and the results are shown in Figure 6; Table 5 shows the corresponding RMSE and MAE values.It was found that the timing residuals obtained by DE405/421 are completely consistent with those obtained by using the DE421 standard profile.Compared with the DE200 standard profile, the timing residuals of DE405/200 increased, the maximum value was more than 150 microseconds, and the values of RMSE and MAE increased by about 30%.The first timing residual of DE405/430 increased significantly, while other values did not change obviously.Affected by the first value, both RMSE and MAE values increased.It can be found that the internal self-consistent ephemeris can generally obtain better results than the mixed ephemerides, which need special attention in the application of pulsar navigation and timing.

Orbit Determination Analysis
Based on the analysis in Section 2.2, the SEPO was used to analyze the influence of different ephemerides on XPNAV-I orbit determination.In this experiment, the realtime photons arrival rate function at the spacecraft was calculated by using the pulsar signal model at the SSB and the spacecraft orbit information, combined with the pulse profile of each ephemeris, and then the photons arrival time series was generated by the scale transformation method [35].The grazing incidence Wolter-I X-ray telescope of XPNAV-I has an area of only 30 cm 2 , but it has a small optical field of view and strong noise suppression ability.For one month's sparse observation data, the amount of signal simulation calculation is large.
For each ephemeris, take the GPS coordinate values [X, Y, Z] recorded by the satellite in the J-2000 coordinate system as the nominal orbit, add errors on the coordinate values in 1 km steps to obtain [X + ∆L, Y + ∆L, Z + ∆L], and the ∆L range is [−10, 10] km.Calculate the χ 2 value under each coordinate and determine the range of the optimal orbit according to the size of the χ 2 value.The χ 2 values obtained from different ephemerides are shown in Figure 7 and Table 6 shows the specific values.
sar navigation and timing.

Orbit Determination Analysis
Based on the analysis in Section 2.2, the SEPO was used to analyze the influence of different ephemerides on XPNAV-I orbit determination.In this experiment, the real-time photons arrival rate function at the spacecraft was calculated by using the pulsar signal model at the SSB and the spacecraft orbit information, combined with the pulse profile of each ephemeris, and then the photons arrival time series was generated by the scale transformation method [35].The grazing incidence Wolter-I X-ray telescope of XPNAV-I has an area of only 30 cm 2 , but it has a small optical field of view and strong noise suppression ability.For one month's sparse observation data, the amount of signal simulation calculation is large.
For each ephemeris, take the GPS coordinate values [, , ] recorded by the satellite in the J-2000 coordinate system as the nominal orbit, add errors on the coordinate values in 1 km steps to obtain [ + ∆,  + ∆,  + ∆], and the ∆ range is [−10,10]km.Calculate the  value under each coordinate and determine the range of the optimal orbit according to the size of the  value.The  values obtained from different ephemerides are shown in Figure 7 and Table 6 shows the specific values.It was found that for the four ephemerides used, the χ 2 value of the simulation profiles can show a certain regularity rather than random distribution.The peak value represents the deviation between the simulated orbit and the actual orbit.The difference between the orbit obtained by DE200 and the actual orbit is less than 1 km, while the deviation of the other three ephemerides is between 1 and 2 km.Due to the limited experimental data, it is impossible to draw a conclusion that the orbit accuracy of DE200 is higher, which may be that the timing parameters used are the crab timing parameters of Jodrell Bank Observatory compiled based on DE200.Further experimental analysis is needed to draw a more accurate conclusion.It should be noted here that the "orbit determination" analyzed in this paper is to reverse the pulse profile based on the known nominal orbit, and it cannot be used for actual pulsar navigation or orbit determination, but for the research of this paper, it can still reflect the influence of different ephemerides on orbit determination accuracy to a certain extent.

Discussion
In the pulsar timing calculation, a relatively strict calculation process was adopted, and the results were improved compared with some research results using the same data [7,17].The timing residuals error of multiple groups of data is less than 50 microseconds, and the corresponding ranging error is less than 15 km.If the corresponding timing accuracy is converted into positioning accuracy, it has satisfied the navigation and positioning requirements of the transfer orbit segment of the deep space explorer.At the same time, how to form an effective timing measurement in a short time is the focus of future research.
Compared with the old ephemeris, the newer ephemeris usually has higher spatial accuracy.Most studies also show that the newly released ephemeris often has better timing performance [6].The above experimental results also show that when the systematic deviation is corrected, the MAE values of the timing residuals of the three newer ephemerides are still slightly better than that of DE200, indicating that their internal accuracy is higher.For some historical reasons, pulsar timing parameters obtained from older ephemeris such as DE200 are still widely used and of a certain value.In order to adapt to the development of research as soon as possible, it is necessary to use the newer ephemeris to obtain timing parameters and promote them.
In this paper, the orbit determination based on single pulsar observation data was used to analyze the navigation accuracy, which is an ideal processing method in the case of a lack of observation data.It should be noted that the "orbit determination" analyzed here is to reverse the pulse profile based on the known nominal orbit, and it cannot be used for actual pulsar navigation or orbit determination, but for the research of this paper, it can still reflect the influence of different ephemerides on orbit determination accuracy to a certain extent.Later, the same observation data can be used for rigorous dynamic orbit determination analysis.
Compared with the modeling of timing and positioning error caused by ephemeris error in some papers, this paper only calculates the difference of different ephemerides in timing and navigation solutions [22][23][24].However, through analysis, it was found that ephemeris error belongs to red noises that are difficult to model, and the absolute errors of celestial coordinates or velocity in different ephemerides are also difficult to determine, so it may be difficult to deal with its impact in ephemeris error modeling.According to the appeal analysis, on the one hand, the whole process can be solved based on one newly released ephemeris to reduce the impact of the internal inconsistency of ephemeris.On the other hand, with the development of complex system analysis and multi-body motion theory, it can also be introduced into the study of a pulsar navigation system to establish a more accurate model.
Due to the limited observation data, there is no experimental comparison for most navigation or timing situations.Subsequent analysis can be based on ground radio observation data and long-term observation data from other satellites to obtain more accurate conclusions.

Conclusions
The effects of JPL DE ephemeris on X-ray pulsar navigation and timing were studied and experimentally analyzed.Firstly, the current theoretical models and data processing methods of X-ray pulsar timing and navigation were summarized and analyzed, and the possible differences of ephemeris in the solution of timing and navigation equations were emphatically analyzed.Based on the measured data released by XPNAV-I, the effects of four commonly used JPL DE ephemerides on timing and orbit determination results were compared, and the following conclusions were summarized: (1) For time scale conversion, the difference is less than 1 ns during the data recording period, and the long-term difference is also less than 5 ns.The ephemeris has a significant impact on the calculation of light-travel delay.For Romer delay, the difference value can be greater than 1 ms, but for Shapiro delay and other relativistic effects caused by gravity, the delay difference value is small.At the same time, with the update and standardization of the ephemeris, the difference value was very small; (2) For the final timing results, the profile obtained by DE200 has a certain systematic deviation compared with other ephemerides, which is approximately equal to the delay difference, and the difference between RMSE and MAE values of timing residuals calculated by each ephemeris is less than 2 µs.However, the cross-correlation operation of the pulse profile obtained from different ephemerides increases the timing residuals.Therefore, when using ephemeris in navigation and timing applications, we must pay attention to its internal self-consistency; (3) When the four ephemerides are used for orbit determination, the orbit differences value is less than 1 km.Since the timing parameters are calculated based on DE200, its navigation or orbit determination accuracy is more accurate than the other three ephemerides, which also shows the importance of adopting a unified ephemeris in all links.At the same time, considering the higher spatial accuracy of the new ephemeris, it is the future development trend to use it to compile and popularize the pulsar timing model as soon as possible.

DE421 8 :Figure 1 .
Figure 1.Comparison of Einstein delay during data observation.To clearly reflect the small difference in Einstein delays during observation, subtract DE200 Einstein delay from DE405, 421 and 430 Einstein delays, respectively.The orange curve "DE405−DE200" in the figure represents ∆ − ∆ ; the blue curve represents ∆ − ∆ ; the green curve represents ∆ − ∆ .

Figure 1 .
Figure 1.Comparison of Einstein delay during data observation.To clearly reflect the small difference in Einstein delays during observation, subtract DE200 Einstein delay from DE405, 421 and 430 Einstein delays, respectively.The orange curve "DE405-DE200" in the figure represents ∆E DE405 − ∆E DE200 ; the blue curve represents ∆E DE421 − ∆E DE200 ; the green curve represents ∆E DE430 − ∆E DE200 .

Table 2 . 19 Figure 2 .
Figure 2. Comparison of time delay during data observation.The meaning of each curve is the same as that in Figure 1, and "DE405-DE200" represents the difference between the light-travel delay calculated by DE405 and DE200, and the rest are analogous.By comparing the results, it was found that Roemer delay is the main component of light-travel delay and the delay with the largest difference in calculation results between different ephemerides.The difference between the Roemer delay calculated based on DE200 and the values of the other three ephemerides is greater than 1 ms.The difference between Roemer delay calculated by DE405, DE421 and DE430 is less than 10 µs.The difference between the two latest ephemeris Roemer delay is 0.1 microsecond.The ephemeris has little effect on parallax delay and Shapiro delay, and the difference between the

Figure 2 .
Figure 2. Comparison of time delay during data observation.The meaning of each curve is the same as that in Figure 1, and "DE405-DE200" represents the difference between the light-travel delay calculated by DE405 and DE200, and the rest are analogous.
(a) Integral period profiles (b) Partial period profiles.

Figure 3 .
Figure 3.Comparison of standard pulse profiles.The number of profile bins is 256.The profiles of DE405, DE421 and DE430 in (a) basically coincide.(b) shows the profile difference during the period of phase value 0.5-0.6 (it can more clearly reflect the photon number differences of the profiles in this interval).

Figure 3 .
Figure 3.Comparison of standard pulse profiles.The number of profile bins is 256.The profiles of DE405, DE421 and DE430 in (a) basically coincide.(b) shows the profile difference during the period of phase value 0.5-0.6 (it can more clearly reflect the photon number differences of the profiles in this interval).

Figure 4 .
Figure 4. Comparison of standard pulse profiles after correcting system deviation.(a) Integral period profiles; (b) Partial period profiles.

Figure 4 .Table 3 .
Figure 4. Comparison of standard pulse profiles after correcting system deviation.(a) Integral period profiles; (b) Partial period profiles.
period profiles; (b) Partial period profiles.

Figure 6 .
Figure 6.Timing residual using DE405 standard pulse profile.In (a), the red curve represents the timing residuals of the original DE200, the purple curve represents the timing residual mixed with the DE200 and standard pulse profile of DE405, and the rest are analogous.(a) DE200; (b) DE421; (c) DE430.

Table 3 .
Pearson correlation coefficient between pulse profiles.P_200_405 represents the Pearson correlation coefficient between DE200 profile and DE405 profile, and the rest are analogous.

Table 3 .
Pearson correlation coefficient between pulse profiles.P_200_405 represents the Pearson correlation coefficient between DE200 profile and DE405 profile, and the rest are analogous.