X射线脉冲星自主导航的光传播时间方程

详细讨论了X射线脉冲星自主导航系统的广义相对论效应,采用DSX体系的后牛顿近似方法,计算出1PN度规下的光线弯曲轨迹和时间延缓以及2PN度规下的引力延缓,得到X射线脉冲自主导航的高精度测量方程.     

二阶后牛顿光线轨迹方程

根据DSX体系的后牛顿近似理论,直接由Lagrange方程导出了轴对称稳态时空中光子的二阶后牛顿轨迹方程,并求得在赤道平面内传播的光线偏转角.在可测量精度范围内,得到的结论与Schwarzschild和Kerr度规下的情况相吻合.     

相相对论性天文参考系研究进展

国际天文联合会(IAU)在2000年的决议上以两个互相等价的相对论N体多参考系理论——Brumberg-Kopeikin体系和Damour-Soffel-Xu体系为基础,构造了严格且自洽的一阶后牛顿(1PN)的局部参考系和全局参考系,并给出了相应的坐标变换规则。回顾了IAU2000决议关于参考系理论的核心内容,并指出该参考系理论的主要优点和理论不足。结合决议发表前后国际上对相对论参考系理论的一系列扩展研究,详细总结了在二阶后牛顿推广、参数后牛顿化以及太阳系非孤立引力系统等方面对参考系理论的研究进展,并讨论了未来的理论研究发展以及对实际天体测量的影响。     

太阳系天体地面VLBI观测的相对论时延模型

以太阳系质心参考系为基础，根据太阳系的质心参考系和非旋转地球质心参考系的坐标转换关系，推导了太阳系天体地面VLBI观测的相对论时间延迟模型，给出了一个通用的解析表达式．根据这一公式可以得到平劲松博士所采用的公式，以及当地心与源的距离无限大时，可得河外射电源VLBI观测的Zhu—Groten模型、Shapiro模型和IERS(92，96)推荐模型．所推导的公式严格解析且无误差，在实际应用中建议采用这一公式．同时详细地讨论了所推导公式的实用范围和各种舍掉项的量级估计，并详细给出了时间延迟理论模型的计算步骤．     

粒子在标量张量理论的二阶后闵可夫斯基近似弱引力场中的引力偏转

标量张量理论是目前为止最成功的相对论性引力理论之一,经受住了已有实验和观测的检验.随着实验和观测精度的提高,理论将得到进一步的检验.为此将该理论的可观测效应计算到二阶后牛顿近似非常必要.标量张量理论的二阶后牛顿近似结果已经由Xie等人给出.为得到有质量粒子在弱引力场中的偏转情况,利用二阶后闵可夫斯基近似下的度规解得到了粒子的测地线方程,并通过迭代的方法求得在非束缚条件下,粒子在弱引力场中的轨迹的解析解.进一步利用此结果计算了粒子在弱引力场中的偏转角(与光线的偏转角不同)并与前人的结果进行了比较.     

地球四极矩定义及参考系转换的探讨

文[1]对地球四极矩在太阳系质心参考系（BRS）中影响人造卫星运动的相对论效应作了仔细的研究，引入相对经模型轨道的厘米级上符合。本文发现，利用文[2]得到的结果与此有矛盾。经过研究，我们认为其原因在于文[2]定义的均匀极矩及其转换上存在缺陷，在实际应用中应考虑其局了性。本文还将文[2]的方法与相对论参考系的最新结果DSX理论进行了比较。     

关于X射线脉冲星自主导航的位置观测方程

在X射线脉冲星自主导航(XNAV)中,位置观测方程表达了X光子到达航天器的时刻(TOA)和航天器位置的关系.具体讨论时,一般用TOA和"时间基准"的差值代替TOA,用太阳系质心系中的位置矢量表示航天器的位置."时间基准"可以取X光子到达太阳系质心(SSB)的真实时间,也可以取X光子到达SSB的"等效时间".讨论了基于这两种时间基准的位置观测方程,给出了时间精度为0.1 ns的位置观测方程,分析了其中各项的物理意义.     

相对论天文参考系的回顾与展望

国际天文联合会(IAU)在1991年和2000年分别通过了关于相对论天文参考系的重要决议。特别是2000年的决议,以两个互相等价的相对论N体多参考系理论:Brumberg-Kopeikin体系和Damour-Soffel-Xu体系为基础,构造了严格且自洽的一阶后牛顿(1PN)的局部参考系和全局参考系,并给出了相应的坐标变换规则。IAU2000决议发表后的10多年里,已经开始应用在一些高精度的天体测量数据处理模型中,但是其工程化应用还没有得到广泛实现,特别是目前为止我国没有将其应用到具体的天文观测或者空间探测计划中。因此有必要对IAU的相对论天文参考系理论做一个系统的解读。首先介绍了IAU在1991年给出的一个简单的相对论参考系决议;随后在第3章和第4章详细讨论了Brumberg-Kopeikin理论和Damour-Soffel-Xu体系;接着,详细给出了IAU2000年大会关于相对论参考系的决议内容;最后,讨论了IAU2000决议的工程化应用、近10多年来的理论发展以及对未来的展望。     

考虑后牛顿效应的二体问题的解和星历计算

本文给出考虑后牛顿(PN)效应的二体问题解所对应的基本关系式,并仿照开普勒(Kepler)运动,给出星历表计算方法和相应的计算公式以及适用于数值研究中的简单形式。     

双星系中两子星的自转对双星轨道变化的后牛顿效应

本文利用解摄动方程的平均值法求得在PPN框架中二体自转对轨道要素产生的后牛顿效应的长期变化影响．利用这一理论对CWCep和DRVul两颗双星中两子星的自转对轨道近星点和平近点角的长期摄动的后牛顿效应做了计算。结果表明：对于两个质量较大快速自转的子星,由此所产生的后牛顿效应的摄动量是不能忽视的。     

The Time Equation of Light Propagation in XNAV

We discuss in detail the general relativistic effect in the X-ray sourcebased navigation for autonomous position determination program (XNAV). By using the post-Newtonian approximate method of the DSX scheme, we calculate the bending of light and the gravitational time delay under the 1PN metric, as well as the gravitational time delay under the 2PN metric, and finally obtain the high-accuracy time equation of light propagation in XNAV.     

On derivation of EIH (Einstein–Infeld–Hoffman) equations of motion from the linearized metric of general relativity theory

The half-century old idea of Infeld to use the variational principle of the general relativity field equations is reminded to show that the commonly employed EIH (Einstein–Infeld–Hoffman) equations of motion may be derived from the linearized (weak-field) metric alone. Based on it, the linearized metric might be sufficient for post-Newtonian celestial mechanics and astrometry enabling one to derive the post-Newtonian equations of motion and rotation of celestial bodies as well as the post-Newtonian equations of light propagation within the general relativity framework.     

On the energy integral for first post-Newtonian approximation

The post-Newtonian approximation for general relativity is widely adopted by the geodesy and astronomy communities. It has been successfully exploited for the inclusion of relativistic effects in practically all geodetic applications and techniques such as satellite/lunar laser ranging and very long baseline interferometry. Presently, the levels of accuracy required in geodetic techniques require that reference frames, planetary and satellite orbits and signal propagation be treated within the post-Newtonian regime. For arbitrary scalar W and vector gravitational potentials \(W^j (j=1,2,3)\), we present a novel derivation of the energy associated with a test particle in the post-Newtonian regime. The integral so obtained appears not to have been given previously in the literature and is deduced through algebraic manipulation on seeking a Jacobi-like integral associated with the standard post-Newtonian equations of motion. The new integral is independently verified through a variational formulation using the post-Newtonian metric components and is subsequently verified by numerical integration of the post-Newtonian equations of motion.     

Post-Newtonian satellite orbits

The first post-Newtonian approximation of general relativity is used to account for the motion of solar system bodies and near-Earth objects which are slow moving and produce weak gravitational fields. The \(n\)-body relativistic equations of motion are given by the Einstein-Infeld-Hoffmann equations. For \(n=2\), we investigate the associated dynamics of two-body systems in the first post-Newtonian approximation. By direct integration of the associated planar equations of motion, we deduce a new expression that characterises the orbit of test particles in the first post-Newtonian regime generalising the well-known Binet equation for Newtonian mechanics. The expression so obtained does not appear to have been given in the literature and is consistent with classical orbiting theory in the Newtonian limit. Further, the accuracy of the post-Newtonian Binet equation is numerically verified by comparing secular variations of known expression with the full general relativistic orbit equation.     

Relativistic effects and solar oblateness from radar observations of planets and spacecraft

We used more than 250 000 high-precision American and Russian radar observations of the inner planets and spacecraft obtained in the period 1961–2003 to test the relativistic parameters and to estimate the solar oblateness. Our analysis of the observations was based on the EPM ephemerides of the Institute of Applied Astronomy, Russian Academy of Sciences, constructed by the simultaneous numerical integration of the equations of motion for the nine major planets, the Sun, and the Moon in the post-Newtonian approximation. The gravitational noise introduced by asteroids into the orbits of the inner planets was reduced significantly by including 301 large asteroids and the perturbations from the massive ring of small asteroids in the simultaneous integration of the equations of motion. Since the post-Newtonian parameters and the solar oblateness produce various secular and periodic effects in the orbital elements of all planets, these were estimated from the simultaneous solution: the post-Newtonian parameters are β = 1.0000 ± 0.0001 and γ = 0.9999 ± 0.0002, the gravitational quadrupole moment of the Sun is J₂ = (1.9 ± 0.3) × 10^?7, and the variation of the gravitational constant is ?/G = (?2 ± 5) × 10^?14 yr^?1. The results obtained show a remarkable correspondence of the planetary motions and the propagation of light to General Relativity and narrow significantly the range of possible values for alternative theories of gravitation.     

R-N场中雷达回波延迟实验的后后牛顿修正

从球对称引力场中光子在赤道面上的运动微分方程出发,采用后牛顿近似方法,讨论了Reissner-Nordstr(o|¨)m度规场中雷达回波延迟实验的后后牛顿修正.通过计算给出了后牛顿以及后后牛顿的修正结果,当荷电量Q为零时,该结果适用于电中性天体,其后牛顿部分的修正与经典文献一致.文中采用的近似展开方法同样适用于讨论更高阶的修正.雷达回波的高阶修正对目前展开的高精度空间引力试验有着重要的参考价值.     

The General Method for Fixing the Gauges of Relativistic Astronomical Reference Systems

This article applies a new scheme of the first post-Newtonian theory (Damour et al., 1991–1994) to the problem of gauge in relativistic reference systems. Choosing and fixing gauge are necessary when the precision of time measurement and application needs to reach the 2PN level (10⁻¹⁶ or better). We present a general method for fixing the gauges of both the global and local coordinate systems, and for determining the expressions of gravitational potentials and coordinate transformations. The results relevant are consistent with the newest IAU resolutions, therefore they can be applied to astronomical practice.     

The linear stability of the post-Newtonian triangular equilibrium in the three-body problem

Continuing a work initiated in an earlier publication (Yamada et al. in Phys Rev D 91:124016, 2015), we reexamine the linear stability of the triangular solution in the relativistic three-body problem for general masses by the standard linear algebraic analysis. In this paper, we start with the Einstein–Infeld–Hoffmann form of equations of motion for N-body systems in the uniformly rotating frame. As an extension of the previous work, we consider general perturbations to the equilibrium, i.e., we take account of perturbations orthogonal to the orbital plane, as well as perturbations lying on it. It is found that the orthogonal perturbations depend on each other by the first post-Newtonian (1PN) three-body interactions, though these are independent of the lying ones likewise the Newtonian case. We also show that the orthogonal perturbations do not affect the condition of stability. This is because these do not grow with time, but always precess with two frequency modes, namely, the same with the orbital frequency and the slightly different one due to the 1PN effect. The condition of stability, which is identical to that obtained by the previous work (Yamada et al. 2015) and is valid for the general perturbations, is obtained from the lying perturbations.     

Coordinate gauge conditions in theN-body equations of motion in the first-order post-Newtonian approximation of general relativity

The explicit forms of the metric as well as the equations of motion in the first-order post-Newtonian approximation are worked out under several gauge conditions. It is noted that the so-called EIH (Einstein, Infeld, and Hoffman) equation of motion for an assembly ofN finite mass points mutually interacting via gravitation is identically obtained under three different gauge conditions, namely the harmonic gauge, Chandrasekhar gauge and a composite Chandrasekhar gauge used by Misneret al. (1970), even though the solutions for the metric are found to be all different. In one case the metric has a component apparently diverging, but finally generates regular affine connections so that the equations of motions become free from any singularity. By use of the Chandrasekhar gauge and his formulation, the second-order contribution to the acceleration of planets in the limit of test particle motion around the Sun has been calculated, the inclusion of which in the EIH set of the equations of motion would extend the relative accuracy of computing the total acceleration of any planet to better than one part in 10¹⁷.     

The 2nd-order Post-Newtonian Orbit Equation of Light

Based on the 2nd-order post-Newtonian approximation under the DSX frame of the general relativity theory, the 2nd-order post-Newtonian orbital equation of light in the axis-symmetrical stationary spacetime is derived, and from this, the angle of deflection of light propagating in the equatorial plane is derived. The obtained results are consistent with those of the Schwarzchild and Kerr metrics within the limits of measuring precision.