Secular effects in the orbits of the Galilean satellites

This paper presents the results of an investigation into the secular behavior of the orbits of the Galilean satellites of Jupiter. Kamel's perturbation method is used to remove all the explicitly periodic variables from the differential equations that describe the long period behavior of the orbits to third order in the masses, and the resulting differential equations for the secular behavior are then solved. Several numerical examples are given to illustrate the sensitivity of the solution to variations in the masses of the satellites.     

The long period behavior of the orbits of the Galilean satellites of Jupiter

Some of the results of an investigation into the long period behavior of the orbits of the Galilean satellites of Jupiter are presented. Special purpose computer programs were used to perform all the algebraic manipulations and series expansions that are necessary to describe the mutual interactions among the satellites.The disturbing function was expanded as a Poisson series in the modified Keplerian elements referred to a Jovicentric coordinate system. The differential equations for the modified Keplerian elements were then formed, and all short period perturbations were removed using Kamel's perturbation method. Approximate analytical solutions for these differential equations are derived, and the general form of the solutions are given.     

Canonical planetary perturbation equations for velocity-dependent forces,and the Lense-Thirring precession

A form of planetary perturbation theory based on canonical equations of motion, rather than on the use of osculating orbital elements, is developed and applied to several problems of interest. It is proved that, with appropriately selected initial conditions on the orbital elements, the two forms of perturbation theory give rise to identical predictions for the observable coordinates and velocities, while the orbital elements themselves may be strikingly different. Differences between the canonical form of perturbation theory and the classical Lagrange planetary perturbation equations are discussed. The canonical form of perturbation theory in some cases has advantages when the perturbing forces are velocity-dependent, but the two forms of perturbation theory are equivalent if the perturbing forces are dependent only on position and not on velocity. The canonical form of the planetary perturbation equations are derived and applied to the Lense Thirring precession of a test body in a Keplerian orbit around a rotating mass source.     

Structure-distortion of polytropic stars by differential rotation: Third-order perturbation theory

In order to specify the structure of a rapidly and differentially rotating gaseous polytrope, we extend Chandrasekhar's perturbation theory to include third-order terms in the perturbation parameter. In the present paper, the theory developed is required for a subsequent numerical treatment of the structure-determination.This research was supported by the Research Development Project of the University of Patras, Greece.     

The gravitational perturbation spectrum in linear satellite theory

A new method for calculating the perturbation spectrum in the framework of Kaula's linear satellite theory (LST) is introduced. The novelty of this approach consists in using recent results on the spectral decomposition of the perturbation frequencies in LST to provide a closed formulation for the amplitude and the phase of each line in the perturbation spectrum. The theory presented here can be applied to perturbations in the elements or in the radial and transverse directions due to the geopotential or to the tides. Separate algorithms are developed for application to orbits with circulating or frozen perigee.     

The oscillation frequencies of magnetic stellar models

A first-order perturbation theory method developed by Goossens to determine the perturbation to the eigenfrequencies of stellar models caused by the presence of a magnetic field is modified slightly, and applied to models with toroidal and poloidal fields. Some limitations of the analysis are pointed out.     

Nonlinear perturbation theory for cosmic ray propagation in random magnetic fields

A nonlinear perturbation theory is applied to the problem of pitch angle diffusion of energetic particles in random magnetic fields. To keep the analysis simple, the discussion is restricted to fluctuation fields, consisting of Alfvén waves. It is shown that the failure of quasilinear theory at small particle velocities parallel to the average field can be overcome by a statistically exact treatment of the particle orbits in the first order fields. In fact, for spherical power spectra which, in addition, do not fall off too steeply with increasing frequency, the conventional perturbation theory also leads to formally convergent expressions for the scattering mean free path. These results are shown to be quite satisfactory, even in a quantitative sense. For more general physically realistic power spectra, however, a divergence-free diffusion theory is indispensible. A simple representation for the resulting pitch-angle diffusion coefficient is suggested.     

Expansion formulae in canonical transformations depending on a small parameter

The theory of perturbation based on Lie transforms is considered. Deprit's equation is reduced to a form which enables us to generate simplified general recursion formulae. These expansions are then modified to speed up the implementation of such perturbation theory in the computerized symbolic manipulation.     

Modulation of ion acoustic waves

The multiple scales perturbation theory is applied to a system of dispersive waves including plasma waves (ion-acoustic waves). Assuming the amplitudes of waves are slowly varying function of space and time, we find that long-time slow modulation of the complex amplitude can be described by the non-linear Schrödinger equation. This result agrees with that obtained by Shimizu and Ichikawa using the reductive perturbation theory, and agrees exactly with the nonrelativistic limit of the nonlinear Schrödinger equation obtained by Nejoh using the stretching method.     

Travel Time Sensitivity Kernels

We derive, following the standard first Born approximation approach used in the geophysics literature, an expression for the travel time perturbation caused by a perturbation to sound speed. In our simple model we employ a point source at one point and calculate the time taken for a wave packet created at the source to move to a second point. In the first Born approximation the travel time delay caused by a perturbation to the background model can be expressed as the integral over the whole sun of some function, called the travel time sensitivity kernel, multiplied by the perturbation. The sensitivity kernels are zero along the geometrical ray connecting the two points and have maximum weight in a tube around the ray; they are the solar equivalent of `the banana-doughnut' kernels discussed in the geophysics literature. Calculating sensitivity kernels that are more accurate than those derived from ray theory is important for the future of inversions done with time-distance helioseismology data as they will allow greater confidence in the results as well as increased resolution.     

Asymptotic expansions in the perturbed two-body problem with application to systems with variable mass

The main theorems of the theory of averaging are formulated for slowly varying standard systems and we show that it is possible to extend the class of perturbation problems where averaging might be used. The application of the averaging method to the perturbed two-body problem is possible but involves many technical difficulties which in the case of the two-body problem with variable mass are avoided by deriving new and more suitable equations for these perturbation problems. Application of the averaging method to these perturbation problems yields asymptotic approximations which are valid on a long time-scale. It is shown by comparison with results obtained earlier that in the case of the two-body problem with slow decrease of mass the averaging method cannot be applied if the initial conditions are nearly parabolic. In studying the two-body problem with quick decrease of mass it is shown that the new formulation of the perturbation problem can be used to obtain matched asymptotic approximations.     

An extended canonical perturbation method

In this investigation, a procedure is described for extending the application of canonical perturbation theories, which have been applied previously to the study of conservative systems only, to the study of non-conservative dynamical systems. The extension is obtained by imbedding then-dimensional non-conservative motion in a 2n-dimensional space can always be specified in canonical form, and, consequently, the motion can be studied by direct application of any canonical perturbation method. The disadvantage of determining a solution to the 2n-dimensional problem instead of the originaln-dimensional problem is minimized if the canonical transformation theory is used to develop the perturbation solution. As examples to illustrate the application of the method, Duffing's equation, the equation for a linear oscillator with cubic damping and the van der Pol equation are solved using the Lie-Hori perturbation algorithm.This research was supported by the Office of Naval Research under Contract N00014-67-a-0126-0013.     

太阳多方模型对行星轨道要素变化的长期摄动影响

本文利用天体力学中的摄动理论和天体物理学中的气体星多方模型理论研究了太阳多方模型对行星轨道要素变化的长期摄动影响。文中给出了太阳日多方指数ｎ＝３的模型由于自转、扁度和内部密度分布等因素对行星轨道要素变化的长期摄动效应的理论结果。研究结果表明：行星轨道要素除长轴、偏心率和轨道倾角不受长期摄动外升交点经度、近日点经度以及平近点角均受长期摄动的影响。最后利用理论结果对行星轨道要素的长期投动效应做了数值计算，数值结果在表１中给出。     

The transformational behaviour of perturbation theories

The transformational behaviour of Hori's noncanonical perturbation theory (Hori 1971) as well as that of the theory of Krylov-Bogoliubof-Mitropolsky is studied. An integration procedure of the perturbation equations is based on the transformation properties that have been established.     

Virtual singularities in the artificial satellite theory

When the coordinate system used in perturbation theory presents a geometrical singularity and when the perturbation technique fails to take account of this, the solution developed may present singularities which are no longer easily explained by purely geometrical means. These singularities have been calledvirtual singularities by Deprit and Rom (1970). We propose to demonstrate that virtual singularities can in general be avoided by the use of Lie transforms. In general, it is sufficient to recognize that the original Hamiltonian function presents the d'Alembert characteristic with respect to pairs of action-angle variables and that the averaging operations preserve this characteristic. We then apply this criterion to the artificial satellite theory (for small to moderate eccentricity) showing that all of three possible virtual singularities can be avoided at the same time. A new set of elliptic elements, well suited to the problem at hand, is introduced.     

Tidal Evolution of the Kepler Candidate Two-planet Systems

The tidal evolution of ten Kepler candidate two-planet systems is investigated by using the general secular perturbation theory, and then a general picture of tidal evolution for these systems is described. Taking the KOI 1239 system as an example, the tidal effect of the system is studied in detail, the results indicate that the dissipative term of the tidal effect causes the attenuation of planets’ orbital eccentricities, and it plays a dominant role in the process of orbital evolution, however, the conservative term of the tidal effect and the relativistic effect may damp the excitation of the eccentricity of the inner planet under the secular perturbation of the outer planet. In addition, the process of tidal evolution is also affected by both the initial eccentricity of the outer planet and the planet's tidal dissipation coeffcient. At the same time, the numerical simulation on the tidal evolution of the KOI 1239 system is also made, and the numerical results are consistent with those of general secular perturbation theory.     

Large-scale bias in the Universe — II. Redshift-space bispectrum

The determination of the density parameter Ω₀ from the large-scale distribution of galaxies is one of the major goals of modern cosmology. However, if galaxies are biased tracers of the underlying mass distribution, linear perturbation theory leads to a degeneracy between Ω₀ and the linear bias parameter b , and the density parameter cannot be estimated. In Matarrese, Verde &38; Heavens we developed a method based on second-order perturbation theory to use the bispectrum to lift this degeneracy by measuring the bias parameter in an Ω₀-independent way. The formalism was developed assuming that one has perfect information on the positions of galaxies in three dimensions. In galaxy redshift surveys, the three-dimensional information is imperfect, because of the contaminating effects of peculiar velocities, and the resulting clustering pattern in redshift space is distorted. In this paper we combine second-order perturbation theory with a model for collapsed, virialized structures, to extend the method to redshift space, and demonstrate that the method should be successful in determining with reasonable accuracy the bias parameter from state-of-the-art surveys such as the Anglo-Australian 2 degree Field Survey and the Sloan Digital Sky Survey.     

Equivalence of the perturbation theories of Hori and Deprit

A comparison is made between the perturbation theories of Hori and Deprit which are based on the use of Poisson brackets. General recurrence formulae are presented for Hori's theory which are analogous to those in Deprit's theory. Explicit relations between the determining functions for the two theories are indicated through the sixth order, these results having been obtained by a novel computer program. A general argument for the equivalence of the theories to all orders is given.     

Construction of a consistent semianalytic theory of a planetary or moon orbiter perturbed by a third body

The third body perturbation of an orbiter of a planet or moon is considered. A very convenient form of the Lagrange equations is given allowing an easy derivation of the various terms of the expansion of the perturbed elements. A careful analysis of the order of magnitude of these terms indicates which ones are required for a consistent theory. It follows that in many practical cases the main term of the disturbing function has to be carried to the second order of the perturbation theory.Paper presented at the 1981 Oberwolfach Conference on Mathematical Methods in Celestial Mechanics.Dedicated to V. Szebehely on the occasion of his 60th birthday     

Restricted quantum-mechanical three-body problems

A new method is presented in a general form to solve the Schrödinger equation of helium-like ions. The wave function is expanded in terms of the eigenfunctions of a moving electron in the field of two Coulombic ions which are fixed in space. This makes the method similar to the Dirac perturbation theory (perturbation theory for time-dependent problems). In the present method an infinitely coupled system of infinitely many second-order ordinary differential equations must be solved instead of one second-order partial differential equation of three variables. The nature of the singular points and boundary conditions are discussed and some general relations are given which are useful for the numerical treatment.