Source author record

Euaggelos E. Zotos

Euaggelos E. Zotos 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

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

37 published item(s)

preprint2021arXiv

Numerical investigation on the Hill's type lunar problem with homogeneous potential

We consider the planar Hill's lunar problem with a homogeneous gravitational potential. The investigation of the system is twofold. First, the starting conditions of the trajectories are classified into three classes, that is bounded, escaping, and collisional. Second, we study the no-return property of the Lagrange point $L_2$ and we observe that the escaping trajectories are scattered exponentially. Moreover, it is seen that in the supercritical case, with $α\geq 2$, the basin boundaries are smooth. On the other hand, in the subcritical case, with $α< 2$ the boundaries between the different types of basins exhibit fractal properties.

preprint2020arXiv

Short-term stability of particles in the WD J0914+1914 white dwarf planetary system

Nearly all known white dwarf planetary systems contain detectable rocky debris in the stellar photosphere. A glaring exception is the young and still evolving white dwarf WD J0914+1914, which instead harbours a giant planet and a disc of pure gas. The stability boundaries of this disc and the future prospects for this white dwarf to be polluted with rocks depend upon the mass and orbit of the planet, which are only weakly constrained. Here we combine an ensemble of plausible planet orbits and masses to determine where observers should currently expect to find the outer boundary of the gas disc. We do so by performing a sweep of the entire plausible phase space with short-term numerical integrations. We also demonstrate that particle-star collisional trajectories, which would lead to the (unseen) signature of rocky metal pollution, occupy only a small fraction of the phase space, mostly limited to particle eccentricities above 0.75. Our analysis reveals that a highly inflated planet on a near-circular orbit is the type of planet which is most consistent with the current observations.

preprint2020arXiv

The grain size survival threshold in one-planet post-main-sequence exoplanetary systems

The size distribution and orbital architecture of dust, grains, boulders, asteroids, and major planets during the giant branch phases of evolution dictate the preponderance and observability of the eventual debris, which have been found to surround white dwarfs and pollute their atmospheres with metals. Here, we utilize the photogravitational planar restricted three-body problem in one-planet giant branch systems in order to characterize the orbits of grains as the parent star luminosity and mass undergo drastic changes. We perform a detailed dynamical analysis of the character of grain orbits (collisional, escape, or bounded) as a function of location and energy throughout giant branch evolution. We find that for stars with main-sequence masses of $2.0M_{\odot}$, giant branch evolution, combined with the presence of a planet, ubiquitously triggers escape in grains smaller than about 1 mm, while leaving grains larger than about 5 cm bound to the star. This result is applicable for systems with either a terrestrial or giant planet, is largely independent of the location of the planet, and helps establish a radiative size threshold for escape of small particles in giant branch planetary systems.

preprint2016arXiv

Determining the type of orbits in the central regions of barred galaxies

We use a simple dynamical model which consists of a harmonic oscillator and a spherical component, in order to investigate the regular or chaotic character of orbits in a barred galaxy with a central spherically symmetric nucleus. Our aim is to explore how the basic parameters of the galactic system influence the nature of orbits, by computing in each case the percentage of chaotic orbits, as well as the percentages of different types of regular orbits. We also give emphasis to the types of regular orbits that support either the formation of nuclear rings or the barred structure of the galaxy. We provide evidence that the traditional x1 orbital family does not always dominate in barred galaxy models since we found several other types of resonant orbits which can also support the barred structure. We also found that sparse enough nuclei, fast rotating bars and high energy models can support the galactic bars. On the other hand, weak bars, dense central nuclei, slowly rotating bars and low energy models favor the formation of nuclear rings. We also compare our results with previous related work.

preprint2016arXiv

Fractal basins of attraction in the planar circular restricted three-body problem with oblateness and radiation pressure

In this paper we use the planar circular restricted three-body problem where one of the primary bodies is an oblate spheroid or an emitter of radiation in order to determine the basins of attraction associated with the equilibrium points. The evolution of the position of the five Lagrange points is monitored when the values of the mass ratio $μ$, the oblateness coefficient $A_1$, and the radiation pressure factor $q$ vary in predefined intervals. The regions on the configuration $(x,y)$ plane occupied by the basins of attraction are revealed using the multivariate version of the Newton-Raphson method. The correlations between the basins of convergence of the equilibrium points and the corresponding number of iterations needed in order to obtain the desired accuracy are also illustrated. We conduct a thorough and systematic numerical investigation demonstrating how the dynamical quantities $μ$, $A_1$, and $q$ influence the basins of attractions. Our results suggest that the mass ratio and the radiation pressure factor are the most influential parameters, while on the other hand the structure of the basins of convergence are much less affected by the oblateness coefficient.

preprint2016arXiv

Orbital and escape dynamics in barred galaxies - I. The 2D system

In this paper we use the two-dimensional (2D) version of a new analytical gravitational model in order to explore the orbital as well as the escape dynamics of the stars in a barred galaxy composed of a spherically symmetric central nucleus, a bar, a flat disk and a dark matter halo component. A thorough numerical investigation is conducted for distinguishing between bounded and escaping motion. Furthermore bounded orbits are further classified into non-escaping regular and trapped chaotic using the Smaller ALingment Index (SALI) method. Our aim is to determine the basins of escape through the two symmetrical escape channels around the Lagrange points $L_2$ and $L_3$ and also to relate them with the corresponding distribution of the escape rates of the orbits. We integrate initial conditions of orbits in several types of planes so as to obtain a more complete view of the overall orbital properties of the dynamical system. We also present evidence that the unstable manifolds which guide the orbits in and out the interior region are directly related with the formation of spiral and ring stellar structures observed in barred galaxies. In particular, we examine how the bar's semi-major axis determines the resulting morphologies. Our numerical simulations indicate that weak barred structures favour the formation of $R_1$ rings or $R_1'$ pseudo-rings, while strong bars on the other hand, give rise to $R_1R_2$ and open spiral morphologies. Our results are compared with earlier related work. The escape dynamics and the properties of the manifolds of the full three-dimensional (3D) galactic system will be given in an accompanying paper.

preprint2016arXiv

Orbital and escape dynamics in barred galaxies - II. The 3D system: Exploring the role of the normally hyperbolic invariant manifolds

A three degrees of freedom (3-dof) barred galaxy model composed of a spherically symmetric nucleus, a bar, a flat disc and a spherically symmetric dark matter halo is used for investigating the dynamics of the system. We use colour-coded plots to demonstrate how the value of the semi-major axis of the bar influences the regular or chaotic dynamics of the 3-dof system. For distinguishing between ordered and chaotic motion we use the Smaller ALingment Index (SALI) method, a fast yet very accurate tool. Undoubtedly, the most important elements of the dynamics are the normally hyperbolic invariant manifolds (NHIMs) located in the vicinity of the index 1 Lagrange points $L_2$ and $L_3$. These manifolds direct the flow of stars over the saddle points, while they also trigger the formation of rings and spirals. The dynamics in the neighbourhood of the saddle points is visualized by bifurcation diagrams of the Lyapunov orbits as well as by the restriction of the Poincaré map to the NHIMs. In addition, we reveal how the semi-major axis of the bar influences the structure of these manifolds which determine the final stellar structure (rings or spirals). Our numerical simulations suggest that in galaxies with weak bars the formation of $R_1$ rings or $R_1'$ pseudo-rings is favoured. In the case of galaxies with intermediate and strong bars the invariant manifolds seem to give rise to $R_1R_2$ rings and twin spiral formations, respectively. We also compare our numerical outcomes with earlier related work and with observational data.

preprint2015arXiv

Introducing a new 3D dynamical model for barred galaxies

The regular or chaotic dynamics of an analytical realistic three dimensional model composed of a spherically symmetric central nucleus, a bar and a flat disk is investigated. For describing the properties of the bar we introduce a new simple dynamical model and we explore the influence on the character of orbits of all the involved parameters of it, such as the mass and the scale length of the bar, the major semi-axis and the angular velocity of the bar as well as the energy. Regions of phase space with ordered and chaotic motion are identified in dependence on these parameters and for breaking the rotational symmetry. First we study in detail the dynamics in the invariant plane $z = p_z = 0$ using the Poincaré map as a basic tool and then we study the full 3 dimensional case using the SALI method as principal tool for distinguishing between order and chaos. We also present strong evidence obtained through the numerical simulations that our new bar model can realistically describe the formation and the evolution of the observed twin spiral structure in barred galaxies.

preprint2015arXiv

Order and chaos in a three dimensional galaxy model

We explore the orbital dynamics of a realistic three dimensional model describing the properties of a disk galaxy with a spherically symmetric central dense nucleus and a triaxial dark matter halo component. Regions of phase space with regular and chaotic motion are identified depending on the parameter values for triaxiality of the dark matter halo and for breaking the rotational symmetry. The four dimensional Poincaré map of the three degrees of freedom system is analyzed by a study of its restriction to various two dimensional invariant subsets of its domain.

preprint2014arXiv

A Hamiltonian system of three degrees of freedom with eight channels of escape: The Great Escape

In this work, we try to shed some light to the nature of orbits in a three-dimensional potential of a perturbed harmonic oscillator with eight possible channels of escape, which was chosen as an interesting example of open three-dimensional Hamiltonian systems. In particular, we conduct a thorough numerical investigation distinguishing between regular and chaotic orbits as well as between trapped and escaping orbits, considering unbounded motion for several values of the energy. In an attempt to discriminate safely and with certainty between ordered and chaotic motion, we use the Smaller ALingment Index (SALI) detector, computed by integrating numerically the basic equations of motion as well as the variational equations. Of particular interest, is to locate the basins of escape towards the different escape channels and connect them with the corresponding escape periods of the orbits. We split our study into three different cases depending on the initial value of the $z$ coordinate which was used for launching the test particles. We found, that when the orbits are started very close to the primary $(x,y)$ plane the respective grids exhibit a high degree of fractalization, while on the other hand for orbits with relatively high values of $z_0$ several well-formed basins of escape emerge thus, reducing significantly the fractalization of the grids. It was also observed, that for values of energy very close to the escape energy the escape times of orbits are large, while for energy levels much higher than the escape energy the vast majority of orbits escape extremely fast or even immediately to infinity. We hope our outcomes to be useful for a further understanding of the escape process in open 3D Hamiltonian systems.

preprint2014arXiv

Classifying orbits in galaxy models with a prolate or an oblate dark matter halo component

We explore the nature of orbits of stars moving in the meridional plane $(R,z)$ of an axially symmetric galactic model with a disk, a spherical nucleus, and a flat biaxial dark matter halo component. In particular, we study the influence of all the involved parameters of the dynamical system, by computing both the percentage of chaotic orbits and the percentages of orbits of the main regular resonant families in each case. To distinguish between ordered and chaotic motion, we use the smaller alignment index (SALI) method to extensive samples of orbits by numerically integrating the equations of motion as well as the variational equations. Moreover, a method based on the concept of spectral dynamics that utilizes the Fourier transform of the time series of each coordinate is used to identify the various families of regular orbits and also to recognize the secondary resonances that bifurcate from them. Two cases are studied for every parameter: (i) the case where the halo component is prolate and (ii) the case where an oblate dark halo is present. Our numerical investigation indicates that all the dynamical quantities affect, more or less, the overall orbital structure. It was observed that the mass of the nucleus, the halo flattening parameter, the scale length of the halo, the angular momentum, and the orbital energy are the most influential quantities, while the effect of all the other parameters is much weaker. It was also found that all the parameters corresponding to the disk only have a minor influence on the nature of orbits. Furthermore, some other quantities, such as the minimum distance to the origin, the horizontal, and the vertical force, were tested as potential chaos detectors. Our analysis revealed that only general information can be obtained from these quantities. We also compared our results with early related work.

preprint2014arXiv

Determining the nature of orbits in a three-dimensional galaxy model hosting a BL Lacertae object

A three-dimensional dynamical model for a galaxy hosting a BL Lacertae object is constructed. The model consists of a logarithmic potential representing an elliptical host galaxy with a bulge of radius $c_b$ and a dense massive nucleus. Using numerical experiments, we try to distinguish between regular and chaotic motion in both 2D and 3D system. In particular, we investigate how the basic parameters of our model, such as the mass of the nucleus, the internal perturbation and the flattening parameters influence the amount and the degree of chaos. Interesting correlations are presented for both 2D and 3D dynamical models. Our numerical results are explained and supported using elementary theoretical arguments and analytical calculations. Of particular interest, is the local integral of motion which have been found to exist in the vicinity of stable periodic points. The obtained numerical outcomes of the present research, are linked and also compared with several data derived from observations.

preprint2014arXiv

Exploring the origin, the nature and the dynamical behaviour of distant stars in galaxy models

We explore the regular or chaotic nature of orbits moving in the meridional plane of an axially symmetric galactic gravitational model with a disk, a dense spherical nucleus and some additional perturbing terms corresponding to influence from nearby galaxies. In order to obtain this we use the Smaller ALingment Index (SALI) technique integrating extensive samples of orbits. Of particular interest is the study of distant, remote stars moving in large galactocentric orbits. Our extensive numerical experiments indicate that the majority of the distant stars perform chaotic orbits. However, there are also distant stars displaying regular motion as well. Most distant stars are ejected into the galactic halo on approaching the dense and massive nucleus. We study the influence of some important parameters of the dynamical system, such as the mass of the nucleus and the angular momentum, by computing in each case the percentage of regular and chaotic orbits. A second order polynomial relationship connects the mass of the nucleus and the critical angular momentum of the distant star. Some heuristic semi-theoretical arguments to explain and justify the numerically derived outcomes are also given. Our numerical calculations suggest that the majority of distant stars spend their orbital time in the halo where it is easy to be observed. We present evidence that the main cause for driving stars to distant orbits is the presence of the dense nucleus combined with the perturbation caused by nearby galaxies. The origin of young O and B stars observed in the halo is also discussed.

preprint2014arXiv

Order and chaos in a new 3D dynamical model describing motion in non axially symmetric galaxies

We present a new dynamical model describing 3D motion in non axially symmetric galaxies. The model covers a wide range of galaxies from a disk system to an elliptical galaxy by suitably choosing the dynamical parameters. We study the regular and chaotic character of orbits in the model and try to connect the degree of chaos with the parameter describing the deviation of the system from axial symmetry. In order to obtain this, we use the Smaller ALingment Index (SALI) technique by numerically integrating the basic equations of motion, as well as the variational equations for extensive samples of orbits. Our results suggest, that the influence of the deviation parameter on the portion of chaotic orbits strongly depends on the vertical distance $z$ from the galactic plane of the orbits. Using different sets of initial conditions, we show that the chaotic motion is dominant in galaxy models with low values of $z$, while in the case of stars with large values of $z$ the regular motion is more abundant, both in elliptical and disk galaxy models.

preprint2014arXiv

Order and Chaos in a Three-dimensional Binary System of Interacting Galaxies

In the present article, we present a galactic gravitational model of three degrees of freedom, in order to investigate and reveal the behavior of orbits in a binary quasar system. The two quasars are hosted in a pair of interacting disk galaxies. We study in detail the regular or chaotic character of motion, in two different cases: the time independent model in both 2D and 3D dynamical systems and the evolving 3D model. Our numerical calculations, indicate that a large fraction of orbits in the 2D system are chaotic in the time independent case. A careful analysis suggest that several Lindblad resonances are also responsible for the chaotic motion of stars in both hosts galaxies. In the time dependent system, we follow the evolution of 3D orbits in our dynamical model, as the two interacting host galaxies develop dense and massive quasars in their cores, by mass transportation from the disks to their nuclei. In this interesting case, there are orbits that change their orbital character from regular to chaotic and vise versa and also orbits that maintain their character during the galactic evolution. These results strongly indicate that the ordered or chaotic nature of 3D orbits depends, not only in the galactic interaction but also in the presence of the quasars in the galactic cores of the host galaxies. The outcomes derived from our dynamical model are compared with observational data. Some theoretical arguments to support the numerically derived outcomes are presented, both in 2D and 3D systems, while a comparison with similar earlier work is also made.

preprint2013arXiv

A rational galactic potential with accurate periodic orbits and quasi-integrals of motion

The motion in a simple, time independent rational galactic potential is studied. The potential is a generalization of a two dimensional harmonic oscillator potential and can be considered to describe plane motion in the central parts of a galaxy. There are cases, where the potential displays accurate periodic orbits, together with large chaotic regions. The starting position of these orbits can be found using analytical arguments. Using a quasi-integral of motion, we study the stability of the central periodic orbit and explain the degree of chaos in the inner regions of the galaxy. Our outcomes are compared with results coming from polynomial potentials with accurate periodic orbits and large chaotic regions.

preprint2013arXiv

Orbit classification in the meridional plane of a disk galaxy model with a spherical nucleus

We investigate the regular or chaotic nature of star orbits moving in the meridional plane of an axially symmetric galactic model with a disk and a spherical nucleus. We study the influence of some important parameters of the dynamical system, such as the mass and the scale length of the nucleus, the angular momentum or the energy, by computing in each case the percentage of chaotic orbits, as well as the percentages of orbits of the main regular resonant families. Some heuristic arguments to explain and justify the numerically derived outcomes are also given. Furthermore, we present a new method to find the threshold between chaos and regularity for both Lyapunov Characteristic Numbers and SALI, by using them simultaneously.

preprint2013arXiv

Revealing the evolution, the stability and the escapes of families of resonant periodic orbits in Hamiltonian systems

We investigate the evolution of families of periodic orbits in a bisymmetrical potential made up of a two-dimensional harmonic oscillator with only one quartic perturbing term, in a number of resonant cases. Our main objective is to compute sufficiently and accurately the position and the period of the periodic orbits. For the derivation of the above quantities (position and period) we deploy in each resonance case semi-numerical methods. The comparison of our semi-numerical results with those obtained by the numerical integration of the equations of motion indicates that, in every case the relative error is always less than 1% and therefore, the agreement is more than sufficient. Thus, we claim that semi-numerical methods are very effective tools for computing periodic orbits. We also study in detail, the case when the energy of the orbits is larger than the escape energy. In this case, the periodic orbits in almost all resonance families become unstable and eventually escape from the system. Our target is to calculate the escape period and the escape position of the periodic orbits and also monitor their evolution with respect to the value of the energy.

preprint2013arXiv

Using new dynamical indicators to distinguish between order and chaos in a galactic potential producing exact periodic orbits and chaotic components

The theory of the inverse problem is used in order to find a two dimensional galactic potential generating a mono-parametric family of elliptic periodic orbits. The potential is made up of a two-dimensional harmonic oscillator with perturbing terms of third and fourth degree and can be considered to describe local motion in the central parts of a barred galaxy. Our numerical calculations indicate that the potential produces also other families of orbits and there are cases, where distinct chaotic components are also observed. The numerical experiments suggest that there are regular as well as chaotic orbits supporting the barred structure. We use a variety of dynamical indicators, in order to determine the character of motion. Of significant interest, is the S(g) dynamical spectrum, which is used in order to identify the islandic motion of resonant orbits and also proves to be a very useful and fast indicator, in order to distinguish between regular and chaotic motion. Comparison with other methods for detecting chaos is also discussed. The present results are compared with outcomes from earlier work.

preprint2012arXiv

A new dynamical indicator for chaos detection in galactic Hamiltonian systems

A new dynamical parameter, the f-indicator, is introduced and used in order to distinguish between regular and chaotic motion in galactic Hamiltonian systems. Two kinds of galactic potentials are used: (i) a global potential, which describes the whole galaxy and (ii) a local potential, which is made up of perturbed harmonic oscillators and describes motion near an equilibrium point. The new indicator is based on the energies of the separable system along the x, y and z axis. Comparison between the outcomes obtained using the new dynamical parameter and other methods, such as the maximum Lyapunov Characteristic Exponent (L.C.E), or the S(c) dynamical spectrum, shows that the new dynamical indicator gives fast and reliable results concerning the regular or chaotic character of the orbits. The new indicator was tested in several Hamiltonian systems of two (2D) degrees and three (3D) degrees of freedom.

preprint2012arXiv

A new dynamical model for the study of galactic structure

In the present article, we present a new gravitational galactic model, describing motion in elliptical as well as in disk galaxies, by suitably choosing the dynamical parameters. Moreover, a new dynamical parameter, the S(g) spectrum, is introduced and used, in order to detect islandic motion of resonant orbits and the evolution of the sticky regions. We investigate the regular or chaotic character of motion, with emphasis in the different dynamical models and make an extensive study of the sticky regions of the system. We use the classical method of the Poincare (r-pr) phase plane and the new dynamical parameter of the S(g) spectrum. The LCE is used, in order to make an estimation of the degree of chaos in our galactic model. In both cases, the numerical calculations, suggest that our new model, displays a wide variety of families of regular orbits, compared to other galactic models. In addition to the regular motion, this new model displays also chaotic regions. Furthermore, the extent of the chaotic regions increases, as the value of the flatness parameter b of the model increases. Moreover, our simulations indicate, that the degree of chaos in elliptical galaxies, is much smaller than that in dense disk galaxies. In both cases numerical calculations show, that the degree of chaos increases linearly, as the flatness parameter b increases. In addition, a linear relationship between the critical value of angular momentum and the b parameter if found, in both cases (elliptical and disk galaxies). Some theoretical arguments to support the numerical outcomes are presented. Comparison with earlier work is also made.

preprint2012arXiv

A semi-numerical method for periodic orbits in a bisymmetrical potential

We use a semi-numerical method to find the position and period of periodic orbits in a bisymmetrical potential, made up of a two dimensional harmonic oscillator, with an additional term of a Plummer potential, in a number of resonant cases. The results are compared with the outcomes obtained by the numerical integration of the equations of motion and the agreement is good. This indicates that the semi-numerical method gives general and reliable results. Comparison with other methods of locating periodic orbits is also made.

preprint2012arXiv

Application of new dynamical spectra of orbits in Hamiltonian systems

In the present article, we investigate the properties of motion in Hamiltonian systems of two and three degrees of freedom, using the distribution of the values of two new dynamical parameters. The distribution functions of the new parameters, define the S(g) and the S(w) dynamical spectra. The first spectrum definition, that is the S(g) spectrum, will be applied in a Hamiltonian system of two degrees of freedom (2D), while the S(w) dynamical spectrum will be deployed in a Hamiltonian system of three degrees of freedom (3D). Both Hamiltonian systems, describe a very interesting dynamical system which displays a large variety of resonant orbits, different chaotic components and also several sticky regions. We test and prove the efficiency and the reliability of these new dynamical spectra, in detecting tiny ordered domains embedded in the chaotic sea, corresponding to complicated resonant orbits of higher multiplicity. The results of our extensive numerical calculations, suggest that both dynamical spectra are fast and reliable discriminants between different types of orbits in Hamiltonian systems, while requiring very short computation time in order to provide solid and conclusive evidence regarding the nature of an orbit. Furthermore, we establish numerical criteria in order to quantify the results obtained from our new dynamical spectra. A comparison to other previously used dynamical indicators, reveals the leading role of the new spectra.

preprint2012arXiv

Are external perturbations responsible for chaotic motion in galaxies?

We study the nature of motion in a logarithmic galactic dynamical model, with an additional external perturbation. Two different cases are investigated. In the first case the external perturbation is fixed, while in the second case it is varying with the time. Numerical experiments suggest, that responsible for the chaotic phenomena is the external perturbation, combined with the dense nucleus. Linear relationships are found to exist, between the critical value of the angular momentum and the dynamical parameters of the galactic system that is, the strength of the external perturbation, the flattening parameter and the radius of the nucleus. Moreover, the extent of the chaotic regions in the phase plane, increases linearly as the strength of the external perturbation and the flattening parameter increases. On the contrary, we observe that the percentage covered by chaotic orbits in the phase plane, decreases linearly, as the scale length of the nucleus increases, becoming less dense. Theoretical arguments are used to support and explain the numerically obtained outcomes. A comparison of the present outcomes with earlier results is also presented.

preprint2012arXiv

Are semi-numerical methods an effective tool for locating periodic orbits in 3D potentials?

A semi-numerical method is used in order to locate the position and calculate the period of periodic orbits in a 3D composite bisymmetrical potential, in a number of resonant cases. The potential consists of a 3D harmonic oscillator and a Plummer sphere. The outcomes are compared with results found using the numerical integration of the equations of motion and the agreement is very good. This agreement strongly suggests, that semi-numerical methods can be used in order to obtain fast and reliable results regarding the position and period of the periodic orbits in 3D composite potentials with a harmonic oscillator part and different kinds of perturbations. Comparison with other methods of obtaining 3D periodic orbits is discussed.

preprint2012arXiv

Chaos in a disk galaxy model induced by asymmetries in the dark halo

We study the regular or chaotic nature of motion in a disk galaxy with a dense nucleus and an asymmetric dark halo. Two cases, the 2D model and the 3D model, are investigated. In the 2D model, a considerable fraction of the phase plane is covered by chaotic orbits. Two factors seem to be responsible for the chaotic motion: (i) the dense nucleus and (ii) the asymmetries in the dark halo. Our numerical experiments suggest, that there are several chaotic components on the Poincare phase plane. Different chaotic components are induced by the asymmetries in the halo. Each chaotic component seems to have a different value of the Lyapunov Characteristic Exponent, for small values of the asymmetry parameter and a unique LCE for larger values of the asymmetry parameter. A comparison of the present results with outcomes from previous work is also presented.

preprint2012arXiv

Chaotic orbits in a 3D galactic dynamical model with a double nucleus

A 3D-dynamical model is constructed for the study of motion in the central regions of a disk galaxy with a double nucleus. Using the results of the 2D-model, we find the regions of initial conditions in the (x,px,z,py)=EJ, (y=pz=0) phase space producing regular or chaotic orbits. The majority of stars are on chaotic orbits. All chaotic orbits come arbitrary close to one or to both nuclei. Regular orbits are those starting near the stable periodic orbits of the 2D-system with small values of z0. All regular orbits circulate around only one of the two nuclei.

preprint2012arXiv

Comparing the behavior of orbits in different 3D dynamical models for elliptical galaxies

We study the behavior of orbits in two different galactic dynamical models, describing the motion in the central parts of a triaxial elliptical galaxy with a dense nucleus. Numerical experiments show that both models display regular motion together with extended chaotic regions. A detailed investigation of the properties of motion is made for the 2D and 3D Hamiltonian systems, using a number of different dynamical parameters, such as the Poincare surface of section, the maximal Lyapunov Characteristic Exponent, the S(c) spectrum, the S(w) spectrum and the P(f) indicator. The numerical calculations suggest that the properties of motion in both potentials are very similar. Our results show that one may use different kinds of gravitational potentials in order to describe the motion in triaxial galaxies while obtaining quantitatively similar results.

preprint2012arXiv

Dark halos acting as chaos controllers in asymmetric triaxial galaxy models

We study the regular or chaotic character of orbits in a 3D dynamical model, describing a triaxial galaxy surrounded by a spherical dark halo component. Our numerical experiments suggest that the percentage of chaotic orbits decreases exponentially as the mass of the dark halo increases. A linear increase of the percentage of the chaotic orbits was observed as the scale length of the halo component increases. In order to distinguish between regular and chaotic motion, we chose to use the total angular momentum Ltot of the 3D orbits as a new indicator. Comparison with other, previously used, dynamical indicators, such as the Lyapunov Characteristic Exponent or the P(f) spectral method, shows that the Ltot indicator gives very fast and reliable results for characterizing the nature of orbits in galactic dynamical models.

preprint2012arXiv

Exploring the nature of orbits in a galactic model with a massive nucleus

In the present article, we use an axially symmetric galactic gravitational model with a disk-halo and a spherical nucleus, in order to investigate the transition from regular to chaotic motion for stars moving in the meridian (r,z) plane. We study in detail the transition from regular to chaotic motion, in two different cases: the time independent model and the time evolving model. In the time dependent model, we follow the evolution of orbits as the galaxy develops a dense and massive nucleus in its core, as mass is transported exponentially from the disk to the galactic center. In addition, we construct some numerical diagrams in which we present the correlations between the main parameters of our galactic model. Our numerical calculations indicate, that stars with values of angular momentum Lz less than or equal to a critical value Lzc, moving near to the galactic plane, are scattered to the halo upon encountering the nuclear region and subsequently display chaotic motion. A linear relationship exists between the critical value of the angular momentum Lzc and the mass of the nucleus Mn. Furthermore, the extent of the chaotic region increases as the value of the mass of the nucleus increases. Moreover, our simulations indicate that the degree of chaos increases linearly, as the mass of the nucleus increases. These results strongly indicate that the ordered or chaotic nature of orbits, depends on the presence of massive objects in the galactic cores of the galaxies. Our results suggest, that for disk galaxies with massive and prominent nuclei, the low angular momentum stars in the associated central regions of the galaxy, must be in predominantly chaotic orbits. Some theoretical arguments to support the numerically derived outcomes are presented. Comparison with similar previous works is also made.

preprint2012arXiv

Investigating the nature of motion in 3D perturbed elliptic oscillators displaying exact periodic orbits

We study the nature of motion in a 3D potential composed of perturbed elliptic oscillators. Our technique is to use the results obtained from the 2D potential in order to find the initial conditions generating regular or chaotic orbits in the 3D potential. Both 2D and 3D potentials display exact periodic orbits together with extended chaotic regions. Numerical experiments suggest, that the degree of chaos increases rapidly, as the energy of the test particle increases. About 97% of the phase plane of the 2D system is covered by chaotic orbits for large energies. The regular or chaotic character of the 2D orbits is checked using the S(c) dynamical spectrum, while for the 3D potential we use the S(c) spectrum, along with the P(f) spectral method. Comparison with other dynamical indicators shows that the S(c) spectrum gives fast and reliable information about the character of motion.

preprint2012arXiv

Order and chaos in a galactic model with a strong nuclear bar

We use a composite gravitational galactic model consisting of a disk, a halo, a massive nucleus and a strong nuclear bar, in order to study the connections between global and local parameters in a realistic dynamical system. The local model is constructed from a two-dimensional perturbed harmonic oscillator and can be derived by expanding the global model in the vicinity of the central stable Lagrange equilibrium point. The frequencies of oscillations are not arbitrary, but they are connected with all the parameters involved with the global model. Moreover, the value of the local energy is also connected with the value of the global energy. Low and high energy stars in the global model display chaotic motion. Comparison with previous research reveals that the presence of the massive nucleus is responsible for the chaotic motion of the low energy stars. In the local motion, the low energy stars show interesting resonance phenomena, but the chaotic motion, if any, is negligible. On the contrary, the high energy stars do not show bounded motion in the local model. This is an indication of particular activity near the center of galaxies possessing massive nuclei.

preprint2012arXiv

Order and chaos in a triaxial galaxy model with a dark halo component

We study the regular or chaotic nature of orbits in a 3D potential describing a triaxial galaxy surrounded by a spherical dark halo component. Our numerical calculations show, that the percentage of chaotic orbits decreases exponentially, as the mass of the dark halo increases. A linear increase of the percentage of the chaotic orbits was observed, as the scale length of the dark halo component increases. In order to distinguish between regular and chaotic character of orbits, we use the total angular momentum Ltot, as a new indicator. Comparison of this new dynamical parameter, with other, previously used chaos indicators, shows that the Ltot gives very fast and reliable results in order to detect the character of orbits in galactic potentials.

preprint2012arXiv

The evolution of chaos in active galaxy models with an oblate or a prolate dark halo component

The evolution of chaotic motion in a galactic dynamical model with a disk, a dense nucleus and a flat biaxial dark halo component is investigated. Two cases are studied: (i) the case where the halo component is oblate and (ii) the case where a prolate halo is present. In both cases, numerical calculations show that the extent of the chaotic regions decreases exponentially as the scale-length of the dark halo increases. On the other hand, a linear relationship exists between the extent of the chaotic regions and the flatness parameter of the halo component. A linear relationship between the critical value of the angular momentum and the flatness parameter is also found. Some theoretical arguments to support the numerical outcomes are presented. An estimation of the degree of chaos is made by computing the Lyapunov Characteristic Exponents. Comparison with earlier work is also made.

preprint2012arXiv

The nature of orbits in a prolate elliptical galaxy model with a bulge and a dense nucleus

We study the transition from regular to chaotic motion in a prolate elliptical galaxy dynamical model with a bulge and a dense nucleus. Our numerical investigation shows that stars with angular momentum Lz less than or equal to a critical value Lzc, moving near the galactic plane, are scattered to higher z, when reaching the central region of the galaxy, thus displaying chaotic motion. An inverse square law relationship was found to exist between the radius of the bulge and the critical value Lzc of the angular momentum. On the other hand, a linear relationship exists between the mass of the nucleus and Lzc. The numerically obtained results are explained using theoretical arguments. Our study shows that there are connections between regular or chaotic motion and the physical parameters of the system, such as the star's angular momentum and mass, the scale length of the nucleus and the radius of the bulge. The results are compared with the outcomes of previous work.

preprint2012arXiv

Transition between order and chaos in a composite disk galaxy model with a massive nucleus and a dark matter halo

We investigate the transition from regular to chaotic motion in a composite galaxy model with a disk-halo, a massive dense nucleus and a dark halo component. We obtain relationships connecting the critical value of the mass of the nucleus or the critical value of the angular momentum, with the mass of the dark halo, where the transition from regular motion to chaos occurs. We also present 3D diagrams connecting the mass of nucleus the energy and the percentage of stars that can show chaotic motion. The fraction of the chaotic orbits observed in the phase plane, as a function of the mass of the dark halo is also computed. We use a semi-numerical method, that is a combination of theoretical and numerical procedure. The theoretical results obtained using the version 8.0 of the Mathematica package, while all the numerical calculations were made using a Bulirsch-Stoer FORTRAN routine in double precision. The results can be obtained in semi-numerical or numerical form and give good description for the connection of the physical quantities entering the model and the transition between regular and chaotic motion. We observe that the mass of the dark halo, the mass of the dense nucleus and the L_z component of the angular momentum, are important physical quantities, as they are linked to the regular or chaotic character of orbits in disk galaxies described by the model. Our numerical experiments suggest, that the amount of the dark matter plays an important role in disk galaxies represented by the model, as the mass of the halo affects, not only the regular or chaotic nature of motion but it is also connected with the existence of the different families of regular orbits. Comparison of the present results with earlier work is also presented.

preprint2012arXiv

Using the S(c) spectrum to distinguish between order and chaos in a 3D galactic potential

The regular and chaotic character of orbits is investigated in a 3D potential describing motion in the central parts of a barred galaxy. This potential is an extension in the 3D space of a 2D potential based on a family of figure-eight orbits, which was produced using the theory of the inverse problem. Starting from the results obtained from the 2D system, we proceed to locate the regions in the phase space of the 3D potential, producing regular or chaotic orbits. In order to obtain this we use a new dynamical parameter, the S(c) spectrum, which proves to be a useful and fast indicator in order to distinguish regular motion from chaos in 3D potentials. Comparison with other methods for detecting chaos is also discussed.