Effects on the nutation of C4,m and S4,m Harmonics

In this paper, we make a study about the influence of the coefficients of the geopotential C_4,m and S_4,m, (m=1,2,3,4) on the nutation, starting from the Hamiltonian theory as developed by Kinoshita (1977).We obtain ten coefficients larger than 0.05 μ as for the nutation in longitude and six for the nutation in obliquity. The present results are included in the reconstruction of the theory of nutation (REN‐2000) at the level of truncation of 0.1 μ as (Souchay et al., 1997). This revised version was published online in July 2006 with corrections to the Cover Date.     

RDAN97: An Analytical Development of Rigid Earth Nutation Series Using the Torque Approach

New series of rigid Earth nutations for the angular momemtum axis, the rotation axis and the figure axis, named RDAN97, are computed using the torque approach. Besides the classical J2 terms coming from the Moon and the Sun, we also consider several additional effects: terms coming from J3 and J4 in the case of the Moon, direct and indirect planetary effects, lunar inequality, J2 tilt, planetary‐tilt, effects of the precession and nutations on the nutations, secular variations of the amplitudes, effects due to the triaxiality of the Earth, new additional out‐of‐phase terms coming from second order effect and relativistic effects. Finally, we obtain rigid Earth nutation series of 1529 terms in longitude and 984 terms in obliquity with a truncation level of 0.1 μ (microarcsecond) and 8 significant digits. The value of the dynamical flattening used in this theory is HD=(C-A)/C=0.0032737674 computed from the initial value pa=50′.2877/yr for the precession rate. These new rigid Earth nutation series are then compared with the most recent models (Hartmann et al., 1998; Souchay and Kinoshita, 1996, 1997; Bretagnon et al., 1997, 1998. We also compute a benchmark series (RDNN97) from the numerical ephemerides DE403/LE403 (Standish et al., 1995) in order to test our model. The comparison between our model (RDAN97) and the benchmark series (RDNN97) shows a maximum difference, in the time domain, of 69 μas in longitude and 29 μas in obliquity. In the frequency domain, the maximum differences are 6 μas in longitude and 4 μ as in obliquity which is below the level of precision of the most recent observations (0.2 mas in time domain (temporal resolution of 1 day) and 0.02 mas in frequency domain). This revised version was published online in July 2006 with corrections to the Cover Date.     

On the effect of the mantle elasticity on the earth's rotation

Starting from the Hamiltonian model for a solid Earth with an elastic mantle previously developped by the authors, analytical expressions are derived which give the nutation series corresponding to the plane perpendicular to the angular momentum vector, to the plane perpendicular to the rotational axis and to the equator of figure, as well as the series that give the polar motion. The effects of the different perturbations — solid Earth, centrifugal and tidal potentials — are calculated separately. The corrections due to the elasticity of the mantle, which mostly correspond to the Oppolzer terms, are calculated with an accuracy of 10^–6 arc sec., given that the intrinsic observational accuracy has reached 0.01 mas.     

Fast computation of satellite gravitational gradient

The computation of the Earth's potential function at high order and degree with the method of reference [1], causes overflow most of the time. The normalized method [2–6] can eliminate the overflows, but leads to formulae much more involved than those in reference [1]; besides, the programming is more complex and the computer time required larger. The method presented in this paper has the following features: each component of the satellite gravitational gradient can be computed; the formulae are short and easy to be programme; the method is much quicker than the normalization method and can be carried out with a micro‐computer, without overflow even in the case of Earth's spherical harmonics of order and degree as high as 1025 or higher. This method satisfies the present demand to compute satellite gravitational gradient with high accuracy. Furthermore, we present formulae for the fast computation, without overflow, of the gravitational gradient corresponding to Earth's spherical harmonics up to order and degree of 3170 × 3170 or higher. This revised version was published online in July 2006 with corrections to the Cover Date.     

Harmonics of high-amplitude wave trains in cosmic ray intensity

Using the ground based neutron monitor data of Deep River, the high-amplitude anisotropic wave train events (HAE) in cosmic ray intensity have been investigated during the period 1991-1994. It has been observed that the phase of diurnal anisotropy for majority of HAE shifts towards later hours; whereas it remains in the corotational/18-h direction for some of the HAE cases. Further, for majority of HAE cases the amplitude of diurnal and semi-diurnal anisotropy significantly deviates from the annual average values. The phase of semi-diurnal and tri-diurnal anisotropy for all HAE cases has shifted to later hours. Furthermore, for tri-diurnal anisotropy the amplitude remains statistically the same. The occurrence of HAE is unaffected by the nature of the B_z component of IMF polarity.     

An Abridged Model of the Precession–Nutation of the Celestial Pole

Precise astrometric observations show that significant systematic differences of the order of 10 milliarcseconds (mas) exist between the observed position of the celestial pole in the International Celestial Reference Frame (ICRF) and the position determined using the International Astronomical Union (IAU) 1976 Precession (Lieske et al., 1977) and the IAU 1980 Nutation Theory (Seidelmann, 1982). The International Earth Rotation Service routinely publishes these 'celestial pole offsets', and the IERS Conventions (McCarthy, 1996) recommends a procedure to account for these errors. The IAU, at its General Assembly in 2000, adopted a new precession/nutation model (Mathews et al., 2002). This model, designated IAU2000A, which includes nearly 1400 terms, provides the direction of the celestial pole in the ICRF with an accuracy of ±0.1 mas. Users requiring accuracy no better than 1 mas, however, may not require the full model, particularly if computational time or storage are issues. Consequently, the IAU also adopted an abridged procedure designated IAU2000B to model the celestial pole motion with an accuracy that does not result in a difference greater than 1 mas with respect to that of the IAU2000A model. That IAU2000B model, presented here, is shown to have the required accuracy for a period of more than 50 years from 1995 to 2050.     

On the inclination functions and a rapid stable procedure for their evaluation together with derivatives

The authors’ individual work on the inclination functions over a period of more than 30 years has led to the need for a joint paper. Intervening papers by other authors have demonstrated misunderstandings needing to be corrected, in particular concerning the key recurrence relation published by the present first author in 1971. This relation is remarkably stable, though this has not always been recognized. The real source of error with the specific functions that are involved in the recurrence relation arises from the possibilities for underflow (as well as overflow) in the computation. The problem exists even with normalized versions of the functions, and is carefully addressed. Very important, for both academic and practical reasons, is a general invariance relation that had been found earlier by the second author, for which a proof is given here for the first time. Some numerical results from our new (and highly efficient) procedure for computing the inclination functions are tabulated, and comparisons made with the results of other authors. Finally, Fortran code for an optimized implementation of this procedure is in supplementary material. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Third Integral and the Dynamics of the Galaxy in the Neighborhood of the Sun

Observational data on the dynamics of stars in the neighborhood of the sun indicate the existence of a third integral besides the integrals of the angular momentum and energy. The Poincaré integral is proposed as a third integral. The consequences of this assumption are derived and compared with available astrophysical data.     

Long period variations of the motion of a satellite due to non-resonant tesseral harmonics of a gravity potential

The main effects of tesseral harmonics of a gravity potential expansion on the motion of a satellite, are short period variations as well as long period variations due to resonances. However, other smaller long period and secular variations can arise from interactions between tesseral terms of the same order. The analytical integration of these effects is developed, using numerical evaluation of Kaula eccentricity and inclination functions. Examples for some Earth's geodetic satellites show that secular effects can reach a few decameters per year. The secular variations can even reach several hundred of meters per year for the Mars natural satellite Phobos.     

Gravitational potential harmonics from the shape of an homogeneous body

The spherical harmonic coefficients of the gravitational potential of an homogeneous body are analytically derived from the harmonics describing its shape. General formulas are given as well as detailed expressions up to the fifth order of the topography harmonics. The volume, surface and inertia tensor of the body are obtained as by-products. The case of a triaxial ellipsoid is given as example and used for numerical checking. Another numerical scheme for verification is provided. The application to Phobos is made and the convergence of the expressions for the harmonics is numerically established.

---

Résumé Les harmoniques du champ de gravitation d'un corps homogène de forme donnée sont calculés analytiquement à partir des harmoniques du développement en série du rayon vecteur exprimant la forme de la surface du corps. Outre la formule générale, des expressions détaillées, au cinquième ordre des harmoniques du rayon vecteur, sont données sous une forme bien adaptée à la programmation. Le volume, la surface et le tenseur d'inertie du corps sont calculés analytiquement `a partir des formules générales. Le cas de l'ellipsoide triaxial est pris comme test des formules établies. Un autre test numérique est fourni dans le cas le plus général. Ceci est appliqué à Phobos, et la convergence des expressions fournissant les harmoniques est numériquement démontrée.

     

An Analytical Solution of the Lagrange Equations Valid also for Very Low Eccentricities: Influence of a Central Potential

This paper presents an analytic solution of the equations of motion of an artificial satellite, obtained using non singular elements for eccentricity. The satellite is under the influence of the gravity field of a central body, expanded in spherical harmonics up to an arbitrary degree and order. We discuss in details the solution we give for the components of the eccentricity vector. For each element, we have divided the Lagrange equations into two parts: the first part is integrated exactly, and the second part is integrated with a perturbation method. The complete solution is the sum of the so-called “main” solution and of the so-called “complementary” solution. To test the accuracy of our method, we compare it to numerical integration and to the method developed in Kaula (Theory of Satellite Geodesy, Blaisdell publ. Co., New York. 1966), expressed in classical orbital elements. For eccentricities which are not very small, the two analytical methods are almost equivalent. For low eccentricities, our method is much more accurate.     

Long-Period Variations of the Eccentricity Vector Valid also for Near Circular Orbits around a Non-Spherical Body

This paper studies the long period variations of the eccentricity vector of the orbit of an artificial satellite, under the influence of the gravity field of a central body. We use modified orbital elements which are non-singular at zero eccentricity. We expand the long periodic part of the corresponding Lagrange equations as power series of the eccentricity. The coefficients characterizing the differential system depend on the zonal coefficients of the geopotential, and on initial semi-major axis, inclination, and eccentricity. The differential equations for the components of the eccentricity vector are then integrated analytically, with a definition of the period of the perigee based on the notion of “free eccentricity”, and which is also valid for circular orbits. The analytical solution is compared to a numerical integration. This study is a generalization of (Cook, Planet. Space Sci., 14, 1966): first, the coefficients involved in the differential equations depend on all zonal coefficients (and not only on the very first ones); second, our method applies to nearly circular orbits as well as to not too eccentric orbits. Except for the critical inclination, our solution is valid for all kinds of long period motions of the perigee, i.e., circulations or librations around an equilibrium point.     

On the Dynamics and Topology of the Elliptic Rectilinear Restricted 3-Body Problem

We study a symmetric collinear restricted 3-body problem, where the equal mass primaries perform elliptic collisions, while a third massless body moves in the line between the primaries, during the time between two consecutive elliptic collisions. After desingularizing binary and triple collisions, we prove the existence of a transversal heteroclinic orbit beginning and ending in triple collision. This orbit is the unique homothetic orbit that the problem possess. Finally, we describe the topology of the compact extended phase space. This revised version was published online in July 2006 with corrections to the Cover Date.     

On the period changes of RR Lyrae variables in the globular cluster M3

O – C diagrams of 78 RR Lyrae variables of the years 1963 to 1978 were used to investigate their comparability with the diagrams of 1895 to 1963. Only 5% have been proved to show a fundamentally deviating trend of the period change.     

DCMT光电扫描主测微器的设计与研究

DCMT主测微器不同于其它类型的子午环测微器,它具有自校准测定仪器参数的功能。该测微器采用了活动光栅的方案,其优点是能观测近极星和各类准直像;活动光栅另一个显著优点是不同赤纬星几乎可用相同的观测时间.对“V”形光栅的工作原理和误差进行了详细讨论,并给出了一组严格的公式。其系统误差来源有:光栅形状改正、光栅驱动方向相对于光栅的倾斜、光栅驱动方向相对于赤径方向的倾斜、星径曲率改正。     

On the internal structure of relativistic jets

A magnetohydrodynamic model is constructed for a cylindrical jet embedded in an external uniform magnetic field. It is shown that, as in the force-free case, the total electric current within the jet can be zero. The particle energetics and the magnetic-field structure are determined in a self-consistent way; all jet parameters depend markedly on the physical conditions in the external medium. In particular, we show that a region with subsonic flow can exist in the central jet regions. In real relativistic jets, most of the energy is transferred by the electromagnetic field only at a sufficiently large magnetization parameter σ>10⁶. We also show that, in general, the well-known solution with a central core B _Z =B ₀/(1+?²/? _c ² ) cannot be realized in the presence of an external medium.     

On the evolution of the cosmic supernova rates

The Galactic radio-emitting X-ray binary Cygnus X-3 is known to be a source of large-scale radio jets associated with periods of intense radio flaring. These jets have been found to have an expansion velocity of ∼0.3 c and are believed (on kinematic grounds) to lie close to the plane of the sky. We present new observations of Cygnus X-3 using the VLBA at 15 GHz. These observations, which included the detection of two small flares, show an additional kind of behaviour with apparent superluminal expansion along both major and minor axes. Evidence for superluminal activity has been found in a number of X-ray binary systems such as GRS 1915+105 and GRO J1655−40 with their superluminal radio jets. Apparently similar morphologies of the Galactic and extragalactic jet sources have led to the X-ray binaries being described as 'micro-quasars'. The superluminal expansion seen in our results appears to be different in nature from these other two sources, and a number of mechanisms are presented and discussed.     

On the variational approach to the periodic n-body problem

This expository paper gathers some of the results obtained by the author in recent works in collaboration with Davide Ferrario and Vivina Barutello, focusing on the periodic n-body problem from the perspective of the calculus of variations and minimax theory. These researches were aimed at developing a systematic variational approach to the equivariant periodic n-body problem in the two and three-dimensional space. The purpose of this paper is to expose the main problems and achievements of this approach. The material here was exposed in the talk that given at the Meeting CELMEC IV promoted by SIMCA (Società italiana di Meccanica Celeste).     

On the instability of D-type ionization fronts

We investigate the stability of D-type ionization fronts (IF). We find the dispersion relation for surface modes of weak-D IF illuminated by a plane-parallel radiation field, including both the inclination of the IF to the ionizing radiation field and the effects of recombination in the ionized gas. It is well known that IF are unstable in the absence of recombination in the ionized gas, but that recombination tends to stabilize them. For finite inclination, however, IF are unstable at essentially all wavelengths in the limit of very cold upstream gas.
We present numerical hydrodynamic simulations of the development of the instabilities. These confirm the presence of long-wavelength instabilities of IF when the sound speed in the neutral gas is a finite fraction of that in the ionized gas and the ionization cross section is finite. They also show the development of the instabilities to non-linear amplitude.     

彗星动力学研究综述

彗星是太阳系中天文学的重要研究对象．它涉及到天体物理、化学、天体力学等多种领域。近年的研究表明,彗星可能是揭开太阳系起源和演化之谜的突破口,还可能与生命的起源有关．本文旨在对彗星动力学的发展情况作一简单的回顾和评述．