We show that the Lyapunov exponents of a periodic orbit can be easily obtained from the eigenvalues of the monodromy matrix. It turns out that the Lyapunov exponents of simply stable periodic orbits are all zero, simply unstable periodic orbits have only one positive Lyapunov exponent, doubly unstable periodic orbits have two different positive Lyapunov exponents and the two positive Lyapunov exponents of complex unstable periodic orbits are equal. We present a numerical example for periodic orbits in a realistic galactic potential. Moreover, the center manifold theorem allowed us to show that stable, simply unstable and doubly unstable periodic orbits are the mothers of families of, respectively, regular, partially and fully chaotic orbits in their neighbourhood.

Observations show that stellar streams originating in satellite dwarf galaxies are frequent in the Universe. While such events are predicted by theory, it is not clear how many of the streams that are generated are washed out afterwards to the point in which it is imposible to detect them. Here we study how these diffusion times are affected by the fact that typical gravitational potentials of the host galaxies can sustain chaotic orbits. We do this by comparing the behaviour of simulated stellar streams that reside in chaotic or non-chaotic regions of the phase-space. We find that chaos does reduce the time interval in which streams can be detected. By analyzing detectability criteria in configuration and velocity space, we find that the impact of these results on the observations depends on the quality of both the data and the underlying stellar halo model. For all the stellar streams, we obtain a similar upper limit to the detectable mass.

The mid-transit times of an exoplanet may be non-periodic. The variations in the timing of the transits with respect to a single period, that is, the transit timing variations (TTVs), can sometimes be attributed to perturbations by other exoplanets present in the system, which may or may not transit the star. Our aim is to compute the mass and the six orbital elements of an non-transiting exoplanet, given only the central times of transit of the transiting body. We also aim to recover the mass of the star and the mass and orbital elements of the transiting exoplanet, suitably modified in order to decrease the deviation between the observed and the computed transit times by as much as possible. We have applied our method, based on a genetic algorithm, to the Kepler-419 system. We were able to compute all fourteen free parameters of the system, which, when integrated in time, give transits within the observational errors. We also studied the dynamics and the long-term orbital evolution of the Kepler-419 planetary system as defined by the orbital elements computed by us, in order to determine its stability.

We present the results of a study of different statistical methods currently used in the literature to analyse the (micro)variability of active galactic nuclei (AGNs) from ground-based optical observations. In particular, we focus on the comparison between the results obtained by applying the so-called C and F statistics, which are based on the ratio of standard deviations and variances, respectively. The motivation for this is that the implementation of these methods leads to different and contradictory results, making the variability classification of the light curves of a certain source dependent on the statistics implemented. For this purpose, we re-analyse the results on an AGN sample observed along several sessions with the 2.15m "Jorge Sahade" telescope (CASLEO), San Juan, Argentina. For each AGN we constructed the nightly differential light curves. We thus obtained a total of 78 light curves for 39 AGNs, and we then applied the statistical tests mentioned above, in order to re-classify the variability state of these light curves and in an attempt to find the suitable statistical methodology to study photometric (micro)variations. We conclude that, although the $C$ criterion is not proper a statistical test, it could still be a suitable parameter to detect variability and that its application allows us to get more reliable variability results, in contrast with the F test.

We built two self-consistent models of triaxial, cuspy, rotating stellar systems adding rotation to non-rotating models presented in previous papers of this series. The final angular velocity of the material is not constant and varies with the distance to the center and with the height over the equator of the systems, but the figure rotation is very uniform in both cases. Even though the addition of rotation to the models modifies their original semiaxes ratios, the final rotating models are considerably flattened and triaxial. An analysis of the orbital content of the models shows that about two thirds of their orbits are chaotic yet the models are very stable over intervals of the order of one Hubble time. The bulk of regular orbits are short axis tubes, while long axis tubes are replaced by tubes whose axes lie on the short-long axes plane, but do not coincide with the major axis. Other types of regular orbits that do not appear in non-rotating systems, like horseshoes and orbits that cross themselves, are also found in the present models. Finally, our frequency maps show empty regions where studies of orbits on fixed potentials found orbits, a likely consequence of the self-consistency of our models that excludes them.

Motivated by recent studies suggesting that the Large Magellanic Cloud (LMC) could be significantly more massive than previously thought, we explore whether the approximation of an inertial Galactocentric reference frame is still valid in the presence of such a massive LMC. We find that previous estimates of the LMC's orbital period and apocentric distance derived assuming a fixed Milky Way are significantly shortened for models where the Milky Way is allowed to move freely in response to the gravitational pull of the LMC. Holding other parameters fixed, the fraction of models favoring first infall is reduced. Due to this interaction, the Milky Way center of mass within the inner 50 kpc can be significantly displaced in phase-space in a very short period of time that ranges from 0.3 to 0.5 Gyr by as much as 30 kpc and 75 km/s. Furthermore, we show that the gravitational pull of the LMC and response of the Milky Way are likely to significantly affect the orbit and phase space distribution of tidal debris from the Sagittarius dwarf galaxy (Sgr). Such effects are larger than previous estimates based on the torque of the LMC alone. As a result, Sgr deposits debris in regions of the sky that are not aligned with the present-day Sgr orbital plane. In addition, we find that properly accounting for the movement of the Milky Way around its common center of mass with the LMC significantly modifies the angular distance between apocenters and tilts its orbital pole, alleviating tensions between previous models and observations. While these models are preliminary in nature, they highlight the central importance of accounting for the mutual gravitational interaction between the MW and LMC when modeling the kinematics of objects in the Milky Way and Local Group.

In several previous investigations we presented models of triaxial stellar systems, both cuspy and non cuspy, that were highly stable and harboured large fractions of chaotic orbits. All our models had been obtained through cold collapses of initially spherical $N$--body systems, a method that necessarily results in models with strongly radial velocity distributions. Here we investigate a different method that was reported to yield cuspy triaxial models with virtually no chaos. We show that such result was probably due to the use of an inadequate chaos detection technique and that, in fact, models with significant fractions of chaotic orbits result also from that method. Besides, starting with one of the models from the first paper in this series, we obtained three different models by rendering its velocity distribution much less radially biased (i.e., more isotropic) and by modifying its axial ratios through adiabatic compression. All three models yielded much higher fractions of regular orbits than most of those from our previous work. We conclude that it is possible to obtain stable cuspy triaxial models of stellar systems whose velocity distribution is more isotropic than that of the models obtained from cold collapses. Those models still harbour large fractions of chaotic orbits and, although it is difficult to compare the results from different models, we can tentatively conclude that chaoticity is reduced by velocity isotropy.

We model the dynamical interaction between the Small and Large Magellanic Clouds and their corresponding stellar cluster populations. Our goal is to explore whether the lack of old clusters ($\gtrsim 7$ Gyr) in the Small Magellanic Cloud (SMC) can be the result of the capture of clusters by the Large Magellanic Cloud (LMC), as well as their ejection due to the tidal interaction between the two galaxies. For this purpose we perform a suite of numerical simulations probing a wide range of parameters for the orbit of the SMC about the LMC. We find that, for orbital eccentricities $e \geq 0.4$, approximately 15 per cent of the SMC clusters are captured by the LMC. In addition, another 20 to 50 per cent of its clusters are ejected into the intergalactic medium. In general, the clusters lost by the SMC are the less tightly bound cluster population. The final LMC cluster distribution shows a spatial segregation between clusters that originally belonged to the LMC and those that were captured from the SMC. Clusters that originally belonged to the SMC are more likely to be found in the outskirts of the LMC. Within this scenario it is possible to interpret the difference observed between the star field and cluster SMC Age-Metallicity Relationships for ages $\gtrsim 7$ Gyr.

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.

Orbits in the principal planes of triaxial potentials are known to be prone to unstable motion normal to those planes, so that three dimensional investigations of those orbits are needed even though they are two dimensional. We present here an investigation of such orbits in the well known logarithmic potential which shows that the third dimension must be taken into account when studying them and that the instability worsens for lower values of the forces normal to the plane. Partially chaotic orbits are present around resonances, but also in other regions. The action normal to the plane seems to be related to the isolating integral that distinguishes regular form partially chaotic orbits, but not to the integral that distinguishes partially from fully chaotic orbits.

Numerical experiments conducted by Fellhauer et al. (MNRAS, 372, 338, 2006) suggest that a supercluster may capture up to about 40 per cent of its mass from the galaxy where it belongs. Nevertheless, in those experiments the cluster was created making appear its mass out of nothing, rather than from mass already present in the galaxy. Here we use a thought experiment, plus a few simple computations, to show that the difference between the dynamical effects of these two scenarios (i.e., mass creation vs. mass concentration) is actually very important. We also present the results of new numerical experiments, simulating the formation of the cluster through mass concentration, that show that trapping depends critically on the process of cluster formation and that the amounts of gained mass are substantially smaller than those obtained from mass creation.

The correlation dimension, that is, the dimension obtained by computing the correlation function of pairs of points of a trajectory in phase space, is a numerical technique introduced in the field of nonlinear dynamics in order to compute the dimension of the manifold in which an orbit moves, without the need of knowing the actual equations of motion that give rise to the trajectory. This technique has been proposed in the past as a method to measure the dimension of stellar orbits in astronomical potentials, i.e., the number of isolating integrals of motion the orbits obey. Although the algorithm can in principle yield that number, some care has to be taken in order to obtain good results. We studied the relevant parameters of the technique, found their optimal values, and tested the validity of the method on a number of potentials previously studied in the literature, using the SALI, Lyapunov exponents and spectral dynamics as gauges.

As a first step towards a comprehensive investigation of stellar motions within globular clusters, we present here the results of a study of stellar orbits in a mildly triaxial globular cluster that follows a circular orbit inside a galaxy. The stellar orbits were classified using the frequency analysis code of Carpintero and Aguilar and, as a check, the Liapunov characteristic exponents were also computed in some cases. The orbit families were obtained using different start spaces. Chaotic orbits turn out to be very common and while, as could be expected, they are particularly abundant in the outer parts of the cluster, they are still significant in the innermost regions. Their relevance for the structure of the cluster is discussed.