The generic gravitational three-body problem is non-integrable and the motion is susceptible to disintegration. For negative total energy the decay outcome is a free body flying apart from a binary. Computer simulations have shown that initial conditions within the chaotic region follow a rather random motion. Decades ago a statistical solution was marked as a goal, yet despite considerable progress, all extant approaches display two flaws. First, probability was equated with phase space volume, ignoring the fact that a significant part of phase space describes regular motion, including post-decay motion. Secondly and relatedly, an adjustable parameter, the strong interaction region, which is a sort of cutoff, was a central ingredient of the theory. Based on an analogy with a particle moving in a leaky container, this paper presents an exact reduction of the statistical prediction, or more precisely, the differential decay rate is shown to factorize. One factor is the outgoing phase-space volume flux, rather than the volume itself, and it is given in a cutoff-independent closed-form. The other factors are the chaotic absorptivity and the regularized phase space volume, which serves as a normalization. The situation is analogous to Kirchhoff's law of thermal radiation, also known as greybody radiation. The statistical prediction is completed by including for the first time not only the outcome probabilities but also the decay rate and its statistics through the introduction of an equation system for the time evolution of the statistical distribution, which accounts for sub-escape excursions and can be described as a differential equation system with memory.

This paper describes and proves a canonical procedure to decouple perturbations and optimize their gauge around backgrounds with one non-homogeneous dimension, namely of co-homogeneity 1, while preserving locality in this dimension. Derivatively-gauged fields are shown to have a purely algebraic action; they can be decoupled from the other fields through gauge-invariant re-definitions; a potential for the other fields is generated in the process; in the remaining action each gauge function eliminates one field without residual gauge. The procedure applies to spherically symmetric and to cosmological backgrounds in either General Relativity or gauge theories. The widely used ``gauge invariant perturbation theory'' is closely related. The supplied general proof elucidates the algebraic mechanism behind it as well as the method's domain of validity and its assumptions.