Note on the Generalized Hansen and Laplace Coefficients

Recently, Breiter et al. [Celest. Mech. Dyn. Astron., 2004, 88, 153–161] reported the computation of Hansen coefficients X _k ^{γ ,m} for non-integer values of γ. In fact, the Hansen coefficients are closely related to the Laplace b _s ^(m), and generalized Laplace coefficients b _s,r ^(m) [Laskar and Robutel, 1995, Celest. Mech. Dyn. Astron., 62, 193–217] that do not require s,r to be integers. In particular, the coefficients X ₀ ^{γ ,m} have very simple expressions in terms of the usual Laplace coefficients b _{γ +2} ^(m), and all their properties derive easily from the known properties of the Laplace coefficients.     

Some Recurrence Relations Between Hansen Coefficients

A new system of recurrence relations for Hansen coefficients is obtained. This system gives a connection between only those coefficients which are included in the disturbing function of planetary or satellite motion and allows to compute efficiently the Hansen coefficients for perturbations both from internal and external bodies. The recurrence process can be realized both from high to low and from low to high harmonical terms of the disturbing function. The corresponding algorithms of evaluation of Hansen coefficients are presented. The efficiency of the obtained system of recurrence relations is discussed.     

Analytic properties of Hansen coefficients

Hansen’s coefficients in the theory of elliptic motion with eccentricity e are studied as functions of the parameter η = (1 − e ²)^1/2. Their analytic behavior in the complex η plane is described and some symmetry relations are derived. In particular, for every Hansen coefficient, multiplication by suitable powers of e and η results in an entire analytic function of η. Consequently, Hansen’s coefficients can be in principle computed by means of rapidly convergent series in powers of η. A representation of Hansen’s coefficients in terms of two entire functions of e ² follows.      

Computation of elliptic Hansen coefficients and their derivatives

The problem of computation of elliptic Hansen coefficients and their derivatives is considered for constructing a motion theory of an artificial Earth satellite with large eccentricity. An algorithm for analytical and numerical computation of these coefficients and their derivatives is described. The recurrence relations for derivatives of the first and second order and initial values for recurrences are obtained. As an example, numerical values of some elliptic Hansen coefficients are given for the orbit with eccentricityk=0.74.     

Numerical Computation of Hansen-like Expansions

Fourier expansions of elliptic motion functions in multiples of the true, eccentric, elliptic and mean anomalies are computed numerically by means of the fast Fourier transform. Both Hansen-like coefficients and their derivatives with respect to eccentricity of the orbit are considered. General behavior of the coefficients and the efficiency (compactness) of the expansions are investigated for various values of eccentricity of the orbit. This revised version was published online in July 2006 with corrections to the Cover Date.     

两个新的Hansen系数的递推公式

推导给出了两个新的Hansen系数X_k~(-(n+1),m)的递推关系:■其中,n、m和k是Hansen系数X_k~(-(n+1),m)的3个指标,e为轨道偏心率.递推公式(R5)可以执行普通Hansen系数的向后递推,需要一行初值,公式简单.递推公式(R6)可以执行偏心率函数的向前递推,需要两行初值,比Vakhidov给出的递推公式明显简单.算例说明,这两种递推是有效的.     

Recursive calculation of hansen coefficients

Hansen coefficients are used in expansions of the elliptic motion. Three methods for calculating the coefficients are studied: Tisserand's method, the Von Zeipel-Andoyer (VZA) method with explicit representation of the polynomials required to compute the Hansen coefficients, and the VZA method with the values of the polynomials calculated recursively. The VZA method with explicit polynomials is by far the most rapid, but the tabulation of the polynomials only extends to 12th order in powers of the eccentricity, and unless one has access to the polynomials in machine-readable form their entry is laborious and error-prone. The recursive calculation of the VZA polynomials, needed to compute the Hansen coefficients, while slower, is faster than the calculation of the Hansen coefficients by Tisserand's method, up to 10th order in the eccentricity and is still relatively efficient for higher orders. The main advantages of the recursive calculation are the simplicity of the program and one's being able to extend the expansions to any order of the eccentricity with ease. Because FORTRAN does not implement recursive procedures, this paper used C for all of the calculations. The most important conclusion is recursion's genuine usefulness in scientific computing.     

Hansen系数及其导数的直接计算方法

回顾总结了7种Hansen系数及其导数的直接计算方法,比较分析了这些方法的计算效率和计算稳定性.研究表明:Hansen系数的递推关系可以用来判别计算结果的稳定性.最后指出, Wnuk方法(双精度计算)和McClain方法(4精度计算)是稳定的,可以用来计算人造卫星轨道摄动.由于大多数人造卫星采用小偏心率轨道,需要计算无奇点摄动,推荐使用McClain方法1 (4精度计算).     

Direct Calculation Methods of Hansen Coefficients and Their Derivatives

Seven direct calculation methods of Hansen coefficients and their derivatives are reviewed. The computational efficiencies of these methods are compared, and their computational stabilities are analyzed. We show that the recursion relations of Hansen coefficients can be used to determine the stabilities of calculation results. Finally, it is pointed out that Wnuk's method (double precision computation) and McClain's methods (quadruple precision computation) are stable, which can be used to calculate orbit perturbations. Because of small orbital eccentricities of most satellites, the perturbation calculations without singularities are required, and McClain's first method (quadruple precision computation) is recommended.     

Decomposition of functions for elliptic orbits

In the algebraA of functions periodic in the mean anomaly we relate the problem of integrating over the mean anomaly with that of decomposing an element ofA as the direct sum of two functions, one in the kernel of the Lie derivative in the Keplerian flow and one in the image of this Lie derivative. We propose recursive rules amenable to general purpose symbolic processors for accomplishing such decomposition in a wide subclass ofA. We introduce the dilogarithmic function to express in exact terms quadratures involving the equation of the center.     

A New Method to Determine the Derivatives of the Laplace Coefficients

The determination of the secular variations of the orbital elements of objects in N-body systems is based on the literal development of the perturbing function. The development makes use of the Laplace coefficients and their derivatives. In this paper a new method is described for the analytical computation of the derivatives of the Laplace coefficients. It is an explicit formula in the sense that it only contains the Laplace coefficients and the parameter on which the Laplace coefficients depend. The advantage of this method is that it is unnecessary to calculate all the derivatives up to the desired order. It is enough to calculate the Laplace coefficients. Easy coding is a further benefit of the method and it provides more accurate numerical results. The paper describes in detail the application of the method through an example and gives comparison with former methods.This revised version was published online in October 2005 with corrections to the Cover Date.     

The significance of the Hansen Ideal space frame

Known and unknown properties of Hansen Ideal coordinates are summarized. It is shown that the ideal space frame is a general and necessary component of basic celestial mechanics and astrodynamics, as well as of any theory of motion. A typical consequence is the intimate correlation of the Hansen frame with the Lagrange constraint within the method of the variation of the parameters. The use of observations in the ideal frame may allow conclusions on the intergalactic fundamental coordinate system (© 2012 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Generalized cumulant correlators and hierarchical clustering

The cumulant correlators, C _pq , are statistical quantities that generalize the better-known S _p parameters; the former are obtained from the two-point probability distribution function of the density fluctuations while the latter describe only the one-point distribution. If galaxy clustering develops from Gaussian initial fluctuations and a small-angle approximation is adopted, standard perturbative methods suggest a particular hierarchical relationship of the C _pq for projected clustering data, such as that obtained from the Automatic Plate Measuring (APM) survey. We establish the usefulness of the two-point cumulants for describing hierarchical clustering by comparing such calculations against available measurements from projected catalogues, finding very good agreement. We extend the idea of cumulant correlators to multipoint generalized cumulant correlators (related to the higher-order correlation functions). We extend previous studies in the highly non-linear regime to express the generalized cumulant correlators in terms of the underlying 'tree amplitudes' of hierarchical scaling models. Such considerations lead to a technique for determining these hierarchical amplitudes, to arbitrary order, from galaxy catalogues and numerical simulations. Knowledge of these amplitudes yields important clues about the phenomenology of gravitational clustering. For instance, we show that a three-point cumulant correlator can be used to separate the tree amplitudes up to sixth order. We also combine the particular hierarchical Ansatz of Bernardeau & Schaeffer with extended and hyper-extended perturbation theory to infer values of the tree amplitudes in the highly non-linear regime.     

Secular Motion in a 2nd Degree and Order-Gravity Field with no Rotation

The motion of a particle about a non-rotating 2nd degree and order-gravity field is investigated. Averaging conditions are applied to the particle motion and a qualitative analysis which reveals the general character of motion in this system is given. It is shown that the orbit plane will either be stationary or precess about the body's axis of minimum or maximum moment of inertia. It is also shown that the secular equations for this system can be integrated in terms of trigonometric, hyperbolic or elliptic functions. The explicit solutions are derived in all cases of interest.This revised version was published online in October 2005 with corrections to the Cover Date.     

Generalized Laplace coefficients and Newcomb derivatives

It is shown how the generalized Laplace coefficients can be employed to deduce explicit formulas for ordinary and Newcomb derivatives of the Laplace coefficients.      

Hansen系数递推的效率

探讨了偏心率函数的递推方法的效率, 给出了一种成批递推方法, 其计算效率明显优于直接计算方法, 而且递推是正向的. 采用该成批递推方法递推时, 偏心率函数的量级从小到大变化, 可以保证递推的精度.     

The Relegation Algorithm

The relegation algorithm extends the method of normalization by Lie transformations. Given a Hamiltonian that is a power series = ₀+ ₁+ ... of a small parameter , normalization constructs a map which converts the principal part ₀into an integral of the transformed system — relegation does the same for an arbitrary function [G]. If the Lie derivative induced by [G] is semi-simple, a double recursion produces the generator of the relegating transformation. The relegation algorithm is illustrated with an elementary example borrowed from galactic dynamics; the exercise serves as a standard against which to test software implementations. Relegation is also applied to the more substantial example of a Keplerian system perturbed by radiation pressure emanating from a rotating source.This revised version was published online in October 2005 with corrections to the Cover Date.     

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. 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. 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.