Comprehensive Study Guide on Uniform Circular Motion and Universal Gravitation

Kinematics of Uniform Circular Motion

Uniform circular motion occurs when an object travels along a circular path at a constant speed. While the scalar speed remains unchanging, the direction of the velocity vector continuously rotates, meaning the object constantly undergoes acceleration directed toward the center of rotation.

The time needed to execute one complete revolution is the period, denoted by TT. The number of complete cycles per unit of time is defined as the rotational frequency ff, which is the reciprocal of the period:

f=1Tf = \frac{1}{T}

The orbital speed (or tangential velocity) vv represents the arc length traveled during one period divided by the period itself:

v=2π⋅rT=2π⋅r⋅fv = \frac{2\pi \cdot r}{T} = 2\pi \cdot r \cdot f

The rate at which the angular position changes over time is known as the angular velocity or angular frequency, denoted by ω\omega:

ω=ΔφΔt=2πT=2π⋅f\omega = \frac{\Delta \varphi}{\Delta t} = \frac{2\pi}{T} = 2\pi \cdot f

The SI unit of angular velocity is expressed as:

[ω]=rads[\omega] = \frac{\text{rad}}{\text{s}}

By substituting ω\omega into the equation for tangential speed, the direct coupling between linear velocity and angular velocity is established:

v=ω⋅rv = \omega \cdot r

Because the direction of the velocity vector changes continuously toward the rotation axis, an object in uniform circular motion experiences centripetal acceleration, denoted by aza_z. Its magnitude can be expressed through angular velocity, tangential velocity, and trajectory radius:

az=ω⋅v=ω2⋅r=v2ra_z = \omega \cdot v = \omega^2 \cdot r = \frac{v^2}{r}

Applications and Quantitative Analysis in Circular Kinematics

The kinematics of circular motion govern phenomena ranging from planetary orbits to laboratory separation techniques and atmospheric events.

The Earth travels practically uniformly on an approximately circular orbit around the Sun at an average orbital distance of r=1.496×1011 mr = 1.496 \times 10^{11}\,\text{m}. The distance ss traversed by the Earth during one full revolution around the Sun is the circumference of this orbit:

s=2π⋅r=2π⋅(1.496×1011 m)=9.400×1011 ms = 2\pi \cdot r = 2\pi \cdot (1.496 \times 10^{11}\,\text{m}) = 9.400 \times 10^{11}\,\text{m}

Given an orbital period of one year, defined as T=365.25 d=3.156×107 sT = 365.25\,\text{d} = 3.156 \times 10^7\,\text{s}, the orbital speed of the Earth is calculated as:

v=2π⋅rT=2π⋅(1.496×1011 m)3.156×107 s=29.79 km/sv = \frac{2\pi \cdot r}{T} = \frac{2\pi \cdot (1.496 \times 10^{11}\,\text{m})}{3.156 \times 10^7\,\text{s}} = 29.79\,\text{km/s}

The resulting centripetal acceleration holding the Earth in its path is:

a=v2r=(29.79×103 m/s)21.496×1011 m=5.930×10−3 m/s2a = \frac{v^2}{r} = \frac{(29.79 \times 10^3\,\text{m/s})^2}{1.496 \times 10^{11}\,\text{m}} = 5.930 \times 10^{-3}\,\text{m/s}^2

Astronaut undergoing high acceleration simulation in a training centrifuge

In human aerospace physiology, large centrifuges accommodating a person are deployed to simulate the intense accelerations experienced during rocket launches. To generate a simulated acceleration of az=3 ga_z = 3\,g (where standard gravitational acceleration is taken as g=9.81 m/s2g = 9.81\,\text{m/s}^2) at an operational arm radius of r=10 mr = 10\,\text{m}, the required frequency is derived from the acceleration relation:

az=v2r=4π2⋅r⋅f2a_z = \frac{v^2}{r} = 4\pi^2 \cdot r \cdot f^2

f=az4π2⋅r=3⋅9.81 m/s24π2⋅10 m=0.27 Hzf = \sqrt{\frac{a_z}{4\pi^2 \cdot r}} = \sqrt{\frac{3 \cdot 9.81\,\text{m/s}^2}{4\pi^2 \cdot 10\,\text{m}}} = 0.27\,\text{Hz}

Converting this frequency to revolutions per minute yields:

0.27 s−1×60 s/min=16 min−10.27\,\text{s}^{-1} \times 60\,\text{s/min} = 16\,\text{min}^{-1}

Theodor Svedberg working with an analytical ultracentrifuge

In 1925, Theodor Svedberg invented the analytical ultracentrifuge, an achievement recognized by the Nobel Prize in Chemistry in 1926. Modern ultracentrifuges spin rapidly enough to subject biological samples to centripetal accelerations reaching up to a=106 ga = 10^6\,g. For a test sample positioned at a radius of r=5.0 cm=0.05 mr = 5.0\,\text{cm} = 0.05\,\text{m}, the operational rotational frequency must be:

f=a4π2⋅r=106⋅9.81 m/s24π2⋅0.05 m=2.2 kHzf = \sqrt{\frac{a}{4\pi^2 \cdot r}} = \sqrt{\frac{10^6 \cdot 9.81\,\text{m/s}^2}{4\pi^2 \cdot 0.05\,\text{m}}} = 2.2\,\text{kHz}

Because centripetal acceleration depends on the parameters according to a=4π2⋅r⋅f2a = 4\pi^2 \cdot r \cdot f^2, doubling the rotation frequency causes a fourfold (4×4\times) increase in acceleration due to the quadratic dependency (a∝f2a \propto f^2), whereas doubling the radius results in a twofold (2×2\times) increase due to linear scaling (a∝ra \propto r).

Large-scale atmospheric phenomena also exhibit rotational dynamics. Consider a tropical cyclone modeled as a rigid, rotating air mass possessing a uniform angular velocity ω\omega across its entire structure. A meteorological station positioned r1=500 km=500×103 mr_1 = 500\,\text{km} = 500 \times 10^3\,\text{m} from the eye records a peripheral wind velocity of v1=230 km/h=63.9 m/sv_1 = 230\,\text{km/h} = 63.9\,\text{m/s}. The centripetal acceleration experienced by the air at this outer boundary is:

a=v12r1=(63.9 m/s)2500×103 m=8.16×10−3 m/s2a = \frac{v_1^2}{r_1} = \frac{(63.9\,\text{m/s})^2}{500 \times 10^3\,\text{m}} = 8.16 \times 10^{-3}\,\text{m/s}^2

The rotational period of the entire cyclone mass is:

T=2π⋅r1v1=2π⋅(500×103 m)63.9 m/s=4.92×104 s=13.7 hT = \frac{2\pi \cdot r_1}{v_1} = \frac{2\pi \cdot (500 \times 10^3\,\text{m})}{63.9\,\text{m/s}} = 4.92 \times 10^4\,\text{s} = 13.7\,\text{h}

Because angular velocity is constant throughout the rigid body (ω=const\omega = \text{const}), tangential speed scales directly with radial distance. At an intermediate distance of r2=250 km=250×103 mr_2 = 250\,\text{km} = 250 \times 10^3\,\text{m}, the local wind speed is:

v2=2π⋅r2T=r2⋅v1r1=(250×103 m)⋅(63.9 m/s)500×103 m=31.9 m/s=115 km/hv_2 = \frac{2\pi \cdot r_2}{T} = \frac{r_2 \cdot v_1}{r_1} = \frac{(250 \times 10^3\,\text{m}) \cdot (63.9\,\text{m/s})}{500 \times 10^3\,\text{m}} = 31.9\,\text{m/s} = 115\,\text{km/h}

In vehicular safety, a critical problem occurs when a driver unexpectedly encounters an obstacle ahead and must choose between emergency straight-line braking and initiating an evasive 180∘180^\circ turn. Assuming the tire-road interface provides the same maximum deceleration in braking as it does centripetal acceleration in cornering (a=aza = a_z):

The stopping distance under linear braking from initial speed v0v_0 to complete rest (v=0v = 0) is determined from kinematics:

v2=2⋅a⋅s+v02=0  ⟹  s=v022⋅av^2 = 2 \cdot a \cdot s + v_0^2 = 0 \implies s = \frac{v_0^2}{2 \cdot a}

Conversely, executing an evasive circular maneuver requires a turning radius of:

az=v02r  ⟹  r=v02aza_z = \frac{v_0^2}{r} \implies r = \frac{v_0^2}{a_z}

Comparing the two expressions reveals that s=r2s = \frac{r}{2}. Under identical acceleration limits, the straight-line braking distance is exactly half the radius of the turn. Consequently, emergency braking affords a significantly greater margin of safety for avoiding an impact than attempting to steer around the hazard.

Centripetal Force, Inertial Frames, and Dynamics

Newton's second law dictates that any acceleration requires an unbalanced net force acting in the direction of the acceleration vector. In uniform circular motion, this force is the centripetal force FzF_z:

Fz=m⋅az=m⋅ω⋅v=m⋅ω2⋅r=m⋅v2rF_z = m \cdot a_z = m \cdot \omega \cdot v = m \cdot \omega^2 \cdot r = m \cdot \frac{v^2}{r}

The centripetal force is directed toward the instantaneous center of the circular trajectory, perpendicular to the instantaneous velocity vector. It is the net physical inward force responsible for continuously deflecting the object's path into a curve.

An observer situated within the rotating system occupies a non-inertial (accelerated) reference frame. Within this rotating frame, the observer experiences the sensation of an outward inertial force termed the centrifugal force, denoted F⃗ZF′\vec{F}'_{ZF}. Its apparent magnitude equals that of the centripetal force:

FZF=m⋅ω2⋅rF_{ZF} = m \cdot \omega^2 \cdot r

The centrifugal force is not a real physical interaction but a fictitious or apparent force arising solely from the inertia of the body resisting rotational acceleration within the non-inertial frame of reference.

When a passenger car drives through a horizontal curve, the real physical centripetal force is supplied entirely by the static friction between the tires and the asphalt. If the car drives at excessive speed, the required centripetal force surpasses the maximum available static frictional force. As a consequence, the circular trajectory can no longer be sustained; the radius of curvature expands uncontrollably and the car slides outward off the roadway.

Consider a driver of mass m=70 kgm = 70\,\text{kg} negotiating a horizontal left-hand curve at v=120 km/h=33.3 m/sv = 120\,\text{km/h} = 33.3\,\text{m/s}. If the side bolster of the vehicle's seat exerts an inward normal force of Fz=600 NF_z = 600\,\text{N} on the driver toward the center of curvature, the radius of the curve is:

r=m⋅v2Fz=70 kg⋅(33.3 m/s)2600 N=130 mr = \frac{m \cdot v^2}{F_z} = \frac{70\,\text{kg} \cdot (33.3\,\text{m/s})^2}{600\,\text{N}} = 130\,\text{m}

In accident reconstruction, a driver named Claudio skidded out of a bend of radius r=50 mr = 50\,\text{m} and collided with a guardrail. The road conditions were dry, corresponding to a static friction coefficient of μHR=0.75\mu_{HR} = 0.75. While Claudio claimed he was not speeding, dynamic equilibrium dictates the threshold cornering speed before loss of control:

Fz=FHR  ⟹  m⋅v2r=μHR⋅m⋅gF_z = F_{HR} \implies \frac{m \cdot v^2}{r} = \mu_{HR} \cdot m \cdot g

v=μHR⋅r⋅g=0.75⋅50 m⋅9.81 m/s2=19.2 m/s=69 km/hv = \sqrt{\mu_{HR} \cdot r \cdot g} = \sqrt{0.75 \cdot 50\,\text{m} \cdot 9.81\,\text{m/s}^2} = 19.2\,\text{m/s} = 69\,\text{km/h}

Claudio must have been traveling at a minimum speed of 69 km/h69\,\text{km/h} at the moment static friction was overcome.

A person swinging a water bucket overhead in a vertical circle

In vertical circular motion, gravity contributes dynamically to the net force along the radial axis. Consider swinging a bucket filled with water (m=4.5 kgm = 4.5\,\text{kg}) on a 40 cm40\,\text{cm} rope, where the center of gravity of the bucket describes a vertical circle of radius r=1.2 mr = 1.2\,\text{m} at a constant speed of v=3.5 m/sv = 3.5\,\text{m/s}.

The constant centripetal force required to sustain this motion is:

Fz=m⋅v2r=4.5 kg⋅(3.5 m/s)21.2 m=46 NF_z = \frac{m \cdot v^2}{r} = \frac{4.5\,\text{kg} \cdot (3.5\,\text{m/s})^2}{1.2\,\text{m}} = 46\,\text{N}

The downward gravitational force acting on the bucket is:

Fg=m⋅g=4.5 kg⋅9.81 m/s2=44.1 NF_g = m \cdot g = 4.5\,\text{kg} \cdot 9.81\,\text{m/s}^2 = 44.1\,\text{N}

At the zenith (top of the loop), gravity points toward the center and helps provide the centripetal acceleration. The tension force Fs,topF_{s,\text{top}} exerted by the rope is therefore:

Fs,top=Fz−Fg=46 N−44.1 N=1.9 NF_{s,\text{top}} = F_z - F_g = 46\,\text{N} - 44.1\,\text{N} = 1.9\,\text{N}

At the nadir (bottom of the loop), the rope tension must support the entire weight of the bucket in addition to supplying the inward centripetal acceleration:

Fs,bottom=Fz+Fg=46 N+44.1 N=100.1 NF_{s,\text{bottom}} = F_z + F_g = 46\,\text{N} + 44.1\,\text{N} = 100.1\,\text{N}

Experimental Determination of Smartphone Accelerometer Position

Modern smartphones contain integrated micro-electro-mechanical systems (MEMS) accelerometers. In gaming and sensor applications (such as detecting device orientation and movement in Pokémon Go), knowing the exact coordinate position of the accelerometer chip inside the phone casing is vital. Because manufacturers do not publish these internal hardware coordinates, the sensor's physical position can be determined experimentally using a rotating turntable and centripetal acceleration.

Geometric configuration of smartphone on turntable measuring edge distance x

The measurement procedure utilizes the sensor logging application "PhyPhox":

  1. Download and open "PhyPhox". Navigate to "Acceleration (without g)" and select the "Amount" (Betrag) tab.

  2. Verify calibration by ensuring the displayed absolute acceleration reads approximately 0 m/s20\,\text{m/s}^2 while the smartphone lies flat and motionless on a level table.

  3. Position the smartphone on the turntable aligned adjacent to its outer edge. Use a ruler to measure the distance xx from the central axis of rotation to the nearest edge of the phone casing. In this experimental setup, x=46 mmx = 46\,\text{mm}.

  4. Begin data recording by pressing the play icon.

  5. Engage the turntable at a preset rotation speed, selected here as 45 rpm45\,\text{rpm} (alternative standard speeds include 33 rpm33\,\text{rpm} and 78 rpm78\,\text{rpm}).

  6. Once the rotation stabilizes and adequate data is gathered, stop the turntable, retrieve the smartphone, pause the recording, and review the acceleration plot.

Acceleration vs time graph recorded in PhyPhox showing plateau at 2.0 m/s^2

Inspection of the graph during steady-state rotation yields a measured centripetal acceleration plateau of az=2.0 m/s2a_z = 2.0\,\text{m/s}^2.

The rotational frequency of the 45 rpm45\,\text{rpm} setting is:

f=4560 s=0.75 Hzf = \frac{45}{60\,\text{s}} = 0.75\,\text{Hz}

The corresponding angular velocity is:

ω=2π⋅f=2π⋅0.75 Hz=4.71 rad/s\omega = 2\pi \cdot f = 2\pi \cdot 0.75\,\text{Hz} = 4.71\,\text{rad/s}

Applying the centripetal acceleration formula gives the total distance rr from the turntable center to the internal accelerometer sensor:

az=ω2⋅r  ⟹  r=azω2=2.0 m/s2(4.71 rad/s)2=0.090 m=90 mma_z = \omega^2 \cdot r \implies r = \frac{a_z}{\omega^2} = \frac{2.0\,\text{m/s}^2}{(4.71\,\text{rad/s})^2} = 0.090\,\text{m} = 90\,\text{mm}

Subtracting the measured spacing xx gives the internal depth dd of the sensor relative to the smartphone edge:

d=r−x=90 mm−46 mm=44 mmd = r - x = 90\,\text{mm} - 46\,\text{mm} = 44\,\text{mm}

Repeating the procedure with the smartphone rotated by 90∘90^\circ provides the perpendicular coordinate, thereby identifying the sensor's precise two-dimensional location on the device's internal circuit board.

Mechanics and Experimental Verification of Bicycle Cornering

When a cyclist navigates a curve, the bicycle must lean inward at an angle α\alpha relative to the vertical. Treating the bicycle and rider as a single lumped mass in uniform circular motion, stable cornering requires the net contact force from the ground to align along the axis connecting the tire contact patch to the system's center of mass. Any misalignment generates an unbalanced torque that causes the cyclist to tip over.

Cyclist leaning into a turn with dashed line connecting ground contact to center of massFree-body diagram of leaning cyclist showing normal force, gravity, friction, and centripetal force

The force equilibrium is resolved into vertical and horizontal components:

Along the vertical (Y) axis, the upward normal force FNF_N balances the downward gravitational force FgF_g:

FN=Fg  ⟹  Fres,y=0F_N = F_g \implies F_{\text{res},y} = 0

Along the horizontal (X) axis, the inward static friction force FHRF_{HR} provides the centripetal acceleration:

FHR=Fres,x=FzF_{HR} = F_{\text{res},x} = F_z

The static friction force FHRF_{HR} adjusts to match the required centripetal demand; it is not automatically equal to the maximum static friction limit given by μHR⋅FN\mu_{HR} \cdot F_N.

From the vector geometry, the tangent of the lean angle α\alpha equals the ratio of the horizontal friction force to the vertical normal force:

tan⁡(α)=FHRFN=FzFg=m⋅v2rm⋅g=v2r⋅g\tan(\alpha) = \frac{F_{HR}}{F_N} = \frac{F_z}{F_g} = \frac{\frac{m \cdot v^2}{r}}{m \cdot g} = \frac{v^2}{r \cdot g}

Solving for the lean angle yields:

α=arctan⁡(v2r⋅g)\alpha = \arctan\left(\frac{v^2}{r \cdot g}\right)

Bird's eye view sketch of experimental circular track showing 8m radius, 5m timing segment, and camera position

An outdoor experiment empirically validates this theoretical model using the following materials and procedures:

  • Measuring tape (minimum length 8 m8\,\text{m})

  • Chalk

  • Smartphone camera with high-speed burst mode

  • Bicycle or electric scooter equipped with a digital speedometer

  • Stopwatch or timing application

The protocol is conducted by student teams of 4 (with 5 being an exception):

  1. Using chalk, a quarter-circle curve of radius r=8 mr = 8\,\text{m} is drawn on asphalt, leaving enough straight runway before and after for acceleration and deceleration. Within this circular arc, a dedicated segment of arc length Δs=5 m\Delta s = 5\,\text{m} is measured and marked.

  2. One observer stands with a clear view of the 5 m5\,\text{m} segment with a stopwatch. Another positions themselves facing the middle of the quarter circle with a camera to photograph the cyclist straight-on at the curve midpoint.

  3. The coordinator signals when the course is clear. The cyclist accelerates to a target speed of approximately 12 km/h12\,\text{km/h} (v1v_1) prior to entering the bend and maintains this speed steadily throughout the arc.

  4. The timer records the elapsed transit time Δt\Delta t over the marked 5 m5\,\text{m} path to calculate the true experimental velocity:

v=ΔsΔtv = \frac{\Delta s}{\Delta t}

  1. The photographer uses burst shooting to capture the cyclist at the exact midpoint, ensuring the lean angle α\alpha can be measured accurately against the vertical.

  2. The run is repeated at higher target speeds of approximately 18 km/h18\,\text{km/h} (v2v_2) and 24 km/h24\,\text{km/h} (v3v_3).

  3. The theoretical angles α1,α2,α3\alpha_1, \alpha_2, \alpha_3 are calculated from the measured speeds using the theoretical formula, and the actual angles are extracted from the photos using angle-measurement software or a triangle ruler.

  4. The absolute angular error is evaluated across all three velocities:

Δα=∣αcalc−αmeas∣\Delta \alpha = |\alpha_{\text{calc}} - \alpha_{\text{meas}}|

Experimental quality is classified based on angular deviation thresholds: an error Δα<4∘\Delta \alpha < 4^\circ indicates high measurement precision; an error between 4∘4^\circ and 8∘8^\circ represents acceptable accuracy with minor procedural errors; and an error Δα>8∘\Delta \alpha > 8^\circ indicates unacceptable discrepancies caused by speed variations, parallax, or improper angular evaluation.

Universal Gravitation

Newton's Law of Universal Gravitation states that any two bodies with masses m1m_1 and m2m_2 attract each other with a gravitational force FGF_G directed along the line connecting their centers of mass:

FG=G⋅m1⋅m2r2F_G = G \cdot \frac{m_1 \cdot m_2}{r^2}

The universal gravitational constant is established as:

G=6.6743×10−11 m3kg⋅s2=6.67×10−11 N⋅m2kg2G = 6.6743 \times 10^{-11}\,\frac{\text{m}^3}{\text{kg}\cdot\text{s}^2} = 6.67 \times 10^{-11}\,\frac{\text{N}\cdot\text{m}^2}{\text{kg}^2}

Although gravitational attraction operates between all matter over arbitrarily large distances, its intensity drops off with the square of the separation distance (∝1/r2\propto 1/r^2). Because GG is very small, gravitational forces between everyday objects are imperceptible and require massive bodies or sensitive instrumentation to detect.

Celestial Mechanics and Gravitational Phenomena

Gravitational theory enables the calculation of surface gravities, celestial masses, and laboratory interactions across diverse physical scales.

The Moon has an average radius of rM=1738 km=1.738×106 mr_M = 1738\,\text{km} = 1.738 \times 10^6\,\text{m} and a total mass of mM=7.35×1022 kgm_M = 7.35 \times 10^{22}\,\text{kg}. The surface gravitational acceleration gMg_M is derived by equating gravitational force to weight:

m⋅gM=G⋅m⋅mMrM2  ⟹  gM=G⋅mMrM2m \cdot g_M = G \cdot \frac{m \cdot m_M}{r_M^2} \implies g_M = \frac{G \cdot m_M}{r_M^2}

gM=6.67×10−11 N⋅m2/kg2⋅(7.35×1022 kg)(1.738×106 m)2=1.62 m/s2g_M = \frac{6.67 \times 10^{-11}\,\text{N}\cdot\text{m}^2/\text{kg}^2 \cdot (7.35 \times 10^{22}\,\text{kg})}{(1.738 \times 10^6\,\text{m})^2} = 1.62\,\text{m/s}^2

When astronaut Neil Armstrong jumped vertically on Earth (gE=9.81 m/s2g_E = 9.81\,\text{m/s}^2), he reached an apex of hE=0.60 mh_E = 0.60\,\text{m}. Assuming identical takeoff velocity and mass, the initial kinetic energy is identical in both environments. At the peak of the jump, kinetic energy has converted entirely into potential energy (Epot=m⋅g⋅hE_{\text{pot}} = m \cdot g \cdot h):

m⋅gE⋅hE=m⋅gM⋅hM  ⟹  hM=hE⋅gEgM=0.60 m⋅9.81 m/s21.62 m/s2=3.6 mm \cdot g_E \cdot h_E = m \cdot g_M \cdot h_M \implies h_M = \frac{h_E \cdot g_E}{g_M} = \frac{0.60\,\text{m} \cdot 9.81\,\text{m/s}^2}{1.62\,\text{m/s}^2} = 3.6\,\text{m}

Portrait sketch of Henry Cavendish

In 1787, Henry Cavendish first measured the weak gravitational attraction between laboratory masses using a torsion balance. Lead has a density of ρ=11.34×103 kg/m3\rho = 11.34 \times 10^3\,\text{kg/m}^3.

To find the diameter dd of two identical lead spheres that, when touching (r=dr = d), attract each other with a target force of FG=0.010 NF_G = 0.010\,\text{N}, the mass of each sphere is expressed as m=ρ⋅V=ρ⋅4π3(d2)3m = \rho \cdot V = \rho \cdot \frac{4\pi}{3}\left(\frac{d}{2}\right)^3:

FG=G⋅m2d2=G⋅ρ2⋅[4π3(d2)3]2d2=G⋅π29⋅ρ2⋅d4F_G = G \cdot \frac{m^2}{d^2} = G \cdot \frac{\rho^2 \cdot \left[\frac{4\pi}{3}\left(\frac{d}{2}\right)^3\right]^2}{d^2} = G \cdot \frac{\pi^2}{9} \cdot \rho^2 \cdot d^4

Solving for the diameter:

d=36⋅FGG⋅π2⋅ρ24=36⋅0.010 N6.67×10−11 N⋅m2/kg2⋅π2⋅(11.34×103 kg/m3)24=1.4 md = \sqrt[4]{\frac{36 \cdot F_G}{G \cdot \pi^2 \cdot \rho^2}} = \sqrt[4]{\frac{36 \cdot 0.010\,\text{N}}{6.67 \times 10^{-11}\,\text{N}\cdot\text{m}^2/\text{kg}^2 \cdot \pi^2 \cdot (11.34 \times 10^3\,\text{kg/m}^3)^2}} = 1.4\,\text{m}

The lead spheres used in Cavendish's actual experiment had diameters of d1=5.1 cm=0.051 md_1 = 5.1\,\text{cm} = 0.051\,\text{m} and d2=20 cm=0.20 md_2 = 20\,\text{cm} = 0.20\,\text{m}. When placed in contact, their center-to-center distance is r=d1+d22r = \frac{d_1 + d_2}{2}. The mutual gravitational force between them is:

FG=G⋅m1⋅m2(d1+d22)2=G⋅π29⋅d13⋅d23(d1+d2)2⋅ρ2F_G = G \cdot \frac{m_1 \cdot m_2}{\left(\frac{d_1 + d_2}{2}\right)^2} = G \cdot \frac{\pi^2}{9} \cdot \frac{d_1^3 \cdot d_2^3}{(d_1 + d_2)^2} \cdot \rho^2

FG=6.67×10−11⋅π29⋅(0.051 m)3⋅(0.20 m)3(0.051 m+0.20 m)2⋅(11.34×103 kg/m3)2=1.6×10−7 NF_G = 6.67 \times 10^{-11} \cdot \frac{\pi^2}{9} \cdot \frac{(0.051\,\text{m})^3 \cdot (0.20\,\text{m})^3}{(0.051\,\text{m} + 0.20\,\text{m})^2} \cdot (11.34 \times 10^3\,\text{kg/m}^3)^2 = 1.6 \times 10^{-7}\,\text{N}

Certain neutron stars rotate with a frequency of f=1.0 Hzf = 1.0\,\text{Hz}. For material on the equatorial surface to remain gravitationally bound rather than being flung off by centrifugal effects, the inward gravitational attraction must equal or exceed the required centripetal force (FG≥FzF_G \ge F_z):

G⋅m⋅Mr2=m⋅v2r=m⋅(2π⋅r⋅f)2rG \cdot \frac{m \cdot M}{r^2} = m \cdot \frac{v^2}{r} = m \cdot \frac{(2\pi \cdot r \cdot f)^2}{r}

Expressing the stellar mass as M=43π⋅r3⋅ρM = \frac{4}{3}\pi \cdot r^3 \cdot \rho:

G⋅43π⋅r3⋅ρr2=4π2⋅r⋅f2G \cdot \frac{\frac{4}{3}\pi \cdot r^3 \cdot \rho}{r^2} = 4\pi^2 \cdot r \cdot f^2

The radius rr cancels completely out of the equation:

G⋅43π⋅ρ=4π2⋅f2  ⟹  ρ=3π⋅f2GG \cdot \frac{4}{3}\pi \cdot \rho = 4\pi^2 \cdot f^2 \implies \rho = \frac{3\pi \cdot f^2}{G}

ρ=3π⋅(1.0 Hz)26.67×10−11 N⋅m2/kg2=1.4×1011 kg/m3\rho = \frac{3\pi \cdot (1.0\,\text{Hz})^2}{6.67 \times 10^{-11}\,\text{N}\cdot\text{m}^2/\text{kg}^2} = 1.4 \times 10^{11}\,\text{kg/m}^3

Asteroid 243 Ida and its small moon Dactyl imaged by the Galileo spacecraft

In May 1999, reports documented that the minor planet Ida (dimensions: length l=56×103 ml = 56 \times 10^3\,\text{m}, width b=24×103 mb = 24 \times 10^3\,\text{m}, height h=21×103 mh = 21 \times 10^3\,\text{m}) is orbited by a small spherical moon named Dactyl (diameter 1.5 km1.5\,\text{km}) at an orbital radius of r=100 km=100×103 mr = 100\,\text{km} = 100 \times 10^3\,\text{m} with an orbital period of T=22 h=79200 sT = 22\,\text{h} = 79200\,\text{s}.

Equating gravitational force to centripetal force yields the mass of Ida:

G⋅m⋅Mr2=m⋅v2r=m⋅(2π⋅rT)2r  ⟹  M=4π2⋅r3G⋅T2G \cdot \frac{m \cdot M}{r^2} = m \cdot \frac{v^2}{r} = m \cdot \frac{\left(\frac{2\pi \cdot r}{T}\right)^2}{r} \implies M = \frac{4\pi^2 \cdot r^3}{G \cdot T^2}

Approximating Ida's volume as a rectangular cuboid (V=l⋅b⋅hV = l \cdot b \cdot h), its density is calculated as:

ρ=MV=4π2⋅r3G⋅T2⋅l⋅b⋅h\rho = \frac{M}{V} = \frac{4\pi^2 \cdot r^3}{G \cdot T^2 \cdot l \cdot b \cdot h}

ρ=4π2⋅(100×103 m)36.67×10−11 N⋅m2/kg2⋅(79200 s)2⋅(56×103 m)⋅(24×103 m)⋅(21×103 m)=3.3×103 kg/m3\rho = \frac{4\pi^2 \cdot (100 \times 10^3\,\text{m})^3}{6.67 \times 10^{-11}\,\text{N}\cdot\text{m}^2/\text{kg}^2 \cdot (79200\,\text{s})^2 \cdot (56 \times 10^3\,\text{m}) \cdot (24 \times 10^3\,\text{m}) \cdot (21 \times 10^3\,\text{m})} = 3.3 \times 10^3\,\text{kg/m}^3

This confirms that Ida has a density similar to that of the Moon, consisting predominantly of silicate rock.

At the center of the Milky Way, a supermassive black hole is orbited by a star with an observed period of T=15 a=15×365.25×86400 s=4.7×108 sT = 15\,\text{a} = 15 \times 365.25 \times 86400\,\text{s} = 4.7 \times 10^8\,\text{s} at a mean orbital radius of r=0.013 LJ=1.2×1014 mr = 0.013\,\text{LJ} = 1.2 \times 10^{14}\,\text{m}.

Balancing gravitational force and centripetal force:

G⋅m⋅Mr2=m⋅v2r  ⟹  G⋅Mr=(2π⋅rT)2G \cdot \frac{m \cdot M}{r^2} = m \cdot \frac{v^2}{r} \implies \frac{G \cdot M}{r} = \left(\frac{2\pi \cdot r}{T}\right)^2

M=4π2⋅r3T2⋅G=4π2⋅(1.2×1014 m)3(4.7×108 s)2⋅6.67×10−11 N⋅m2/kg2=4.9×1036 kgM = \frac{4\pi^2 \cdot r^3}{T^2 \cdot G} = \frac{4\pi^2 \cdot (1.2 \times 10^{14}\,\text{m})^3}{(4.7 \times 10^8\,\text{s})^2 \cdot 6.67 \times 10^{-11}\,\text{N}\cdot\text{m}^2/\text{kg}^2} = 4.9 \times 10^{36}\,\text{kg}

Comparing this mass to the solar mass (MS=1.9884×1030 kgM_S = 1.9884 \times 10^{30}\,\text{kg}):

MMS=4.9×1036 kg1.9884×1030 kg=2.5×106\frac{M}{M_S} = \frac{4.9 \times 10^{36}\,\text{kg}}{1.9884 \times 10^{30}\,\text{kg}} = 2.5 \times 10^6

The supermassive black hole contains approximately 2.5×1062.5 \times 10^6 solar masses.