Comparison between Deprit and Dragt-Finn perturbation methods

In this paper, the relationship between the Dragt-Finn transform and the classical Lie transform introduced by Deprit is discussed. The relative performance of the algorithms used for the computations of the transformed functions is compared, and the relation between their generators is given. These generators produce the same transform which insures the construction of the same invariants.     

Attitude dynamics of a rigid body on a Keplerian orbit: A simplification

An infinitestimal contact transformation is proposed to simplify at first order the Hamiltonian representing the attitude of a triaxial rigid body on a Keplerian orbit around a mass point. The simplified problem reduces to the Euler-Poinsot model, but with moments of inertia depending on time through the longitude in orbit. Should the orbit be circular, the moments of inertia would be constant.     

Keplerian expansions in terms of Henrard's practical variables

Formulae for the Keplerian expansions in terms of Henrard's practical variables are given. Two different methods were applied: one using the Bessel functions and one based on the Lie transforms. The former involves less series products, but the latter is more flexible and universal.     

An Addendum to the Damour and Deruelle's Solution of the Post-Newtonian Two-Body Problem

A complement to the the analytical solution given by Damour and Deruelle (D&D) in order to solve the differential equations that describe the conic-like and rectilinear post-Newtonian two-body problem is presented. We show the geometric relation between the eccentric and true anomalies for the elliptic-like and hyperbolic-like cases and expressions for their velocities are obtained. A post-Newtonian version of the Lambert's theorem is shown. This revised version was published online in July 2006 with corrections to the Cover Date.     

Kepler Equation solver

Kepler's Equation is solved over the entire range of elliptic motion by a fifth-order refinement of the solution of a cubic equation. This method is not iterative, and requires only four transcendental function evaluations: a square root, a cube root, and two trigonometric functions. The maximum relative error of the algorithm is less than one part in 10¹⁸, exceeding the capability of double-precision computer arithmetic. Roundoff errors in double-precision implementation of the algorithm are addressed, and procedures to avoid them are developed.     

Keplerian periodic orbits in the isosceles problem

The planar isosceles three-body problem where the two symmetric bodies have small masses is considered as a perturbation of the Kepler problem. We prove that the circular orbits can be continued to saddle orbits of the Isosceles problem. This continuation is not possible in the elliptic case. Their perturbed orbits tend to a continued circular one or approach a triple collision. The basic tool used is the study of the Poincaré maps associated with the periodic solutions. This revised version was published online in July 2006 with corrections to the Cover Date.     

Practical Symplectic Methods with Time Transformation for the Few-Body Problem

The use of the extended phase space and time transformations for constructing efficient symplectic algorithms for the investigation of long term behavior of hierarchical few-body systems is discussed. Numerical experiments suggest that the time-transformed generalized leap-frog, combined with symplectic correctors, is one of the most efficient methods for such studies. Applications extend from perturbed two-body motion to hierarchical many-body systems with large eccentricities. This revised version was published online in July 2006 with corrections to the Cover Date.     

The restricted two-body problem in constant curvature spaces

We perform the bifurcation analysis of the Kepler problem on and . An analog of the Delaunay variables is introduced. We investigate the motion of a point mass in the field of a Newtonian center moving along a geodesic on and (the restricted two-body problem). For the case of a small curvature, the pericenter shift is computed using the perturbation theory. We also present the results of numerical analysis based on an analogy with the motion of a rigid body.     

Relativistic effects in the motion of artificial satellites: The oblateness of the central body II

The motion of artificial satellites in the gravitational field of an oblate body is discussed in the post — Newtonian framework using the technique of canonical Lie transformations. Two Lie transformations are used to derive explicit results for the longperiodic and secular perturbations for satellite orbits in the Einstein case.     

Perturbation theory for asteroid motion in the gravitational field of a migrating planet

A new perturbation method for the determination of proper elements of an asteroid in the gravitational field of a migrating planet is developed. The article is published in the original.     

Position and velocity perturbations in the orbital frame in terms of classical element perturbations

The transformation of classical orbit element perturbations to perturbations in position and velocity in the radial, transverse and normal directions of the orbital frame is developed. The formulation is given for the case of mean anomaly perturbations as well as for eccentric and true anomaly perturbations. Approximate formulas are also developed for the case of nearly circular orbits and compared with those found in the literature.     

On the numerical transformation of variables in perturbation theory

Once the generating function of a Lie-type transformation is known, canonical variables can be transformed numerically by application of a Runge-Kutta type integration method or any other appropriate numerical integration algorithm. The proposed approach avails itself of the fact, that the transformation is defined by a system of differential equations with a small parameter as the independent variable. The integration of such systems arising in the perturbation theories of Hori and Deprit is discussed. The method allows to compute numerical values of periodic perturbations without deriving explicitly the perturbation series. This saving of an algebraic work is achieved at the expense of multiple evaluations of the generator's derivatives.     

Radiative Transfer in Inhomogeneous Atmospheres. III

The approach proposed in the previous parts of this series of papers is used to solve the radiative transfer problem in scattering and absorbing multicomponent atmospheres. Linear recurrence relations are obtained for both the reflectance and transmittance of these kinds of atmospheres, as well as for the emerging intensities when the atmosphere contains energy sources. Spectral line formation in a one-dimensional inhomogeneous atmosphere is examined as an illustration of the possibility of generalizing our approach to the matrix case. It is shown that, in this case as well, the question reduces to solving an initial value problem for linear differential equations. Some numerical calculations are presented.     

Radiative Transfer in Inhomogeneous Atmospheres. I

This series of papers is devoted to multiple scattering of light in plane parallel, inhomogeneous atmospheres. The approach proposed here is based on Ambartsumyan's method of adding layers. The main purpose is to show that one can avoid difficulties with solving various boundary value problems in the theory of radiative transfer, including some standard problems, by reducing them to initial value problems. In this paper the simplest one dimensional problem of diffuse reflection and transmission of radiation in inhomogeneous atmospheres with finite optical thicknesses is considered as an example. This approach essentially involves first determining the reflection and transmission coefficients of the atmosphere, which, as is known, are a solution of the Cauchy problem for a system of nonlinear differential equations. In particular, it is shown that this system can be replaced with a system of linear equations by introducing auxiliary functions P and S. After the reflectivity and transmissivity of the atmosphere are determined, the radiation field in it is found directly without solving any new equations. We note that this approach can be used to obtain the required intensities simultaneously for a family of atmospheres with different optical thicknesses. Two special cases of the functional dependence of the scattering coefficient on the optical thickness, for which the solutions of the corresponding equations can be expressed in terms of elementary functions, are examined in detail. Some numerical calculations are presented and interpreted physically to illustrate specific features of radiative transport in inhomogeneous atmospheres.     

Radiative Transfer in Inhomogeneous Atmospheres. II

This paper is a continuation of a study of radiative transfer in one-dimensional inhomogeneous atmospheres. Two of the most important characteristics of multiple scattering in these media are calculated: the photon escape probability and the average number of scattering events. The latter is determined separately for photons leaving the medium and for photons that have undergone thermalization in the medium. The problem of finding the radiation field in an inhomogeneous atmosphere containing energy sources is also examined. It is assumed that the power of these sources, as well as the scattering coefficient, can vary arbitrarily with depth. It is shown that knowledge of the reflection and transmission coefficients of the atmosphere makes it possible to reduce all these problems to solving some first order linear differential equations with specified initial conditions. A series of new analytic results are obtained. Numerical calculations are done for two types of atmosphere with different depth dependences for the scattering coefficient. These are interpreted physically.     

The Main Problem in Satellite Theory Revisited

Using the elimination of the parallax followed by the Delaunay normalization, we present a procedure for calculating a normal form of the main problem (J ₂ perturbation only) in satellite theory. This procedure is outlined in such a way that an object-oriented automatic symbolic manipulator based on a hierarchy of algebras can perform this computation. The Hamiltonian after the Delaunay normalization is presented to order six explicitly in closed form, that is, in which there is no expansion in the eccentricity. The corresponding generating function and transformation of coordinates, too lengthy to present here to the same order; the generator is given through order four. This revised version was published online in July 2006 with corrections to the Cover Date.     

Asteroid Mean Elements: Higher Order and Iterative Theories

Mean orbital elements are obtained from osculating ones by removing the short periodic perturbations. Large catalogues of asteroid mean elements need to be computed, as a first step in the computation of proper elements, used to study asteroid families. The algorithms for this purpose available so far are only accurate to first order in the masses of the perturbing planet; the mean elements have satisfactory accuracy for most of the asteroid belt, but degraded accuracy in the neighbourhoods of the main mean motion resonances, especially the 2:1. We investigate a number of algorithms capable of improving this approximation; they belong to the two classes of Breiter-type methods and iterative methods. The former are obtained by applying some higher order numerical integration scheme, such as Runge–Kutta, to the differential equation whose solution is a transformation removing the fast angular variables from the equations; they can be used to compute a full second order theory, however, only if the full second order determining function is explicitly computed, and this is computationally too cumbersome for a complicated problem such as the N-body. The latter are fixed point iterative schemes, with the first order theory as an iteration step, used to compute the inverse map from mean to osculating elements; formally the method is first order, but because they implement a fixed frequency perturbation theory, they are more accurate than conventional single iteration methods; a similar method is already in use in our computation of proper from mean elements. Many of these methods are tested on a sample of asteroid orbits taken from the Themis family, up to the edge of the 2:1 resonance, and the dispersion of the values of the computed mean semimajor axis over 100 000 years is used as quality control. The results of these tests indicate that the iterative methods are superior, in this specific application, to the Breiter methods, in accuracy and reliability. This is understood as the result of the cancellations occurring between second order perturbation terms: the incomplete second order theory, resulting from the use of a Breiter method with the first order determining function only, can be less accurate than complete, fixed frequency theories of the first order. We have therefore computed new catalogues of asteroid mean and proper elements, incorporating an iterative algorithm in both steps (osculating to mean and mean to proper elements). This new data set, significantly more reliable even in the previously degraded regions of Themis and Cybele, is in the public domain. This revised version was published online in July 2006 with corrections to the Cover Date.     

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.     

任意倾角的共轨运动研究

共轨运动天体与摄动天体的半长径相同,处于1:1平运动共振中.太阳系内多个行星的特洛伊天体即为处于蝌蚪形轨道的共轨运动天体,其中一些高轨道倾角特洛伊天体的轨道运动与来源仍未被完全理解.利用一个新发展的适用于处理1:1平运动共振的摄动函数展开方式,对三维空间中的共轨运动进行考察,计算不同初始轨道根数情况下共轨轨道的共振中心、共振宽度,分析轨道类型与初始轨道根数的关系.并将分析方法所得结果与数值方法的结果相互比较验证,得到了广阔初始轨道根数空间内共轨运动的全局图景.     

The 3s 2S–4p 2Po transition probability in Mg ii

The superposition of configurations (SOC) method has been used to calculate f -values for the Mg  ii doublet at 1240 Å, which has been observed in the interstellar medium near stars such as ζ Ophiuchi. Our best value for the multiplet oscillator strength is 0.00083. SOC calculations have also been undertaken for the stronger (3s−3p) doublet at 2800 Å giving a multiplet f -value of 0.92.