Diffusion and scaling in escapes from two-degree-of-freedom Hamiltonian systems

Abstract

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.