Source author record

C. Efthymiopoulos

C. 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

15works
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

15 published item(s)

preprint2019arXiv

Manifold spirals in barred galaxies with multiple pattern speeds

In the manifold theory of spiral structure in barred galaxies, the usual assumption is that the spirals rotate with the same pattern speed as the bar. Here we generalize the manifold theory under the assumption that the spirals rotate with different pattern speed than the bar. More generally, we consider the case when one or more modes, represented by the potentials V_2, V_3, \ldots, co-exist in the galactic disc in addition to the bar's mode V_{bar}, but rotate with pattern speeds Ω_2, Ω_3, \ldots incommensurable between themselves and with Ω_{bar}. Through a perturbative treatment (assuming that V_2,V_3... are small with respect to V_{bar}) we then show that the unstable Lagrangian points L_1, L_2 of the pure bar model (V_{bar},Ω_{bar}) are `continued' in the full model as periodic orbits, when we have one extra pattern speed different from Ω_{bar}, or as epicyclic `Lissajous-like' unstable orbits, when we have more than one extra pattern speeds. As an example we compute the generalized orbits GL_1, GL_2 and their manifolds in a Milky-way like model with bar and spiral pattern speeds assumed different. We find that the manifolds produce a time-varying morphology consisting of segments of spirals or `pseudorings'. These structures are repeated after a period equal to half the relative period of the imposed spirals with respect to the bar. Along one period, the manifold-induced time-varying structures are found to continuously support at least some part of the imposed spirals, except at short intervals around those times at which the relative phase of the imposed spirals with respect to the bar becomes equal to \pmπ/2. A connection of these effects to the phenomenon of recurrent spirals is discussed.

preprint2016arXiv

Analytical study of chaos and applications

We summarize various cases where chaotic orbits can be described analytically. First we consider the case of a magnetic bottle where we have non-resonant and resonant ordered and chaotic orbits. In the sequence we consider the hyperbolic Henon map, where chaos appears mainly around the origin, which is an unstable periodic orbit. In this case the chaotic orbits around the origin are represented by analytic series (Moser series). We find the domain of convergence of these Moser series and of similar series around other unstable periodic orbits. The asymptotic manifolds from the various unstable periodic orbits intersect at homoclinic and heteroclinic orbits that are given analytically. Then we consider some Hamiltonian systems and we find their homoclinic orbits by using a new method of analytic prolongation. An application of astronomical interest is the domain of convergence of the analytical series that determine the spiral structure of barred-spiral galaxies.

preprint2016arXiv

Structures induced by companions in galactic discs

Using N-body simulations we study the structures induced on a galactic disc by repeated flybys of a companion in decaying eccentric orbit around the disc. Our system is composed by a stellar disc, bulge and live dark matter halo, and we study the system's dynamical response to a sequence of a companion's flybys, when we vary i) the disc's temperature (parameterized by Toomre's Q-parameter) and ii) the companion's mass and initial orbit. We use a new 3D Cartesian grid code: MAIN (Mesh-adaptive Approximate Inverse N-body solver). The main features of MAIN are reviewed, with emphasis on the use of a new Symmetric Factored Approximate Sparse Inverse (SFASI) matrix in conjunction with the multigrid method that allows the efficient solution of Poisson's equation in three space variables. We find that: i) companions need to be assigned initial masses in a rather narrow window of values in order to produce significant and more long-standing non-axisymmetric structures (bars and spirals) in the main galaxy's disc by the repeated flyby mechanism. ii) a crucial phenomenon is the antagonism between companion-excited and self-excited modes on the disc. Values of $Q >1.5$ are needed in order to allow for the growth of the companion-excited modes to prevail over the the growth of the disc's self-excited modes. iii) We give evidence that the companion-induced spiral structure is best represented by a density wave with pattern speed nearly constant in a region extending from the ILR to a radius close to, but inside, corotation.

preprint2015arXiv

Analytical description of the structure of chaos

We consider analytical formulae that describe the chaotic regions around the main periodic orbit $(x=y=0)$ of the Hénon map. Following our previous paper (Efthymiopoulos, Contopoulos, Katsanikas $2014$) we introduce new variables $(ξ, η)$ in which the product $ξη=c$ (constant) gives hyperbolic invariant curves. These hyperbolae are mapped by a canonical transformation $Φ$ to the plane $(x,y)$, giving "Moser invariant curves". We find that the series $Φ$ are convergent up to a maximum value of $c=c_{max}$. We give estimates of the errors due to the finite truncation of the series and discuss how these errors affect the applicability of analytical computations. For values of the basic parameter $κ$ of the Hénon map smaller than a critical value, there is an island of stability, around a stable periodic orbit $S$, containing KAM invariant curves. The Moser curves for $c \leq 0.32$ are completely outside the last KAM curve around $S$, the curves with $0.32<c<0.41$ intersect the last KAM curve and the curves with $0.41\leq c< c_{max} \simeq 0.49$ are completely inside the last KAM curve. All orbits in the chaotic region around the periodic orbit $(x=y=0)$, although they seem random, belong to Moser invariant curves, which, therefore define a "structure of chaos". Orbits starting close and outside the last KAM curve remain close to it for a stickiness time that is estimated analytically using the series $Φ$. We finally calculate the periodic orbits that accumulate close to the homoclinic points, i.e. the points of intersection of the asymptotic curves from $x=y=0$, exploiting a method based on the self-intersections of the invariant Moser curves. We find that all the computed periodic orbits are generated from the stable orbit $S$ for smaller values of the Hénon parameter $κ$, i.e. they are all regular periodic orbits.

preprint2015arXiv

Effective power-law dependence of Lyapunov exponents on the central mass in galaxies

Using both numerical and analytical approaches, we demonstrate the existence of an effective power-law relation $L\propto m^p$ between the mean Lyapunov exponent $L$ of stellar orbits chaotically scattered by a supermassive black hole in the center of a galaxy and the mass parameter $m$, i.e. ratio of the mass of the black hole over the mass of the galaxy. The exponent $p$ is found numerically to obtain values in the range $p \approx 0.3$--$0.5$. We propose a theoretical interpretation of these exponents, based on estimates of local `stretching numbers', i.e. local Lyapunov exponents at successive transits of the orbits through the black hole's sphere of influence. We thus predict $p=2/3-q$ with $q\approx 0.1$--$0.2$. Our basic model refers to elliptical galaxy models with a central core. However, we find numerically that an effective power law scaling of $L$ with $m$ holds also in models with central cusp, beyond a mass scale up to which chaos is dominated by the influence of the cusp itself. We finally show numerically that an analogous law exists also in disc galaxies with rotating bars. In the latter case, chaotic scattering by the black hole affects mainly populations of thick tube-like orbits surrounding some low-order branches of the $x_1$ family of periodic orbits, as well as its bifurcations at low-order resonances, mainly the Inner Lindbland resonance and the 4/1 resonance. Implications of the correlations between $L$ and $m$ to determining the rate of secular evolution of galaxies are discussed.

preprint2015arXiv

Resonant normal form and asymptotic normal form behavior in magnetic bottle Hamiltonians

We consider normal forms in `magnetic bottle' type Hamiltonians of the form $H=\frac{1}{2}(ρ^2_ρ+ω^2_1ρ^2) +\frac{1}{2}p^2_z+hot$ (second frequency $ω_2$ equal to zero in the lowest order). Our main results are: i) a novel method to construct the normal form in cases of resonance, and ii) a study of the asymptotic behavior of both the non-resonant and the resonant series. We find that, if we truncate the normal form series at order $r$, the series remainder in both constructions decreases with increasing $r$ down to a minimum, and then it increases with $r$. The computed minimum remainder turns to be exponentially small in $\frac{1}{ΔE}$, where $ΔE$ is the mirror oscillation energy, while the optimal order scales as an inverse power of $ΔE$. We estimate numerically the exponents associated with the optimal order and the remainder's exponential asymptotic behavior. In the resonant case, our novel method allows to compute a `quasi-integral' (i.e. truncated formal integral) valid both for each particular resonance as well as away from all resonances. We applied these results to a specific magnetic bottle Hamiltonian. The non resonant normal form yields theorerical invariant curves on a surface of section which fit well the empirical curves away from resonances. On the other hand the resonant normal form fits very well both the invariant curves inside the islands of a particular resonance as well as the non-resonant invariant curves. Finally, we discuss how normal forms allow to compute a critical threshold for the onset of global chaos in the magnetic bottle.

preprint2014arXiv

Analytical invariant manifolds near unstable points and the structure of chaos

It is known that the asymptotic invariant manifolds around an unstable periodic orbit in conservative systems can be represented by convergent series (Cherry 1926, Moser 1956, 1958, Giorgilli 2001). The unstable and stable manifolds intersect at an infinity of homoclinic points, generating a complicated homoclinic tangle. In the case of simple mappings it was found (Da Silva Ritter et al. 1987) that the domain of convergence of the formal series extends to infinity along the invariant manifolds. This allows in practice to study the homoclinic tangle using only series. However in the case of Hamiltonian systems, or mappings with a finite analyticity domain,the convergence of the series along the asymptotic manifolds is also finite. Here, we provide numerical indications that the convergence does not reach any homoclinic points. We discuss in detail the convergence problem in various cases and we find the degree of approximation of the analytical invariant manifolds to the real (numerical) manifolds as i) the order of truncation of the series increases, and ii) we use higher numerical precision in computing the coefficients of the series. Then we introduce a new method of series composition, by using action-angle variables, that allows the calculation of the asymptotic manifolds up to an a arbitrarily large extent. This is the first case of an analytic development that allows the computation of the invariant manifolds and their intersections in a Hamiltonian system for an extent long enough to allow the study of homoclinic chaos by analytical means.

preprint2014arXiv

Analytical study of the structure of chaos near unstable points

In a 2D conservative Hamiltonian system there is a formal integral $Φ$ besides the energy H. This is not convergent near a stable periodic orbit, but it is convergent near an unstable periodic orbit. We explain this difference and we find the convergence radius along the asymptotic curves. In simple mappings this radius is infinite. This allows the theoretical calculation of the asymptotic curves and their intersections at homoclinic points. However in more complex mappings and in Hamiltonian systems the radius of convergence is in general finite and does not allow the theoretical calculation of any homoclinic point. Then we develop a method similar to analytic continuation, applicable in systems expressed in action-angle variables, that allows the calculation of the asymptotic curves to an arbitrary length. In this way we can study analytically the chaotic regions near the unstable periodic orbit and near its homoclinic points.

preprint2014arXiv

Modeling resonant trojan motion in planetary systems

We consider the dynamics of a small trojan companion of a hypothetical giant exoplanet under the secular perturbations of additional planets. By a suitable choice of action-angle variables, the problem is amenable to the study of the slow modulation, induced by secular perturbations, to the dynamics of an otherwise called `basic' Hamiltonian model of two degrees of freedom (planar case). We present this Hamiltonian decomposition, which implies that the slow chaotic diffusion at resonances is best described by the paradigm of modulational diffusion.

preprint2013arXiv

Combining particle acceleration and coronal heating via data-constrained calculations of nanoflares in coronal loops

We model nanoflare heating of extrapolated active-region coronal loops via the acceleration of electrons and protons in Harris-type current sheets. The kinetic energy of the accelerated particles is estimated using semi-analytical and test-particle-tracing approaches. Vector magnetograms and photospheric Doppler velocity maps of NOAA active region 09114, recorded by the Imaging Vector Magnetograph (IVM), were used for this analysis. A current-free field extrapolation of the active-region corona was first constructed. The corresponding Poynting fluxes at the footpoints of 5000 extrapolated coronal loops were then calculated. Assuming that reconnecting current sheets develop along these loops, we utilized previous results to estimate the kinetic-energy gain of the accelerated particles and we related this energy to nanoflare heating and macroscopic loop characteristics. Kinetic energies of 0.1 to 8 keV (for electrons) and 0.3 to 470 keV (for protons) were found to cause heating rates ranging from $10^{-6}$ to 1 $\mathrm{erg\, s^{-1} cm^{-3}}$. Hydrodynamic simulations show that such heating rates can sustain plasma in coronal conditions inside the loops and generate plasma thermal distributions which are consistent with active region observations. We concluded the analysis by computing the form of X-ray spectra generated by the accelerated electrons using the thick target approach that were found to be in agreement with observed X-ray spectra, thus supporting the plausibility of our nanoflare-heating scenario.

preprint2012arXiv

Order in de Broglie - Bohm quantum mechanics

A usual assumption in the so-called {\it de Broglie - Bohm} approach to quantum dynamics is that the quantum trajectories subject to typical `guiding' wavefunctions turn to be quite irregular, i.e. {\it chaotic} (in the dynamical systems' sense). In the present paper, we consider mainly cases in which the quantum trajectories are {\it ordered}, i.e. they have zero Lyapunov characteristic numbers. We use perturbative methods to establish the existence of such trajectories from a theoretical point of view, while we analyze their properties via numerical experiments. Using a 2D harmonic oscillator system, we first establish conditions under which a trajectory can be shown to avoid close encounters with a moving nodal point, thus avoiding the source of chaos in this system. We then consider series expansions for trajectories both in the interior and the exterior of the domain covered by nodal lines, probing the domain of convergence as well as how successful the series are in comparison with numerical computations or regular trajectories. We then examine a Hénon - Heiles system possessing regular trajectories, thus generalizing previous results. Finally, we explore a key issue of physical interest in the context of the de Broglie - Bohm formalism, namely the influence of order in the so-called {\it quantum relaxation} effect. We show that the existence of regular trajectories poses restrictions to the quantum relaxation process, and we give examples in which the relaxation is suppressed even when we consider initial ensembles of only chaotic trajectories, provided, however, that the system as a whole is characterized by a certain degree of order.

preprint2011arXiv

Quantum vortices and trajectories in particle diffraction

We investigate the phenomenon of the diffraction of charged particles by thin material targets using the method of the de Broglie-Bohm quantum trajectories. The particle wave function can be modeled as a sum of two terms $ψ=ψ_{ingoing}+ψ_{outgoing}$. A thin separator exists between the domains of prevalence of the ingoing and outgoing wavefunction terms. The structure of the quantum-mechanical currents in the neighborhood of the separator implies the formation of an array of \emph{quantum vortices}. The flow structure around each vortex displays a characteristic pattern called `nodal point - X point complex'. The X point gives rise to stable and unstable manifolds. We find the scaling laws characterizing a nodal point-X point complex by a local perturbation theory around the nodal point. We then analyze the dynamical role of vortices in the emergence of the diffraction pattern. In particular, we demonstrate the abrupt deflections, along the direction of the unstable manifold, of the quantum trajectories approaching an X-point along its stable manifold. Theoretical results are compared to numerical simulations of quantum trajectories. We finally calculate the {\it times of flight} of particles following quantum trajectories from the source to detectors placed at various scattering angles $θ$, and thereby propose an experimental test of the de Broglie - Bohm formalism.

preprint2011arXiv

Wavepacket approach to particle diffraction by thin targets: Quantum trajectories and arrival times

We develop a wavepacket approach to the diffraction of charged particles by a thin material target and we use the de Broglie-Bohm quantum trajectories to study various phenomena in this context. We find the form of the separator, i.e.the limit between the domains of prevalence of the ingoing and outgoing quantum flow. The structure of the quantum-mechanical currents in the neighborhood of the separator implies the formation of an array of \emph{quantum vortices} (nodal point - X point complexes). We show how the deformation of the separatior near Bragg angles explains the emergence of a diffraction pattern by the de Broglie - Bohm trajectories. We calculate the arrival time distributions for particles scattered at different angles. The predictions of the de Broglie - Bohm theory for $ΔT$ turn to be different from estimates of the same quantity using other theories on time observables like the sum-over-histories or the Kijowski approach. We propose an experimental setup aiming to test such predictions. Finally, we explore the semiclassical limit of short wavelength and short quantum coherence lengths, and demonstrate how, in this case, results with the de Broglie - Bohm trajectories are similar to the classical results of Rutherford scattering.

preprint2009arXiv

Energy localization on q-tori, long term stability and the interpretation of FPU recurrences

We focus on two approaches that have been proposed in recent years for the explanation of the so-called FPU paradox, i.e. the persistence of energy localization in the `low-q' Fourier modes of Fermi-Pasta-Ulam nonlinear lattices, preventing equipartition among all modes at low energies. In the first approach, a low-frequency fraction of the spectrum is initially excited leading to the formation of `natural packets' exhibiting exponential stability, while in the second, emphasis is placed on the existence of `q-breathers', i.e periodic continuations of the linear modes of the lattice, which are exponentially localized in Fourier space. Following ideas of the latter, we introduce in this paper the concept of `q-tori' representing exponentially localized solutions on low-dimensional tori and use their stability properties to reconcile these two approaches and provide a more complete explanation of the FPU paradox.

preprint2008arXiv

The production of Tsallis entropy in the limit of weak chaos and a new indicator of chaoticity

We study the connection between the appearance of a `metastable' behavior of weakly chaotic orbits, characterized by a constant rate of increase of the Tsallis q-entropy (Tsallis 1988), and the solutions of the variational equations of motion for the same orbits. We demonstrate that the variational equations yield transient solutions, lasting for long time intervals, during which the length of deviation vectors of nearby orbits grows in time almost as a power-law. The associated power exponent can be simply related to the entropic exponent for which the q-entropy exhibits a constant rate of increase. This analysis leads to the definition of a new sensitive indicator distinguishing regular from weakly chaotic orbits, that we call `Average Power Law Exponent' (APLE). We compare the APLE with other established indicators of the literature. In particular, we give examples of application of the APLE in a) a thin separatrix layer of the standard map, b) the stickiness region around an island of stability in the same map, and c) the web of resonances of a 4D symplectic map. In all these cases we identify weakly chaotic orbits exhibiting the `metastable' behavior associated with the Tsallis q-entropy.