Orbit Propagation with Lie Transfer Maps in the Perturbed Kepler Problem

The Lie transfer map method may be applied to orbit propagation problems in celestial mechanics. This method, described in another paper, is a perturbation method applicable to Hamiltonian systems. In this paper, it is used to calculate orbits for zonal perturbations to the Kepler (two-body) problem, in both expansion in the eccentricity and closed form. In contrast with a normal form method like that of Deprit, the Lie transformations here are used to effect a propagation of phase space in time, and not to transform one Hamiltonian into another.     

Research career of an astronomer who has studied celestial mechanics

Celestial mechanics has been a classical field of astronomy. Only a few astronomers were in this field and not so many papers on this subject had been published during the first half of the 20 thcentury.However, as the beauty of classical dynamics and celestial mechanics attracted me very much, I decided to take celestial mechanics as my research subject and entered university, where a very famous professor of celestial mechanics was a member of the faculty. Then as artificial satellites were launched starting from October 1958, new topics were investigated in the field of celestial mechanics. Moreover, planetary rings,asteroids with moderate values of eccentricity, inclination and so on have become new fields of celestial mechanics. In fact I have tried to solve such problems in an analytical way. Finally, to understand what gravitation is I joined the TAMA300 gravitational wave detector group.     

Limits of predictability of gravitational systems

The concept of finite predictability of gravitational many-body systems is related to the non-deterministic nature of celestial mechanics and of dynamics, in general. The basic, fundamental reasons for the uncertainty of predictions are as follows: (1) the initial conditions are known only approximately since they are obtained either from observations or from approximate computations; (2) the equations of motion given by a selected model describe the actual system only approximately; (3) the physical constants of the dynamical system have error limits; (4) the differential equations of motion are non-integrable and numerical integration methods must be used for solution, generating errors in the final result at every integration step.In addition to these reasons, mostly depending on our techniques, there are some more fundamental reasons depending on the nature of the dynamical system investigated. These are the appearance of regions of instability, non-integrability and chaotic motion.Details, effects and controls of these regions for finite predictability are discussed for various dynamical systems of importance in celestial mechanics with special emphasis on planetary systems.     

Specialized celestial mechanics systems for symbolic manipulation

Construction and application of the current high accuracy analytical theories of motion of celestial bodies necessitates the development of specialized software for the implementation of analytical algorithms of celestial mechanics. This paper describes a typical software package of this kind. This package includes a universal Poisson processor for the rational functions of many variables, a tensorial processor for purposes of relativistic celestial mechanics, a Keplerian processor valid for the solutions of the two body problem in the form of a Poisson series, Taylor expansions in powers of time and closed expressions, and an analytical generator of celestial mechanics functions, facilitating the immediate implementation of the present analytical methods of celestial mechanics. The package is completed with a numerical-analytical interface designed, in particular, for the fast evaluation of the long Poisson series.     

Long-Term Stability Analysis of Quasi Integrable Degenerate Systems Through the Spectral Formulation of the Nekhoroshev Theorem

We describe a numerical application of the Nekhoroshev theorem to investigate the long-term stability of quasi-integrable systems. We extend the results of a previous paper to a class of degenerate systems, which are typical in celestial mechanics.     

Celestial mechanics: Past, present, future

Over all steps of its development celestial mechanics has played a key role in solar system researches and verification of the physical theories of gravitation, space and time. This is particularly characteristic for celestial mechanics of the second half of the 20th century with its various physical applications and sophisticated mathematical techniques. This paper is attempted to analyze, in a simple form (without mathematical formulas), the celestial mechanics problems already solved, the problems that can be and should be solved more completely, and the problems still waiting to be solved.     

First Post-Newtonian Angular Momentum In Relativistic Celestial Mechanics

With a new theory on the 1PN celestial mechanics recently developed by Damour, Soffel and Xu (1991,1992,1993,1994), definitions and expressions of the 1PN spin angular momentum are investigated and analysed. The total spin angular momentum of a system of extended bodies such as the solar system is calculated and expressed as the function of local parameters and observables under reasonable assumptions, which would find its application in the evolution and dynamics of systems of celestial bodies. This revised version was published online in July 2006 with corrections to the Cover Date.     

Computing with certainty individual members of families of periodic orbits of a given period

The accurate computation of families of periodic orbits is very important in the analysis of various celestial mechanics systems. The main difficulty for the computation of a family of periodic orbits of a given period is the determination within a given region of an individual member of this family which corresponds to a periodic orbit. To compute with certainty accurate individual members of a specific family we apply an efficient method using the Poincaré map on a surface of section of the considered problem. This method converges rapidly, within relatively large regions of the initial conditions. It is also independent of the local dynamics near periodic orbits which is especially useful in the case of conservative dynamical systems that possess many periodic orbits, often of the same period, close to each other in phase space. The only computable information required by this method is the signs of various function evaluations carried out during the integration of the equations of motion. This method can be applied to any system of celestial mechanics. In this contribution we apply it to the photogravitational problem.     

The Long-term Error Estimation Method for the Numerical Integrations of Celestial Orbits

Numerical methods have become a very important type of tool for celestial mechanics, especially in the study of planetary ephemerides. The errors generated during the computation are hard to know beforehand when applying a certain numerical integrator to solve a certain orbit. In that case, it is not easy to design a certain integrator for a certain celestial case when the requirement of accuracy were extremely high or the time-span of the integration were extremely large. Especially when a fixed-step method is applied, the caution and effort it takes would always be tremendous in finding a suitable time-step, because it is about whether the accuracy and time-cost of the final result are acceptable. Thus, finding the best balance between efficiency and accuracy with the least time cost appeared to be a major obstruction in the face of both numerical integrator designers and their users. To solve this problem, we investigate the variation pattern of truncation error and the pattern of rounding error distributions with time-step and time-span of the integration. According to those patterns, we promote an error estimation method that could predict the distribution of rounding errors and the total truncation errors with any time-step at any time-spot with little experimental cost, and test it with the Adams-Cowell method in the calculation of circular periodic orbits. This error estimation method is expected to be applied to the comparison of the performance of different numerical integrators, and also it can be of great help for finding the best solution to certain cases of complex celestial orbits calculations.     

天体轨道长期数值积分的误差估计方法

数值积分方法是进行天体力学研究的重要工具, 尤其对于行星历表的研究工作而言. 由于在使用数值方法计算天体轨道时, 最终误差通常是难以预知的, 所以在面对精度要求较高或者积分时间较长的工作时具体积分方案的设计---尤其是当使用定步长方法时的步长选择---需要十分谨慎, 因为这将意味着是否能在时间成本可以被接受的范围内使解的精度达到要求. 因此, 在使用数值方法解决实际问题时如何快速寻找效率与精度之间的最佳平衡点是每一个数值积分方法的设计者与使用者都会面临的难题. 为解决这一问题, 在定步长条件下对数值积分方法的舍入误差概率分布函数以及截断误差积累量对步长的依赖关系和随时间的增长关系进行了深入研究. 基于所得结论, 提出了一种仅需较少的数值实验资料即可对选择任意时间步长积分至任意积分时刻时的舍入误差概率分布函数与截断误差积累量进行准确估计的方法, 并使用Adams-Cowell方法对该误差估计方法在圆周期轨道条件下进行了验证. 该误差估计方法在未来有望用于不同数值算法的性能对比研究, 同时也可以对数值积分方法求解实际轨道问题时的决策工作带来重要帮助.     

相对论天体力学和后牛顿运动方程

叙述了与Astrod工程有关的相对论天体力学基础内容。包括相对论天体力学、广义相对论基本原理、PPN方法体系、PPN多体问题、PPN二体问题。高阶PN二体问题等     

Conditionally periodic solutions of the canonical systems of differential equations in non-autonomous resonant case

A formal method of constructing of conditionally periodic solutions of canonical systems of differential equations in the vicinity of a commemsurability of frequencies is proposed. The method is a union of the rapid convergence method and (well-known in celestial mechanics) Delaunay-Zeipel's method of canonical transformations. For a successful application of the method an existence of stationary resonant solutions of an averaged system of the differential equations is necessary.     

Onset of secular chaos in planetary systems: period doubling and strange attractors

As a result of resonance overlap, planetary systems can exhibit chaotic motion. Planetary chaos has been studied extensively in the Hamiltonian framework, however, the presence of chaotic motion in systems where dissipative effects are important, has not been thoroughly investigated. Here, we study the onset of stochastic motion in presence of dissipation, in the context of classical perturbation theory, and show that planetary systems approach chaos via a period-doubling route as dissipation is gradually reduced. Furthermore, we demonstrate that chaotic strange attractors can exist in mildly damped systems. The results presented here are of interest for understanding the early dynamical evolution of chaotic planetary systems, as they may have transitioned to chaos from a quasi-periodic state, dominated by dissipative interactions with the birth nebula.     

Time scaling as an infinitesimal canonical transformation

We show that time scaling transformations for Hamiltonian systems are infinitesimal canonical transformations in a suitable extended phase space constructed from geometrical considerations. We compute its infinitesimal generating function in some examples: regularization and blow up in celestial mechanics, classical mechanical systems with homogeneous potentials and Scheifele theory of satellite motion.Research partially supported by CONACYT (México), Grant PCCBBNA 022553 and CICYT (Spain).     

Virial oscillations of celestial bodies

A general approach to the solution of the perturbed oscillation problem for celestial bodies is considered. The solution sought describes unperturbed virial oscillations (zero approximation) affected by external perturbing effects. In the general case, these perturbations can be expressed by an arbitrary given function of time, Jacobi's function and its first derivative. Standard methods and modes of perturbation theory are used for solution of the problem.It is shown that while studying the evolution of a celestial body as a dissipative system in the framework of perturbed virial oscillations, the analytical expression for perturbing function can be derived, assuming the celestial body to be an oscillating electrical dipole emitting electromagnetic energy.The general covariant form of Jacobi's equation is derived and its spur is examined. It is shown that the scalar form of Jacobi's equation appears to be more universal than Newton's laws of motion from which it is derived.     

Semi-analytical integration algorithms based on the use of several kinds of anomalies as temporal variable

One of the main problems in celestial mechanics is the management of long developments in Fourier or Poisson series used to describe the perturbed motion in the planetary system.In this work we shall develop a software package that is suitable for managing these objects. This package includes algorithms to obtain the inverse of the distance based on an iterative method, a set of integration algorithms according to several sets of temporal variables.This paper contains a comparative study on the use of the true, eccentric, and elliptic anomalies in semi-analytical methods on celestial mechanics.     

Analysis of interval constants in calendars affiliated with the Shoushili

We study interval constants that are related to motions of the Sun and Moon,i.e., the Qi, Intercalation, Revolution and Crossing interval, in calendars affiliated with the Shoushi calendar(Shoushili), such as Datongli and Chiljeongsannaepyeon. It is known that these interval constants were newly introduced in the Shoushili calendar and revised afterward, except for the Qi interval constant, and the revised values were adopted in later calendars affiliated with the Shoushili. We first investigate the accuracy of these interval constants and then the accuracy of calendars affiliated with the Shoushili in terms of these constants by comparing times for the new moon and the maximum solar eclipse calculated by each calendar with modern methods of calculation. During our study, we found that the Qi and Intercalation interval constants used in the early Shoushili were well determined, whereas the Revolution and Crossing interval constants were relatively poorly measured. We also found that the interval constants used by the early Shoushili were better than those of the later one, and hence better than those of Datongli and Chiljeongsannaepyeon. On the other hand, we found that the early Shoushili is, in general, a worse calendar than Datongli for use in China but a better one than Chiljeongsannaepyeon for use in Korea in terms of times for the new moon and when a solar eclipse occurs, at least for the period 1281 – 1644.Finally, we verified that the times for sunrise and sunset in the Shoushili-Li-Cheng and Mingshi are those at Beijing and Nanjing, respectively.     

双星轨道拟合的研究进展

双星轨道拟合是天文学的一项基础性研究工作。其主要目的是给出双星系统的二体轨道参数,这些参数不仅是高精度、高网格密度星表参考架的必要组成部分,而且也为理解各种有关观测现象提供了必要的动力学基础;更重要的是,双星轨道拟合可以直接估计恒星物理和星系天文学等领域极有应用价值的恒星质量参数。因此,长期以来双星轨道拟合工作一直受到研究者的广泛关注。近年来,随着高精度的恒星运动学观测资料的大量积累,双星轨道拟合更成为天体测量和天体力学的一个共同的热点课题,有关研究也有了长足的进展。综述了双星轨道拟合的历史及现状,其中着重介绍了目前所用的主要观测资料和各种具体的拟合模型、拟合方法;简要描述了几种主要的双星星表;展望了今后双星轨道拟合工作的发展趋势。     

A Geometric Interpretation of Integrable Motions

Integrability, one of the classic issues in galactic dynamics and in general in celestial mechanics, is here revisited in a Riemannian geometric framework, where Newtonian motions are seen as geodesics of suitable -mechanical- manifolds. The existence of constants of motion that entail integrability is associated with the existence of Killing tensor fields on the mechanical manifolds. Such tensor fields correspond to hidden symmetries of non-Noetherian kind. Explicit expressions for Killing tensor fields are given for the N = 2 Toda model, and for a modified Hénon-Heiles model, recovering the already known analytic expressions of the second conserved quantity besides energy for each model respectively.     

Minimum energy configurations in the N-body problem and the celestial mechanics of granular systems

Minimum energy configurations in Celestial Mechanics are investigated. It is shown that this is not a well defined problem for point-mass celestial mechanics but well-posed for finite density distributions. This naturally leads to a granular mechanics extension of usual Celestial Mechanics questions such as relative equilibria and stability. This paper specifically studies and finds all relative equilibria and minimum energy configurations for N?=?1, 2, 3 and develops hypotheses on the relative equilibria and minimum energy configurations for N ? 1 bodies.