In an axisymmetric cylindrical coordinate system, the steady Collisionless Boltzmann equation is
Let and . Multiply by and integrate over velocities. The spatial terms give . Integration by parts in supplies , and the derivative supplies . The force term vanishes because does not depend on . This assumes sufficient decay of the galactic distribution function to discard all velocity boundary terms. The radial Jeans equation is consequently
Reflection symmetry about makes even and odd. At the midplane the latter vanishes, and . Multiplying by proves
All quadratic velocities here are population moments, not individual stellar velocities. In particular, symmetry does not make the midplane tilt contribution to asymmetric drift disappear: an odd mixed moment can have a nonzero derivative at zero.
For application to the Milky Way, assume negligible mean radial and vertical flows, and define the circular speed by . Write and . Rearranging gives the exact stress-support relation
Thus the stellar asymmetric drift is the right-hand side divided by . For , this denominator is approximately . If and locally, the bracket becomes , making the density and dispersion-gradient contributions explicit.
A radially declining random-motion stress supplies some of the force needed to support a galactic disk, so its stars need less ordered azimuthal motion than a population on circular orbits. Hotter, often older, stellar populations therefore generally rotate more slowly than cold populations in the same gravitational potential; this is the observed stellar asymmetric drift. It is a stress-support effect, not a frictional slowing of each orbit. For a velocity ellipsoid aligned approximately with spherical coordinates, , and the midplane tilt term contributes rather than zero.
Local measurements of mean rotation and all components of the velocity dispersion, combined with tracer-density and dispersion gradients, can therefore constrain and the radial Newtonian gravitational field. A hotter sample's mean rotation must not be identified directly with the circular speed. The tracer number density need not be the gravitating mass density. Nonaxisymmetric streaming, a nonstationary disk, poorly measured radial gradients, or a neglected tilt term can bias this inference; these are assumptions of the derivation rather than extra free corrections to the circular force.
Divide the velocity by the speed calculated above. The unit tangent vector is
Its derivative is
Because , the principal normal vector has the direction of , giving
This is a unit vector, and for every real . The curvature of a space curve is , so the principal normal vector is defined everywhere on this curve. The same perpendicularity follows from constant tangent component forces perpendicular principal normals: is constant, and differentiating this identity forces the normal component along to vanish.
Differentiate the parametrized curve to obtain . Its speed obeys
Since , the arc length is
At infinity, the meteorite's incoming asymptote has speed and perpendicular distance from the Earth’s centre. The magnitude of its specific angular momentum is the magnitude of . On the asymptote, this is , so conservation of angular momentum gives . Its mechanical energy at infinity is . Comparing with part (ii) gives
For this Kepler orbit is a hyperbola. The smallest occurs at periapsis, where . Thus
The middle equality follows by substituting . A trajectory entering the Earth satisfies , equivalently
Therefore the sufficient collision condition is
This is gravitational focusing: Newtonian gravity increases the collision radius measured at infinity. Equality is the ideal grazing trajectory, which touches the surface; a strict inequality penetrates it. For the incoming motion is radial and collides, with as the limiting value of the above expression.
Let be the uniform mass density. In spherical coordinates with the chosen diameter as polar axis, the squared perpendicular distance to that axis is . The moment of inertia is therefore
This derives the moment of inertia of a uniform solid sphere about any central axis by spherical symmetry.
Put . The center of mass follows the circle of radius , with position . Let and . Its velocity is , and the contact point is displaced from the center of mass by . With positive spin angular velocity counterclockwise, its contact-point velocity is
The cylinder is fixed, so rolling without slipping requires this velocity to vanish. Consequently . If instead the spin angle is positive clockwise, opposite to increasing , then the printed positive-coefficient convention is
The spin and orbital rotation have opposite directions when both are measured about the same oriented axis; the angular speed is . This specifies the sign convention in the printed expression.
The instantaneous contact point is at rest, so the normal reaction and static friction do no work. Conservation of energy equates the loss of gravitational potential energy from release to the total kinetic energy:
On the descent from , decreases, so
Integrating from release to the bottom gives
The integral is finite: at the release point its denominator is asymptotic to . Contact is maintained, since the inward radial equation of motion gives the positive normal reaction throughout this range.
For frictionless sliding, the normal reaction passes through the sphere's centre and produces no torque; Newtonian gravity also exerts no central torque. With initial zero spin, the sphere remains nonrotating. All the lost gravitational potential energy becomes translational kinetic energy, so
A hollow spherical shell of the same mass and radius places its material farther from the centre, giving a larger moment of inertia. Under rolling without slipping, a given centre speed therefore requires more rotational kinetic energy; the same gravitational drop yields a smaller speed. Replacing the solid sphere by a hollow shell increases and leaves unchanged. The sliding case remains nonrotating and has the same center of mass path and gravitational drop. This is the general comparison for rolling inside a fixed circular cylinder.
In polar coordinates, the curve is the logarithmic spiral , . It starts at and winds counterclockwise outwards through one and a half turns. It crosses the horizontal axis at
The following original sketch uses representative positive values; changing rescales it and changing changes its rate of expansion.
Figure 1.
Logarithmic spiral over one and a half counterclockwise turns
. Details.
Differentiate the coordinates:
Hence . Integrating the arc length element gives, for ,
For a negative endpoint parameter the length between and zero uses the absolute value of the final factor.
A second differentiation gives
The numerator is positive, so the absolute value in the PDF's curvature of a plane curve formula leaves it unchanged. Consequently
The speed is , and therefore . Thus the requested elapsed time is
This travel time at reciprocal-curvature speed is independent of because the growing arc length per unit angle is exactly matched by the growing speed.
In the rotating-hoop bead problem, use a signed angle in the hoop plane from the downward vertical. The gravitational potential energy is , the rotating-frame speed is , and the distance from the rotation axis is . The centrifugal potential gives
Gravity is stationary in these rotating coordinates, and the frictionless constraint does no work. Projecting the rotating equation along the hoop, or differentiating its effective potential in the one-dimensional Euler-Lagrange equation, gives
If , the two off-axis positions are
Their unsigned angle to the downward vertical is the same; they lie on opposite sides of the hoop. The second derivative of the effective potential there is . Equivalently, writing gives to first order
The same holds at . Thus both equilibria are stable, with small-oscillation angular frequency .
For , the Lorentz transformation with coincident origins is
For two simultaneous events in , , so . Distinct spatial positions therefore generally give different times in the boosted frame: this is relativity of simultaneity.
Choose the emission events at and . After emission, the photon world lines are and . Applying the Lorentz transformation to either gives
Thus the two photons move at speed in on parallel lines whose intercepts differ by . Measuring their positions at the same after both emissions gives the photon separation under a collinear Lorentz boost:
This separation is constant. It is not the contracted distance between the stationary sources: the two emission events are not simultaneous in , and one photon has already moved when the other is emitted.
Set initially and choose its positive sense to agree with the initial tangential motion. The inward radial component is , while the tangential speed is . Hence
Substitution in the Binet equation gives , or . This is an orbit equation for combined inverse-square and inverse-cube attraction. Solving the constant-coefficient equation and applying both initial conditions gives
On , the sine is nonnegative, so stays finite and positive. At , , hence and the particle is back at its original spatial position after one revolution. But , so its radial velocity is now , rather than the original . The return is not a periodic return in position and velocity.
After that revolution, the first zero of is at , where . The physical branch therefore has as ; it must not be continued into negative . Since and has a simple zero, reaching infinity takes infinite time. As a check, the specific potential energy is and the conserved energy per unit mass is , giving escape speed at infinity, also obtained from .
Differentiation gives . The squared speed is , so the arc length is
The speed never vanishes, so this is a regular curve throughout its parameter range.
Use the arch geometry of two hinged rods sliding on a floor, with and on the other side of . As long as remains in contact with the floor,
The centres of mass of the two rods are and . Their squared speeds are respectively and . Each uniform rod has moment of inertia about its midpoint, and their angular velocities are and . Applying the kinetic energy decomposition about the center of mass to each rod gives
The gravitational potential energy is . The fixed pivot does no work; neither does the frictionless floor reaction while slides horizontally. Conservation of energy from rest therefore gives
On the falling branch is the negative square root. Just before the hinge reaches the floor, , and its speed is . Thus, under the maintained-contact interpretation,
Figure 1.
Geometry and centers of mass of two hinged rods sliding on a floor
. Details.
There is a contact hypothesis behind this calculation. A floor can push upward but cannot pull downward. For and , the required normal reaction at , found from the second rod's linear and angular equations, is
For the usual acute arch, this stays nonnegative throughout exactly when . Indeed, for , its nonnegativity is , whose minimum is at . For larger initial angles, the endpoint lifts before the hinge hits the floor, and the one-coordinate formula no longer governs the whole motion. The displayed impact speed is the intended floor-constrained answer; for unrestricted initial angles, maintained contact is an additional assumption. Subsequent free motion and possible reimpact would require a different model and contact/impact conditions.
The PDF gives initial position and initial rotating-frame velocity zero. Use those values; the TeX erroneously places dots over the position coordinates.
For rotation about the third axis, the supplied acceleration relation and Newton's second law give the equations of motion in a rotating frame
Thus . Set to combine the planar equations:
The two characteristic roots are and . From and ,
The triple-angle formulas therefore give the particularly simple astroid motion of a harmonic oscillator in a rotating frame
The orbit satisfies ; the zero speed at each cusp is consistent with the smooth time-parametrized motion.
Its rotating-frame speed is
Hence, for , the first maximum occurs at , and
The maximum repeats periodically. For a signed negative angular frequency, use and ; at the particle remains at rest.
Figure 1.
Astroid orbit and the first maximum speed in the rotating frame
. Details.
The total relativistic energy and momentum in the original frame are
A Lorentz boost of signed velocity makes the total spatial momentum . Setting it to zero gives the centre-of-momentum frame velocity
It is a positive-weighted average of subluminal velocities, so .
In the rest mass of a coalescing relativistic collision, four-momentum conservation gives provided no energy or momentum escapes in additional radiation. For the single outgoing particle, , so its speed in the original frame is . Its rest mass is the invariant norm of the total four-momentum:
Expanding and using gives
The positive square root is taken for . The factor multiplying is the relative Lorentz factor, at least one. Therefore , with equality when ; otherwise the initial relative kinetic energy becomes internal rest energy.
Use zero-relative-velocity mass shedding: the released sand initially shares the balloon's vertical velocity, rather than being actively ejected. Let and let a positive mass be shed in a short interval . The momentum of the balloon and that newly released sand changes by
The leading external impulse is . Newton's second law applied to this combined material system therefore gives
Writing for the remaining balloon alone would omit the outgoing momentum flux. A nonzero relative ejection velocity would instead produce a rocket equation thrust term.
Initial mechanical equilibrium gives . Write , so up to depletion. Then
Integrating from rest,
Since and , the final upward speed is
The acceleration is nonnegative throughout release, consistent with the upward motion.
A four-vector is a collection of four components that transforms by the same Lorentz transformation as . Use the metric signature . The Lorentzian inner product is
A nonzero four-vector is a timelike vector, null vector or spacelike vector according as is negative, zero or positive. Reversing the metric signature reverses the signs used to name these three classes, without changing their geometric meaning.
For a frame moving at speed along the positive axis, write and use the Lorentz factor . The component Lorentz transformation is
Expanding the first two squares gives
The other two components are unchanged, so . Consequently a timelike vector remains timelike under a Lorentz transformation.
For a nonzero null vector, and . Set and . Then and . If the zero four-vector is included among null vectors, take and any unit vector.
For two future-pointing null vectors, write and with . The sum of future-pointing null vectors satisfies
because the ordinary inner product of two unit vectors is at most one. The sum is null if their spatial directions coincide, and timelike otherwise. Its positive time component makes the sum nonzero and future-pointing.
First apply momentum conservation to the rocket and the infinitesimal amount of expelled gas. During , the exhaust mass is and its inertial velocity, to first order, is . With no dust, the momentum balance is
It gives the rocket equation , while .
For a dust-collecting rocket, the incoming dust mass is with zero initial momentum. The same balance has rocket mass after the interval, so
Consequently
The second term in the acceleration equation accounts for the momentum required to bring collected dust up to rocket speed.
Set and , so . As long as the engines operate and , the velocity increases and time can be eliminated:
Integrating from , , gives
Thus the required mass–velocity relation is
If , this formula would give as . In this interval, . A real rocket retains positive body mass even after its fuel is gone, and collected dust only increases that residual mass. It therefore cannot follow this branch all the way to zero mass: its fuel must run out while . Indeed the hypothetical time integral is integrable at when , so even the formal zero-mass endpoint would occur at finite time.
If , the same relation instead gives as . The mass first decreases, reaches its minimum at , and then increases because dust arrives faster than fuel is expelled. The time integral now diverges at : under indefinitely maintained thrust, approaches asymptotically while stored dust grows without bound. At that limiting speed, thrust balances the dust-loading term , and . This is the physical reason the total rocket mass no longer signals fuel exhaustion. A rocket with a finite initial fuel supply still runs out at finite time and before reaching ; increasing total mass does not create new fuel.
Write for the genuine surface density of a disk, and for its three-dimensional mass density. The delta function in the printed surface density formula belongs to , not to . Introduce a reference length to make logarithms dimensionless. The intended model is the infinite, self-gravitating, scale-free Mestel disk, with no extra source or imposed gravitational field.
The circular speed satisfies . Thus the flat galaxy rotation curve gives
A reflection-symmetric harmonic function with this midplane boundary value is
Verify this continuation using the Poisson equation for Newtonian gravity. For , set ; then
Hence off the astrophysical disk, and reflection gives the lower-half-space solution. Across the astrophysical disk the derivative jumps by . The distributional Poisson equation for Newtonian gravity therefore gives
At the origin the enclosed disk mass tends to zero linearly with radius, so there is no additional central point mass. This verifies the Mestel disk potential-density pair. The astrophysical disk has infinite total mass and a logarithmic Newtonian gravitational potential, so it is not an isolated finite-mass model with Newtonian gravitational potential zero at infinity. The usual scale-free boundary condition is important: the midplane rotation curve alone would also permit an added term , representing an extra uniform sheet without changing the radial circular force. That contribution is excluded in the intended Mestel disk model.
Use a mass-weighted planar galactic distribution function, so is the stellar mass in a small planar phase space element. A number-weighted function instead needs the stellar mass factor when computing . Assume a steady collisionless stellar system and isotropy in the two in-plane velocity components. Put and write . The stationary Collisionless Boltzmann equation becomes
Since this holds in every velocity direction, . In coordinates with , this is exactly . Therefore the planar isotropic distribution is , where is the specific orbital energy. Stationarity is essential; instantaneous isotropy alone would not imply this result.
Integrating over the two-dimensional velocity plane gives the planar isotropic distribution inversion:
On the astrophysical disk , so . Differentiate the integral with respect to its lower limit:
The boundary value as verifies the integrated equation as well as its derivative. Choosing the implicit length unit recovers the printed normalization. Changing the additive energy zero changes this prefactor accordingly.
At a fixed radius, the normalized velocity density is
It is a product of centered Gaussian distributions. Differentiating the supplied Gaussian integral with respect to its coefficient gives the second moments, and odd moments vanish. Thus the in-plane velocity dispersions are
An exactly planar astrophysical disk has and . The nonzero circular speed is a property of the force field, not a statement that this hot stellar distribution has net rotation.
Reversing every retrograde star folds the azimuthal Gaussian to a half-normal distribution. Equivalently the new steady galactic distribution function is , because and are integrals of the motion. Its density and even velocity moments are unchanged. Its streaming velocity is
This maximally prograde stellar distribution still has radial motion and a spread of azimuthal speeds; it does not place every star on a circular orbit. In particular while . The mean speed is smaller than the root-mean-square speed , which explains why it is not the circular speed.
Synchronous radius and its limits. Kepler's third law for a circumplanetary orbit gives
Here is the gravitational constant. A planet-synchronous orbit repeats after one spin period. A planet-stationary orbit additionally requires a circular orbit, zero orbital inclination to the equator and prograde motion. Then an antenna fixed on the ground points continuously at the same satellite; merely matching the orbital period does not give this property.
Write the planetary radius as . Requiring the semi-major axis to exceed gives . The stellar tidal force restricts the circumplanetary orbit to the Hill sphere, of radius . Thus
These are the surface and Hill-radius constraints in the idealized spherical, small- model. Long-term prograde stability generally requires a radius appreciably inside the Hill sphere; for nonzero orbital eccentricity, the surface constraint applies to the periapsis, not just to the semi-major axis.
Collision time and the launch population. A phase-mixed isotropic swarm occupies a shell of radial scale . Its number density scales as , its relative speed as , and its geometric collision cross-section as . Consequently the total collision rate scales as . To give the numerical normalization used here, adopt an effective shell volume . At a fixed position the velocity directions are uniformly distributed in the tangent plane, so their mean relative speed is . For an unordered pair, the collision cross-section is , giving
Thus the mean interval between collisions is in this large-population kinetic model. The effective shell width fixes an order-one coefficient: small orbital eccentricity and random planes alone do not specify a unique radial probability density. Phase mixing, negligible gravitational focusing, , and uncorrelated encounters are implicit in this estimate. Exactly identical orbital periods with perfectly fixed phases do not themselves produce a memoryless collision process.
For nearly planet-stationary orbits with , the swarm volume is smaller by a factor of order , while the relative speed is smaller by the same factor, since vertical motion of scale dominates the eccentric motion. The two changes cancel in the rate : there is no parametric factor in the collision time within the same phase-mixed kinetic approximation. Numerical factors and phase correlations can differ. This is not an argument that bringing all satellites into one nearly circular plane makes their phases random.
With , an Inhomogeneous Poisson process has cumulative hazard function
The survival function of the first collision is . Setting the expected number of collisions to one gives
This is a characteristic first-event population, with probability of an earlier event. It is not a median: the median has an additional factor .
Which population collides next? Immediately after the first disruption let and . The latter follows from mass conservation for equal-density spherical pieces. Using the same unordered-pair counting and geometric collision cross-sections as above, define ; then
Here stand for the usual large-population approximations to . The factor two distinguishing identical and different species is essential. The probability that the next collision is fragment–fragment is
It is the most likely type if and . It has probability greater than one half if
For , a strongly fragment-dominated next event therefore requires ; the fragment–satellite comparison is more restrictive than the satellite–satellite comparison. If , the next event must be fragment–fragment, provided fragments remain.
Population equations and the normalization discrepancy. Every satellite–satellite event produces fragments and destroys two satellites. Every satellite–fragment event produces a net fragments and destroys one satellite; every fragment–fragment event destroys two fragments. Consistent collision counting therefore gives
For , and . The last two terms in are then twice those in the printed equation. This is a genuine factor-of-two inconsistency: equal-size fragment pairs have a collision cross-section smaller by , so their event rate must be if the satellite event rate is ; destroying both fragments necessarily gives the sink .
If the printed equation is taken as a prescribed approximate rate model instead, its implicit event rates are , and . Under precisely that mixed normalization, its corresponding satellite equation is
The printed source term also neglects the consumed fragment in a satellite–fragment event, a legitimate relative approximation. That approximation does not repair the pair-counting discrepancy.
The ensuing cascade. This is a two-size fragmentation cascade. First use the printed approximate model, dropping satellite–satellite events as requested. Put , , and take . Then
Integrating this linear differential equation, with and , gives
Initially the collisional cascade grows if , with approximate early exponential growth time when the satellite population is nearly fixed. The fragment population reaches its maximum at
Afterwards satellites are depleted and fragment–fragment losses dominate. For the continuum solution has , , with
When no satellites remain initially, directly. Small integer populations eventually invalidate these deterministic differential equations.
The consistently counted model has, to leading order in , exactly the same curve and peak, but both retained time derivatives are twice as large. Its growth time is , and . Keeping and the consumed fragment replaces by and in the curve by . This explicitly separates the physical collision bookkeeping from the printed normalization while giving the evolution under both conventions.
For equal in the two cases, their closest distances are conjugate roots expressed with the same . Consequently
At either central-force radial turning point, the velocity is tangential and conservation of angular momentum gives . Thus
These attractive and repulsive inverse-square closest approaches pair a shorter, faster attractive passage with a longer, slower repulsive one. The speed-product conclusion assumes ; at , the attractive orbit hits the singular centre with unbounded speed, so the formal product with the repulsive zero speed is not defined.
For , the same asymptotic data give and . At the central-force radial turning point,
The positive root and the tangential speed, again with , are
For , and . The repulsive inverse-square force bends the trajectory away from the origin, as in the right-hand sketch. Its polar equation of a repulsive inverse-square orbit is , with and measured from closest approach.
For , the regular radial turning point is and its speed is zero; the particle reverses direction there. The formula still gives zero, while its rationalized form containing should be interpreted only by a limit.