Short times characterisations of stochasticity in nonintegrable galactic potentials

This paper proposes a new, potentially useful, way in which to characterise the degree of stochasticity exhibited by orbits in a fixed time-independent galactic potential. This approach differs from earlier work involving the computation of Liapounov exponents in two ways, namely (1) by focusing on the statistical properties of ensembles of trajectories, rather than individual orbits, and (2) 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 will evolve relatively quickly towards a time-independent invariant measure {GAMMA}, which is arguably the natural unit to consider if one is interested in self-consistent equilibria. The basic idea proposed here is to compute short time Liapounov characteristic numbers χ({DELTA}t) over time intervals {DELTA}t for orbits in an ensemble that samples this invariant measure, and to analyse the overall distribution of these χ's. This is done in detail for one model potential, namely the sixth order truncation of the Toda lattice potential. One especially significant conclusion is that time averages and ensemble averages coincide, so that the form of the distribution of short time χ({DELTA}t)'s for such an ensemble is actually encoded in the calculation of χ(t) for a single orbit over long times t>>t_H_. The distribution of short time χ's is analysed as a function of the energy E of orbits in the ensemble and the length of the short time sampling interval {DELTA}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.