We investigate regular and chaotic two-dimensional (2D) and three-dimensional (3D) orbits of stars in models of a galactic potential consisting in a disk, a halo and a bar, to find the origin of boxy components, which are part of the bar or (almost) the bar itself. Our models originate in snapshots of an N-body simulation, which develops a strong bar. We consider three snapshots of the simulation and for the orbital study we treat each snapshot independently, as an autonomous Hamiltonian system. The calculated corotation-to-bar-length ratios indicate that in all three cases the bar rotates slowly, while the orientation of the orbits of the main family of periodic orbits changes along its characteristic. We characterize the orbits as regular, sticky, or chaotic after integrating them for a 10 Gyr period by using the GALI$_2$ index. Boxiness in the equatorial plane is associated either with quasi-periodic orbits in the outer parts of stability islands, or with sticky orbits around them, which can be found in a large range of energies. We indicate the location of such orbits in diagrams, which include the characteristic of the main family. They are always found about the transition region from order to chaos. By perturbing such orbits in the vertical direction we find a class of 3D non-periodic orbits, which have boxy projections both in their face-on and side-on views.

Studies of dynamical stability (chaotic versus regular motion) in galactic dynamics often rely on static analytical models of the total gravitational potential. Potentials based upon self-consistent N-body simulations offer more realistic models, fully incorporating the time-dependent nature of the systems. Here we aim at analysing the fractions of chaotic motion within different morphological components of the galaxy. We wish to investigate how the presence of chaotic orbits evolves with time, and how their spatial distribution is associated with morphological features of the galaxy. We employ a time-dependent analytical potential model that was derived from an N-body simulation of a strongly barred galaxy. With this analytical potential we may follow the dynamical evolution of ensembles of orbits. Using the Generalized Alignment Index (GALI) chaos detection method, we study the fraction of chaotic orbits, sampling the dynamics of both the stellar disc and of the dark matter halo. Within the stellar disc, the global trend is for chaotic motion to decrease in time, specially in the region of the bar. We scrutinized the different changes of regime during the evolution (orbits that are permanently chaotic, permanently regular, those that begin regular and end chaotic, and those that begin chaotic and end regular), tracing the types of orbits back to their common origins. Within the dark matter halo, chaotic motion also decreases globally in time. The inner halo (r < 5 kpc) is where most chaotic orbits are found and it is the only region where chaotic orbits outnumber regular orbits, in the early evolution.

We provide a concise presentation of the Smaller (SALI) and the Generalized Alignment Index (GALI) methods of chaos detection. These are efficient chaos indicators based on the evolution of two or more, initially distinct, deviation vectors from the studied orbit. After explaining the motivation behind the introduction of these indices, we sum up the behaviors they exhibit for regular and chaotic motion, as well as for stable and unstable periodic orbits, focusing mainly on finite-dimensional conservative systems: autonomous Hamiltonian models and symplectic maps. We emphasize the advantages of these methods in studying the global dynamics of a system, as well as their ability to identify regular motion on low dimensional tori. Finally we discuss several applications of these indices to problems originating from different scientific fields like celestial mechanics, galactic dynamics, accelerator physics and condensed matter physics.

Self-consistent N-body simulations are efficient tools to study galactic dynamics. However, using them to study individual trajectories (or ensembles) in detail can be challenging. Such orbital studies are important to shed light on global phase space properties, which are the underlying cause of observed structures. The potentials needed to describe self-consistent models are time-dependent. Here, we aim to investigate dynamical properties (regular/chaotic motion) of a non-autonomous galactic system, whose time-dependent potential adequately mimics certain realistic trends arising from N-body barred galaxy simulations. We construct a fully time-dependent analytical potential, modeling the gravitational potentials of disc, bar and dark matter halo, whose time-dependent parameters are derived from a simulation. We study the dynamical stability of its reduced time-independent 2-degrees of freedom model, charting the different islands of stability associated with certain orbital morphologies and detecting the chaotic and regular regions. In the full 3-degrees of freedom time-dependent case, we show representative trajectories experiencing typical dynamical behaviours, i.e., interplay between regular and chaotic motion for different epochs. Finally, we study its underlying global dynamical transitions, estimating fractions of (un)stable motion of an ensemble of initial conditions taken from the simulation. For such an ensemble, the fraction of regular motion increases with time.

We study the distinction and quantification of chaotic and regular motion in a time-dependent Hamiltonian barred galaxy model. Recently, a strong correlation was found between the strength of the bar and the presence of chaotic motion in this system, as models with relatively strong bars were shown to exhibit stronger chaotic behavior compared to those having a weaker bar component. Here, we attempt to further explore this connection by studying the interplay between chaotic and regular behavior of star orbits when the parameters of the model evolve in time. This happens for example when one introduces linear time dependence in the mass parameters of the model to mimic, in some general sense, the effect of self-consistent interactions of the actual N-body problem. We thus observe, in this simple time-dependent model also, that the increase of the bar's mass leads to an increase of the system's chaoticity. We propose a new way of using the Generalized Alignment Index (GALI) method as a reliable criterion to estimate the relative fraction of chaotic vs. regular orbits in such time-dependent potentials, which proves to be much more efficient than the computation of Lyapunov exponents. In particular, GALI is able to capture subtle changes in the nature of an orbit (or ensemble of orbits) even for relatively small time intervals, which makes it ideal for detecting dynamical transitions in time-dependent systems.

We use probability density functions (pdfs) of sums of orbit coordinates, over time intervals of the order of one Hubble time, to distinguish weakly from strongly chaotic orbits in a barred galaxy model. We find that, in the weakly chaotic case, quasi-stationary states arise, whose pdfs are well approximated by $q$-Gaussian functions (with $1<q<3$), while strong chaos is identified by pdfs which quickly tend to Gaussians ($q=1$). Typical examples of weakly chaotic orbits are those that "stick" to islands of ordered motion. Their presence in rotating galaxy models has been investigated thoroughly in recent years due of their ability to support galaxy structures for relatively long time scales. In this paper, we demonstrate, on specific orbits of 2 and 3 degree of freedom barred galaxy models, that the proposed statistical approach can distinguish weakly from strongly chaotic motion accurately and efficiently, especially in cases where Lyapunov exponents and other local dynamic indicators appear to be inconclusive.

As originally formulated, the Generalized Alignment Index (GALI) method of chaos detection has so far been applied to distinguish quasiperiodic from chaotic motion in conservative nonlinear dynamical systems. In this paper we extend its realm of applicability by using it to investigate the local dynamics of periodic orbits. We show theoretically and verify numerically that for stable periodic orbits the GALIs tend to zero following particular power laws for Hamiltonian flows, while they fluctuate around non-zero values for symplectic maps. By comparison, the GALIs of unstable periodic orbits tend exponentially to zero, both for flows and maps. We also apply the GALIs for investigating the dynamics in the neighborhood of periodic orbits, and show that for chaotic solutions influenced by the homoclinic tangle of unstable periodic orbits, the GALIs can exhibit a remarkable oscillatory behavior during which their amplitudes change by many orders of magnitude. Finally, we use the GALI method to elucidate further the connection between the dynamics of Hamiltonian flows and symplectic maps. In particular, we show that, using for the computation of GALIs the components of deviation vectors orthogonal to the direction of motion, the indices of stable periodic orbits behave for flows as they do for maps.

The distinction between chaotic and regular behavior of orbits in galactic models is an important issue and can help our understanding of galactic dynamical evolution. In this paper, we deal with this issue by applying the techniques of the Smaller (and Generalized) ALingment Indices, SALI (and GALI), to extensive samples of orbits obtained by integrating numerically the equations of motion in a barred galaxy potential. We estimate first the fraction of chaotic and regular orbits for the two-degree-of-freedom (DOF) case (where the galaxy extends only in the (x,y)-space) and show that it is a non-monotonic function of the energy. For the three DOF extension of this model (in the z-direction), we give similar estimates, both by exploring different sets of initial conditions and by varying the model parameters, like the mass, size and pattern speed of the bar. We find that regular motion is more abundant at small radial distances from the center of the galaxy, where the relative non-axisymmetric forcing is relatively weak, and at small distances from the equatorial plane, where trapping around the stable periodic orbits is important. We also find that the variation of the bar pattern speed, within a realistic range of values, does not affect much the phase space's fraction of regular and chaotic motions. Using different sets of initial conditions, we show that chaotic motion is dominant in galaxy models whose bar component is more massive, while models with a fatter or thicker bar present generally more regular behavior. Finally, we find that the fraction of orbits that are chaotic correlates strongly with the bar strength.

We study the dynamics of 2D and 3D barred galaxy analytical models, focusing on the distinction between regular and chaotic orbits with the help of the Smaller ALigment Index (SALI), a very powerful tool for this kind of problems. We present briefly the method and we calculate the fraction of chaotic and regular orbits in several cases. In the 2D model, taking initial conditions on a Poincar\'{e} $(y,p_y)$ surface of section, we determine the fraction of regular and chaotic orbits. In the 3D model, choosing initial conditions on a cartesian grid in a region of the $(x, z, p_y)$ space, which in coordinate space covers the inner disc, we find how the fraction of regular orbits changes as a function of the Jacobi constant. Finally, we outline that regions near the $(x,y)$ plane are populated mainly by regular orbits. The same is true for regions that lie either near to the galactic center, or at larger relatively distances from it.

We use the Smaller ALignment Index (SALI) method of chaos detection, to study the global dynamics of conservative dynamical systems described by differential or difference equations. In particular, we consider the well--known 2 and 4--dimensional symplectic standard map, as well as an autonomous Hamiltonian system of 2 and 3 degrees of freedom describing the motion of a star in models of barred galaxies. The application of SALI helped us to compute rapidly and accurately the percentage of regular and chaotic motion for particular values of the parameters of these systems. We were also able to perform a computationally efficient determination of the dependence of these percentages on the variation of several parameters of the studied models.

The chaotic or ordered character of orbits in galactic models is an important issue, since it can influence dynamical evolution. This distinction can be achieved with the help of the Smaller Alingment Index - (SALI). We describe here briefly this method and its advantages. Then we apply it to that case of 2D and 3D barred galaxy potentials. In particular, we find the fraction of chaotic and ordered orbits in such potentials and present how this fraction changes when the main parameters of the model are varied. For this, we consider models with different bar mass, bar thickness or pattern speed. Varying only one parameter at a time, we find that bars that are more massive, or thinner, or faster, have a larger fraction of chaotic orbits.

A very important issue in the area of galactic dynamics is the detection of chaotic and ordered motion inside galaxies. In order to achieve this target, we use the Smaller ALignment Index (SALI) method, which is a very suitable tool for this kind of problems. Here, we apply this index to 3D barred galaxy potentials and we present some results on the chaotic behavior of the model when its main parameters vary.

The dynamical evolution of barred galaxies depends crucially on the fraction and their spacial distribution of chaotic orbits in them. In order to distinguish between the two kinds of orbits, we use the Smaller Alignment Index (SALI) method, a very powerful method which can be applied to problems of galactic dynamics. Using model potentials, and taking into account the full 3D distribution of matter, we discuss how the distribution of chaotic orbits depends on the main model parameters, like the mass of the various components and the bar axial ratio.