The characteristic motion of a disk as it spins and rolls without slipping on a surface before coming to rest is a daily-life phenomenon that exhibits interesting dynamics. In this work, we explore an analogy between the disk coming to rest and the merging of two black holes. By comparing the analytical descriptions of both systems, we highlight how the precession frequency of the disk and the orbital frequency of black holes exhibit a similar increase over time. Experimentally, we investigate the motion of the disk coming to rest by measuring the vertical vibrations induced in the supporting platform. These vibrations are related to the precession angular velocity of the disk. Notably, these vibrations share common features with the gravitational waves emitted by merging black holes, including their frequency and amplitude increasing with time. We conclude by examining both the strengths and limitations of this analogy and by presenting an open challenge for students to further explore these connections.

Many of us have observed a spinning coin coming to rest on a table and have probably been captivated by its relatively long-lasting dynamics accompanied by a ringing with progressively increasing frequency. A disk coming to rest on a solid surface typically exhibits both precessional and rotational motion. The precession frequency progressively increases as the angle the disk forms with the platform decreases until the disk reaches full contact with the platform. This phenomenon has received increased attention since the promotion of Euler's disk, an educational toy with minimal dissipation that can roll and spin for several minutes. In the last few decades, the dynamics of a spinning disk have been the object of theoretical1,2 and experimental3–5 research works.

The motion of a spinning disk coming to rest presents some analogies with another phenomenon that occurs on a completely different scale: the merging of two compact astrophysical bodies.6 The analogy between spinning disks and merging black holes was first mentioned by Bildsten2 in his reply to Moffat's model1 of Euler's disk. The scientific community's interest in merging events has significantly increased following the first detection of the gravitational waves generated by two merging black holes.7 This historic measurement confirmed once again Einstein's theory of general relativity by detecting the gravitational wave produced during the final stages of the merger.

While previous educational and research works investigated either the disk dynamics or the emitted sound, here we measure and analyze the vertical vibration of the platform on which the disk spins and rolls. The experimental configuration proposed here is especially simple, and relevant signals can be acquired with a smartphone. This configuration has been developed as a class project by a group of master's students in the context of an advanced course in experimental physics. The manuscript is organized as follows. We first lay out analytical descriptions of the spinning disk coming to rest and of black hole mergers and highlight the similarities and differences between the systems (Sec. II). We experimentally show how the platform vibration is related to the precession of the disk (Sec. III) and compare it to the gravitational signal created by the merging of two black holes (Sec. IV). This work could serve as a reference for physics instructors in mechanics and astrophysics and help the students develop the ability to recognize analogies.

This section presents the analytical descriptions of the mechanics of a spinning disk coming to rest and that of two merging black holes. The schematics of the two phenomena are presented in Fig. 1. A side-by-side comparison of the main quantities derived in this section can be found in Table I.

Fig. 1.

Schematics of a spinning disk (a), where Ω and ωr are the precession and rotation angular velocities, respectively, and the merger of two black holes (b).

Fig. 1.

Schematics of a spinning disk (a), where Ω and ωr are the precession and rotation angular velocities, respectively, and the merger of two black holes (b).

Close modal
Table I.

Comparison of the relevant quantities of a disk coming to rest and two black holes merging.

Coming to rest of a spinning disk Merger of two orbiting black holes
Energy conservation, dEdisk/dt+Pdiss=0  Energy conservation, dEorb/dt+Prad=0 
Disk energy, Edisk=32mgRα  Orbital energy, Eorb=12GμM/l 
Dissipated power, Pdiss=kαγ  Gravitational radiation, Prad=325G4μ2M3/c5l5 
Coming to rest time, tc=3α0γ+1mgR/2k(γ+1)  Coalescence time, tc=5256c5l04/G3μM2 
Contact angle, α=α0(1t/tc)1/γ+1  Orbital separation, l=l0(1t/tc)1/4 
Precession angular velocity,  Orbital angular velocity, 
Ω=Ω0(1t/tc)1/2(γ+1)  ωK=ωK0(1t/tc)3/8 
Coming to rest of a spinning disk Merger of two orbiting black holes
Energy conservation, dEdisk/dt+Pdiss=0  Energy conservation, dEorb/dt+Prad=0 
Disk energy, Edisk=32mgRα  Orbital energy, Eorb=12GμM/l 
Dissipated power, Pdiss=kαγ  Gravitational radiation, Prad=325G4μ2M3/c5l5 
Coming to rest time, tc=3α0γ+1mgR/2k(γ+1)  Coalescence time, tc=5256c5l04/G3μM2 
Contact angle, α=α0(1t/tc)1/γ+1  Orbital separation, l=l0(1t/tc)1/4 
Precession angular velocity,  Orbital angular velocity, 
Ω=Ω0(1t/tc)1/2(γ+1)  ωK=ωK0(1t/tc)3/8 

We describe how to derive expressions for the disk's precession angular velocity Ω and for the black holes' orbital angular velocity ωK as a function of time, which are at the core of our analogy. The descriptions are developed in terms of the variation of the systems' total energy as a function of time. The total energy of the disk during its motion, which combines rotation about its symmetry axis and precession, depends on the angle, α(t), it forms with the underlying surface, which decreases over time due to dissipation. This angle can be related to the disk's precession angular velocity. On the contrary, the total energy of the black hole system depends on the distance between them, l(t), which decreases over time due to the radiation of gravitational waves and is related to the system's orbital angular velocity. The power laws describing Ω(t) and ωK(t) exhibit singular behaviors corresponding to the moments when the disk comes to rest and the two black holes merge, respectively. The power laws are characterized by their exponents that depend on the kind of dissipation involved: for the case of the black hole merger it is very well described by the quadrupole mass moment approximation (see Ferrari et al.6) resulting in a tiny range of variation for the associated exponent, whereas different dissipation mechanisms could be involved with the disk dynamics, including friction at the contact point and air resistance, resulting in a much wider range of variation for the measured exponent. While the black holes emit gravitational waves, the disk produces vibrations on the underlying platform, the measurement of which we present in Sec. III.

To study the motion of a spinning disk coming to rest (Fig. 1(a)), we consider an infinitely thin disk that spins on a platform without slipping. More detailed descriptions and general conditions of the disk's motion can be found in the literature.1,2 We consider the case where the disk precesses with angular velocity Ω around a vertical axis passing through its center of mass, and rotates with angular velocity ωr about its symmetry axis normal to the disk. It thus has total angular velocity ω=Ω+ωr. The equation accounting for the energy variation is
(1)
where Edisk is the mechanical energy of the disk and Pdiss is the dissipated power.
The mechanical energy is the sum of the potential and kinetic energies,
(2)
where m is the disk's mass, g is the gravitational acceleration, R is the disk radius, α is the contact angle with the surface, and I=mR2/4 is the disk's moment of inertia about its diameter. The total angular momentum is
(3)
where j is the unit vector colinear with ω, as defined in Fig. 1(a). The torque at the contact point due to the weight of the disk can be written as
(4)
where z is the vertical unit vector, as represented in Fig. 1(a). From Euler's second law of motion,
(5)
Since ω=Ωsinα, we obtain ω2=(4g/R)sinα and we can re-write the disk's total energy as
(6)
As documented in previous works,1,2 the dissipated power is expected to have the form
(7)
in which k is a constant that depends on the disk's geometry and on the properties of the surrounding fluid, and γ>0. Specifically, Moffat1 suggested γ=2 by considering the viscous dissipation in the boundary layer between the disk and the surface in the final stages of the motion. Subsequently, Bildsten2 argued that this model is valid only for a viscous boundary layer extending along the entire lower disk's surface, therefore for α<αc=4.8×103 rad, whereas if it does not, α>αc, and the dissipation is governed by γ=5/4.
We now consider small angles and thus approximate sinαα. From Eqs. (1), (6), and (7), we can write
(8)
which yields the solution
(9)
where α0=α(t=0) is the initial contact angle, while
(10)
The contact angle goes to zero when ttc, which corresponds to the disk coming to rest.
Considering that ωΩα, we obtain
(11)
Therefore,
(12)
in which Ω0=2g/Rα0. The disk's precession frequency thus diverges as the disk comes to rest.
We here report a summary of the dynamics of two orbiting black holes merging together.6 For our description, we locate the origin of the reference frame at the center of mass of the system (Fig. 1(b)). The energy variation during orbiting is
(13)
where Eorb is the mechanical energy of the system and Prad is the power radiated in the form of gravitational waves. The gravitational waves emitted are perturbations of the space-time, i.e., perturbations of the metric tensor in Einstein's theory of relativity, and are represented by the gravitational wave strain.8 An expression for the radiated power can be derived under the adiabatic assumption, which means that the system has the time to compensate the energy lost in gravitational waves by changing its orbital energy. By considering the quadrupole mass moment approximation, the following expression for the radiated power is obtained:6 
(14)
in which G is the gravitational constant, c is the speed of light, l=r1+r2 is the distance between the two black holes, and M=m1+m2 and μ=m1m2/M are the total and reduced masses, respectively. The kinetic energy is
(15)
where r12=m22l2/M2 and r22=m12l2/M2. The orbital frequency is derived from Kepler's law,
(16)
The potential energy is
(17)
Therefore,
(18)
and
(19)
Substituting Eq. (19) in Eq. (13) and using Eqs. (18) and (14), we obtain
(20)
Denoting with l0=l(t=0) the initial orbital separation between the two black holes, the solution of Eq. (20) reads
(21)
which can be rewritten as
(22)
where tc=5c5l04/(256G3μM2). By substituting Eq. (22) into Eq. (16), we obtain the orbital angular velocity
(23)
where ωK0=GM/l03 is the initial orbital angular velocity.

The orbital angular velocity diverges when ttc. Since gravitational waves in vacuum are not dispersive, the time evolution of their strain measured on Earth reflects the time variation of the orbital separation distance during the black hole merger. It is worth noting that, since the emission frequency is twice the orbital frequency,6 one gravitational wavelength is emitted in half an orbiting period.

The comparison between the main physical parameters for both the spinning disk and the merger of two black holes is reported in Table I. The disk's precession angular velocity and the orbital angular velocity of the black holes exhibit similar scaling with singularities at t=tc. The frequencies of both systems increase as the singularity is approached. While the scaling of the black hole merger is characterized by a well-known exponent (−3/8), the scaling of the disk generally depends on the type of dissipation considered via the parameter γ. This exponent will be experimentally determined in our system.

Our experiments focused on the vertical vibration of the platform on which the disk spins. The experimental setup is represented in Fig. 2(a). A custom-made aluminum disk with mass m=154.2 ± 0.1 g, thickness 12.7 ± 0.1 mm, and diameter 76.1 ± 0.1 mm is released by hand on the central area of a black acrylic platform with thickness 0.51 cm, length 60.0 cm, and width 29.0 cm (see the Figshare repository of supplementary materialsupplementary material for a sample video). The platform lies on two planar supports placed underneath its short sides. The platform's vertical acceleration is measured with a smartphone via its accelerometer with a sampling frequency of 405 Hz, and the application phyphox.9 The smartphone lies on the platform as depicted in Fig. 2(a). Prior to our experiments, we measured the free oscillations of the platform by hitting it at different positions and consistently found damped oscillations with frequency fp11 Hz, which is significantly lower than the frequency range relevant for our study, f[20,36] Hz.

Fig. 2.

Experimental setup and acquisition. (a) Schematic representation of the experimental setup, which includes a smartphone for measuring the platform vertical acceleration and a tablet for acquiring side-view videos of the spinning disk. (b) Comparison between contact point and platform accelerations. (c) Three snapshots (top row) and corresponding zoomed-in processed images (bottom row) from a sample side-view video (see the Figshare repository of supplementary materialsupplementary material for the video).

Fig. 2.

Experimental setup and acquisition. (a) Schematic representation of the experimental setup, which includes a smartphone for measuring the platform vertical acceleration and a tablet for acquiring side-view videos of the spinning disk. (b) Comparison between contact point and platform accelerations. (c) Three snapshots (top row) and corresponding zoomed-in processed images (bottom row) from a sample side-view video (see the Figshare repository of supplementary materialsupplementary material for the video).

Close modal

A typical acceleration signal acquired while the disk was spinning on the platform is reported in Fig. 2(b).

To check if the platform acceleration is a good proxy for the disk motion, we simultaneously acquired a side-view video of the disk's dynamics and tracked the motion of the disk's contact point with the platform. This requires appropriate contrast between the disk, the platform, and the background, obtained by illuminating a white panel behind the disk. The video was acquired with a tablet at 240 fps. A custom-made code for image analysis was then used to process each frame and track the position of the contact point. In this code, the contact point is defined as the center of the contact region, which is the black region that is imaged as connecting the disk to the platform. The width of the contact region ranges from 1 mm to 1 cm depending on the stage of the disk motion. In Fig. 2(c), we present three snapshots of a sample video and the corresponding processed images. A comparison between the contact point and platform accelerations (Fig. 2(b)) shows that they exhibit the same periodicity and similar amplitude trends. A complete precession of the disk thus corresponds to one platform oscillation as a result of the variation in time of the excitation point's position, which induces a vibration of the platform. Furthermore, the amplitudes of both signals increase as the disk approaches rest. While measuring the motion of the contact point requires the appropriate recording settings and image analysis we described, measuring the platform vibration only requires the smartphone to be placed on the platform itself. The latter allows for faster data acquisition over a consistent number of runs and is also more adapted to classroom demonstrations.

In Fig. 3, we show a side-by-side comparison between the vibration of the platform on which the spinning disk comes to rest (a)–(c) and the gravitational waves emitted during the merger of a binary black hole system7 (d)–(f). Despite these systems involving substantially different spatial and energy scales, we note that the timescales are accidentally comparable. Moreover, the time evolutions of the platform acceleration (a) and the gravitational wave strain (d) exhibit interesting similarities. Initially, both the amplitude and the oscillation frequency of the signal increase toward the singularity. For the disk, this phase corresponds to spinning with decreasing contact angle until full contact with the platform is reached. For the binary system, this phase includes the spiraling of the black holes toward the merger. In the final part of the evolution, both systems exhibit a ringdown characterized by amplitude decreasing over time, a feature that is not described by the models in Sec. II. For the disk system, the platform still vibrates after the disk achieves full contact because the oscillations generated previously require some time to be damped out.

Fig. 3.

Signal comparison. (a) Vibrational acceleration of the platform as a function of time. t=0 corresponds to the time at which the disk starts spinning. tc corresponds to the disk's settling time. (b) Precession angular velocity of the disk Ω as a function of time t. Points represent experimental data, while the line is a fit from Eq. (12). The measurement time tc considered is represented in (a). (c) Wavelet coefficients of the signal in (a). (d) Gravitational wave strain from Ref. 7. (e) Orbital angular velocity of the black holes ωK as a function of time t. Points represent experimental data, while the line is a fit from Eq. (23). The measurement time tc considered is represented in (d). (f) Wavelet coefficients of the signal in (d).

Fig. 3.

Signal comparison. (a) Vibrational acceleration of the platform as a function of time. t=0 corresponds to the time at which the disk starts spinning. tc corresponds to the disk's settling time. (b) Precession angular velocity of the disk Ω as a function of time t. Points represent experimental data, while the line is a fit from Eq. (12). The measurement time tc considered is represented in (a). (c) Wavelet coefficients of the signal in (a). (d) Gravitational wave strain from Ref. 7. (e) Orbital angular velocity of the black holes ωK as a function of time t. Points represent experimental data, while the line is a fit from Eq. (23). The measurement time tc considered is represented in (d). (f) Wavelet coefficients of the signal in (d).

Close modal

Scaling laws are presented in Figs. 3(b) and 3(e) for the disk and black hole system, respectively. The time interval that we analyze for each measurement is denoted tc, an example of which is represented in Figs. 3(a) and 3(d). In general, tc is the characteristic time for the disk to achieve full contact with the platform or for the black holes to merge, and depends on the values of contact angle α0 and separation distance l0, respectively, which we do not measure. For the gravitational signal, the fit of the data yields the exponent 0.35 ± 0.03 with 95% confidence bounds, in agreement with the exponent 3/8 expected from Eq. (23). For the disk, the fit of the data in Fig. 3(b) yields the exponent 0.13 ± 0.01 with 95% confidence bounds. We repeated the experiment 30 times and obtained the mean exponent 0.13 ± 0.03 with 95% confidence bounds. From this value and Eq. (12), we obtain γ=2.9 ± 0.9. This is in the same order of magnitude but generally larger than the values theoretically estimated by Moffatt1 (γ=2) and Bildsten2 (γ=5/4), which are derived in ideal conditions, i.e., no friction at the contact point is considered, which is not respected by our setup.

We have decomposed the energy content of the signal in time and frequency using a wavelet transform,10 a procedure commonly employed to analyze gravitational wave signals.7,11 Compared to classical Fourier analysis, the wavelet transform retains the time information, allowing for a precise frequency characterization of the signal at the various stages of the dynamics. In Figs. 3(c) and 3(f), we show the wavelet coefficients of the signals in Figs. 3(a) and 3(d), respectively.

Here, we use a mother wavelet known as “Morlet,” defined as a plane wave modulated by a Gaussian function.12 The wavelet coefficients represent the energy of the processed signals. Both wavelet transforms exhibit energy transfer to high frequencies as the singularity is approached. This energy transfer is more evident for the black hole merger, as it is expected from the larger exponent characterizing its scaling law. We note that noise appears in the wavelet transform of the disk, while it is absent for the gravitational signal. Indeed, background noise was removed from the raw gravitational signal,7 while for the disk, we considered the signal as acquired by the accelerometer.

We investigated an analogy between a disk spinning on a surface and a phenomenon occurring at a very different scale, i.e., the merging of two orbiting black holes. We presented a simple experimental configuration for measuring the vibration of the platform on which the disk spins. Our setup allows for straightforward signal acquisitions that do not require any post-processing. By simultaneously measuring the platform vibration and the disk contact point, we showed that the platform frequency corresponds to the precession frequency of the disk. As highlighted by the theoretical expressions presented in Sec. II, the disk precession frequency is expected to have a similar scaling to the frequency of the gravitational waves emitted during the merger. The platform perturbation was thus compared to the space-time perturbation measured by the LIGO detector. Both signals exhibit increasing frequency and amplitude and energy transfer to high frequency. The corresponding scaling laws are characterized by different exponents that are 0.13 in our system and 0.35 from the black hole data due to the different mechanics and dissipation laws involved.

This work originates from a class project developed by master's students in a Physics Department and can be proposed in advanced courses of experimental physics also by varying the experimental conditions. For the disk, the exponent depends upon the type of dissipation involved in the dynamics via the parameter γ, which appears in the expression of the dissipated power Pdiss=kαγ. The scaling laws of the two systems would be equivalent if γ=1/3. Students could thus be challenged to vary the experimental conditions to obtain an exponent closer to the scaling law of the black hole merger.

Please click on Figshare repository of supplementary materialthis link to access the supplementary material at the reviewer website. Video acquired at 240 fps and reproduced here at 30 fps (8 times slower). Print readers can see the supplementary material at DOI for figshare repositoryhttps://doi.org/10.60893/figshare.ajp.c.7829717.

The authors are grateful to Andrea De Luca, who accompanied us during a significant part of this project and for developing the code that tracks the motion of the contact point.

The authors have no conflicts to disclose.

1.
H. K.
Moffat
, “
Euler's disk and its finite-time singularity
,”
Nature
404
,
833
834
(
2000
).
2.
L.
Bildsten
, “
Viscous dissipation for Euler's disk
,”
Phys. Rev. E
66
(
5
),
056309
(
2002
).
3.
G.
Van den Engh
,
P.
Nelson
, and
J.
Roach
, “
Numismatic gyrations
,”
Nature
408
(
6812
),
540
540
(
2000
).
4.
K.
Easwar
,
F.
Rouyer
, and
N.
Menon
, “
Speeding to a stop: The finite-time singularity of a spinning disk
,”
Phys. Rev. E
66
(
4
),
045102
(
2002
).
5.
H.
Caps
,
S.
Dorbolo
,
S.
Ponte
,
H.
Croisier
, and
N.
Vandewalle
, “
Rolling and slipping motion of Euler's disk
,”
Phys. Rev. E
69
(
5
),
056610
(
2004
).
6.
V.
Ferrari
,
L.
Gualtieri
, and
P.
Pani
,
General Relativity and Its Applications. Black Holes, Compact Stars and Gravitational Waves
(
Taylor and Francis Group
,
LLC
., Milton Park, Abingdon, UK,
2021
).
7.
B. P.
Abbott
et al, “
Observation of gravitational waves from a binary black hole merger
,”
Phys. Rev. Lett.
116
(
6
),
061102
(
2016
).
8.
The gravitational wave strain, as measured on the earth surface, may be understood as (ss0)/s0, in which s0 is the distance between two generic points in space in unperturbed conditions and s is the same distance when the gravitational wave passes between these two points. For a detailed description on how the gravitational wave strain is measured, see the book by Ferrari, Gualtieri, and Pani.6 
10.
I.
Daubechies
,
Ten Lectures on Wavelets
(
Society for Industrial and Applied Mathematics
, Philadelphia, PA,
1992
).
11.
M.
Imbrogno
,
C.
Meringolo
, and
S.
Servidio
, “
Extreme gravitational interactions in the problem of three black holes in general relativity
,”
Class. Quantum Grav.
40
(
7
),
075008
(
2023
).
12.
R.
Büssow
, “
An algorithm for the continuous Morlet wavelet transform
,”
Mech. Syst. Sig. Process.
21
(
8
),
2970
2979
(
2007
).