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

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.

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

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.

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.

We investigate the properties of motion in a map model derived from a galactic Hamiltonian made up of perturbed elliptic oscillators. The phase space portrait is obtained in all three different cases using the map and numerical integration of the equations of motion. Our numerical calculations suggest that the map describes very well the properties of regular motion in all three cases. There are cases where the map fails to describe chaotic motion and cases where the map describes satisfactorily the chaotic phase plane. We also compare the Lyapunov characteristic number and the spectra of orbits derived by the map and numerical integration in each case. The agreement of the outcomes is satisfactory.