Sławomir Breiter

The long-term stability of extrasolar system HD 37124. Numerical study of resonance effects

We describe numerical tools for the stability analysis of extrasolar planetary systems. In par- ticular, we consider the relative Poincare variables and symplectic integration of the equations of motion. We apply the tangent map to derive a numerically efficient algorithm of the fast indicator Mean Exponential Growth factor of Nearby Orbits (MEGNO), a measure of the maximal Lyapunov exponent, that helps to distinguish chaotic and regular configurations. The results concerning the three-planet extrasolar system HD 37124 are presented and discussed. The best-fitting solutions found in earlier works are studied more closely. The system involves Jovian planets with similar masses. The orbits have moderate eccentricities, nevertheless the best-fitting solutions are found in dynamically active region of the phase space. The long-term stability of the system is determined by a net of low-order two-body and three-body mean mo- tion resonances. In particular, the three-body resonances may induce strong chaos that leads to self-destruction of the system after Myr of apparently stable and bounded evolution. In such a case, numerically efficient dynamical maps are useful to resolve the fine structure of the phase space and to identify the sources of unstable behaviour.

The YORP effect on 25143 Itokawa

Context. The asteroid 25143 Itokawa is one of the candidates for the detection of the Yarkovsky-O’Keefe-Radzievskii-Paddack (YORP) effect in the rotation period. Previous studies were carried out up to the 196 608 facets triangulation model and were not able to provide a good theoretical estimate of this effect, raising questions about the influence of the mesh resolution and the centre of mass location on the evolution the rotation period. Aims. The YORP effect on Itokawa is computed for different topography models up to the highest resolution Gaskell mesh of 3 145 728 triangular faces in an attempt to find the best possible YORP estimate. Other, lower resolution models are also studied and the question of the dependence of the rotation period drift on the density distribution inhomogeneities is reexamined. A comparison is made with 433 Eros models possessing a similar resolution. Methods. The Rubincam approximation (zero conductivity) is assumed in the numerical simulation of the YORP effect in rotation period. The mean thermal radiation torques are summed over triangular facets assuming Keplerian heliocentric motion and uniform rotation around a body-fixed axis. Results. There is no evidence of YORP convergence in Gaskell model family. Differently simplified meshes may converge quickly to their parent models, but this does not prove the quality of YORP computed from the latter. We confirm the high sensitivity of the YORP effect to the fine details of the surface for 25 143 Itokawa and 433 Eros. The sensitivity of the Itokawa YORP to the centre of mass shift is weaker than in earlier works, but instead the results prove to be sensitive to the spin axis orientation in the body frame. Conclusions. Either the sensitivity of the YORP effect is a physical phenomenon and all present predictions are questionable, or the present thermal models are too simplified.

Analysis of the rotation period of asteroids (1865) Cerberus, (2100) Ra-Shalom, and (3103) Eger - search for the YORP effect

Context. The spin state of small asteroids can change on a long timescale by the Yarkovsky-O’Keefe-Radzievskii-Paddack (YORP) effect, the net torque that arises from anisotropically scattered sunlight and proper thermal radiation from an irregularly-shaped asteroid. The secular change in the rotation period caused by the YORP effect can be detected by analysis of asteroid photometric lightcurves. Aims. We analyzed photometric lightcurves of near-Earth asteroids (1865) Cerberus, (2100) Ra-Shalom, and (3103) Eger with the aim to detect possible deviations from the constant rotation caused by the YORP effect. Methods. We carried out new photometric observations of the three asteroids, combined the new lightcurves with archived data, and used the lightcurve inversion method to model the asteroid shape, pole direction, and rotation rate. The YORP effect was modeled as a linear change in the rotation rate in time dω/dt .V alues of dω/dt derived from observations were compared with the values predicted by theory. Results. We derived physical models for all three asteroids. We had to model Eger as a nonconvex body because the convex model failed to fit the lightcurves observed at high phase angles. We probably detected the acceleration of the rotation rate of Eger dω/dt = (1.4 ± 0.6) × 10 −8 rad d −2 (3σ error), which corresponds to a decrease in the rotation period by 4. 2m s yr −1 . The photometry of Cerberus and Ra-Shalom was consistent with a constant-period model, and no secular change in the spin rate was detected. We could only constrain maximum values of |dω/dt| < 8 × 10 −9 rad d −2 for Cerberus, and |dω/dt| < 3 × 10 −8 rad d −2 for Ra-Shalom.

LUNISOLAR RESONANCES REVISITED

Lunisolar resonances arise in the artificial satellite problem without short-periodic terms. The basic model including the Earths J2and a Hill-type model for the Sun or the Moon admits 20 different periodic terms which may lead to a resonance involving the satellites perigee, node and the longitude of the perturbing body. Some of the resonances have been studied separately since 1960s. The present paper reviews all single resonances, attaching an appropriate fundamental model to each case. Only a part of resonances match known fundamental models. An extended fundamental model is proposed to account for some complicated phenomena. Most of the double resonance cases still remain unexplored.

Efficient Lie-Poisson Integrator for Secular Spin Dynamics of Rigid Bodies

A fast and efficient numerical integration algorithm is presented for the problem of the secular evolution of the spin axis. Under the assumption that a celestial body rotates around its maximum moment of inertia, the equations of motion are reduced to the Hamiltonian form with a Lie-Poisson bracket. The integration method is based on the splitting of the Hamiltonian function, and so it conserves the Lie-Poisson structure. Two alternative partitions of the Hamiltonian are investigated, and second-order leapfrog integrators are provided for both cases. Non-Hamiltonian torques can be incorporated into the integrators with a combination of Euler and Lie-Euler approximations. Numerical tests of the methods confirm their useful properties of short computation time and reliability on long integration intervals.

Stress field and spin axis relaxation for inelastic triaxial ellipsoids

A compact formula for the stress tensor inside a self-gravitating, triaxial ellipsoid in an arbitrary rotation state is given. It contains no singularity in the incompressible medium limit. The stress tensor and the quality factor model are used to derive a solution for the energy dissipation resulting in the damping (short-axis mode) or excitation (long axis) of wobbling. In the limit of an ellipsoid of revolution, we compare our solution with earlier ones and show that, with appropriate corrections, the differences in damping times estimates are much smaller than it has been claimed.

YORP torques with 1D thermal model

A numerical model of the Yarkovsky–O’Keefe–Radzievskii–Paddack (YORP) effect for objects defined in terms of a triangular mesh is described. The algorithm requires that each surface triangle can be handled independently, which implies the use of a 1D thermal model. Insolation of each triangle is determined by an optimized ray–triangle intersection search. Surface temperature is modelled with a spectral approach; imposing a quasi-periodic solution we replace heat conduction equation by the Helmholtz equation. Non-linear boundary conditions are handled by an iterative, fast Fourier transform based solver. The results resolve the question of the YORP effect in rotation rate independence on conductivity within the non-linear 1D thermal model regardless of the accuracy issues and homogeneity assumptions. A seasonal YORP effect in attitude is revealed for objects moving on elliptic orbits when a non-linear thermal model is used.

Tesseral harmonic perturbations in radial transverse and binormal components

The formulae for the perturbations in radial, transverse and binormal components of the Earth artificial satellite motion have been derived. Perturbations due to the tesseral part of the geopotential are considered. The geopotential expressed in terms of the orbital elements has the form proposed by Wnuk (1988). The formulae for the perturbations have been obtained using the Hori (1966) method. They can be effectively applied in calculation of the perturbations in the components including the coefficients of the high order and degree tesseral harmonics. The derived formulae reveal no singularities at zero eccentricity.

The Prograde C7 Resonance for Earth and Mars Satellite Orbits

The resonance C7 is a 1:1 eccentricity (apsidal) resonance between the longitude of a satellites pericentre and the mean longitude of the Sun. A previous paper by the author (Breiter, 1999) identified it as the strongest of the lunisolar apsidal resonances. After the reduction to a single degree of freedom, the problem is studied qualitatively for the prograde orbits around the Earth and Mars. Pitchfork, saddle-node, and saddle connection bifurcations give rise to a complicated phase flow, which may involve up to nine critical points.

Explicit Symplectic Integrator for Highly Eccentric Orbits

An explicit symplectic integrator is constructed for perturbed elliptic orbits of an arbitrary eccentricity. The perturbation should be Hamiltonian, but it may depend on time explicitly. The main feature of the integrator is the use of KS variables in the ten-dimensional extended phase space. As an example of its application the motion of an Earth satellite under the action of the planets oblateness and of lunar perturbations is studied. The results confirm the superiority of the method over a classical Wisdom–Holman algorithm in both accuracy and computation time.