This article aims to investigate the points of equilibrium and the associated convergence basins in a seventh-order generalized H\'enon-Heiles potential. Using the well-known Newton-Raphson iterator we numerically locate the position of the points of equilibrium, while we also obtain their linear stability. Furthermore, we demonstrate how the two variable parameters, entering the generalized H\'enon-Heiles potential, affect the convergence dynamics of the system as well as the fractal degree of the basin diagrams. The fractal degree is derived by computing the (boundary) basin entropy as well as the uncertainty dimension.

Continuing the series of papers on a new model for a barred galaxy, we investigate the heteroclinic connections between the two normally hyperbolic invariant manifolds sitting over the two index-1 saddle points of the effective potential. The heteroclinic trajectories and the nearby periodic orbits of similar shape populate the bar region of the galaxy and a neighbourhood of its nucleus. Thereby we see a direct relation between the important structures of the interior region of the galaxy and the projection of the heteroclinic tangle into the position space. As a side result, we obtain a detailed picture of the primary heteroclinic intersection surface in the phase space.

In this paper, the near-optimal capture problem is numerically investigated to locate the optimal capture points around Charon using the Pluto-Charon planar circular restricted three-body problem (PCRTBP). Capture orbits are divided into the pre- and post-maneuver portions. In the pre-maneuver portion, the gravitational capture conditions are discussed by backward numerical integration. In the post-maneuver portion, the smaller alignment index (SALI) is applied to numerically investigate the orbital characters of long-term capture orbits around Charon. The initial conditions corresponding to three types of motion: (i) bounded (regular or chaotic), (ii) escaping and (iii) collisional are classified and studied. Combining results of the pre- and post-maneuver portions, the near-optimal capture method is presented to find optimal capture points by overlay figures, and then the corresponding capture orbits are constructed in the PCRTBP. The different results between the Pluto-Charon system and the Earth-Moon system are compared and analyzed. Our results obtained in this paper, including the optimal capture points and the corresponding capture orbits, could be applied in future space mission design.

The dynamics of a test particle in a non-spinning binary black hole system of equal masses is numerically investigated. The binary system is modeled in the context of the pseudo-Newtonian circular restricted three-body problem, such that the primaries are separated by a fixed distance and move in a circular orbit around each other. In particular, the Paczy\'{n}ski-Wiita potential is used for describing the gravitational field of the two non-Newtonian primaries. The orbital properties of the test particle are determined through the classification of the initial conditions of the orbits, using several values of the Jacobi constant, in the Hill's regions of possible motion. The initial conditions are classified into three main categories: (i) bounded, (ii) escaping and (iii) displaying close encounters. Using the smaller alignment index (SALI) chaos indicator, we further classify bounded orbits into regular, sticky or chaotic. To gain a complete view of the dynamics of the system, we define grids of initial conditions on different types of two-dimensional planes. The orbital structure of the configuration plane, along with the corresponding distributions of the escape and collision/close encounter times, allow us to observe the transition from the classical Newtonian to the pseudo-Newtonian regime. Our numerical results reveal a strong dependence of the properties of the considered basins with the Jacobi constant as well as with the Schwarzschild radius of the black holes.

The theory of the post-Newtonian (PN) planar circular restricted three-body problem is used for numerically investigating the orbital dynamics of a test particle (e.g., a comet, asteroid, meteor or spacecraft) in the planar Sun-Jupiter system with a scattering region around Jupiter. For determining the orbital properties of the test particle, we classify large sets of initial conditions of orbits for several values of the Jacobi constant in all possible Hill region configurations. The initial conditions are classified into three main categories: (i) bounded, (ii) escaping and (iii) collisional. Using the smaller alignment index chaos indicator (SALI), we further classify bounded orbits into regular, sticky or chaotic. In order to get a spherical view of the dynamics of the system, the grids of the initial conditions of the orbits are defined on different types of two-dimensional planes. We locate the different types of basins and we also relate them with the corresponding spatial distributions of the escape and collision time. Our thorough analysis exposes the high complexity of the orbital dynamics and exhibits an appreciable difference between the final states of the orbits in the classical and PN approaches. Furthermore, our numerical results reveal a strong dependence of the properties of the considered basins with the Jacobi constant, along with a remarkable presence of fractal basin boundaries. Our outcomes are compared with earlier ones, regarding other planetary systems.

We elucidate the escape properties of stars moving in the combined gravitational field of a binary system of two interacting dwarf spheroidal galaxies. A galaxy model of three degrees of freedom is adopted for describing the dynamical properties of the Hamiltonian system. All the numerical values of the involved parameters are chosen having in mind the real binary system of the dwarf spheroidal galaxies NGC 147 and NGC 185. We distinguish between bounded (regular, sticky or chaotic) and escaping motion by classifying initial conditions of orbits in several types of two dimensional planes, considering only unbounded motion for several energy levels. We analyze the orbital structure of all types of two dimensional planes of initial conditions by locating the basins of escape and also by measuring the corresponding escape time of the orbits. Furthermore, the properties of the normally hyperbolic invariant manifolds (NHIMs), located in the vicinity of the index-1 saddle points $L_1$, $L_2$, and $L_3$, are also investigated. These manifolds are of great importance, as they control the flow of stars (between the two galaxies and toward the exterior region) over the different saddle points. In addition, bifurcation diagrams of the Lyapunov periodic orbits as well as restrictions of the Poincar\'e map to the NHIMs are presented for revealing the dynamics in the neighbourhood of the saddle points. Comparison between the current outcomes and previous related results is also made.

The H\'enon-Heiles potential was first proposed as a simplified version of the gravitational potential experimented by a star in the presence of a galactic center. Currently, this system is considered a paradigm in dynamical systems because despite its simplicity exhibits a very complex dynamical behavior. In the present paper, we perform a series expansion up to the fifth-order of a potential with axial and reflection symmetries, which after some transformations, leads to a generalized H\'enon-Heiles potential. Such new system is analyzed qualitatively in both regimes of bounded and unbounded motion via the Poincar\'e sections method and plotting the exit basins. On the other hand, the quantitative analysis is performed through the Lyapunov exponents and the basin entropy, respectively. We find that in both regimes the chaoticity of the system decreases as long as the test particle energy gets far from the critical energy. Additionally, we may conclude that despite the inclusion of higher order terms in the series expansion, the new system shows wider zones of regularity (islands) than the ones present in the H\'enon-Heiles system.

The escape dynamics of the stars in a barred galaxy composed of a spherically symmetric central nucleus, a bar, a flat thin disk and a dark matter halo component is investigated by using a realistic three degrees of freedom (3-dof) dynamical model. Modern colour-coded diagrams are used for distinguishing between bounded and escaping motion. In addition, the smaller alignment index (SALI) method is deployed for determining the regular, sticky or chaotic nature of bounded orbits. We reveal the basins of escape corresponding to the escape through the two symmetrical escape channels around the Lagrange points $L_2$ and $L_3$ and also we relate them with the corresponding distribution of the escape times of the orbits. Furthermore, we demonstrate how the stable manifolds, around the index-1 saddle points, accurately define the fractal basin boundaries observed in the colour-coded diagrams. The development scenario of the fundamental vertical Lyapunov periodic orbit is thoroughly explored for obtaining a more complete view of the unfolding of the singular behaviour of the dynamics at the cusp values of the parameters. Finally, we examine how the combination of the most important parameters of the bar (such as the semi-major axis and the angular velocity) influences the observed stellar structures (rings and spirals) which are formed by escaping stars guided by the invariant manifolds near the saddle points.

We use a simple dynamical model in order to investigate the regular or chaotic character of orbits in a barred galaxy with a central, spherically symmetric, dense nucleus and a flat disk. In particular, we explore how the total orbital energy influences the overall orbital structure of the system, by computing in each case the percentage of regular, sticky and chaotic orbits. In an attempt to distinguish safely and with certainty between ordered and chaotic motion, we apply the Smaller ALingment Index (SALI) as a chaos detector to extensive samples of orbits obtained by integrating numerically the basic equations of motion as well as the variational equations. We integrate large sets of initial conditions of orbits in several types of two dimensional planes for better understanding of the orbital properties. Our numerical calculations suggest, that the value of the energy has a huge impact on the percentages of the orbits, thus indicating that a rotating barred galaxy is indeed a very interesting stellar quantity.

The escape mechanism of orbits in a star cluster rotating around its parent galaxy in a circular orbit is investigated. A three degrees of freedom model is used for describing the dynamical properties of the Hamiltonian system. The gravitational field of the star cluster is represented by a smooth and spherically symmetric Plummer potential. We distinguish between ordered and chaotic orbits as well as between trapped and escaping orbits, considering only unbounded motion for several energy levels. The Smaller Alignment Index (SALI) method is used for determining the regular or chaotic nature of the orbits. The basins of escape are located and they are also correlated with the corresponding escape time of the orbits. Areas of bounded regular or chaotic motion and basins of escape were found to coexist in the $(x,z)$ plane. The properties of the normally hyperbolic invariant manifolds (NHIMs), located in the vicinity of the index-1 Lagrange points $L_1$ and $L_2$, are also explored. These manifolds are of paramount importance as they control the flow of stars over the saddle points, while they also trigger the formation of tidal tails observed in star clusters. Bifurcation diagrams of the Lyapunov periodic orbits as well as restrictions of the Poincar\'e map to the NHIMs are deployed for elucidating the dynamics in the neighbourhood of the saddle points. The extended tidal tails, or tidal arms, formed by stars with low velocity which escape through the Lagrange points are monitored. The numerical results of this work are also compared with previous related work.

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\'e 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.

Locating the position of periodic orbits in galaxies is undoubtedly an issue of paramount importance. We reveal the position and the stability of periodic orbits of stars moving in the meridional plane $(R,z)$ of an axially symmetric galactic model with a disk, a spherical nucleus, and a biaxial dark matter halo component. In particular, we study how all the involved parameters of the dynamical system influence the position and the stability of all resonant families. To locate the position and measure the stability of periodic orbits we use a highly sensitive numerical code which is able to identify resonant periodic orbits of the type $n:m$. Two cases are studied for every parameter: (i) the case where the dark matter halo component is prolate and (ii) the case where an oblate dark matter halo is present. Our numerical exploration reveals that all the dynamical quantities affect, more or less, the position and the stability of the periodic orbits. It is shown that the mass of the nucleus, the mass of the disk, the halo flattening parameter, the scale length of the halo, the angular momentum, and the total orbital energy are the most influential quantities, while the effect of all other parameters is much weaker.

We consider the planar circular equilateral restricted four body-problem where a test particle of infinitesimal mass is moving under the gravitational attraction of three primary bodies which move on circular orbits around their common center of gravity, such that their configuration is always an equilateral triangle. The case where all three primaries have equal masses is numerically investigated. A thorough numerical analysis takes place in the configuration $(x,y)$ as well as in the $(x,C)$ space in which we classify initial conditions of orbits into four main categories: (i) bounded regular orbits, (ii) trapped chaotic orbits, (iii) escaping orbits and (iv) collision orbits. Interpreting the collision motion as leaking in the phase space we related our results to both chaotic scattering and the theory of leaking Hamiltonian systems. We successfully located the escape and the collision basins and we managed to correlate them with the corresponding escape and collision times of orbits. We hope our contribution to be useful for a further understanding of the escape and collision properties of motion in this interesting dynamical system.

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 $\mu$, 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 $\mu$, $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.

The ordered or chaotic character of orbits of stars moving in the meridional $(R,z)$ plane of an analytic axisymmetric time-independent disk galaxy model with an additional spherically symmetric central nucleus is investigated. Our aim is to determine how the total energy influences the orbital structure of the galaxy. For this purpose we monitor how the percentage of chaotic orbits as well as the rates of orbits composing the main regular families evolve as a function of the value of the energy. In order to distinguish with certainty between chaotic and ordered motion we use the SALI method in extensive sets of initial conditions of orbits. Moreover, a spectral method is applied for identifying the various regular families and also for recognizing the secondary resonances that bifurcate from them. Our numerical computations suggest that for low energy levels the observed amount of chaos is high and the orbital content is rather poor, while for high energy levels, corresponding to global motion, regular motion dominates and many secondary higher resonances emerge. We also compared our results with previous related work.

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.

The orbital dynamics of a spacecraft, or a comet, or an asteroid in the Earth-Moon system in a scattering region around the Moon using the three dimensional version of the circular restricted three-body problem is numerically investigated. The test particle can move in bounded orbits around the Moon or escape through the openings around the Lagrange points $L_1$ and $L_2$ or even collide with the surface of the Moon. We explore in detail the first four of the five possible Hill's regions configurations depending on the value of the Jacobi constant which is of course related with the total orbital energy. We conduct a thorough numerical analysis on the phase space mixing by classifying initial conditions of orbits in several two-dimensional types of planes and distinguishing between four types of motion: (i) ordered bounded, (ii) trapped chaotic, (iii) escaping and (iv) collisional. In particular, we locate the different basins and we relate them with the corresponding spatial distributions of the escape and collision times. Our outcomes reveal the high complexity of this planetary system. Furthermore, the numerical analysis suggests a strong dependence of the properties of the considered basins with both the total orbital energy and the initial value of the $z$ coordinate, with a remarkable presence of fractal basin boundaries along all the regimes. Our results are compared with earlier ones regarding the planar version of the Earth-Moon system.

The escape dynamics in an analytical gravitational model which describes the motion of stars in a binary system of interacting dwarf spheroidal galaxies is investigated in detail. We conduct a numerical analysis distinguishing between regular and chaotic orbits as well as between trapped and escaping orbits, considering only unbounded motion for several energy levels. In order to distinguish safely and with certainty between ordered and chaotic motion, we apply the Smaller ALingment Index (SALI) method. It is of particular interest to locate the escape basins through the openings around the collinear Lagrangian points $L_1$ and $L_2$ and relate them with the corresponding spatial distribution of the escape times of the orbits. Our exploration takes place both in the configuration $(x,y)$ and in the phase $(x,\dot{x})$ space in order to elucidate the escape process as well as the overall orbital properties of the galactic system. Our numerical analysis reveals the strong dependence of the properties of the considered escape basins with the total orbital energy, with a remarkable presence of fractal basin boundaries along all the escape regimes. It was also observed, that for energy levels close to the critical escape energy the escape rates of orbits are large, while for much higher values of energy most of the orbits have low escape periods or they escape immediately to infinity. We hope our outcomes to be useful for a further understanding of the escape mechanism in binary galaxy models.

We investigate in detail the escape dynamics in an analytical gravitational model which describes the motion of stars in a quasar galaxy with a disk and a massive nucleus. We conduct a thorough numerical analysis distinguishing between regular and chaotic orbits as well as between trapped and escaping orbits, considering only unbounded motion for several energy levels. In order to distinguish safely and with certainty between ordered and chaotic motion we apply the Smaller ALingment Index (SALI) method. It is of particular interest to locate the escape basins through the openings around the collinear Lagrangian points $L_1$ and $L_2$ and relate them with the corresponding spatial distribution of the escape times of the orbits. Our exploration takes place both in the configuration $(x,y)$ and in the phase $(x,\dot{x})$ space in order to elucidate the escape process as well as the overall orbital properties of the galactic system. Our numerical analysis reveals the strong dependence of the properties of the considered escape basins with the total orbital energy, with a remarkable presence of fractal basin boundaries along all the escape regimes. We hope our outcomes to be useful for a further understanding of the escape mechanism in active galaxy models.

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.

The case of the planar circular photogravitational restricted three-body problem where the more massive primary is an emitter of radiation is numerically investigated. A thorough numerical analysis takes place in the configuration $(x,y)$ and the $(x,C)$ space in which we classify initial conditions of orbits into three main categories: (i) bounded, (ii) escaping and (iii) collisional. Our results reveal that the radiation pressure factor has a huge impact on the character of orbits. Interpreting the collisional motion as leaking in the phase space we related our results to both chaotic scattering and the theory of leaking Hamiltonian systems. We successfully located the escape as well as the collisional basins and we managed to correlate them with the corresponding escape and collision times. We hope our contribution to be useful for a further understanding of the escape and collision properties of motion in this interesting version of the restricted three-body problem.

We numerically investigate the orbital dynamics of a spacecraft, or a comet, or an asteroid in the Pluto-Charon system in a scattering region around Charon using the planar circular restricted three-body problem. The test particle can move in bounded orbits around Charon or escape through the necks around the Lagrangian points $L_1$ and $L_2$ or even collide with the surface of Charon. We explore four of the five possible Hill's regions configurations depending on the value of the Jacobi constant which is of course related with the total orbital energy. We conduct a thorough numerical analysis on the phase space mixing by classifying initial conditions of orbits and distinguishing between three types of motion: (i) bounded, (ii) escaping and (iii) collisional. In particular, we locate the different basins and we relate them with the corresponding spatial distributions of the escape and collision times. Our results reveal the high complexity of this planetary system. Furthermore, the numerical analysis shows a strong dependence of the properties of the considered basins with the total orbital energy, with a remarkable presence of fractal basin boundaries along all the regimes. Our results are compared with earlier ones regarding the Saturn-Titan planetary system.

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\'e 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.

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\'e 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.

We use the planar circular restricted three-body problem in order to numerically investigate the orbital dynamics of orbits of a spacecraft, or a comet, or an asteroid in the Saturn-Titan system in a scattering region around the Titan. The orbits can escape through the necks around the Lagrangian points $L_1$ and $L_2$ or collide with the surface of the Titan. We explore all the four possible Hill's regions depending on the value of the Jacobi constant. We conduct a thorough numerical analysis on the phase space mixing by classifying initial conditions of orbits and distinguishing between three types of motion: (i) bounded, (ii) escaping and (iii) collisional. In particular, we locate the different basins and we relate them with the corresponding spatial distributions of the escape and crash times. Our results reveal the high complexity of this planetary system. Furthermore, the numerical analysis shows a strong dependence of the properties of the considered basins with the total orbital energy, with a remarkable presence of fractal basin boundaries along all the regimes. We hope our contribution to be useful in both space mission design and study of planetary systems.

The aim of this work is to compare the orbital dynamics in three different models describing the properties of a star cluster rotating around its parent galaxy in a circular orbit. In particular, we use the isochrone and the Hernquist potentials to model the spherically symmetric star cluster and we compare our results with the corresponding ones of a previous work in which the Plummer model was applied for the same purpose. Our analysis takes place both in the configuration $(x,y)$ and in the phase $(x,\dot{x})$ space in order to elucidate the escape process as well as the overall orbital properties of the tidally limited star cluster. We restrict our investigation into two dimensions and we conduct a thorough numerical analysis distinguishing between ordered and chaotic orbits as well as between trapped and escaping orbits, considering only unbounded motion for several energy levels above the critical escape energy. It is of particular interest to determine the escape basins towards the two exit channels (near the Lagrangian points $L_1$ and $L_2$) and relate them with the corresponding escape times of the orbits.

We numerically investigate the case of the planar circular restricted three-body problem where the more massive primary is an oblate spheroid. A thorough numerical analysis takes place in the configuration $(x,y)$ and the $(x,E)$ space in which we classify initial conditions of orbits into three categories: (i) bounded, (ii) escaping and (iii) collisional. Our results reveal that the oblateness coefficient has a huge impact on the character of orbits. Interpreting the collisional motion as leaking in the phase space we related our results to both chaotic scattering and the theory of leaking Hamiltonian systems. We successfully located the escape as well as the collisional basins and we managed to correlate them with the corresponding escape and collision times. We hope our contribution to be useful for a further understanding of the escape and collision properties of motion in this interesting version of the restricted three-body problem.

The regular or chaotic character of orbits of stars moving in the meridional plane (R,z) of an axially symmetric elliptical galaxy with a dense, massive spherical nucleus and a dark matter halo component is under investigation. In particular, we explore how the flattening of an elliptical galaxy influences the overall orbital structure of the system, by computing in each case the percentage of chaotic orbits, as well as the percentages of orbits composing the main regular families. In an attempt to discriminate safely and with certainty between regular and chaotic motion, we use the Smaller ALingment Index (SALI) method to extensive samples of orbits obtained by integrating numerically the basic equations of motion as well as the variational equations. In addition, a technique which is based mainly on the field of spectral dynamics that utilizes the Fourier transform of the time series of each coordinate is used for classifying the regular orbits into different families and also to recognize the secondary resonances that usually bifurcate from them. Three cases are considered in our work: (i) the case where the elliptical galaxy is prolate (ii) the case where a spherically symmetric elliptical galaxy is present and (iii) the case where the elliptical galaxy has an oblate shape. Comparison between the current results and early related work is also made.

The escape dynamics in a simple analytical gravitational model which describes the motion of stars in a Seyfert galaxy is investigated in detail. We conduct a thorough numerical analysis distinguishing between regular and chaotic orbits as well as between trapped and escaping orbits, considering only unbounded motion for several energy levels. In order to distinguish safely and with certainty between ordered and chaotic motion, we apply the Smaller ALingment Index (SALI) method. It is of particular interest to locate the escape basins through the openings around the collinear Lagrangian points $L_1$ and $L_2$ and relate them with the corresponding spatial distribution of the escape times of the orbits. Our exploration takes place both in the physical $(x,y)$ and in the phase $(x,\dot{x})$ space in order to elucidate the escape process as well as the overall orbital properties of the galactic system. Our numerical analysis reveals the strong dependence of the properties of the considered escape basins with the total orbital energy, with a remarkable presence of fractal basin boundaries along all the escape regimes. It was also observed, that for energy levels close to the critical escape energy the escape rates of orbits are large, while for much higher values of energy most of the orbits have low escape periods or they escape immediately to infinity. We also present evidence obtained through numerical simulations that our model can describe the formation and the evolution of the observed spiral structure in Seyfert galaxies. We hope our outcomes to be useful for a further understanding of the escape mechanism in galaxies with active nuclei.

We explore the regular or chaotic nature of orbits of stars moving in the meridional (R,z) plane of an axially symmetric time-dependent disk galaxy model with a central, spherically symmetric nucleus. In particular, mass is linearly transported from the disk to the galactic nucleus, in order to mimic, in a way, the case of self-consistent interactions of an actual N-body simulation. We thus try to unveil the influence of this mass transportation on the different families of orbits of stars by monitoring how the percentage of chaotic orbits, as well as the percentages of orbits of the main regular resonant families, evolve as the galaxy develops a dense and massive nucleus in its core. The SALI method is applied to samples of orbits in order to distinguish safely between ordered and chaotic motion. In addition, a method based on the concept of spectral dynamics is used for identifying the various families of regular orbits and also for recognizing the secondary resonances that bifurcate from them. Our computations strongly suggest that the amount of the observed chaos is substantially increased as the nucleus becomes more massive. Furthermore, extensive numerical calculations indicate that there are orbits which change their nature from regular to chaotic and vice versa and also orbits which maintain their orbital character during the galactic evolution. The present outcomes are compared to earlier related work.

The aim of this work is to explore the escape process of three-dimensional orbits in a star cluster rotating around its parent galaxy in a circular orbit. The gravitational field of the cluster is represented by a smooth, spherically symmetric Plummer potential, while the tidal approximation was used to model the steady tidal field of the galaxy. We conduct a thorough numerical analysis distinguishing between regular and chaotic orbits as well as between trapped and escaping orbits, considering only unbounded motion for several energy levels. It is of particular interest to locate the escape basins towards the two exit channels and relate them with the corresponding escape times of the orbits. For this purpose, we split our investigation into three cases depending on the initial value of the $z$ coordinate which was used for launching the stars. The most noticeable finding is that the majority of stars initiated very close to the primary $(x,y)$ plane move in chaotic orbits and they remain trapped for vast time intervals, while orbits with relatively high values of $z_0$ on the other hand, form well-defined basins of escape. It was also observed, that for energy levels close to the critical escape energy the escape rates of orbits are large, while for much higher values of energy most of the orbits have low escape periods or they escape immediately to infinity. We hope our outcomes to be useful for a further understanding of the dissolution process and the escape mechanism in open star clusters.

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.

We determine the character of orbits of stars moving in the meridional plane $(R,z)$ of an axially symmetric time-independent disk galaxy model with a spherical central nucleus. In particular, we try to reveal the influence of the value of the angular momentum on the different families of orbits of stars, by monitoring how the percentage of chaotic orbits, as well as the percentages of orbits of the main regular resonant families evolve when angular momentum varies. The smaller alignment index (SALI) was computed by numerically integrating the equations of motion as well as the variational equations to extensive samples of orbits in order to distinguish safely between ordered and chaotic motion. In addition, 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. Our investigation takes place both in the physical $(R,z)$ and the phase $(R,\dot{R})$ space for a better understanding of the orbital properties of the system. Our numerical computations reveal that low angular momentum stars are most likely to move in chaotic orbits, while on the other hand, the vast majority of high angular momentum stars perform regular orbits.

The main objective of this work is to determine the character of orbits of stars moving in the meridional $(R,z)$ plane of an axially symmetric time-independent disk galaxy model with a central massive nucleus and an additional spherical dark matter halo component. In particular, we try to reveal the influence of the scale length of the dark matter halo on the different families of orbits of stars, by monitoring how the percentage of chaotic orbits, as well as the percentages of orbits of the main regular resonant families evolve when this parameter varies. The smaller alignment index (SALI) was computed by numerically integrating the equations of motion as well as the variational equations to extensive samples of orbits in order to distinguish safely between ordered and chaotic motion. In addition, 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. Our numerical computations reveal that when the dark matter halo is highly concentrated, that is when the scale length has low values the vast majority of star orbits move in regular orbits, while on the other hand in less concentrated dark matter halos the percentage of chaos increases significantly. We also compared our results with early related work.

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.

We investigate the regular or chaotic nature of orbits of stars moving in the meridional plane $(R,z)$ of an axially symmetric galactic model with a flat disk and a central, non spherical and massive nucleus. In particular, we study the influence of the flattening parameter of the central nucleus on the nature of orbits, by computing in each case the percentage of chaotic orbits, as well as the percentages of orbits of the main regular families. In an attempt to maximize the accuracy of our results upon distinguishing between regular and chaotic motion, we use both the Fast Lyapunov Indicator (FLI) and the Smaller ALingment Index (SALI) methods to extensive samples of orbits obtained by integrating numerically the equations of motion as well as the variational equations. Moreover, a technique which is based mainly on the field of spectral dynamics that utilizes the Fourier transform of the time series of each coordinate is used for identifying the various families of regular orbits and also to recognize the secondary resonances that bifurcate from them. Varying the value of the flattening parameter, we study three different cases: (i) the case where we have a prolate nucleus (ii) the case where the central nucleus is spherical and (iii) the case where an oblate massive nucleus is present. Furthermore, we present some additional findings regarding the reliability of short time (fast) chaos indicators, as well as a new method to define the threshold between chaos and regularity for both FLI and SALI, by using them simultaneously. Comparison with early related work is also made.

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.

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.

In this article, we present a galactic gravitational model of three degrees of freedom (3D), in order to study and reveal the character of the orbits of the stars, in a binary stellar system composed of a primary quiet or active galaxy and a small satellite companion galaxy. Our main dynamical analysis will be focused on the behavior of the primary galaxy. We investigate in detail the regular or chaotic nature of motion, in two different cases: (i) the time-independent model in both 2D and 3D dynamical systems and (ii) the time-evolving 3D model. Our numerical calculations indicate that the percentage of the chaotic orbits increases when the primary galaxy has a dense and massive nucleus. The presence of the dense galactic core also increases the stellar velocities near the center of the galaxy. Moreover, for small values of the distance R between the two bodies, low-energy stars display chaotic motion, near the central region of the galaxy, while for larger values of the distance R, the motion in active galaxies is entirely regular for low-energy stars. Our simulations suggest that in galaxies with a satellite companion, the chaotic nature of motion is not only a result of the galactic interaction between the primary galaxy and its companion, but also a result caused by the presence of the dense nucleus in the core of the primary galaxy. Theoretical arguments are presented in order to support and interpret the numerically derived outcomes. Furthermore, we follow the 3D evolution of the primary galaxy, when mass is transported adiabatically from the disk to the nucleus. Our numerical results are in satisfactory agreement with observational data obtained from the M51-type binary stellar systems. A comparison between the present research and similar and earlier work is also made.

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.

An axially symmetric galactic gravitational model composed of a dense, massive and spherical nucleus with an additional dark matter halo component was used, to distinguish between the regular and chaotic character of orbits of stars that move in the meridional plane (R,z). We investigated two different cases: (i) a flat-disk galaxy (ii) an elliptical galaxy. It is of particular interest to reveal how the portion of the dark matter inside the main body of the galaxy influences the ordered or chaotic nature of motion. Varying the ratio of dark matter to stellar mass, we monitored the evolution not only of the percentage of chaotic orbits, but also of the percentages of orbits that compose the main regular resonant families, by classifying regular orbits into different families. Moreover we tried, to reveal how the starting position of the parent periodic orbits of each regular family changes with respect to the fractional portion of dark matter. We compared our results with previous similar work.

We investigate the regular or chaotic nature of orbits of stars moving in the meridional plane (R,z) of an axially symmetric galactic model with a dense, massive spherical nucleus and a dark matter halo component. In particular, we study the influence of the fractional portion of the dark matter, by computing in each case the percentage of chaotic orbits, as well as the percentages of orbits of the main regular resonant families. In an attempt to distinguish between regular and chaotic motion, we use the Fast Lyapunov Indicator (FLI) method to extensive samples of orbits obtained by integrating numerically the equations of motion as well as the variational equations. Furthermore, a technique which is based mainly on the field of spectral dynamics that utilizes the Fourier transform of the time series of each coordinate is used for identifying the various families of regular orbits and also to recognize the secondary resonances that bifurcate from them. Two cases are studied in our work: (i) the case where we have a disk galaxy model (ii) the case where our model represents an elliptical galaxy. Comparison with early related work is also made.

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.

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.

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.

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.

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.

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.

In the present article, we investigate the behavior of orbits in a time independent axially symmetric galactic type potential. This dynamical model can be considered to describe the motion in the central parts of a galaxy, for values of energies larger than the energy of escape. We use the classical method of the surface of section, in order to visualize and interpret the structure of the phase space of the dynamical system. Moreover, the Lyapunov Characteristic Exponent (LCE), is used in order to make an estimation of the degree of the chaoticity of the orbits in our galactic model. Our numerical calculations suggest that in this galactic type potential, there are two kinds of orbits: (i) escaping orbits and (ii) trapped orbits which do not escape at all. Furthermore, a large number of orbits of the dynamical system, display chaotic motion. Among the chaotic orbits, there are orbits that escape fast and also orbits that remain trapped for vast time intervals. When the value of the test particle's energy exceeds slightly the energy of escape, the amount of the trapped regular orbits increases, as the the value of the angular momentum increases. Therefore, the extent of the chaotic regions observed in the phase plane decreases as the value of the energy increases. Moreover, we calculate the average value of the escape period of the chaotic orbits and we try to correlate it with the value of the energy and also with the maximum value of the z component of the orbits. In addition, we find that the value of the LCE corresponding to each chaotic region, for different values of the energy, increases exponentially as the value of the energy increases. Some theoretical arguments in order to support the numerically obtained outcomes are presented.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

A 3D dynamical model with a quasi-homogeneous core and a disk component is used for the chaos control in the central parts of elliptical galaxy. Numerical experiments in the 2D system show a very complicated phase plane with a large chaotic sea, considerable sticky layers and a large number of islands, produced by secondary resonances. When the mass of the disk increases, the chaotic regions decrease gradually, and, finally, a new phase plane with only regular orbits appears. This evolution indicates that disks in elliptical galaxies can act as the chaos controllers. Starting from the results obtained in the 2D system, we locate the regions in the phase space of the 3D system, producing regular and chaotic orbits. For this we introduce and use a new dynamical parameter, the S(w) spectrum, which proves to be useful as a fast indicator and allows us to distinguish the regular motion from chaos in the 3D potentials. Other methods for detecting chaos are also discussed.

A 3D dynamical model is used to study the motion in the central parts of an elliptical galaxy, hosting a massive and dense nucleus. Our aim is to investigate the regular or chaotic character of the motion, with emphasis in the different chaotic components, as well as the sticky regions of the dynamical system. In order to define the character of the motion in the 2D system, we use the classical method of the Poincare x-px phase plane, the Lyapunov Characteristic Exponent (LCE) and the dynamical parameter - the S(c) spectrum. Then the results obtained from the 2D system are used to investigate the properties of the 3D system. For this, we introduce and use a new dynamical parameter, the S(k) spectrum, which proves to be a very reliable and fast method to detect the islandic motion and the evolution of the sticky orbits in the 3D system. Numerical experiments conducted by the new S(k) spectrum suggest that the different chaotic components in the 3D system do not interact for time intervals much larger than the age of the galaxy. The results indicate that the different sticky regions do not lead to a unified chaotic sea. Thus, the behavior of the 3D sticky orbits differs from that observed in the 2D system. Furthermore, the 3D motion near the center of a triaxial elliptical galaxy seems to be very complicated, displaying several families of resonant orbits, different chaotic components and sticky regions, while only a small fraction of orbits is regular. The comparison with earlier results reveals the importance of the conception of the new dynamical spectrum.