Earth’s precession–nutation motion: the error analysis of the theories IAU 2000 and IAU 2006 applying the VLBI data of the years 1984–2006

The long-term systematic errors of the analytical theories IAU 2000 and IAU 2006 of the Earth’s precession–nutational motion are studied making use of the VLBI data of 1984–2007. Several independent methods give indubitable evidence of the significant quadratic error in the IAU 2000 residuals of the precessional angle while the adopted value of the secular decrease /cy of the Earth’s ellipticity e (derived from Satellite Laser Ranging data) should manifest itself in the residuals of as the negative quadratic trend . The problem with the precession of the IAU 2006 theory adopted as a new international standard and based on the precession model P03 (Capitaine et al., Astron Astrophys 432:355–367, 2005) appears to be even more serious because the above mentioned quadratic term has already been incorporated into the P03 precession. Our analysis of the VLBI data demonstrates that the quadratic trend of the IAU 2006 residuals does amount to the expected value (30.0 ± 3) mas/cy². It means, first, that the theoretical precession rate of IAU 2006 should be augmented by the large secular correction and, second, that the available VLBI data have potentiality of estimating the rate . And indeed, processing these data by the numerical theory ERA of the Earth’s rotation (Krasinsky, Celest Mech Dyn Astron 96:169–217, 2006, Krasinsky and Vasilyev, Celest Mech Dyn Astron 96:219–237, 2006) yields the estimate /cy statistically in accordance with the satellite-based . On the other hand, applying IAU 2000/2006 models, the positive value /cy is found which is incompatible with the SLR estimate and, evidently, has no physical meaning. The large and steadily increasing error of the precession motion of the IAU 2006 theory makes the task of replacing IAU 2006 by a more accurate model be most pressing.     

Completely Integrable Systems Connected with Lie Algebras

In a recent paper Ballersteros and Ragnisco (1998) have proposed a new method of constructing integrable Hamiltonian systems. A new class of integrable systems may be devised using the following sequence: , where A is a Lie algebra is a Lie–Poisson structure on R ₃, C is a Casimir for is a reduced Poisson bracket and (A, ▵) is a bialgebra. We study the relation between a Lie-Poisson stucture Λ and a reduced Poisson bracket , which is a key element in using the Lie algebra A to constructing this sequence. New examples of Lie algebras and their related integrable Hamiltonian systems are given. This revised version was published online in July 2006 with corrections to the Cover Date.     

The Evolution of Dilatonic Universe with a Double Exponential Potential

In this paper, dilaton in Weyl-Scaled induced gravitational theory is regarded as a candidate of dark energy. When the potential of dilaton field is taken as the form of a double exponential , we find that there exist attractor solutions in dilatonic dark energy model, and these attractors correspond to an equations of state and a cosmic density parameter , which are important features for a dark energy model that can meet the current observations. We find out the sufficient condition of the existence of a late time de Sitter attractor.     

Numerical theory of rotation of the deformable Earth with the two-layer fluid core. Part 2: Fitting to VLBI data

VLBI-based offsets of the Celestial Pole positions, as well as the variations of UT (series of Goddard Space Flight Center, 1984–2005) are processed applying the Earth’s rotation theory (ERA) 2005 constructed by the numerical integration of the differential equations of rotation of the deformable Earth. The equations were published earlier (Krasinsky 2006) as the first part of the work. The resulting weighted root mean square (WRMS) errors of the residuals , for the angles of nutation and precession are 0.136 and 0.129 mas, respectively. They are significantly less than the corresponding values 0.172 and 0.165 mas for the IAU 2000 model adopted as the international standard. In ERA 2005, the angles , are related to the inertial ecliptical frame J2000, the angle including the precessional secular motion. As the published observational data are theory-dependent being related to IAU 2000, a procedure to confront the numerical theory to the observed Celestial Pole offsets and UT variations is developed. Processing the VLBI data has shown that beside the well known 435-day FCN mode of the free core nutation, there exits a second mode, FICN, caused by the inner part of the fluid core, with the period of 420 day close to that of the FCN mode. Beatings between the two modes are responsible for the apparent damping and excitation of the free oscillations, and are implicitly modeled by ERA 2005. The nutational and precessional motions in ERA 2005 are proved to be mutually consistent but only in case the relativistic correction for the geodetic precession is applied. Otherwise, the overall WRMS error of the residuals would increase by 35%. Thus, the effect of the geodetic precession in the Earth rotation is confirmed experimentally. The other finding is the reliable estimation δ_c = 3.844 ± 0.028° of the phase lag δ_c of the tides in the fluid core. When processing the UT variations, a simple model of the elastic interaction between the mantle and fluid core at their common boundary made it possible to satisfactory describe the largest observed oscillations of UT with the period of 18.6 year, reducing the WRMS error of the UT residuals to the value 0.18 ms (after removing the secular, annual and semi-annual terms).     

The First DIRECT Distance to a Detached Eclipsing Binary in M33

We present the first direct distance determination to a detached eclipsing binary in M33, which was found by the DIRECT Project. Located in the OB 66 association, it was one of the most suitable detached eclipsing binaries found by DIRECT for distance determination, given its 4.8938 day period. We obtained follow-up BV photometry and spectroscopy from which we determined the parameters of the system. It contains two O7 main sequence stars with masses of and and radii of and , respectively. We derive temperatures of K and K and determine the reddening . Using HST photometry for flux calibration in the V band, we obtain a preliminary distance modulus of mag ( kpc). The photometry and thus distance is subject to revision in the final paper.     

Dynamical dark energy from extra dimensions

We consider multidimensional cosmological model with a higher-dimensional product manifold M = R × × /Γ, where is d _o-dimensional Ricci-flat external (our) space and /Γ is d ₁-dimensional compact hyperbolic internal space. M2-brane solution for this model has the stage of accelerating expansion of the external space. We apply this model to explain the late time acceleration of our Universe. Recent observational data (the Hubble parameter at the present time and the redshift when the deceleration parameter changes its sign) fix fully all free parameters of the model. As a result, we find that considered model has too big size of the internal space at the present time and variation of the effective four-dimensional fine structure constant strongly exceeds the observational limits. The article is published in the original.     

Rotational evolution of exoplanets under the action of gravitational and magnetic perturbations

We investigate the evolution of the rotational axes of exoplanets under the action of gravitational and magnetic perturbations. The planet is assumed to be dynamically symmetrical and to be magnetised along its dynamical-symmetry axis. By qualitative methods of the bifurcation theory of multiparametric PDEs, we have derived a gallery of 69 phase portraits. The portraits illustrate evolutionary trajectories of the angular momentum of a planet for a variety of the initial conditions, for different values of the ratio between parameters describing gravitational and magnetic perturbations, and for different rates of the orbital evolution. We provide examples of the phase portraits, that reveal the differences in topology and the evolutionary track of in the vicinity of an equilibrium state. We determine the bifurcation properties, i.e., the way of reorganisation of phase trajectories in the vicinities of equilibria; and we point out the combinations of parameters’ values that permit ip-overs from a prograde to a retrograde spin mode.     

Motion On the Sphere: Integrability and Families of Orbits

We study the problem of the motion of a unit mass on the unit sphere and examine the relation between integrability and certain monoparametric families of orbits. In particular we show that if the potential is compatible with a family of meridians, it is integrable with an integral linear in the velocities, while a family of parallels guarantees integrability with an integral quadratic in the velocities.     

Analytical theory of a lunar artificial satellite with third body perturbations

We present here the first numerical results of our analytical theory of an artificial satellite of the Moon. The perturbation method used is the Lie Transform for averaging the Hamiltonian of the problem, in canonical variables: short-period terms (linked to l, the mean anomaly) are eliminated first. We achieved a quite complete averaged model with the main four perturbations, which are: the synchronous rotation of the Moon (rate ), the oblateness J ₂ of the Moon, the triaxiality C ₂₂ of the Moon ( ) and the major third body effect of the Earth (ELP2000). The solution is developed in powers of small factors linked to these perturbations up to second-order; the initial perturbations being sorted ( is first-order while the others are second-order). The results are obtained in a closed form, without any series developments in eccentricity nor inclination, so the solution apply for a wide range of values. Numerical integrations are performed in order to validate our analytical theory. The effect of each perturbation is presented progressively and separately as far as possible, in order to achieve a better understanding of the underlying mechanisms. We also highlight the important fact that it is necessary to adapt the initial conditions from averaged to osculating values in order to validate our averaged model dedicated to mission analysis purposes.     

Periodic solutions for resonance 4/3: Application to the construction of an intermediary orbit for Hyperion’s motion

The motion of Hyperion is an almost perfect application of second kind and second genius orbit, according to Poincaré’s classification. In order to construct such an orbit, we suppose that Titan’s motion is an elliptical one and that the observed frequencies are such that 4n _H−3n _T+3n _ω=0, where n _H, n _T are the mean motions of Hyperion and Titan, n _ω is the rate of rotation of Hyperion’s pericenter. We admit that the observed motion of Hyperion is a periodic motion such as . Then, .N _H, N _T, k∈ N ⁺. With that hypothesis we show that Hyperion’s orbit tends to a particular periodic solution among the periodic solutions of the Keplerian problem, when Titan’s mass tends to zero. The condition of periodicity allows us to construct this orbit which represents the real motion with a very good approximation. This revised version was published online in July 2006 with corrections to the Cover Date.     

Tidal effects in the earth–moon system and the earth's rotation

Analytical expressions for tidal torques induced by a tide‐arising planet which perturbs rotation of a nonrigid body are derived. Corresponding expressions both for secular and periodic perturbations of the Euler's angles are given for the case of the earth's rotation. Centennial secular rates of the nutation angle θ and of the earth's angular velocity ω, as well as the centennial logarithmic decrement ν of the Chandler wobble are evaluated:  mas, . In the Universal Time (UT) a large out‐of‐phase (sine) dissipative term with the period 18.6 years and the amplitude 2.3 ms is found. Corrections to nutation coefficients, which presumably have not been taken into account in IAU theory, are given. This revised version was published online in July 2006 with corrections to the Cover Date.     

Secular frequencies of 3-D exoplanetary systems

Using a 12th order expansion of the perturbative potential in powers of the eccentricities and the inclinations, we study the secular effects of two non-coplanar planets which are not in mean–motion resonance. By means of Lie transformations (which introduce an action–angle formulation of the Hamiltonian), we find the four fundamental frequencies of the 3-D secular three-body problem and compute the long-term time evolutions of the Keplerian elements. To find the relations between these elements, the main combinations of the fundamental frequencies common to these evolutions are identified by frequency analysis. This study is performed for two different reference frames: a general one and the Laplace plane. We underline the known limitations of the linear Laplace–Lagrange theory and point out the great sensitivity of the 3-D secular three-body problem to its initial values. This analytical approach is applied to the exoplanetary system Andromedae in order to search whether the eccentricities evolutions and the apsidal configuration (libration of ) observed in the coplanar case are maintained for increasing initial values of the mutual inclination of the two orbital planes. Anne-Sophie Libert is FNRS Research Fellow.     

The large-number hypothesis and the Earth's expansion

The equation of state of the terrestrial material obtained from seismic data is adopted to construct three zone earth models under hypothesis of variable constant of gravityG as proposed by Dirac. Three hypotheses are investigated: variableG without creation, creation such thatm (mass) G ^–1, and multiplicative creation,mG ^–2. It is shown that, with the currently accepted value of the Hubble constant, , and for each hypothesis. On the multiplicative creation, the Earth radius would have been 5100 km, which is in agreement with estimate by some geophysicists.     

Der offene sternhaufen NGC 2362

The very young open star cluster NGC 2362 was investigated by the strip method on charts of two photographs taken with the 1-m Schmidt telescope of the European Southern Observatory. Up to the limiting magnitudeM _v ^* =5 _. ^m 8 the cluster contains 100 stars and can be described by the Gaussian density law (6). Further results are: Mass = 246 , central mass density ₀ = 43.1 = 246 pc^-3 , radiusR2.6 pc, mean velocity of the stars = 0.64 km s^–1.

---

---

Mitteilungen Serie A.     

On the structure and mass of the galactic halo

On the basis of a globular cluster study a crude estimate of the total mass of the galactic halo within 20 kpc from the centre is done. It gives a minimal halo mass of the order of , yielding possibilities for a mass as large as . The content of the interstellar matter in the halo is estimated too. It is found that the gas content is a few percents the minimal mass, the gas temperature is very high — about 1×10⁶ K, the magnetic field weak — about 0.25 nT. A weak nonthermal radio emission might be expected from such a halo.     

Solution of Lambert's problem for short arcs

Approximation formulas are found for and , wherex(t) satisfies ,x(0)=x ₀,x(1)=x ₁. The results are applied to an example of two-body motion.     

Two Sets of Improved Approximate Expressions of the Gyrosynchrotron Radiation

In this paper two sets of improved approximate expressions of emissivity , absorptivity , effective temperature T_eff, and frequency of peak brightness _p of gyrosynchrotron radiation are presented respectively for the ranges from 5 to 10 and 10 to 100 of harmonic numbers s(= /_B). The expressions are designed for the range from 20° to 80° of viewing angle , and the range 2 to 7 of electron energy spectral index . They are expressed by a power-law function in which the indexes are fitted by polynomial expressions of . Their statistical errors are, respectively, 24% and 32% for and for and 28% for . Their accuracies are much better than those of linear fitting of the power-law index.     

Charged particle diffusion in the presence of Alfvén waves: The perpendicular particle flux

The diffusion of charged particles through a weak stochastic electro-magnetic field which is superimposed on a constant background magnetic field is considered. The stochastic electromagnetic fields are assumed to consist of unpolarized Alfvén waves propagating at arbitrary angles to the direction . When the Alfvén waves are propagating in directions other than and the particle gyro-radius,r _g, is sufficiently large (but may be smaller than the correlation length of the stochastic fields) it is shown that the particle flux perpendicular to the direction is , wherev is the particle speed andf the particle density. The expression forK differs from those calculated by previous authors. For small particle gyro-radii the flux S has a different functional form and is identical to that found by Urch (1977) to describe particle diffusion when the Alfvén waves only propagate in the direction .