---
title: 'Kepler Billiards: Dynamics & Integrability'
url: https://www.emergentmind.com/topics/kepler-billiards
type: topic
---

# Kepler Billiards: Dynamics & Integrability

Searching arXiv for recent papers on Kepler billiards and related integrability results.
Kepler billiards are mechanical billiard systems in which a point particle moves between impacts under a Kepler potential and undergoes elastic reflection at a prescribed boundary. In the planar one-center setting, the underlying Hamiltonian is \(H(q,p)=\tfrac12|p|^2-\mu/|q|\) with \(\mu>0\), and the reflection law preserves the tangential component of momentum while reversing the normal component [2507.15879]. The subject combines celestial mechanics, billiard dynamics, projective and conformal regularization, and rigidity theory. Recent work has clarified both the classical integrable examples—primarily conic boundaries focused at the Kepler center—and strong non-integrability results, including a zero-energy classification that identifies ellipses with the force center at a focus as the only smooth integrable planar tables under a natural \(K\)-convexity hypothesis [2507.15879].

## 1. Definition and basic dynamical structure

A Kepler billiard in the plane consists of a smooth, simple, closed curve \(\Gamma\subset \mathbb R^2\) bounding a domain, together with a distinguished Kepler center \(O\), typically placed at the origin, and the Hamiltonian
\[
H(p,q)=\frac{|p|^2}{2}-\frac{\mu}{|q|}, \qquad \mu>0,
\]
defined on \(\mathbb R^2\setminus\{O\}\) [2507.15879]. Between reflections, the particle follows a Keplerian conic arc: ellipse, parabola, or hyperbola depending on the energy. At the boundary, the impact is elastic: the angle of incidence equals the angle of reflection, equivalently the tangential component of velocity is preserved and the normal component is reversed [2507.15879].

For positive-energy formulations used in high-energy rigidity, the center may be written as \(c\in\Omega\) and the Hamiltonian as
\[
H(q,p)=\tfrac12|p|^2-\frac{\mu}{|q-c|},
\]
with the billiard map encoded on \(\partial\Omega\times(0,\pi)\) by arc-length and incidence-angle coordinates; in that formulation the map is area-preserving [2507.08446]. This is the natural Keplerian analogue of the Birkhoff billiard map.

The zero-energy level \(H=0\) is distinguished because all free-flight arcs are parabolas focused at the Kepler center [2507.15879]. In contrast, for fixed negative energy, the free flights are elliptic arcs with a common semi-major axis \(a=\frac{m\mu}{-2E}\) in the focused-conic integrable setting [2312.05542]. This energy dependence is central to the modern theory: zero energy is amenable to exact regularization into free motion, while high positive energies admit asymptotic comparison with punctured Birkhoff billiards [2507.15879; 2507.08446].

## 2. Classical integrable examples with conic boundaries

The standard integrable Kepler billiards arise when the reflection wall is a conic section with one focus at the Kepler center. In Cartesian coordinates, a focused ellipse or hyperbola may be written as
\[
\frac{(x-c_K)^2}{a_K^2}+\frac{y^2}{a_K^2-c_K^2}=1,
\]
with foci \(F=(0,0)\) and \(F'=(2c_K,0)\) in the elliptic case [2312.05542]. In the broader conformal formulation, the wall is a single branch of a conic with one focus at the origin, and every ellipse or hyperbola focused at the origin defines an integrable Kepler billiard [2110.03376].

These systems possess a second first integral beyond energy. One formulation uses angular momentum \(L\) and the Runge–Lenz vector \(A\). For a focused conic wall, the invariant can be written as
\[
D=L^2-2c_KA_1,
\]
where \(A_1\) is the appropriate component of the Runge–Lenz vector [2312.05542]. An equivalent Cartesian expression is
\[
I=(q_1p_2-q_2p_1)^2-2c\Big[(q_1p_2-q_2p_1)p_2-\frac{\kappa q_1}{|q|}\Big],
\]
which is preserved both by the Kepler flow and by elastic reflections at the focused conic [2110.03376].

A geometric hallmark of integrability is the “second-focus circle.” If \(\mathcal K_i\) are the successive Kepler ellipses between reflections, each with one focus at the Kepler center and second focus \(F_i\), then all \(F_i\) lie on a fixed circle centered at the second focus \(F'\) of the reflecting conic:
\[
|F'-F_i|=R=\mathrm{const.}
\]
[2312.05542]. This converts the extra first integral into an explicitly geometric invariant. The same paper shows that the family of flight ellipses is confined between two confocal conics \(E_\pm\), providing caustic-type envelope curves, and analyzes the envelope of the directrices of the flight ellipses [2312.05542].

A plausible implication is that focused conics furnish the canonical Liouville-integrable model for one-center Kepler billiards in much the same way that confocal quadrics do for free Birkhoff billiards. That interpretation is consistent with later rigidity results identifying ellipses as the only smooth integrable planar zero-energy tables under natural convexity assumptions [2507.15879].

## 3. Zero-energy regime and Levi–Civita regularization

At zero energy, the planar Kepler billiard admits a particularly transparent regularization through the Levi–Civita map
\[
q=z^2,\qquad p=\frac{w}{2\bar z},
\]
combined with the rescaled Hamiltonian
\[
\widetilde H(q,p)=|q|\,H(q,p),
\]
which becomes
\[
K(z,w)=\frac{|w|^2}{8},\qquad K=1
\]
on the regularized side [2507.15879]. Under this transformation, focused parabolic Kepler arcs in the \(q\)-plane lift to straight line segments in the \(z\)-plane, and the reflection law is conformally preserved [2507.15879].

The boundary \(\Gamma\) lifts to the two-sheeted preimage
\[
\widetilde\Gamma=\{z: z^2\in\Gamma\},
\]
which is centrally symmetric. If every focused parabola meets \(\Gamma\) in at most two points—the \(K\)-convexity condition—then every line meets \(\widetilde\Gamma\) in at most two points, which implies strict convexity of \(\widetilde\Gamma\) [2507.15879]. Thus the zero-energy Kepler billiard becomes a standard Birkhoff billiard in a centrally symmetric strictly convex domain.

This reduction is the foundation of the zero-energy classification theorem. If the original Kepler billiard admits a nontrivial \(C^1\) first integral independent of the Hamiltonian, then the pulled-back Birkhoff billiard also admits a nontrivial \(C^1\) first integral [2507.15879]. The problem is thereby transferred from a singular central-force billiard to an ordinary free billiard.

The special role of zero energy had already appeared in conformal correspondence results relating Kepler and Hooke dynamics by the square map \(q=z^2\), but the 2025 zero-energy theorem makes this regularization the core of a rigidity argument rather than merely a source of examples [2110.03376; 2507.15879].

## 4. Zero-energy rigidity and classification

The principal rigidity result states that if \(\Gamma\subset\mathbb R^2\) is a smooth \((C^\infty)\), simple, closed curve satisfying the \(K\)-convexity condition that every Kepler parabolic orbit focused at the origin meets \(\Gamma\) in at most two points, and if the zero-energy Kepler billiard inside \(\Gamma\) admits a \(C^1\) first integral not depending only on the Hamiltonian, then \(\Gamma\) must be an ellipse with one focus at the Kepler center [2507.15879].

The proof is a direct translation of a theorem of Bialy–Mironov via the Levi–Civita map. Once the lifted boundary \(\widetilde\Gamma\) is shown to be a strictly convex, centrally symmetric, \(C^\infty\) curve, the Bialy–Mironov rigidity theorem applies: a Birkhoff billiard in \(\mathbb R^2\) bounded by a \(C^1\), strictly convex, centrally symmetric curve and admitting a nontrivial \(C^1\) first integral near the boundary must have elliptical boundary [2507.15879]. Since \(\widetilde\Gamma\) is then an ellipse in the \(z\)-plane, its image under \(z\mapsto z^2\) is an ellipse in the original \(q\)-plane, and one focus is the origin [2507.15879].

This settles the zero-energy case of the conjectural classification that the only integrable planar Kepler billiards are the classical focused conics [2507.15879]. The result is exact rather than perturbative, requires only smoothness and \(K\)-convexity, and does not depend on analyticity. It also sharpens the geometric meaning of zero-energy integrability: any nontrivial first integral forces a hidden free-billiard structure whose only compatible smooth geometry is elliptic.

A common misconception is that the presence of a Kepler potential should enlarge the class of integrable boundaries because Kepler motion itself is maximally superintegrable. The zero-energy classification shows the opposite in the smooth planar setting: the combination of Kepler arcs and elastic reflection is sufficiently rigid that integrability collapses to the focused ellipse [2507.15879].

## 5. High-energy non-integrability and the Keplerian Birkhoff conjecture

A complementary line of work treats analytic integrability at high positive energies. For a strictly convex planar domain \(\Omega\) with real-analytic boundary, there exists \(h_0\) such that for all \(h\ge h_0\), if \(\partial\Omega\) is not an ellipse then for every attractive center \(c\in\Omega\) except at most one point, the Kepler billiard at energy \(h\) does not admit an analytic first integral; if \(\partial\Omega\) is an ellipse, then exactly the two foci are analytically integrable at all energies [2507.08446].

The mechanism is symbolic dynamics. The generating functions for large \(h\) split into “direct” and “indirect” branches,
\[
S_d(\xi,\eta;h)=\sqrt h\,|\gamma(\xi)-\gamma(\eta)|+O(h^{-1/2}),
\]
\[
S_i(\xi,\eta;h)=\sqrt h\bigl(|\gamma(\xi)|+|\gamma(\eta)|\bigr)+O(h^{-1/2}),
\]
with further clockwise and anti-clockwise branches near antipodal configurations [2507.08446]. In the high-energy limit, the Kepler billiard is approximated by a “punctured Birkhoff billiard” in which line segments may pass through the center. From this limit one constructs three kinds of chaotic building blocks: those based on first-kind points, nondegenerate punctured triangles, and degenerate punctured triangles [2507.08446].

In each case, one obtains an invariant set on which the induced dynamics is semi-conjugate to the full two-symbol Bernoulli shift. Positive topological entropy follows, excluding analytic integrability [2507.08446]. This gives a partial affirmative answer to a Keplerian analogue of the classical Birkhoff–Poritsky conjecture.

The same work introduces focal points of the second kind, defined by the property that every billiard ray through \(c\) returns to \(c\) after two elastic bounces. In a real-analytic, strictly convex table there can be at most two such points; if there are exactly two, the boundary must be an ellipse and they are its classical foci; in a non-elliptic analytic table there is at most one [2507.08446]. The existence of an infinite-dimensional family of non-elliptic analytic tables with a single focal point of the second kind explains why the theorem allows at most one exceptional center rather than none [2507.08446].

Taken together with the zero-energy rigidity theorem, these results delineate two complementary regimes. At zero energy, smooth integrability forces an ellipse with a focal center [2507.15879]. At sufficiently high positive energy, analytic integrability is again essentially restricted to ellipses with the center at a focus, up to at most one exceptional center in the non-elliptic case [2507.08446].

## 6. Geometric invariants, projective correspondences, and extensions

The extra first integral \(D=C^2-2hA_1\) reappears in several broader contexts, linking Kepler billiards with projective dynamics and secular celestial mechanics [2504.17645]. In that formulation, \(C\) denotes angular momentum, \(A_1\) a component of the Laplace–Runge–Lenz vector, and \(h\) the focus–wall distance. The same quantity arises in the two-center problem and in a partially averaged system in the secular three-body problem, providing a structural bridge between billiard integrability and classical perturbation theory [2504.17645].

Projective methods also produce integrable Kepler billiards on surfaces of constant curvature. For \(M=\mathbb R^2\), \(S^2\), or \(H^2\), with Hamiltonian
\[
H(x,p)=\tfrac12|p|_g^2-\mu/r,
\]
and reflecting wall taken from a confocal family, the resulting billiard is Liouville-integrable with independent integrals \(H\) and \(D=C^2-2hA_1\) [2504.17645]. This extends the planar focused-conic picture to spherical and hyperbolic geometries.

A parallel generalization concerns higher-dimensional space forms and two-center mechanical billiards. In \(\mathbb R^n\), \(S^n\), and \(H^n\), billiards with two Kepler centers and finite unions of confocal quadrics as walls admit a complete involutive set of integrals \(\{I_1,I_2,C_3,\dots,C_n\}\) [2303.12443]. Although these are two-center rather than one-center systems, they show that the confocal-quadric paradigm persists in higher-dimensional mechanical billiards.

In dimension three, quaternionic regularization via the Kustaanheimo–Stiefel transformation identifies integrable one-center zero-energy Kepler billiards with \(S^1\)-invariant Hooke billiards in four dimensions. The resulting admissible reflection walls in \(\mathbb R^3\) are planes, centered spheres, spheroids with focus at the origin, one sheet of a two-sheeted circular hyperboloid with focus at the origin, and paraboloids with focus at the origin [2311.08920]. This gives a geometric explanation for the rotationally invariant families of known three-dimensional integrable Kepler billiards.

These extensions suggest that the integrable core of the theory is governed less by the specific Newtonian singularity than by regularization-compatible confocal geometry. That is an inference rather than an explicit theorem, but it is consistent with the repeated appearance of conformal, projective, and symplectic reductions across the planar, curved, and higher-dimensional settings [2504.17645; 2303.12443; 2311.08920].

## 7. Poncelet-type geometry, variants, and chaotic regimes

The integrable planar elliptic Kepler billiard has a refined Poncelet-type structure. For non-zero-energy orbits in a billiard bounded by an ellipse whose one focus is the Kepler center, the lines joining consecutive second foci are tangent to a fixed circle [2511.06464]. The dynamics can be linearized on an elliptic curve, with one reflection corresponding to translation by a fixed divisor class, and \(n\)-periodicity is governed by Cayley-type determinant conditions [2511.06464]. This places elliptic Kepler billiards within the algebro-geometric tradition of Poncelet porisms and elliptic-curve linearization.

A distinct variant replaces reflection by refraction. In a two-zone model with an inner Kepler region and outer harmonic region, trajectories alternate Keplerian hyperbolae and harmonic ellipses, with the interface obeying a generalized Snell law
\[
\sqrt{V_I(\bar z)}\,\sin(\alpha_I)=\sqrt{V_E(\bar z)}\,\sin(\alpha_E)
\]
[2108.11159]. For a circular interface the return map is integrable with a global action \(I\), and under small perturbations KAM tori and Aubry–Mather sets persist [2108.11159]. Although these are not reflective Kepler billiards in the strictest sense, they show how Keplerian flight segments remain compatible with twist-map methods once the boundary interaction is modified.

Three-dimensional non-integrable behavior has also been studied directly. For reflective and refractive Keplerian billiards in \(\mathbb R^3\), the existence of two distinct non-antipodal nondegenerate central configurations on the boundary implies topological chaos for sufficiently large inner energy, via a symbolic dynamics built from homothetic radial trajectories [2410.15935]. A triaxial ellipsoid provides such a chaotic example, despite the complete integrability of the zero-potential ellipsoidal billiard [2410.15935].

These developments clarify a broad point. Integrable Kepler billiards are not generic mechanical billiards with a central potential; they are special constructions tied to focused conics, confocal families, or exact regularization schemes. Outside those structures, the dominant phenomena are positive entropy, symbolic dynamics, broken invariant curves, and KAM/Aubry–Mather remnants rather than Liouville integrability [2507.08446; 2410.15935; 2108.11159].

Source: https://www.emergentmind.com/topics/kepler-billiards