Orbits in a three-dimensional potential subjected to periodic driving, V(x^i,t)=[1+m_0 sin(omega t) V_0(x^i), divide naturally into two types, regular and chaotic, between which transitions are seemingly impossible. The chaotic orbits divide in turn into two types, apparently separated by entropy barriers, namely `sticky' orbit segments, which are `locked' to the driving frequency and exhibit little systematic energy diffusion, and `wildly' chaotic segments, which are not so locked and can exhibit significant energy diffusion. Attention focuses on how the relative abundance of these different orbit types and the transition rate between sticky and wildly chaotic orbits depends on amplitude m_0, and on the extent to which these quantities can be altered by weak friction and/or noise and by pseudo-random variations in the driving frequency, idealized as an Ornstein-Uhlenbeck process with (in general) nonzero autocorrelation time t_c. When, in the absence of perturbations, there exist large measures of both regular and chaotic orbits, the primary effect of weak noise is to increase the relative measure of chaotic orbits. Alternatively, when almost all the orbits are already chaotic, noise serves primarily to accelerate transitions from sticky to wildly chaotic behavior. The presence or absence of friction is unimportant and the details of the noise seem largely immaterial. In particular, there is only a weak, roughly logarithmic dependence on amplitude, and additive and multiplicative noise typically have virtually identical effects. By contrast, allowing for high frequency `noisy' variations in omega tends to weaken effects of the driving, decreasing the relative measure of chaotic orbits and suppressing large scale energy diffusion.

This paper focuses on symmetric potentials subjected to periodic driving. Four unperturbed potentials V_0(r) were considered, namely the Plummer potential and Dehnen potentials with \gamma=0.0, 0.5, and 1.0, each subjected to a time-dependence of the form V(r,t)=V_0(r)(1+m_0\sin(\omega t)). In each case, the orbits divide clearly into regular and chaotic, distinctions which appear absolute. In particular, transitions from regularity to chaos are seemingly impossible. Over finite time intervals, chaotic orbits subdivide into what can be termed `sticky' chaotic orbits, which exhibit no large scale secular changes in energy and remain trapped in the phase space region where they started; and `wildly chaotic' orbits, which do exhibit systematic drifts in energy as the orbits diffuse to inhabit different phase space regions. This latter distinction is not absolute, apparently corresponding instead to orbits penetrating a `leaky' phase space barrier. The three different orbit types can be identified simply in terms of the frequencies for which their Fourier spectra have the most power. An examination of the statistical properties of orbit ensembles as a function of driving frequency \omega allows one to identify the specific resonances that determine orbital structure. Attention focuses also on how, for fixed amplitude m_0, such quantities as the mean energy shift, the relative measure of chaotic orbits, and the mean value of the largest Lyapunov exponent vary with driving frequency \omega and how, for fixed \omega, the same quantities depend on m_0.

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 investigates chaos and chaotic phase mixing in triaxial Dehnen potentials which have been proposed to describe realistic ellipticals. Earlier work is extended by exploring the effects of (1) variable axis ratios, (2) `graininess' associated with stars and bound substructures, idealised as friction and white noise, and (3) large-scale organised motions presumed to induce near-random forces idealised as coloured noise with finite autocorrelation time. Three important conclusions are: (1) not all the chaos can be attributed to the cusp; (2) significant chaos can persist even for axisymmetric systems; and (3) introducing a supermassive black hole can increase both the relative number of chaotic orbits and the size of the largest Lyapunov exponent. Sans perturbations, distribution functions associated with initially localised chaotic ensembles evolve exponentially towards a nearly time-independent form at a rate L that correlates with the finite time Lyapunov exponents associated with the evolving orbits. Perturbations accelerate phase space transport by increasing the rate of phase mixing in a given phase space region and by facilitating diffusion along the Arnold web that connects different phase space regions, thus facilitating an approach towards a true equilibrium. The details of the perturbation appear unimportant. All that matters are the amplitude and the autocorrelation time, upon which there is a weak logarithmic dependence. Even comparatively weak perturbations can increase L by a factor of three or more, a fact that has potentially significant implications for violent relaxation.

This note illustrates the possibility in simple loaded string models of trapping most of the system energy in a single degree of freedom for very long times, demonstrating in particular that the robustness of the trapping is enhanced by increasing the `connectance' of the system, i.e., the extent to which many degrees of freedom are coupled directly by the interaction Hamiltonian, and/or the strength of the couplings.

Recent work has shown that statistical arguments, seemingly well-justified in higher dimensions, can also be used to derive reasonable, albeit less accurate, estimates of the largest Lyapunov exponent ${\chi}$ in lower-dimensional Hamiltonian systems. This letter explores the detailed assumptions incorporated into these arguments. The predicted values of ${\chi}$ are insensitive to most of these details, which can in any event be relaxed straightforwardly, but {\em can} depend sensitively on the nongeneric form of the auto-correlation function characterising the time-dependence of an orbit. This dependence on dynamics implies a fundamental limitation to the application of thermodynamic arguments to such lower-dimensional systems.

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.

This paper discusses three new issues that necessarily arise in realistic attempts to apply nonlinear dynamics to galaxy evolution, namely: (i) the meaning of chaos in many-body systems, (ii) the time-dependence of the bulk potential, which can trigger intervals of {\em transient chaos}, and (iii) the self-consistent nature of any bulk chaos, which is generated by the bodies themselves, rather than imposed externally. Simulations and theory both suggest strongly that the physical processes associated with galactic evolution should also act in nonneutral plasmas and charged particle beams. This in turn suggests the possibility of testing this physics in real laboratory experiments, an undertaking currently underway.

This talk summarises several different lines of argument suggesting that one should not expect cuspy nonaxisymmetric galaxies to exist as robust, long-lived collisionless equilibria, i.e., that such objects should not be idealised as time-independent solutions to the collisionless Boltzmann equation.

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 a numerical investigation of phase mixing in time-independent Hamiltonian systems that admit a coexistence of regular and chaotic phase space regions, allowing also for low amplitude perturbations idealised as periodic driving, friction, and/or white and colored noise. The evolution of initially localised ensembles of orbits was probed through lower order moments and coarse-grained distribution functions. In the absence of time-dependent perturbations, regular ensembles disperse initially as a power law in time and only exhibit a coarse-grained approach towards an invariant equilibrium over comparatively long times. Chaotic ensembles generally diverge exponentially fast on a time scale related to a typical finite time Lyapunov exponent, but can exhibit complex behaviour if they are impacted by the effects of cantori or the Arnold web. Viewed over somewhat longer times, chaotic ensembles typical converge exponentially towards an invariant or near-invariant equilibrium. This, however, need not correspond to a true equilibrium, which may only be approached over very long time scales. Time-dependent perturbations can dramatically increase the efficiency of phase mixing, both by accelerating the approach towards a near-equilibrium and by facilitating diffusion through cantori or along the Arnold web so as to accelerate the approach towards a true equilibrium. The efficacy of such perturbations typically scales logarithmically in amplitude, but is comparatively insensitive to most other details, a conclusion which reinforces the interpretation that the perturbations act via a resonant coupling.

Cuspy triaxial potentials admit a large number of chaotic orbits, which moreover exhibit extreme "stickiness" that makes the process of chaotic mixing surprisingly inefficient. Environmental effects, modeled as noise and/or periodic driving, help accelerate phase space transport but probably not as much as in simpler potentials. This could mean that cuspy triaxial ellipticals cannot exist as time-independent systems.

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 paper explores the phenomenon of energy relaxation for stars in a galaxy embedded in a high density environment that is subjected continually to perturbations reflecting the presence of other nearby galaxies and/or incoherent internal pulsations. The analysis is similar to earlier analyses of energy relaxation induced by binary encounters between nearby stars and between stars and giant molecular clouds in that the perturbations are idealised as a sum of near-random events which can be modeled as diffusion and dynamical friction. However, the analysis differs in one important respect: because the time scale associated with these perturbations need not be short compared with the characteristic dynamical time t_D for stars in the original galaxy, the diffusion process cannot be modeled as resulting from a sequence of instantaneous kicks, i.e., white noise. Instead, the diffusion is modeled as resulting from random kicks of finite duration, i.e., coloured noise characterised by a nonzero autocorrelation time t_c. A detailed analysis of coloured noise generated by sampling an Ornstein-Uhlenbeck process leads to a simpling scaling in terms of t_c and an effective diffusion constant D. Interpreting D and t_c following early work by Chandrasekhar (1941) (the `nearest neighbour approximation') implies that, for realistic choices of parameter values, energy relaxation associated with an external environment and/or internal pulsations could be important on times short compared with the age of the Universe.

This talk summarises what is currently understood about the phenomenon that has come to be known as {\it chaotic mixing}. The first part presents a concise statement as to what chaotic mixing actually is, and then explains why it should be important in galactic dynamics. The second discusses in detail what is currently known, describing the manifestations of chaotic mixing both for flows in time-independent Hamiltonian systems and for systems subjected to comparatively weak time-dependent perturbations. The third part discusses what one might expect if one were to allow for a strongly time-dependent potential, including possible implications for violent relaxation.

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.

(Abridged) This paper studies chaotic orbit ensembles evolved in triaxial generalisations of the Dehnen potential which have been proposed to model ellipticals with a strong density cusp that manifest significant deviations from axisymmetry. Allowance is made for a possible supermassive black hole, as well as low amplitude friction, noise, and periodic driving which can mimic irregularities associated with discreteness effects and/or an external environment. The degree of chaos is quantified by determining how (1) the relative number of chaotic orbits and (2) the size of the largest Lyapunov exponent depend on the steepness of the cusp and the black hole mass, and (3) the extent to which Arnold webs significantly impede phase space transport, both with and without perturbations. In the absence of irregularities, chaotic orbits tend to be extremely `sticky,' so that different pieces of the same chaotic orbit can behave very differently for 10000 dynamical times or longer, but even very low amplitude perturbations can prove efficient in erasing many -- albeit not all -- these differences. The implications thereof are discussed both for the structure and evolution of real galaxies and for the possibility of constructing approximate near-equilibrium models using Schwarzschild's method. Much of the observed qualitative behaviour can be reproduced with a toy potential given as the sum of an anisotropic harmonic oscillator and a spherical Plummer potential, which suggests that the results may be generic.

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.

This paper summarises an investigation of the statistical properties of orbits escaping from three different two-degree-of-freedom Hamiltonian systems which exhibit global stochasticity. Each H=H_{0}+eH', with H_{0} integrable and eH' a nonintegrable correction, not necessarily small. For e below a critical e_{0} escapes are impossible energetically. For somewhat higher values, escape is allowed energetically but many orbits never escape and the escape probability P for a generic orbit ensemble decays exponentially. At or near a critical e_{1}>e_{0}, there is an abrupt qualitative change in behaviour. Above e_{1}, P typically exhibits (1) a rapid evolution towards a nonzero P_{0}(e) followed by (2) a much slower subsequent decay towards zero which, in at least one case, is well fit by a power law P=const x t^{m}, with m=0.35-0.40. In all three cases, P_{0} and the time T required to converge towards P_{0} scales in e-e_{1}, i.e., P_{0}=const x (e-e_{1})^{a} and T=const x (e-e_{1})^{b}, and T also scales in the size r of the region sampled for initial conditions, i.e., T=const x r^{-d}. To within statistical uncertainties, the best fit values of the critical exponents are the same for all three potentials, namely: a=0.5, b=0.4, and d=0.1, and satisfy a-b-d=0. The transitional behaviour observed near e_{1} is attributed to the breakdown of some especially significant KAM tori or cantori. The power law behaviour at late times is interpreted as reflecting intrinsic diffusion of chaotic orbits through cantori surrounding islands of regular orbits.

This paper summarises a numerical investigation of how the usual manifestations of chaos and regularity for flows in time-independent Hamiltonians can be alterred by a systematic time-dependence of the form arising naturally in an expanding Universe. If the time-dependence is not too strong, the observed behaviour can be understood in an adiabatic approximation. One still infers sharp distinctions between regular and chaotic behaviour, even though ``regular'' does not mean ``periodic'' and ``chaotic'' will not in general imply strictly exponential sensitivity towards small changes in initial conditions. However, these distinctions are no longer absolute, as it is possible for a single orbit to change from regular to chaotic and/or visa versa. If the time-dependence becomes too strong, the distinction between regular and chaotic can disappear so that no sensitive dependence on initial conditions is manifest.

First passage time experiments were used to explore the effects of low amplitude noise as a source of accelerated phase space diffusion in two-dimensional Hamiltonian systems, and these effects were then compared with the effects of periodic driving. The objective was to quantify and understand the manner in which ``sticky'' chaotic orbits that, in the absence of perturbations, are confined near regular islands for very long times, can become ``unstuck'' much more quickly when subjected to even very weak perturbations. For both noise and periodic driving, the typical escape time scales logarithmically with the amplitude of the perturbation. For white noise, the details seem unimportant: Additive and multiplicative noise typically have very similar effects, and the presence or absence of a friction related to the noise by a Fluctuation-Dissipation Theorem is also largely irrelevant. Allowing for colored noise can significantly decrease the efficacy of the perturbation, but only when the autocorrelation time becomes so large that there is little power at frequencies comparable to the natural frequencies of the unperturbed orbit. Similarly, periodic driving is relatively inefficient when the driving frequency is not comparable to these natural frequencies. This suggests strongly that noise-induced extrinsic diffusion, like modulational diffusion associated with periodic driving, is a resonance phenomenon. The logarithmic dependence of the escape time on amplitude reflects the fact that the time required for perturbed and unperturbed orbits to diverge a given distance scales logarithmically in the amplitude of the perturbation.

This paper summarises a numerical investigation which aimed to identify and characterise regular and chaotic behaviour in time-dependent Hamiltonians H(r,p,t) = p^2/2 + U(r,t), with U=R(t)V(r) or U=V[R(t)r], where V(r) is a polynomial in x, y, and/or z, and R = const * t^p is a time-dependent scale factor. When p is not too negative, one can distinguish between regular and chaotic behaviour by determining whether an orbit segment exhibits a sensitive dependence on initial conditions. However, chaotic segments in these potentials differ from chaotic segments in time-independent potentials in that a small initial perturbation will usually exhibit a sub- or super-exponential growth in time. Although not periodic, regular segments typically exhibit simpler shapes, topologies, and Fourier spectra than do chaotic segments. This distinction between regular and chaotic behaviour is not absolute since a single orbit segment can seemingly change from regular to chaotic and visa versa. All these observed phenomena can be understood in terms of a simple theoretical model.

The objective of the work summarised here has been to exploit and extend ideas from plasma physics and accelerator dynamics to formulate a unified description of collisionless relaxation that views violent relaxation, Landau damping, and phase mixing as (manifestations of) a single phenomenon. This approach embraces the fact that the collisionless Boltzmann equation (CBE), the basic object of the theory, is an infinite-dimensional Hamiltonian system, with the distribution function f playing the role of the fundamental dynamical variable, and that, interpreted appropriately, an evolution described by the {\it CBE} is no different fundamentally from an evolution described by any other Hamiltonian system. Equilibrium solutions correspond to extremal points of the Hamiltonian subject to the constraints associated with Liouville's Theorem. Stable equilibria correspond to energy minima. The evolution of a system out of equilibrium involves (in general nonlinear) phase space oscillations which may -- or may not -- interfere destructively so as to damp away.

Everything you ever wanted to know about what has come to be known as ``chaotic mixing:'' This paper describes the evolution of localised ensembles of initial conditions in 2- and 3-D time-independent potentials which admit both regular and chaotic orbits. The coarse-grained approach towards an invariant, or near-invariant, distribution was probed by tracking (1) phase space moments through order 4 and (2) binned reduced distributions f(Z_a,Z_b,t) for a,b=x,y,z,p_x,p_y,p_z, computed at fixed time intervals. For ``unconfined'' chaotic orbits in 2-D systems not stuck near islands by cantori, the moments evolve exponentially: Quantities like the dispersion in p_x, which start small and eventually asymptote towards a larger value, initially grow exponentially in time at a rate comparable to the largest short time Lyapunov exponent. Quantities like |<p_x>|, that can start large but eventually asymptote towards zero, decrease exponentially. With respect to a discrete L^p norm, reduced distributions f(t) generated from successive decay exponentially towards a near-invariant f_{niv}, although a plot of Df(t)=||f(t)-f_{niv}|| can exhibit considerable structure. Regular ensembles behave very differently, both moments and Df evolving in a fashion better represented by a power law time dependence. ``Confined'' chaotic orbits, initially stuck near regular islands because of cantori, exhibit an intermediate behaviour. The behaviour of ensembles evolved in 3-D potentials is qualitatively similar, except that, in this case, it is relatively likely to find one direction in configuration space which is ``less chaotic'' than the other two, so that quantities like L_{ab} depend more sensitively on which phase space variables one tracks.

This paper analyses the effect of low amplitude friction and noise in accelerating phase space transport in time-independent Hamiltonian systems that exhibit global stochasticity. Numerical experiments reveal that even very weak non-Hamiltonian perturbations can dramatically increase the rate at which an ensemble of orbits penetrates obstructions like cantori or Arnold webs, thus accelerating the approach towards an invariant measure, i.e., a near-microcanonical population of the accessible phase space region. An investigation of first passage times through cantori leads to three conclusions, namely: (i) that, at least for white noise, the detailed form of the perturbation is unimportant, (ii) that the presence or absence of friction is largely irrelevant, and (iii) that, overall, the amplitude of the response to weak noise scales logarithmically in the amplitude of the noise.

This paper discusses the possibility of constructing time-independent solutions to the collisionless Boltzmann equation which depend on quantities other than global isolating integrals such as energy and angular momentum. The key point is that, at least in principle, a self-consistent equilibrium can be constructed from any set of time-independent phase space building blocks which, when combined, generate the mass distribution associated with an assumed time-independent potential. This approach provides a way to justify Schwarzschild's (1979) method for the numerical construction of self-consistent equilibria with arbitrary time-independent potentials, generalising thereby an approach developed by Vandervoort (1984) for integrable potentials. As a simple illustration, Schwarzschild's method is reformulated to allow for a straightforward computation of equilibria which depend only on one or two global integrals and no other quantities, as is reasonable, e.g., for modeling axisymmetric configurations characterised by a nonintegrable potential.

This paper proposes a geometric interpretation of flows generated by the collisionless Boltzmann equation (CBE), focusing on the coarse-grained approach towards equilibrium. The CBE is a noncanonical Hamiltonian system with the distribution function f the fundamental dynamical variable, the mean field energy H[f] playing the role of the Hamiltonian and the natural arena of physics being the infinite-dimensional phase space of distribution functions. Every time-independent equilibrium f_0 is an energy extremal with respect to all perturbations that preserve the constraints associated with Liouville's Theorem, local energy minima corresponding to linearly stable equilibria. If an initial f(t=0) is sufficiently close to some linearly stable lower energy f_0, its evolution involves linear phase space oscillations about f_0 which, in many cases, would be expected to exhibit linear Landau damping. If f(t=0) is far from any stable extremal, the flow will be more complicated but, in general, one would anticipate that the evolution involves nonlinear oscillations about some lower energy f_0. In this picture, the coarse-grained approach towards equilibrium usually termed violent relaxation is interpreted as nonlinear Landau damping. The evolution of a generic initial f(t=0) involves a coherent initial excitation, not necessarily small, being converted into incoherent motion associated with nonlinear oscillations about some equilibrium f_0 which, in general, will exhibit destructive interference.

Time-independent Hamiltonian flows are viewed as geodesic flows in a curved manifold, so that the onset of chaos hinges on properties of the curvature two-form entering into the Jacobi equation. Attention focuses on ensembles of orbit segments evolved in 2-D potentials, examining how various orbital properties correlate with the mean value and dispersion, <K> and k, of the trace K of the curvature. Unlike most analyses, which have attributed chaos to negative curvature, this work exploits the fact that geodesics can be chaotic even if K is everywhere positive, chaos arising as a parameteric instability triggered by regular variations in K along the orbit. For ensembles of fixed energy, with both regular and chaotic segments, simple patterns connect the values of <K> and k for different segments, both with each other and with the short time Lyapunov exponent X. Often, but not always, there is a near one-to- one correlation between <K> and k, a plot of these quantities approximating a simple curve. X varies smoothly along this curve, chaotic segments located furthest from the regular regions tending systematically to have the largest X's. For regular orbits, <K> and k also vary smoothly with ``distance'' from the chaotic phase space regions, as probed, e.g., by the location of the initial condition on a surface of section. Many of these observed properties can be understood qualitatively in terms of a one-dimensional Mathieu equation.

This talk provides a critical assessment of collisionless galactic dynamics, focusing on the interpretation and limitations of the collisionless Boltzmann equation and the physical mechanisms associated with collisionless relaxation. Numerical and theoretical arguments are presented to motivate the idea that the evolution of a system far from equilibrium should be interpreted as involving nonlinear gravitational Landau damping, which implies a greater overall coherence and remembrance of initial conditions than is implicit in the conventional theory of violent relaxation.

This paper summarises a number of new, potentially significant, results, obtained recently by the author and his collaborators, which impact on various issues related to the gravitational N-body problem, both Newtonianly and in the context of general relativity. Topics addressed include: (1) direct N-body simulations and their interpretation, with reference to the observed exponential instability towards small changes in initial conditions and the phenomenon of ``nonviolent relaxation;'' (2) the Hamiltonian structure of the collisionless Boltzmann equation of general relativity, i.e., the Vlasov-Einstein system; (3) ``transient ensemble dynamics,'' i.e.,,the short time statistical characterisation of collections of orbits in nonintegrable mean field potentials; and (4) the structural stability of the smooth potential approximation typically used in galactic dynamics.

This paper proposes an alternative characterisation of the degree of stochasticity exhibited by orbits in a fixed galactic potential. This differs from earlier work involving Liapounov exponents by focusing on the statistical properties of ensembles of trajectories, rather than individual orbits, and by restricting attention to the properties of these ensembles over time scales shorter than the age of the Universe, t_H. For many potentials, generic ensembles of initial conditions corresponding to stochastic orbits exhibit a rapid evolution towards an invariant measure $\Gamma$, the natural unit to consider if one is interested in self-consistent equilibria. The basic idea is to compute short time Liapounov exponents ${\chi}(Dt)$ over time intervals $Dt$ for orbits in an ensemble that samples $\Gamma$, and to analyse the overall distribution of these ${\chi}$'s. One finds that time averages and ensemble averages coincide, so that the form of the distribution of short time ${\chi}(Dt)$'s for such an ensemble is actually encoded in the calculation of ${\chi}(t)$ for a single orbit over long times $t>>t_{H}$. The distribution of short time ${\chi}$'s is analysed as a function of the energy E of orbits in the ensemble and the length of the short time sampling interval ${\D}t$. For relatively high energies, the distribution is essentially Gaussian, the dispersion decreasing with time as $t^{-p}$, with an exponent 0<p<1/2 that depends on the energy E. (to appear in Astronomy and Astrophysics)

ABBREVIATED ABSTRACT: This paper summarises an investigation of the effects of weak friction and noise in time-independent, nonintegrable potentials which admit both regular and stochastic orbits. The aim is to understand the qualitative effects of internal and external irregularities associated, e.g., with discreteness effects or couplings to an external environment, which stars in any real galaxy must experience. The two principal conclusions are: (1) These irregularities can be important on time scales much shorter than the natural relaxation time scale t_R associated with the friction and noise. For stochastic orbits friction and noise induce an average exponential divergence from the unperturbed Hamiltonian trajectory at a rate set by the value of the local Lyapunov exponent. Even weak noise can make a pointwise interpretation of orbits suspect already on time scales much shorter than t_R. (2) The friction and noise can also have significant effects on the statistical properties of ensembles of stochastic orbits, these also occurring on time scales much shorter than t_R. Potential implications for galactic dynamics are discussed, including the problem of shadowing.

Results are reported from a numerical investigation of orbits in a truncated Toda potential which is perturbed by weak friction and noise. Two significant conclusions are shown to emerge: (1) Despite other nontrivial behaviour, configuration, velocity, and energy space moments associated with these perturbations exhibit a simple scaling in the amplitude of the friction and noise. (2) Even very weak friction and noise can induce an extrinsic diffusion through cantori on a time scale much shorter than that associated with intrinsic diffusion in the unperturbed system.