Source author record

Christos Efthymiopoulos

Christos Efthymiopoulos 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

19works
8topics
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

19 published item(s)

preprint2022arXiv

A detailed dynamical model for inclination-only dependent lunisolar resonances. Effect on the "eccentricity growth" mechanism

The focus of this paper is on inclination-only dependent lunisolar resonances, which shape the dynamics of a MEO (Medium Earth Orbit) object over secular time scales (i.e. several decades). Following the formalism of arXiv:2107.14507, we discuss an analytical model yielding the correct form of the separatrices of each one of the major lunisolar resonances in the "action" space $(i, e)$ (inclination, eccentricity) for any given semi-major axis $a$. We then highlight how our method is able to predict and explain the main structures found numerically in Fast Lyapunov Indicator (FLI) cartography. We focus on explaining the dependence of the FLI maps from the initial phase of the argument of perigee $ω$ and of the longitude of the ascending node $Ω$ of the object and of the moon $Ω_L$. In addition, on the basis of our model, we discuss the role played by the $Ω-Ω_L$ and the $2 Ω-Ω_L$ resonances, which overlap with the inclination-only dependent ones as they sweep the region for increasing values of $a$, generating large domains of chaotic motion. Our results provide a framework useful in designing low-cost satellite deployment or space debris mitigation strategies, exploiting the natural dynamics of lunisolar resonances that increase an object's eccentricity up until it reaches a domain where friction leads to atmospheric re-entry.

preprint2022arXiv

Kolmogorov algorithm for isochronous Hamiltonian systems

We present a Kolmogorov-like algorithm for the computation of a normal form in the neighborhood of an invariant torus in `isochronous' Hamiltonian systems, i.e., systems with Hamiltonians of the form $\mathcal{H}=\mathcal{H}_0+\varepsilon \mathcal{H}_1$ where $\mathcal{H}_0$ is the Hamiltonian of $N$ linear oscillators, and $\mathcal{H}_1$ is expandable as a polynomial series in the oscillators' canonical variables. This method can be regarded as a normal form analogue of a corresponding Lindstedt method for coupled oscillators. We comment on the possible use of the Lindstedt method itself under two distinct schemes, i.e., one producing series analogous to those of the Birkhoff normal form scheme, and another, analogous to the Kolomogorov normal form scheme in which we fix in advance the frequency of the torus.

preprint2021arXiv

A deep dive into the $2g+h$ resonance: separatrices, manifolds and phase space structure of navigation satellites

Despite extended past studies, several questions regarding the resonant structure of the medium-Earth orbit (MEO) region remain hitherto unanswered. This work describes in depth the effects of the $2g+h$ lunisolar resonance. In particular, (i) we compute the correct forms of the separatrices of the resonance in the inclination-eccentricity space for fixed semi-major axis. This allows to compute the change in the width of the $2g+h$ resonance as the altitude increases. (ii) We discuss the crucial role played by the value of the inclination of the Laplace plane, $i_{L}$. Since $i_L$ is comparable to the resonance's separatrix width, the parametrization of all resonance bifurcations has to be done in terms of the proper inclination $i_{p}$, instead of the mean one. (iii) The subset of circular orbits constitutes an invariant subspace embedded in the full phase space, the center manifold $\mathcal{C}$. Using $i_p$ as a label, we compute its range of values for which $\mathcal{C}$ becomes a normally hyperbolic invariant manifold (NHIM). The structure of invariant tori in $\mathcal{C}$ allows to explain the role of the initial phase $h$ noticed in several works. (iv) Through Fast Lyapunov Indicator (FLI) cartography, we portray the stable and unstable manifolds of the NHIM as the altitude increases. Manifold oscillations dominate in phase space between $a=24,000$ km and $a=30,000$ km as a result of the sweeping of the $2g+h$ resonance by the $h-Ω_{\rm{Moon}}$ and $2h-Ω_{\rm{Moon}}$ resonances. The noticeable effects of the latter are explained as a consequence of the relative inclination of the Moon's orbit with respect to the ecliptic. The role of the phases $(h,Ω_{\rm{Moon}})$ in the structures observed in the FLI maps is also clarified. Finally,(v) we discuss how the understanding of the manifold dynamics could inspire end-of-life disposal strategies.

preprint2021arXiv

Closed form perturbation theory in the restricted three-body problem without relegation

We propose a closed-form normalization method suitable for the study of the secular dynamics of small bodies in heliocentric orbits perturbed by the tidal potential of a planet with orbit external to the orbit of the small body. The method makes no use of relegation, thus, circumventing all convergence issues related to that technique. The method is based on a convenient use of a book-keeping parameter keeping simultaneously track of all the small quantities in the problem. The book-keeping affects both the Lie series and the Poisson structure employed in successive perturbative steps. In particular, it affects the definition of the normal form remainder at every normalization step. We show the results obtained by assuming Jupiter as perturbing planet and we discuss the validity and limits of the method.

preprint2021arXiv

Closed-form perturbation theory in the Sun-Jupiter restricted three body problem without relegation

We present a closed-form normalization method suitable for the study of the secular dynamics of small bodies inside the trajectory of Jupiter. The method is based on a convenient use of a book-keeping parameter introduced not only in the Lie series organization but also in the Poisson bracket structure employed in all perturbative steps. In particular, we show how the above scheme leads to a redefinition of the remainder of the normal form at every step of the formal solution of the homological equation. An application is given for the semi-analytical representation of the orbits of main-belt asteroids.

preprint2018arXiv

Semi-analytic computations of the speed of Arnold diffusion along single resonances in a priori stable Hamiltonian systems

Cornerstone models of Physics, from the semi-classical mechanics in atomic and molecular physics to planetary systems, are represented by quasi-integrable Hamiltonian systems. Since Arnold's example, the long-term diffusion in Hamiltonian systems with more than two degrees of freedom has been represented as a slow diffusion within the `Arnold web', an intricate web formed by chaotic trajectories. With modern computers it became possible to perform numerical integrations which reveal this phenomenon for moderately small perturbations. Here we provide a semi-analytic model which predicts the extremely slow-time evolution of the action variables along the resonances of multiplicity one. We base our model on two concepts: (i) By considering a (quasi-)stationary phase approach to the analysis of the Nekhoroshev normal form, we demonstrate that only a small fraction of the terms of the associated optimal remainder provide meaningful contributions to the evolution of the action variables. (ii) We provide rigorous analytical approximations to the Melnikov integrals of terms with stationary or quasi-stationary phase. Applying our model to an example of three degrees of freedom steep Hamiltonian provides the speed of Arnold diffusion, as well as a precise representation of the evolution of the action variables, in very good agreement (over several orders of magnitude) with the numerically computed one.

preprint2016arXiv

Analytical forms of chaotic spiral arms

We develop an analytical theory of chaotic spiral arms in galaxies. This is based on the Moser theory of invariant manifolds around unstable periodic orbits. We apply this theory to the chaotic spiral arms, that start from the neighborhood of the Lagrangian points L1 and L2 at the end of the bar in a barred-spiral galaxy. The series representing the invariant manifolds starting at the Lagrangian points L1, L2, or unstable periodic orbits around L1 and L2, yield spiral patterns in the configuration space. These series converge in a domain around every Lagrangian point, called "Moser domain" and represent the orbits that constitute the chaotic spiral arms. In fact, these orbits are not only along the invariant manifolds, but also in a domain surrounding the invariant manifolds. We show further that orbits starting outside the Moser domain but close to it converge to the boundary of the Moser domain, which acts as an attractor. These orbits stay for a long time close to the spiral arms before escaping to infinity.

preprint2016arXiv

Geostationary secular dynamics revisited: application to high area-to-mass ratio objects

The long-term dynamics of the geostationary Earth orbits (GEO) is revisited through the application of canonical perturbation theory. We consider a Hamiltonian model accounting for all major perturbations: geopotential at order and degree two, lunisolar perturbations with a realistic model for the Sun and Moon orbits, and solar radiation pressure. The long-term dynamics of the GEO region has been studied both numerically and analytically, in view of the relevance of such studies to the issue of space debris or to the disposal of GEO satellites. Past studies focused on the orbital evolution of objects around a nominal solution, hereafter called the forced equilibrium solution, which shows a particularly strong dependence on the area-to-mass ratio. Here, we i) give theoretical estimates for the long-term behavior of such orbits, and ii) we examine the nature of the forced equilibrium itself. In the lowest approximation, the forced equilibrium implies motion with a constant non-zero average `forced eccentricity', as well as a constant non-zero average inclination, otherwise known in satellite dynamics as the inclination of the invariant `Laplace plane'. Using a higher order normal form, we demonstrate that this equilibrium actually represents not a point in phase space, but a trajectory taking place on a lower-dimensional torus. We give analytical expressions for this special trajectory, and we compare our results to those found by numerical orbit propagation. We finally discuss the use of proper elements, i.e., approximate integrals of motion for the GEO orbits.

preprint2016arXiv

New Hamiltonian expansions adapted to the Trojan problem

A number of studies, referring to the observed Trojan asteroids of various planets in our Solar System, or to hypothetical Trojan bodies in extrasolar planetary systems, have emphasized the importance of so-called secondary resonances in the problem of the long term stability of Trojan motions. Such resonances describe commensurabilities between the fast, synodic, and secular frequency of the Trojan body, and, possibly, additional slow frequencies produced by more than one perturbing bodies. The presence of secondary resonances sculpts the dynamical structure of the phase space. Hence, identifying their location is a relevant task for theoretical studies. In the present paper we combine the methods introduced in two recent papers (Paez & Efthymiopoulos, 2015, Paez & Locatelli, 2015) in order to analytically predict the location of secondary resonances in the Trojan problem (SEE FILE FOR COMPLETE ABSTRACT)

preprint2016arXiv

Origin of chaos in 3-d Bohmian trajectories

We study the 3-d Bohmian trajectories of a quantum system of three harmonic oscillators. We focus on the mechanism responsible for the generation of chaotic trajectories. We demonstrate the existence of a 3-d analogue of the mechanism found in earlier studies of 2-d systems, based on moving 2-d `nodal point - X-point complexes'. In the 3-d case, we observe a foliation of nodal point - X-point complexes, forming a `3-d structure of nodal and X-points'. Chaos is generated when the Bohmian trajectories are scattered at one or more close encounters with such a structure.

preprint2016arXiv

The theory of secondary resonances in the spin-orbit problem

We study the resonant dynamics in a simple one degree of freedom, time dependent Hamiltonian model describing spin-orbit interactions. The equations of motion admit periodic solutions associated with resonant motions, the most important being the synchronous one in which most evolved satellites of the Solar system, including the Moon, are observed. Such primary resonances can be surrounded by a chain of smaller islands which one refers to as secondary resonances. Here, we propose a novel canonical normalization procedure allowing to obtain a higher order normal form, by which we obtain analytical results on the stability of the primary resonances as well as on the bifurcation thresholds of the secondary resonances. The procedure makes use of the expansion in a parameter, called the detuning, measuring the shift from the exact secondary resonance. Also, we implement the so-called `book-keeping' method, i.e., the introduction of a suitable separation of the terms in orders of smallness in the normal form construction, which deals simultaneously with all the small parameters of the problem. Our analytical computation of the bifurcation curves is in excellent agreement with the results obtained by a numerical integration of the equations of motion, thus providing relevant information on the parameter regions where satellites can be found in a stable configuration.

preprint2015arXiv

Modeling Trojan dynamics: diffusion mechanisms through resonances

In the framework of the ERTBP, we study an example of the influence of secondary resonances over the long term stability of Trojan motions. By the integration of ensembles of orbits, we find various types of chaotic diffusion, slow and fast. We show that the distribution of escape times is bi-modular, corresponding to two populations of short and long escape times. The objects with long escape times produce a power-law tail in the distribution.

preprint2014arXiv

Trojan resonant dynamics, stability, and chaotic diffusion, for parameters relevant to exoplanetary systems

We investigate the dynamics of small trojan exoplanets in domains of secondary resonances within the tadpole domain of motion. We consider the limit of a massless trojan companion of a giant planet. Without other planets, this is a case of the elliptic restricted three body problem (ERTBP). The presence of more planets (the restricted multi-planet problem, RMPP) induces new direct and indirect secular effects on the trojan's dynamics. In the theoretical part of this paper, we develop a Hamiltonian formalism in action-angle variables, which allows to treat in a unified way resonant dynamics and secular effects on the trojan body in both the ERTBP or the RMPP. Our formalism leads to a decomposition of the Hamiltonian in two parts, $H=H_b+H_{sec}$. $H_b$, called the basic model, describes resonant dynamics in the short-period (epicyclic) and synodic (libration) degrees of freedom. $H_{sec}$ contains only terms depending on slow (secular) angles. $H_b$ is formally identical in the ERTBP and the RMPP, apart from a re-definition of angular variables. An important physical consequence is that the slow chaotic diffusion proceeds in both the ERTBP and the RMPP by a qualitatively similar dynamical mechanism better approximated by the paradigm of `modulational diffusion'. In the numerical part, we focus on the ERTBP for making a numerical demonstration of the chaotic diffusion process along resonances. Using color stability maps, we provide a survey of the resonant web for characteristic mass parameters of the primary, in which the secondary resonances from 1:5 to 1:12 (ratio of the short over the synodic period) and their resonant multiplets appear. We give numerical examples of diffusion of weakly chaotic orbits in the resonant web. We make a statistics of the escaping times in the resonant domain, and find power-law tails of the distribution of escaping times for slowly diffusing chaotic orbits.

preprint2013arXiv

Chirikov and Nekhoroshev diffusion estimates: bridging the two sides of the river

We present theoretical and numerical results pointing towards a strong connection between the estimates for the diffusion rate along simple resonances in multidimensional nonlinear Hamiltonian systems that can be obtained using the heuristic theory of Chirikov and a more formal one due to Nekhoroshev. We show that, despite a wide-spread impression, the two theories are complementary rather than antagonist. Indeed, although Chirikov's 1979 review has thousands of citations, almost all of them refer to topics such as the resonance overlap criterion, fast diffusion, the Standard or Whisker Map, and not to the constructive theory providing a formula to measure diffusion along a single resonance. However, as will be demonstrated explicitly below, Chirikov's formula provides values of the diffusion coefficient which are quite well comparable to the numerically computed ones, provided that it is implemented on the so-called optimal normal form derived as in the analytic part of Nekhoroshev's theorem. On the other hand, Chirikov's formula yields unrealistic values of the diffusion coefficient, in particular for very small values of the perturbation, when used in the original Hamiltonian instead of the optimal normal form. In the present paper, we take advantage of this complementarity in order to obtain accurate theoretical predictions for the local value of the diffusion coefficient along a resonance in a specific 3DoF nearly integrable Hamiltonian system. Besides, we compute numerically the diffusion coefficient and a full comparison of all estimates is made for ten values of the perturbation parameter, showing a very satisfactory agreement.

preprint2013arXiv

Low-dimensional q-Tori in FPU Lattices: Dynamics and Localization Properties

This is a continuation of our study concerning q-tori, i.e. tori of low dimensionality in the phase space of nonlinear lattice models like the Fermi-Pasta-Ulam (FPU) model. In our previous work we focused on the beta FPU system, and we showed that the dynamical features of the q-tori serve as an interpretational tool to understand phenomena of energy localization in the FPU space of linear normal modes. In the present paper i) we employ the method of Poincare - Lindstedt series, for a fixed set of frequencies, in order to compute an explicit quasi-periodic representation of the trajectories lying on q-tori in the alpha model, and ii) we consider more general types of initial excitations in both the alpha and beta models. Furthermore we turn into questions of physical interest related to the dynamical features of the q-tori. We focus on particular q-tori solutions describing low-frequency `packets' of modes, and excitations of a small set of modes with an arbitrary distribution in q-space. In the former case, we find formulae yielding an exponential profile of energy localization, following an analysis of the size of the leading order terms in the Poincare - Lindstedt series. In the latter case, we explain the observed localization patterns on the basis of a rigorous result concerning the propagation of non-zero terms in the Poincare - Lindstedt series from zeroth to subsequent orders. Finally, we discuss the extensive (i.e. independent of the number of degrees of freedom) properties of some q-tori solutions.

preprint2013arXiv

The speed of Arnold diffusion

A detailed numerical study is presented of the slow diffusion (Arnold diffusion) taking place around resonance crossings in nearly integrable Hamiltonian systems of three degrees of freedom in the so-called `Nekhoroshev regime'. The aim is to construct estimates regarding the speed of diffusion based on the numerical values of a truncated form of the so-called remainder of a normalized Hamiltonian function, and to compare them with the outcomes of direct numerical experiments using ensembles of orbits. In this comparison we examine, one by one, the main steps of the so-called analytic and geometric parts of the Nekhoroshev theorem. We are led to two main results: i) We construct in our concrete example a convenient set of variables, proposed first by Benettin and Gallavotti (1986), in which the phenomenon of Arnold diffusion in doubly resonant domains can be clearly visualized. ii) We determine, by numerical fitting of our data the dependence of the local diffusion coefficient "D" on the size "||R_{opt}||" of the optimal remainder function, and we compare this with a heuristic argument based on the assumption of normal diffusion. We find a power law "D\propto ||R_{opt}||^{2(1+b)}", where the constant "b" has a small positive value depending also on the multiplicity of the resonance considered.

preprint2008arXiv

Explicit Construction of First Integrals with Quasi-monomial Terms from the Painlevé Series

The Painlevé and weak Painlevé conjectures have been used widely to identify new integrable nonlinear dynamical systems. For a system which passes the Painlevé test, the calculation of the integrals relies on a variety of methods which are independent from Painlevé analysis. The present paper proposes an explicit algorithm to build first integrals of a dynamical system, expressed as `quasi-polynomial' functions, from the information provided solely by the Painlevé - Laurent series solutions of a system of ODEs. Restrictions on the number and form of quasi-monomial terms appearing in a quasi-polynomial integral are obtained by an application of a theorem by Yoshida (1983). The integrals are obtained by a proper balancing of the coefficients in a quasi-polynomial function selected as initial ansatz for the integral, so that all dependence on powers of the time $τ=t-t_0$ is eliminated. Both right and left Painlevé series are useful in the method. Alternatively, the method can be used to show the non-existence of a quasi-polynomial first integral. Examples from specific dynamical systems are given.

preprint2006arXiv

Invariant manifolds, phase correlations of chaotic orbits and the spiral structure of galaxies

In the presence of a strong $m=2$ component in a rotating galaxy, the phase space structure near corotation is shaped to a large extent by the {\it invariant manifolds} of the short period family of unstable periodic orbits terminating at L$_1$ or L$_2$. The main effect of these manifolds is to create robust {\it phase correlations} among a number of chaotic orbits large enough to support a {\it spiral} density wave outside corotation. The phenomenon is described theoretically by soliton-like solutions of a Sine-Gordon equation. Numerical examples are given in an N-Body simulation of a barred spiral galaxy. In these examples, we demonstrate how the projection of unstable manifolds in configuration space reproduces essentially the entire observed bar-spiral pattern.

preprint2006arXiv

Special Features of Galactic Dynamics

The present lecture notes are an introduction to selected topics of {\it Galactic Dynamics}. The focus is on topics that we consider more relevant to the main theme of this workshop, {\it Celestial Mechanics}. This is not intended to be a review article. In fact, any of the topics below could be the subject of a separate review. Only the main ideas and notions are introduced, as well as some important currently open problems in each topic. Some relevant results from our own research are also presented. We discuss topics related mostly to the so-called {\it ellipsoidal components} of galaxies. These are \textbf{a)} the dark halos of both elliptical and disk galaxies, \textbf{b)} the luminous matter in elliptical galaxies, and \textbf{c)} the bulges of disk galaxies. We shall only occasionally refer to the dynamics of disks, bars or spiral structure.