-
Galactic and stellar perturbations of long period comet motion -- Practical considerations
Authors:
Piotr A. Dybczyński,
Sławomir Breiter
Abstract:
Context. Due to our increasing knowledge on the Galactic and stellar neighborhood of the Solar System, modern long-period comet motion studies have to take into account both stellar perturbations and the overall Galactic potential. Aims. Our aim is to propose algorithms and methods to perform numerical integration of a Solar System small body equations of motion much faster and at the same time wi…
▽ More
Context. Due to our increasing knowledge on the Galactic and stellar neighborhood of the Solar System, modern long-period comet motion studies have to take into account both stellar perturbations and the overall Galactic potential. Aims. Our aim is to propose algorithms and methods to perform numerical integration of a Solar System small body equations of motion much faster and at the same time with greater precision. Methods. We propose a new formulation of the equations of motion formulated in the Solar System barycentric frame but accurately accounting for the differential perturbations caused by the Galactic potential. To use these equations effectively we provide numerical ephemerides of the Galactic positions of the Sun and a set of potential stellar perturbers. Results. The proposed methods offer the precision higher by several orders of magnitude and simultaneously greatly reduce the necessary CPU time. The application of this approach is presented with the example of a detailed dynamical study of the past motion of comet C/2015 XY1.
△ Less
Submitted 19 October, 2021; v1 submitted 15 September, 2021;
originally announced September 2021.
-
Analytical solution of the Colombo top problem
Authors:
J. Haponiak,
S. Breiter,
D. Vokrouhlicky
Abstract:
The Colombo top is a basic model in the rotation dynamics of a celestial body moving on a precessing orbit and perturbed by a gravitational torque. The paper presents a detailed study of analytical solution to this problem. By solving algebraic equations of degree 4, we provide the expressions for the extreme points of trajectories as functions of their energy. The location of stationary points (k…
▽ More
The Colombo top is a basic model in the rotation dynamics of a celestial body moving on a precessing orbit and perturbed by a gravitational torque. The paper presents a detailed study of analytical solution to this problem. By solving algebraic equations of degree 4, we provide the expressions for the extreme points of trajectories as functions of their energy. The location of stationary points (known as the Cassini states) is found as the function of the two parameters of the problem. Analytical solution in terms the Weierstrass and the Jacobi elliptic functions is given for regular trajectories. Some trajectories are expressible through elementary functions: not only the homoclinic orbits, as expected, but also a special periodic solution whose energy is equal to that of the first Cassini state (unnoticed in previous studies).
△ Less
Submitted 31 March, 2020;
originally announced March 2020.
-
Analysis of the rotation period of asteroids (1865) Cerberus, (2100) Ra-Shalom, and (3103) Eger - search for the YORP effect
Authors:
J. Durech,
D. Vokrouhlicky,
A. R. Baransky,
S. Breiter,
O. A. Burkhonov,
W. Cooney,
V. Fuller,
N. M. Gaftonyuk,
J. Gross,
R. Ya. Inasaridze,
M. Kaasalainen,
Yu. N. Krugly,
O. I. Kvaratshelia,
E. A. Litvinenko,
B. Macomber,
F. Marchis,
I. E. Molotov,
J. Oey,
D. Polishook,
J. Pollock,
P. Pravec,
K. Sarneczky,
V. G. Shevchenko,
I. Slyusarev,
R. Stephens
, et al. (6 additional authors not shown)
Abstract:
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. We analyz…
▽ More
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. 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. 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. Values of dω/ dt derived from observations were compared with the values predicted by theory. 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) x 10^{-8} rad/d (3σerror), which corresponds to a decrease in the rotation period by 4.2 ms/yr. 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 x 10^{-9} rad/d for Cerberus, and |dω/ dt| < 3 x 10^{-8} rad/d for Ra-Shalom.
△ Less
Submitted 8 October, 2012;
originally announced October 2012.
-
Stress field and spin axis relaxation for inelastic triaxial ellipsoids
Authors:
Slawomir Breiter,
Agata Rozek,
David Vokrouhlicky
Abstract:
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 ellip…
▽ More
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.
This version implements corrections of misprints found in the MNRAS published text.
△ Less
Submitted 19 April, 2013; v1 submitted 21 August, 2012;
originally announced August 2012.
-
YORP effect with anisotropic radiation
Authors:
Slawomir Breiter,
David Vokrouhlicky
Abstract:
The influence of optical scattering and thermal radiation models an the Yarkovsky-O'Keefe-Radzievskii-Paddack (YORP) effect is studied. Lambertian formulation is compared with Hapke scattering and emission laws and Lommel-Seeliger reflection. Although the form of reflectivity function strongly influences mean torques due to scattering or thermal radiation alone, their combined contribution to the…
▽ More
The influence of optical scattering and thermal radiation models an the Yarkovsky-O'Keefe-Radzievskii-Paddack (YORP) effect is studied. Lambertian formulation is compared with Hapke scattering and emission laws and Lommel-Seeliger reflection. Although the form of reflectivity function strongly influences mean torques due to scattering or thermal radiation alone, their combined contribution to the rotation period YORP is not much different from the standard Lambertian values. For higher albedo values the differences between the Hapke and Lambert models become significant for the YORP in attitude.
△ Less
Submitted 8 September, 2010;
originally announced September 2010.
-
YORP torques with 1D thermal model
Authors:
Slawomir Breiter,
Przemyslaw Bartczak,
Maria Czekaj
Abstract:
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 modeled with a spec…
▽ More
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 modeled with a spectral approach; imposing a quasi-periodic solution we replace heat conduction equation by the Helmholtz equation. Nonlinear boundary conditions are handled by an iterative, FFT based solver. The results resolve the question of the YORP effect in rotation rate independence on conductivity within the nonlinear 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 nonlinear thermal model is used.
△ Less
Submitted 14 July, 2010; v1 submitted 21 June, 2010;
originally announced June 2010.
-
YORP torque as the function of shape harmonics
Authors:
S. Breiter,
H. Michalska
Abstract:
The second order analytical approximation of the mean YORP torque components is given as an explicit function of the shape spherical harmonics coefficients for a sufficiently regular minor body. The results are based upon a new expression for the insolation function, significantly simpler than in previous works. Linearized plane parallel model of the temperature distribution derived from the ins…
▽ More
The second order analytical approximation of the mean YORP torque components is given as an explicit function of the shape spherical harmonics coefficients for a sufficiently regular minor body. The results are based upon a new expression for the insolation function, significantly simpler than in previous works. Linearized plane parallel model of the temperature distribution derived from the insolation function allows to take into account a nonzero conductivity. Final expressions for the three average components of the YORP torque related with rotation period, obliquity, and precession are given in a form of Legendre series of the cosine of obliquity. The series have good numerical properties and can be easily truncated according to the degree of Legendre polynomials or associated functions, with first two terms playing the principal role. The present version fixes the errors discovered in the text that appeared in Monthly Notices RAS (388, pp. 297-944).
△ Less
Submitted 14 May, 2009; v1 submitted 16 May, 2008;
originally announced May 2008.
-
The long-term stability of extrasolar system HD 37124. Numerical study of resonance effects
Authors:
Krzysztof Gozdziewski,
Slawomir Breiter,
Wojciech Borczyk
Abstract:
We describe numerical tools for the stability analysis of extrasolar planetary systems. In particular, 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 MEGNO (a measure of the maximal Lyapunov exponent) that helps to distinguish chaotic and regular configur…
▽ More
We describe numerical tools for the stability analysis of extrasolar planetary systems. In particular, 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 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 fit 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 fit 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 motion resonances. In particular, the three-body resonances may induce strong chaos that leads to self-destruction of the system after Myrs 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 behavior.
△ Less
Submitted 29 October, 2007;
originally announced October 2007.