skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: The Dynamics of Digits: Calculating Pi with Galperin’s Billiards
In Galperin billiards, two balls colliding with a hard wall form an analog calculator for the digits of the number π . This classical, one-dimensional three-body system (counting the hard wall) calculates the digits of π in a base determined by the ratio of the masses of the two particles. This base can be any integer, but it can also be an irrational number, or even the base can be π itself. This article reviews previous results for Galperin billiards and then pushes these results farther. We provide a complete explicit solution for the balls’ positions and velocities as a function of the collision number and time. We demonstrate that Galperin billiard can be mapped onto a two-particle Calogero-type model. We identify a second dynamical invariant for any mass ratio that provides integrability for the system, and for a sequence of specific mass ratios we identify a third dynamical invariant that establishes superintegrability. Integrability allows us to derive some new exact results for trajectories, and we apply these solutions to analyze the systematic errors that occur in calculating the digits of π with Galperin billiards, including curious cases with irrational number bases.  more » « less
Award ID(s):
1912542
PAR ID:
10158998
Author(s) / Creator(s):
; ; ; ; ;
Date Published:
Journal Name:
Mathematics
Volume:
8
Issue:
4
ISSN:
2227-7390
Page Range / eLocation ID:
509
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. This work is motivated by an article by Wang, Casati, and Prosen[Phys. Rev. E vol. 89, 042918 (2014)] devoted to a study of ergodicityin two-dimensional irrational right-triangular billiards. Numericalresults presented there suggest that these billiards are generally notergodic. However, they become ergodic when the billiard angle is equalto \pi/2 π / 2 times a Liouvillian irrational, morally a class of irrational numberswhich are well approximated by rationals. In particular, Wang etal. study a special integer counter that reflects the irrationalcontribution to the velocity orientation; they conjecture that thiscounter is localized in the generic case, but grows in the Liouvilliancase. We propose a generalization of the Wang-Casati-Prosen counter:this generalization allows to include rational billiards intoconsideration. We show that in the case of a 45°\!\!:\!45°\!\!:\!90° 45 ° : 45 ° : 90 ° billiard, the counter grows indefinitely, consistent with theLiouvillian scenario suggested by Wang et al. 
    more » « less
  2. Hybrid (and impact) systems are dynamical systems experiencing both continuous and discrete transitions. In this work, we derive necessary and sufficient conditions for when a given differential form is invariant, with special attention paid to the case of the existence of invariant volumes. Particular attention is given to impact systems where the continuous dynamics are Lagrangian and subject to nonholonomic constraints. A celebrated result for volume-preserving dynamical systems is Poincaré recurrence. In order to be recurrent, trajectories need to exist for long periods of time, which can be controlled in continuous-time systems through e.g. compactness. For hybrid systems, an additional mechanism can occur which breaks long-time existence: Zeno (infinitely many discrete transitions in a finite amount of time).We demonstrate that the existence of a smooth invariant volume severely inhibits Zeno behavior; hybrid systems with the “boundary identity property” along with an invariant volume-form have almost no Zeno trajectories (although Zeno trajectories can still exist). This leads to the result that many billiards (e.g. the classical point, the rolling disk, and the rolling ball) are recurrent independent on the shape of the compact table-top. 
    more » « less
  3. A bstract A novel mechanism, “catalyzed baryogenesis”, is proposed to explain the observed baryon asymmetry in our universe. In this mechanism, the motion of a ball-like catalyst provides the necessary out-of-equilibrium condition, its outer wall has CP-violating interactions with the Standard Model particles, and its interior has baryon number violating interactions. We use the electroweak-symmetric ball model as an example of such a catalyst. In this model, electroweak sphalerons inside the ball are active and convert baryons into leptons. The observed baryon number asymmetry can be produced for a light ball mass and a large ball radius. Due to direct detection constraints on relic balls, we consider a scenario in which the balls evaporate, leading to dark radiation at testable levels. 
    more » « less
  4. Abstract We show that for any setDof at least two digits in a given baseb, almost all even integers taking digits only inDwhen written in basebsatisfy the Goldbach conjecture. More formally, if$$\mathcal {A}$$is the set of numbers whose digits basebare exclusively fromD, almost all elements of$$\mathcal {A}$$satisfy the Goldbach conjecture. Moreover, the number of even integers in$$\mathcal {A}$$which are less thanXand not representable as the sum of two primes is less than$$|\mathcal {A}\cap \{1,\ldots ,X\}|^{1-\delta }$$. 
    more » « less
  5. This paper shows how techniques for linear dynamical systems can be used to reason about the behavior of general loops. We present two main results. First, we show that every loop that can be expressed as a transition formula in linear integer arithmetic has a best model as a deterministic affine transition system. Second, we show that for any linear dynamical system f with integer eigenvalues and any integer arithmetic formula G, there is a linear integer arithmetic formula that holds exactly for the states of f for which G is eventually invariant. Combining the two, we develop a monotone conditional termination analysis for general loops. 
    more » « less