Source author record

Christoph Lhotka

Christoph Lhotka appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

13works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

13 published item(s)

preprint2020arXiv

Asteroid migration due to the Yarkovsky effect and the distribution of the Eos family

Based on a linearized model of the Yarkovsky effect, we investigate in this paper the dependence of the semimajor axis drift $Δa$ of a celestial body on its size, spinning obliquity, initial orbit and thermal parameters on its surface. With appropriate simplification and approximation, we obtain the analytical solutions to the perturbation equations for the motion of asteroids influenced by the Yarkovsky effect, and they are then verified by numerical simulations of the full equations of motion. These solutions present explicitly the dependencies of $Δa$ on the thermal and dynamical parameters of the asteroid. With these analytical formulae for $Δa$, we investigate the combined seasonal and diurnal Yarkovsky effects. The critical points where the migration direction reverses are calculated and the consequent selective effects according to the size and rotation state of asteroids are discussed. %Solely the Yarkovsky effect is found to be able to produce some ring structure in the aged circumstellar debris disk. Finally, we apply the analytical formulae to calculate the migration of Eos family members. The space distribution of asteroids is well reproduced. Our calculations suggest that statistically the orientations of spin axes of family members satisfy a random-obliquity distribution, and the rotation rate $ω_{\rm rot}$ of asteroid depends on its size $R$ by $ω_{\rm rot}\propto R^{-1}$.

preprint2020arXiv

Dynamics of charged dust in the orbit of Venus

We study the dynamics of co-orbital dust in the inner solar system, i.e. the role of the solar radiation pressure, Poynting-Robertson effect, solar wind, and the interplanetary magnetic field on the location, width and stability of resonant motion of charged, and micron sized dust grains situated in the 1:1 mean motion resonance with planet Venus. We find deviations and asymmetry between $L_4$ and $L_5$ in locations of libration centers and libration width under the influence of non-gravitational effects via both analytical and numerical methods. The triangular Lagrangian points become unstable once we take into consideration solar radiation pressure, the Poynting-Robertson effect and solar wind drag. The Lorentz force could further destabilize the orbits, especially for small dust particles. We also make a comparison between the circular, elliptic restricted three-body model and a more complete model including all planets.

preprint2019arXiv

Resonances in the Earth's Space Environment

We study the presence of resonances in the region of space around the Earth. We consider a massless body (e.g, a dust particle or a small space debris) subject to different forces: the gravitational attraction of the geopotential, the effects of Sun and Moon. We distinguish different types of resonances: tesseral resonances are due to a commensurability involving the revolution of the particle and the rotation of the Earth, semi-secular resonances include the rates of variation of the mean anomalies of Moon and Sun, while secular resonances just depend on the rates of variation of the arguments of perigee and the longitudes of the ascending nodes of the perturbing bodies. We characterize such resonances, giving precise statements on the regions where the resonances can be found and provide examples of some specific commensurability relations.

preprint2016arXiv

Charged dust grain dynamics subject to solar wind, Poynting-Robertson drag, and the interplanetary magnetic field

We investigate the combined effect of solar wind, Poynting-Robertson drag, and the frozen-in interplanetary magnetic field on the motion of charged dust grains in our solar system. For this reason we derive a secular theory of motion by the means of averaging method and validate it with numerical simulations of the un-averaged equations of motions. The theory predicts that the secular motion of charged particles is mainly affected by the z-component of the solar magnetic axis, or the normal component of the interplanetary magnetic field. The normal component of the interplanetary magnetic field leads to an increase or decrease of semi-major axis depending on its functional form and sign of charge of the dust grain. It is generally accepted that the combined effects of solar wind and photon absorption and re-emmision (Poynting-Robertson drag) lead to a decrease in semi-major axis on secular time scales. On the contrary, we demonstrate that the interplanetary magnetic field may counteract these drag forces under certain circumstances. We derive a simple relation between the parameters of the magnetic field, the physical properties of the dust grain as well as the shape and orientation of the orbital ellipse of the particle, which is a necessary conditions for the stabilization in semi-major axis.

preprint2016arXiv

Gravity field and solar component of the precession rate and nutation coefficients of Comet 67P Churyumov-Gerasimenko

The aim of this study is first to determine the gravity field of the comet 67P Churyumov-Gerasimenko and second to derive the solar component of the precession rate and nutation coefficients of the spin axis of the comet nucleus, i.e. without the direct, usually larger, effect of outgassing. The gravity field, and related moments of inertia, are obtained from two polyhedra, that are provided by the OSIRIS and NAV-CAM experiments on Rosetta, and are based on the assumption of uniform density for the comet nucleus. We also calculate the forced precession rate as well as the nutation coefficients on the basis of Kinoshita's theory of rotation of the rigid Earth and adapted it to be able to indirectly include the effect of outgassing on the rotational parameters. The 2nd degree denormalized Stokes coefficients of comet 67P C-G turn out to be (bracketed numbers refer to second shape model) C20 = -6.74 [-7.93] * 10^-2, C22 = 2.60 [2.71] * 10^-2 consistent with normalized principal moments of inertia A/MR^2 = 0.13 [0.11], B/MR^2 = 0.23 [0.22], with polar moment c = C/MR^2 = 0.25, depending on the choice of the polyhedron model. The obliquity between the rotation axis and the mean orbit normal is 52 degree, and the precession rate only due to solar torques becomes 20-30 "/y. Oscillations in longitude caused by the gravitational pull of the Sun turn out to be of the order of 1', oscillations in obliquity can be estimated to be of the order of 0.5'.

preprint2016arXiv

Poynting-Robertson drag and solar wind in the space debris problem

We analyze the combined effect of Poynting-Robertson and solar wind drag on space debris. We derive a model within Cartesian, Gaussian and Hamiltonian frameworks. We focus on the geosynchronous resonance, although the results can be easily generalized to any resonance. By numerical and analytical techniques, we compute the drift in semi-major axis due to Poynting-Robertson and solar wind drag. After a linear stability analysis of the equilibria, we combine a careful investigation of the regular, resonant, chaotic behavior of the phase space with a long-term propagation of a sample of initial conditions. The results strongly depend on the value of the area-to-mass ratio of the debris, which might show different dynamical behaviors: temporary capture or escape from the geosynchronous resonance, as well as temporary capture or escape from secondary resonances involving the rate of variation of the longitude of the Sun. Such analysis shows that Poynting-Robertson and solar wind drag must be taken into account, when looking at the long-term behavior of space debris. Trapping or escape from the resonance can be used to place the debris in convenient regions of the phase space.

preprint2015arXiv

Effective stability around the Cassini state in the spin-orbit problem

We investigate the long-time stability in the neighborhood of the Cassini state in the conservative spin-orbit problem. Starting with an expansion of the Hamiltonian in the canonical Andoyer-Delaunay variables, we construct a high-order Birkhoff normal form and give an estimate of the effective stability time in the Nekhoroshev sense. By extensively using algebraic manipulations on a computer, we explicitly apply our method to the rotation of Titan. We obtain physical bounds of Titan's latitudinal and longitudinal librations, finding a stability time greatly exceeding the estimated age of the Universe. In addition, we study the dependence of the effective stability time on three relevant physical parameters: the orbital inclination, $i$, the mean precession of the ascending node of Titan orbit, $\dotΩ$, and the polar moment of inertia, $C$.

preprint2014arXiv

The effect of Poynting-Robertson drag on the triangular Lagrangian points

We investigate the stability of motion close to the Lagrangian equilibrium points L4 and L5 in the framework of the spatial, elliptic, restricted three- body problem, subject to the radial component of Poynting-Robertson drag. For this reason we develop a simplified resonant model, that is based on averaging theory, i.e. averaged over the mean anomaly of the perturbing planet. We find temporary stability of particles displaying a tadpole motion in the 1:1 resonance. From the linear stability study of the averaged simplified resonant model, we find that the time of temporary stability is proportional to beta a1 n1 , where beta is the ratio of the solar radiation over the gravitational force, and a1, n1 are the semi-major axis and the mean motion of the perturbing planet, respectively. We extend previous results (Murray (1994)) on the asymmetry of the stability indices of L4 and L5 to a more realistic force model. Our analytical results are supported by means of numerical simulations. We implement our study to Jupiter-like perturbing planets, that are also found in extra-solar planetary systems.

preprint2013arXiv

The influence of orbital dynamics, shape and tides on the obliquity of Mercury

Earth-based radar observations of the rotational dynamics of Mercury (Margot et al. 2012) combined with the determination of its gravity field by MESSENGER (Smith et al. 2012) give clues on the internal structure of Mercury, in particular its polar moment of inertia C, deduced from the obliquity (2.04 +/- 0.08) arcmin. The dynamics of the obliquity of Mercury is a very-long term motion (a few hundreds of kyrs), based on the regressional motion of Mercury's orbital ascending node. This paper, following the study of Noyelles & D'Hoedt (2012), aims at first giving initial conditions at any time and for any values of the internal structure parameters for numerical simulations, and at using them to estimate the influence of usually neglected parameters on the obliquity, like J3, the Love number k2 and the secular variations of the orbital elements. We use for that averaged representations of the orbital and rotational motions of Mercury, suitable for long-term studies. We find that J3 should alter the obliquity by 250 milli-arcsec, the tides by 100 milli-arcsec, and the secular variations of the orbital elements by 10 milli-arcsec over 20 years. The resulting value of C could be at the most changed from 0.346mR^2 to 0.345mR^2.

preprint2012arXiv

Stability estimates of nearly-integrable systems with dissipation and non-resonant frequency

We consider a dissipative vector field which is represented by a nearly-integrable Hamiltonian flow to which a non symplectic force is added, so that the phase space volume is not preserved. The vector field depends upon two parameters, namely the perturbing and dissipative parameters, and by a drift function. We study the general case of an l-dimensional, time-dependent vector field. Assuming to start with non-resonant initial conditions, we prove the stability of the variables which are actions of the conservative system (namely, when the dissipative parameter is set to zero) for exponentially long times. In order to construct the normal form, a suitable choice of the drift function must be performed. We also provide some simple examples in which we construct explicitly the normal form, we make a comparison with a numerical integration and we compute theoretical bounds on the parameters as well as we give explicit stability estimates.

preprint2012arXiv

Stability of nearly-integrable systems with dissipation

We study the stability of a vector field associated to a nearly-integrable Hamiltonian dynamical system to which a dissipation is added. Such a system is governed by two parameters, named the perturbing and dissipative parameters, and it depends on a drift function. Assuming that the frequency of motion satisfies some resonance assumption, we investigate the stability of the dynamics, and precisely the variation of the action variables associated to the conservative model. According to the structure of the vector field, one can find linear and exponential stability times, which are established under smallness con- ditions on the parameters. We also provide some applications to concrete examples, which exhibit a linear or exponential stability behavior.

preprint2012arXiv

The orbit of 2010 TK7. Possible regions of stability for other Earth Trojan asteroids

Recently the first Earth Trojan has been observed (Mainzer et al., ApJ 731) and found to be on an interesting orbit close to the Lagrange point L4 (Connors et al., Nature 475). In the present study we therefore perform a detailed investigation on the stability of its orbit and moreover extend the study to give an idea of the probability to find additional Earth-Trojans. Our results are derived using different approaches: a) we derive an analytical mapping in the spatial elliptic restricted three-body problem to find the phase space structure of the dynamical problem. We explore the stability of the asteroid in the context of the phase space geometry, including the indirect influence of the additional planets of our Solar system. b) We use precise numerical methods to integrate the orbit forward and backward in time in different dynamical models. Based on a set of 400 clone orbits we derive the probability of capture and escape of the Earth Trojan asteroids 2010 TK7. c) To this end we perform an extensive numerical investigation of the stability region of the Earth's Lagrangian points. We present a detailed parameter study in the regime of possible stable tadpole and horseshoe orbits of additional Earth-Trojans, i.e. with respect to the semi-major axes and inclinations of thousands of fictitious Trojans. All three approaches underline that the Earth Trojan asteroid 2010 TK7 finds himself in an unstable region on the edge of a stable zone; additional Earth-Trojan asteroids may be found in this regime of stability.

preprint2009arXiv

An Overview of the 13:8 Mean Motion Resonance between Venus and Earth

It is known since the seminal study of Laskar (1989) that the inner planetary system is chaotic with respect to its orbits and even escapes are not impossible, although in time scales of billions of years. The aim of this investigation is to locate the orbits of Venus and Earth in phase space, respectively to see how close their orbits are to chaotic motion which would lead to unstable orbits for the inner planets on much shorter time scales. Therefore we did numerical experiments in different dynamical models with different initial conditions -- on one hand the couple Venus-Earth was set close to different mean motion resonances (MMR), and on the other hand Venus' orbital eccentricity (or inclination) was set to values as large as e = 0.36 (i = 40deg). The couple Venus-Earth is almost exactly in the 13:8 mean motion resonance. The stronger acting 8:5 MMR inside, and the 5:3 MMR outside the 13:8 resonance are within a small shift in the Earth's semimajor axis (only 1.5 percent). Especially Mercury is strongly affected by relatively small changes in eccentricity and/or inclination of Venus in these resonances. Even escapes for the innermost planet are possible which may happen quite rapidly.