{} {To quantify the amount of chaos that exists in the local phase space.} {A sample of orbits from four different models of the Solar neighbourhood phase space are analysed by a new chaos identification (and quantification) technique. While three of the used models bear the signature of the perturbation due to both the Galactic bar and the spiral pattern, the last of the models is a bar only one. We explore the models by inter-comparing the corresponding values of chaos strength that is induced at the various energy levels .}{(1) We find that of all the viable models that have been demonstrated to successfully reproduce the local phase space structure, i.e. those that include the bar as well as the spiral, bear strong chaoticity, though the model that implies the highest degree of chaos is the one in which overlap of the major resonances of the bar and the spiral occurs. The bar only model is found to display regularity. (2) We advance chaos to be primarily responsible for the splitting of the Hyades-Pleiades mode (the larger mode) of the local velocity distribution}{}

It has been shown that the dissipative gas infall during galaxy formation has the capability to modify the shape of dark matter halos. In this paper we perform the first detailed analysis of particle orbits in a cosmological dark matter halos to understand {\it how} and {\it why} baryons alter its shape. We perform a series of numerical experiments where we grow a baryonic core inside a live dark matter halo extracted from a cosmological simulation. We follow the orbits of more than 300 particles with more that 50000 timesteps. Our results clearly show that the dissipational component is responsible for repeatedly deflecting orbits which visit often the center of the system. Moreover the gravitational potential time dependence associated with the growth of the baryonic mass, shifts the frequencies of the orbits, making them extremely chaotic. This randomization makes the orbits explore a large phase space. When this effect takes place for a significant number of orbits it will be manifested in the density distribution as an approach to a rounder configuration. As a consequence, the influence of the central mass on the shape of the phase space decreases with increasing distance from the center. We discuss the importance of future analysis of controlled experiments (i.e. using analytic potentials instead of live DM halos) to better decipher the dynamics of this phenomenon.

This paper describes how parametric resonances associated with a galactic potential subjected to relatively low amplitude, strictly periodic time-dependent perturbations can be impacted by pseudo-random variations in the pulsation frequency, modeled as coloured noise. One aim thereby is to allow for the effects of a changing oscillation frequency as the density distribution associated with a galaxy evolves during violent relaxation. Another is to mimic the possible effects of internal substructures, satellite galaxies, and or a high density environment. The principal conclusion is that allowing for a variable frequency does not vitiate the effects of parametric resonance; and that, in at least some cases, such variations can increase the overall importance of parametric resonance associated with systematic pulsations.

This paper focuses on the dynamical implications of close supermassive black hole binaries both as an example of resonant phase mixing and as a potential explanation of inversions and other anomalous features observed in the luminosity profiles of some elliptical galaxies. The presence of a binary comprised of black holes executing nearly periodic orbits leads to the possibility of a broad resonant coupling between the black holes and various stars in the galaxy. This can result in efficient chaotic phase mixing and, in many cases, systematic increases in the energies of stars and their consequent transport towards larger radii. Allowing for the presence of a supermassive black hole binary with plausible parameter values near the center of a spherical, or nearly spherical, galaxy characterised initially by a Nuker density profile enables one to reproduce in considerable detail the central surface brightness distributions of such galaxies as NGC 3706.

This paper examines discreteness effects in nearly collisionless N-body systems of charged particles interacting via an unscreened r^-2 force, allowing for bulk potentials admitting both regular and chaotic orbits. Both for ensembles and individual orbits, as N increases there is a smooth convergence towards a continuum limit. Discreteness effects are well modeled by Gaussian white noise with relaxation time t_R = const * (N/log L)t_D, with L the Coulomb logarithm and t_D the dynamical time scale. Discreteness effects accelerate emittance growth for initially localised clumps. However, even allowing for discreteness effects one can distinguish between orbits which, in the continuum limit, feel a regular potential, so that emittance grows as a power law in time, and chaotic orbits, where emittance grows exponentially. For sufficiently large N, one can distinguish two different `kinds' of chaos. Short range microchaos, associated with close encounters between charges, is a generic feature, yielding large positive Lyapunov exponents X_N which do not decrease with increasing N even if the bulk potential is integrable. Alternatively, there is the possibility of larger scale macrochaos, characterised by smaller Lyapunov exponents X_S, which is present only if the bulk potential is chaotic. Conventional computations of Lyapunov exponents probe X_N, leading to the oxymoronic conclusion that N-body orbits which look nearly regular and have sharply peaked Fourier spectra are `very chaotic.' However, the `range' of the microchaos, set by the typical interparticle spacing, decreases as N increases, so that, for large N, this microchaos, albeit very strong, is largely irrelevant macroscopically. A more careful numerical analysis allows one to estimate both X_N and X_S.

This paper explores how orbits in a galactic potential can be impacted by large amplitude time-dependences of the form that one might associate with galaxy or halo formation or strong encounters between pairs of galaxies. A period of time-dependence with a strong, possibly damped, oscillatory component can give rise to large amounts of transient chaos, and it is argued that chaotic phase mixing associated with this transient chaos could play a major role in accounting for the speed and efficiency of violent relaxation. Analysis of simple toy models involving time-dependent perturbations of an integrable Plummer potential indicates that this chaos results from a broad, possibly generic, resonance between the frequencies of the orbits and harmonics thereof and the frequencies of the time-dependent perturbation. Numerical computations of orbits in potentials exhibiting damped oscillations suggest that, within a period of 10 dynamical times t_D or so, one could achieve simultaneously both `near-complete' chaotic phase mixing and a nearly time-independent, integrable end state.

Integrations in fixed N-body realisations of smooth density distributions corresponding to a chaotic galactic potential can be used to derive reliable estimates of the largest (finite time) Lyapunov exponent X_S associated with an orbit in the smooth potential generated from the same initial condition, even though the N-body orbit is typically characterised by an N-body exponent X_N >> X_S. This can be accomplished either by comparing initially nearby orbits in a single N-body system or by tracking orbits with the same initial condition evolved in two different N-body realisations of the same smooth density.

This paper summarises an investigation of chaos in a toy potential which mimics much of the behaviour observed for the more realistic triaxial generalisations of the Dehnen potentials, which have been used to model cuspy triaxial galaxies both with and without a supermassive black hole. The potential is the sum of an anisotropic harmonic oscillator potential, V_o=(1/2)*(a^{2}x^{2}+b^{2}y^{2}+c^{2}z^{2}), and a spherical Plummer potential, V_o=-M_{BH}/(r^{2}+e^{2})^{1/2} with r^{2}=x^{2}+y^{2}+z^{2}. Attention focuses on three issues related to the properties of ensembles of chaotic orbits which impact on chaotic mixing and the possibility of constructing self-consistent equilibria: (1) What fraction of the orbits are chaotic? (2) How sensitive are the chaotic orbits, i.e., how large are their largest (short time) Lyapunov exponents? (3) To what extent is the motion of chaotic orbits impeded by Arnold webs, i.e.,, how `sticky' are the chaotic orbits? These questions are explored as functions of the axis ratio a:b:c, black hole mass M_BH, softening length e, and energy E with the aims of understanding how the manifestations of chaos depend on the shape of the system and why the black hole generates chaos. The simplicity of the model makes it amenable to a perturbative analysis. That it mimics the behaviour of more complicated potentials suggests that much of this behaviour should be generic.

This paper continues a numerical investigation of orbits evolved in `frozen,' time-independent N-body realisations of smooth time-independent density distributions corresponding to both integrable and nonintegrable potentials, allowing for N as large as 300,000. The principal focus is on distinguishing between, and quantifying, the effects of graininess on initial conditions corresponding, in the continuum limit, to regular and chaotic orbits. Ordinary Lyapunov exponents X do not provide a useful diagnostic for distinguishing between regular and chaotic behaviour. Frozen-N orbits corresponding in the continuum limit to both regular and chaotic characteristics have large positive X even though, for large N, the `regular' frozen-N orbits closely resemble regular characteristics in the smooth potential. Viewed macroscopically both `regular' and `chaotic' frozen-N orbits diverge as a power law in time from smooth orbits with the same initial condition. There is, however, an important difference between `regular' and `chaotic' frozen-N orbits: For regular orbits, the time scale associated with this divergence t_G ~ N^{1/2}t_D, with t_D a characteristic dynamical time; for chaotic orbits t_G ~ (ln N) t_D. At least for N>1000 or so, clear distinctions exist between phase mixing of initially localised orbit ensembles which, in the continuum limit, exhibit regular versus chaotic behaviour. For both regular and chaotic ensembles, finite-N effects are well mimicked, both qualitatively and quantitatively, by energy-conserving white noise with amplitude ~ 1/N. This suggests strongly that earlier investigations of the effects of low amplitude noise on phase space transport in smooth potentials are directly relevant to real physical systems.

This paper uses the assumptions of ergodicity and a microcanonical distribution to compute estimates of the largest Lyapunov exponents in lower-dimensional Hamiltonian systems. That the resulting estimates are in reasonable agreement with the actual values computed numerically corroborates the intuition that chaos in such systems can be understood as arising generically from a parametric instability and that this instability can be modeled by a stochastic-oscillator equation (cf. Casetti, Clementi, and Pettini, Phys. Rev. E 54, 5969 (1996)), linearised perturbations of a chaotic orbit satisfying a harmonic-oscillator equation with a randomly varying frequency.

This paper summarises a numerical investigation of the statistical properties of orbits evolved in `frozen,' time-independent N-body realisations of smooth, time-independent density distributions, allowing for 10^(2.5)<N<10^(5.5). Two principal conclusions were reached: (1) In the limit of a nearly `unsoftened' two-body kernel, the value of the largest Lyapunov exponent does NOT appear to decrease systematically with increasing N. (2) Nevertheless, there is a clear quantifiable sense in which, on the average, as N increases chaotic orbits in the frozen-N systems come to more closely approximate characteristics in the smooth potential. When viewed in configuration or velocity space, or as probed by collisionless invariants like angular momentum, frozen-N orbits typically diverge from smooth potential characteristics as a power law in time on a time scale proportional to N^(1/2)t_D, with t_D a characteristic dynamical, or crossing, time.

This talk provides a progress report on an extended collaboration which has aimed to address two basic questions, namely: Should one expect to see cuspy, triaxial galaxies in nature? And can one construct realistic cuspy, triaxial equilibrium models that are robust? Three technical results are described: (1) Unperturbed chaotic orbits in cuspy triaxial potentials can be extraordinarily sticky, much more so than orbits in many other three-dimensional potentials. (2) Even very weak perturbations can be important by drastically reducing, albeit not completely eliminating, this stickiness. (3) A simple toy model facilitates a simple understanding of why black holes and cusps can serve as an effective source of chaos. These results suggest that, when constructing models of galaxies using Schwarzschild's method or any analogue thereof, astronomers would be well advised to use orbital building blocks that have been perturbed by `noise' or other weak irregularities, since such building blocks are likely to be more nearly time-independent than orbits evolved in the absence of all perturbations.

This paper summarises an investigation of the effects of low amplitude noise and periodic driving on phase space transport in 3-D Hamiltonian systems, a problem directly applicable to systems like galaxies, where such perturbations reflect internal irregularities and.or a surrounding environment. A new diagnsotic tool is exploited to quantify how, over long times, different segments of the same chaotic orbit can exhibit very different amounts of chaos. First passage time experiments are used to study how small perturbations of an individual orbit can dramatically accelerate phase space transport, allowing `sticky' chaotic orbits trapped near regular islands to become unstuck on suprisingly short time scales. Small perturbations are also studied in the context of orbit ensembles with the aim of understanding how such irregularities can increase the efficacy of chaotic mixing. For both noise and periodic driving, the effect of the perturbation scales roughly in amplitude. For white noise, the details are unimportant: additive and multiplicative noise tend to have similar effects and the presence or absence of a friction related to the noise by a Fluctuation- Dissipation Theorem is largely irrelevant. Allowing for coloured noise can significantly decrease the efficacy of the perturbation, but only when the autocorrelation time, which vanishes for white noise, becomes so large that t here is little power at frequencies comparable to the natural frequencies of the unperturbed orbit. This suggests strongly that noise-induced extrinsic diffusion, like modulational diffusion associated with periodic driving, is a resonance phenomenon. Potential implications for galaxies are discussed.