Regular and chaotic motion in elliptical galaxies

Here I review recent work, by other authors and by myself, on some particular topics related to the regular and chaotic motion in elliptical galaxies. I show that it is quite possible to build highly stable triaxial stellar systems that include large fractions of chaotic orbits and that partially and fully chaotic orbits fill different regions of space, so that it is important not to group them together under the single denomination of chaotic orbits. Partially chaotic orbits should not be confused with weakly fully chaotic orbits either, and their spatial distributions are also different. Slow figure rotation (i.e., rotation in systems with zero angular momentum) seems to be always present in highly flattened models that result from cold collapses, with the rotational velocity diminishing or becoming negligibly small for less flattened models. Finally, I comment on the usefulness and limitation of the classification of regular orbits via frequency analysis.Comment: 12 pages, 2 figures (both are mosaics of 6 and 4 individual figures, respectively, in eps format). It is an invited talk delivered at the workshop "Chaos in Astronomy 2007", in memory of N. Voglis, held in Athens (Greece), 17 - 20 September 2008, and accepted for publication in the Proccedings of that worksho

On the trapping of stars by a newborn stellar supercluster

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.Comment: 6 pages, 3 figures. Submitted to MNRA

Models of cuspy triaxial stellar systems. II. Regular orbits

In the first paper of this series we used the N--body method to build a dozen cuspy (gamma ~ 1) triaxial models of stellar systems, and we showed that they were highly stable over time intervals of the order of a Hubble time, even though they had very large fractions of chaotic orbits (more than 85 per cent in some cases). The models were grouped in four sets, each one comprising models morphologically resembling E2, E3, E4 and E5 galaxies, respectively. The three models within each set, although different, had the same global properties and were statistically equivalent. In the present paper we use frequency analysis to classify the regular orbits of those models. The bulk of those orbits are short axis tubes (SATs), with a significant fraction of long axis tubes (LATs) in the E2 models that decreases in the E3 and E4 models to become negligibly small in the E5 models. Most of the LATs in the E2 and E3 models are outer LATs, but the situation reverses in the E4 and E5 models where the few LATs are mainly inner LATs. As could be expected for cuspy models, most of the boxes are resonant orbits, i.e., boxlets. Nevertheless, only the (x, y) fishes of models E3 and E4 amount to about 10 per cent of the regular orbits, with most of the fractions of the other boxlets being of the order of 1 per cent or less.Comment: Accepted for publication in the Monthly Notices of the Royal Astronomical Societ

Predicting the progress of diffusively limited chemical reactions in the presence of chaotic advection

The effects of chaotic advection and diffusion on fast chemical reactions in two-dimensional fluid flows are investigated using experimentally measured stretching fields and fluorescent monitoring of the local concentration. Flow symmetry, Reynolds number, and mean path length affect the spatial distribution and time dependence of the reaction product. A single parameter \lambda*N, where \lambda is the mean Lyapunov exponent and N is the number of mixing cycles, can be used to predict the time-dependent total product for flows having different dynamical features.Comment: 4 pages, 4 figures, updated reference

Models of cuspy triaxial stellar systems. III: The effect of velocity anisotropy on chaoticity

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.Comment: 11 pages, 14 figures. Accepted for publication in MNRA

The structure of flame filaments in chaotic flows

The structure of flame filaments resulting from chaotic mixing within a combustion reaction is considered. The transverse profile of the filaments is investigated numerically and analytically based on a one-dimensional model that represents the effect of stirring as a convergent flow. The dependence of the steady solutions on the Damkohler number and Lewis number is treated in detail. It is found that, below a critical Damkohler number Da(crit), the flame is quenched by the flow. The quenching transition appears as a result of a saddle-node bifurcation where the stable steady filament solution collides with an unstable one. The shape of the steady solutions for the concentration and temperature profiles changes with the Lewis number and the value of Da(crit) increases monotonically with the Lewis number. Properties of the solutions are studied analytically in the limit of large Damkohler number and for small and large Lewis number.Comment: 17 pages, 13 figures, to be published in Physica

On the effect of chaotic orbits on dynamical friction

Chaotic orbits suffer significant changes as a result of small perturbations. One can thus wonder whether the dynamical friction suffered by a satellite on a regular orbit, and interacting with the stars of a galaxy, will be different if the bulk of the stars of the galaxy are in regular or chaotic orbits. In order to check that idea, we investigated the orbital decay (caused by dynamical friction) of a rigid satellite moving within a larger stellar system (a galaxy) whose potential is nonintegrable. We performed numerical experiments using two kinds of triaxial galaxy models: (1) the triaxial generalization of Dehnen´s spherical mass model (Dehnen; Merritt & Fridman); (2) a modified Satoh model (Satoh; Carpintero, Muzzio, & Wachlin). The percentages of chaotic orbits present in these models were increased by perturbing them. In the first case, a central compact object (black hole) was introduced; in the second case, the perturbation was produced by allowing the galaxy to move on a circular orbit in a logarithmic potential. The equations of motion were integrated with a non-self-consistent code. Our results show that the presence of chaotic orbits does not affect significantly the orbital decay of the satellite.Fil: Cora, Sofia Alejandra. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - La Plata. Instituto de Astrofísica La Plata. Universidad Nacional de La Plata. Facultad de Ciencias Astronómicas y Geofísicas. Instituto de Astrofísica la Plata; Argentina. Universidad Nacional de la Plata. Facultad de Ciencias Astronómicas y Geofísicas; ArgentinaFil: Vergne, Maria Marcela. Universidad Nacional de la Plata. Facultad de Ciencias Astronómicas y Geofísicas; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - La Plata. Instituto de Astrofísica La Plata. Universidad Nacional de La Plata. Facultad de Ciencias Astronómicas y Geofísicas. Instituto de Astrofísica la Plata; ArgentinaFil: Muzzio, Juan Carlos. Universidad Nacional de la Plata. Facultad de Ciencias Astronómicas y Geofísicas; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - La Plata. Instituto de Astrofísica La Plata. Universidad Nacional de La Plata. Facultad de Ciencias Astronómicas y Geofísicas. Instituto de Astrofísica la Plata; Argentin