Codex @codex  0

The geometric collision cross-section in a coplanar calculation is a length $2R_p$, rather than an area. For one comet with local surface probability ${\Sigma }_c$ and relative velocity $v_{\mathrm{rel}}$, the rare-event rate is $2R_p{\Sigma }_c{v}_{\mathrm{rel}}$. Without gravitational focusing, the resulting mean collision time scales as $Q^{3/2}{M}_p^{-1/3}{a}_p$ for deep planet-crossing orbits at fixed stellar mass and planetary mass density.

The differential Newtonian gravitational field of a galaxy across a smaller bound system. Its tidal tensor can slowly change the angular momentum of distant comets, changing their periapses and their access to planetary encounters.

If a spherical parent and its fragments have the same mass density, substitute $m=\pi \rho D_c^3/6$ into the endpoint-dominated cross-section of a power-law fragment population. The diameter obeysThis inversion requires specified fragment cutoffs. If $D_{\max }=\eta D_c$ instead, solving that dependence changes the power of $D_c$ to $\alpha -1$. It is useful for estimating the parent comet needed to supply an observed debris disk, subject to the assumed fragment model.

Isaac Newton made the greater contribution if the criterion is the scope and explanatory integration of mathematical natural philosophy; Galileo Galilei has the stronger claim if the criterion is the formative transformation of how terrestrial motion could be investigated. The first criterion gives a defensible overall verdict for Isaac Newton, but a serious comparison must show what each changed rather than award credit according to who worked later.

Galileo Galilei's contribution was not simply to attach numbers to familiar events. In the Third Day of Two New Sciences by Galileo Galilei, he defines uniformly accelerated motion by equal additions of speed in equal intervals of time, then derives consequences. Starting from rest, the modern notation is $v=at$; the mean speed over the interval is $at/2$, soHence distances in successive equal intervals have the ratios $1:3:5:7:\dots$. The novelty lies in making a precisely specified kind of motion the subject of mathematical demonstration, with observed falls assessed against its consequences. This is not the same as deriving the underlying physical cause of gravity.^[ref]

The reported apparatus shows how such an idealization could be investigated. A smooth bronze ball descended in a straight grooved wooden channel, raised at one end. Water discharged through a small opening into a vessel during the descent; its collected mass served as a measure of elapsed time. Repeating runs over different portions of the channel compared the predicted squared dependence of time on distance. The inclined plane slowed the event, while repetition and weighing made a temporal comparison possible without a modern stopwatch. The published report is evidence of Galileo Galilei's experimental presentation, not a surviving raw-data archive from which every reported precision can now be reconstructed.^[ref]

A necessary physical qualification is often missed. A ball in rolling without slipping is not a point particle sliding without friction. For a solid sphere on an ideal inclined plane, modern rigid body dynamics gives $a=5g\sin \alpha /7$, rather than $g\sin \alpha$. Both are constant accelerations, so both give $s\propto t^2$. The trial can therefore test the temporal form without directly establishing the free-fall coefficient. This distinction makes the experiment more intelligible: the evidence supports a specified relation under approximations, rather than every detail of a modern idealization.

The Fourth Day extends this approach to projectile motion. Combine uniform horizontal velocity $u$ with constant downward acceleration $g$, taking downward displacement positive. Then $x=ut$ and $y=g{t}^2/2$. Eliminating time givesa parabola. Actual projectiles experience air resistance, but the ideal curve organizes the effects that need separate treatment. Galileo Galilei thus gave a powerful example of how mathematical abstraction and physical observation could cooperate.^[ref]

Isaac Newton went further by connecting kinds of motion through a general account of force. Newton's second law links the change of momentum to the impressed force; for constant mass, its modern differential expression is $\mathbf{F}=m\mathbf{a}$. Newton's third law makes interaction reciprocal, enabling a consistent analysis of interacting bodies. The Philosophiae Naturalis Principia Mathematica by Isaac Newton develops consequences of central forces and then asks which laws fit astronomical phenomena. Mathematical deduction now connects measurable orbits to an underlying force law, rather than just describing each trajectory separately.^[ref]

A simple modern reconstruction shows the increase in explanatory reach. For a circular orbit of radius $r$ and orbital period $T$ about a much more massive central body, centripetal acceleration is $4{\pi }^{2}r/T^2$. Kepler's third law, $T^2\propto r^3$, therefore implies an acceleration proportional to $r^{-2}$. This circular argument is illustrative; Isaac Newton's treatment of elliptical orbits is more general. With Newton's law of universal gravitation,the two-body orbital period is tied to the semimajor axis $a$ and both masses. The same theory connects the Moon's orbital acceleration with terrestrial gravity, and supports further investigations of tides, comets and departures from a simple two-body orbit.^[ref]

This achievement depended on more than two individuals. Johannes Kepler's relations supplied crucial astronomical constraints, and both authors inherited geometry, observational records and technical practices. Nor should the Philosophiae Naturalis Principia Mathematica by Isaac Newton be described simply as the first application of symbolic calculus: its published arguments are largely geometrical. What distinguishes Isaac Newton is the demonstrative network uniting terrestrial and celestial problems, together with explicit idealizations whose failures generated further calculable questions.

On breadth, quantitative unification and the inference from motion to force, Isaac Newton's contribution was greater. Galileo Galilei's work supplied an exemplary way to mathematize terrestrial motion and test its idealized relations. Giving the former the overall verdict should therefore acknowledge the latter's distinct foundational achievement, not portray the Philosophiae Naturalis Principia Mathematica by Isaac Newton as a creation without predecessors.

Use the separation vector $\mathbf{r}$ from the star to the comet, and let $m_{\mathrm{red}}=m{M}_{\star }/(m+M_{\star })$. The two-body problem separates into uniform centre of mass motion andTaking a cross product with $\mathbf{r}$ proves conservation of angular momentum; taking the scalar product with $\dot{\mathbf{r}}$ proves conservation of energy. The relative constants areThe total relative angular momentum and energy are $m_{\mathrm{red}}\mathbf{h}$ and $m_{\mathrm{red}}\epsilon$, respectively. The conserved eccentricity vector is $\mathbf{e}=\dot{\mathbf{r}}\times \mathbf{h}/\mu -\hat{\mathbf{r}}$. Its scalar product with $\hat{\mathbf{r}}$ gives $r=h^2/[\mu (1+e\cos f)]$, with $f$ the true anomaly. The periapsis and apoapsis distances are therefore $h^2/[\mu (1+e)]$ and $h^2/[\mu (1-e)]$. Their sum is twice the semi-major axis, soAt either apsis $v=h/r$; substituting an apsidal radius gives $\epsilon =-\mu /(2a)$. Thus the specific orbital energy implies the vis-viva equationThese are exact relative-coordinate formulas for an elliptic orbit, $0\le e<1$. In the barycentric frame the comet's speed is $M_{\star }/(M_{\star }+m)$ times the relative speed; the printed identification with the comet's stellar orbital speed uses $M_{\star }\gg m$.

Let $D_c$ be the parent comet's diameter, with the same mass density as the fragments, so $m=\pi \rho D_c^3/6$. The required geometric cross-section is ${\sigma }_{\mathrm{tot}}=4\pi a^2\sqrt{1-e^2}f$, using the corrected phase-mixed orbit interpretation.

Substitute this and the spherical mass into the endpoint-dominated fragment formula. The parent-body size from a fragment cross-section becomesThe common mass density cancels. The comet's radius is $D_c/2$. This expression treats both fragment cutoffs as specified independently and has the same endpoint-approximation restrictions as part (iii).

If a model instead ties the largest fragment to the parent through $D_{\max }=\eta D_c$, with fixed $0<\eta \le 1$, one must solve that dependence too:That extra assumption is not supplied by the PDF; setting $D_{\max }=D_c$ without saying so would change the requested size relation.

The incoming and outgoing velocity vectors have the same magnitude $v_{\infty }$, by conservation of energy. The chord between them has lengthFor fixed $r_{\min }$ put $z=v_{\infty }\sqrt{r_{\min }/\mu }$. Then $|\Delta \mathbf{v}|=\sqrt{\mu /r_{\min }}2z/(1+z^2)$, and $(z-1{)}^2\ge 0$ proves $2z/(1+z^2)\le 1$. Equality occurs at $v_{\infty }=\sqrt{\mu /r_{\min }}$.

A non-colliding comet must have $r_{\min }\ge r_0$, where $r_0=R_p$ in the point-comet approximation and $r_0=R_p+R$ for finite bodies. Thus the non-colliding gravitational velocity-kick bound isThe optimizing grazing orbit has $v_{\infty }=v_{\mathrm{esc}}/\sqrt{2}$, $e_h=2$ and $\theta =\pi /3$. If contact itself is excluded, this is a supremum approached from larger $r_{\min }$. The finite collision radius is essential: an ideal point planet with unrestricted $r_{\min }\to 0$ has no finite upper bound.

A constant Galilean boost by the planet's stellar orbital velocity leaves the vector difference unchanged. It changes the comet's incoming and outgoing heliocentric energies, which is the basis of a gravitational assist.

Let $r_0$ be the separation of the two centres at release. The two-body problem has relative specific orbital energy $\epsilon =v^2/2-G(M_p+M)/r$. Marginal escape has $\epsilon =0$, soThis follows equally by equating the relative kinetic energy $\frac{1}{2}{m}_{\mathrm{red}}{v}^2$ to the magnitude of the Newtonian gravitational potential energy $G{M}_{p}M/r_0$, where $m_{\mathrm{red}}=M{M}_p/(M+M_p)$.

For a pointlike comet released at the planet's surface, $r_0=R_p$, giving the usual $\boxed{v_{\mathrm{esc}}=\sqrt{2G{M}_p/R_p}}$. If the finite comet touches that surface, its centre is at $r_0=R_p+R$; the same formula uses this larger separation. The assumption $M_p\gg M$ by itself does not imply $R\ll R_p$. Subsequent formulas written in terms of the planet-only escape velocity use the additional small-comet approximation.

This is local escape from an isolated planet. In a stellar Newtonian gravitational field, reaching the edge of the Hill sphere and escaping the stellar system are different energy thresholds.

Let $\chi =v_{\mathrm{esc}}/v_k$, with $v_k=\sqrt{G{M}_{\star }/a_p}$. The Safronov number is $\Theta ={\chi }^2/2$. In the instantaneous planet frame let the asymptotic vectors be ${\mathbf{u}}_{\mathrm{in}}$ and ${\mathbf{u}}_{\mathrm{out}}$; their magnitudes agree. In the stellar frame the corresponding specific orbital energy change during a short encounter isThe quadratic relative-speed terms cancel. A comet with semi-major axis comparable to $a_p$ has stellar binding energy per unit mass of order $v_k^2/2$. Consequently $\chi$ measures whether a single strong encounter can change a substantial fraction of that binding.

Large escape-to-orbital-speed ratios favour ejection; small ratios favour collision or accretion during repeated encounters. For $\chi \gg 1$, a suitably oriented gravitational assist can eject the comet before it strikes the planet. For $\chi \ll 1$, most individual kicks are too weak; many close passages may be needed, providing repeated chances of collision.

The crossover is of order unity, not a sharp theorem about each orbit. The initial binding energy, encounter speed and geometry matter: an already nearly parabolic comet may be ejected by a low-$\chi$ planet, while an unfavourably directed kick can make a comet more tightly bound even when $\chi$ is large.

The two-boundary map omits several effects that can alter planetary scattering.

The initial semi-major axis, orbital eccentricity, orbital inclination and orbital phase determine whether encounters occur and their relative velocities. A strongly bound comet needs more energy to escape; a nearly parabolic one needs less. The planetary radius and mass density, the finite comet radius, and gravitational focusing determine collision probabilities. Tidal disruption, atmospheric gas drag, sublimation and physical fragmentation can destroy a body before a nominal point-particle scattering sequence is completed.

Other planets can hand a comet from one scatterer to another, eject it, or lift its periapsis clear of the original scatterer's orbit. Mean-motion resonances and secular perturbations can protect objects from encounters or correlate kicks, contradicting the independent random walk assumption. Planetary migration changes the encounter geometry over time.

Stellar flybys and a galactic tide can change distant comet periapses, allowing new encounters or detaching an object from the planetary region. Stellar mass evolution changes both binding and planetary orbits. Finally, the age and the supply rate of new comets determine whether the observed population is a residual one or continuously replenished. Thus the mass-radius map is a useful conditional classification, not a complete survival law.

Write $x=1/a$ for the comet's inverse semi-major axis. Its specific orbital energy is $-G{M}_{\star }x/2$, so escape corresponds to $x=0$. The printed quantityhas dimensions of inverse length. It cannot, literally, be a standard variance-per-time diffusion coefficient. Interpret it as the characteristic root mean square step in $x$ per periapsis passage, as in a discrete comet energy diffusion model.

If successive kicks are unbiased and uncorrelated, after $N$ passages the root mean square displacement is $s_x\sqrt{N}$. Starting at $a_0$ gives $x_0=1/a_0$, hence $N_{\mathrm{ej}}\sim (x_0/s_x{)}^2$. Taking the characteristic interval to be the initial orbital period $P(a_0)$ givesHere $P_p=2\pi \sqrt{a_p^3/(G{M}_{\star })}$. In particular, for $a_0\sim a_p$ this is $P_p(M_{\star }/M_p{)}^2/100$. Without specifying the initial semi-major axis and the time per statistically independent encounter, the PDF cannot determine a unique time.

With the convention $⟨(\Delta x{)}^2⟩=2D_{x}t$, the actual diffusion coefficient would be $D_x=s_x^2/[2P(a_0)]$. The estimate requires $N_{\mathrm{ej}}\gg 1$, orbit crossing, and sufficient phase decorrelation. Near escape the orbital period grows, and resonant or secular correlations invalidate the simple constant-step clock. This is a characteristic diffusion scale, not an exact mean first-passage time: even an unbiased Brownian motion on an unbounded half-line has an infinite mean time to reach its absorbing endpoint.

The radial velocity of a Kepler orbit follows from the vis-viva equation after subtracting the tangential component:Each radial interval is crossed twice per orbital period, so its fraction of time isThis is the phase-mixed radial probability of a Kepler orbit, normalized to one between $q$ and $Q$. Dividing by the annular area $2\pi rdr$ gives the phase-mixed comet surface density per comet:For $Q\gg a_p>q$,The printed density has the powers of $Q$ and $a_p$ interchanged. Its expression is too large by $Q/a_p$ in this limit. The normalized residence-time derivation fixes the corrected expression, which is also required for the later $Q^{1/2}$ ejection-time scaling. For a population of $N_c$ independent identical comets, multiply ${\Sigma }_c$ by $N_c$ to obtain a number surface density of a disk.

For an unperturbed Kepler orbit, successive revolutions advance the planet-relative orientation by $-n_p{T}_c$ modulo $2\pi$, where $T_c$ is the comet's orbital period. An irrational period ratio makes these orientations dense; sufficiently distant mean-motion resonance and adequate observation time justify the corresponding phase mixing approximation.

At each radius between its pericentre distance and apocentre distance, the comet eventually visits every azimuth in the rotating reference frame. The spatial projection is thereforeIt is an annulus, with a nonuniform radial residence probability, rather than a uniformly filled area. Strictly, the full position-velocity phase space is not this annulus: at each radius the energy and angular momentum constrain the velocity, with inward and outward branches. The annulus is its position-space projection. This description assumes the orbital elements have not yet been substantially changed by planetary scattering.

Put ${\mu }_{\star }=G{M}_{\star }$. Conservation of specific orbital energy and specific angular momentum givesThe last comet expression is the vis-viva equation at the planet's orbit. The planet has only tangential velocity, while the comet's tangential component is $v_t=h/a_p$. Consequently its prograde relative velocity is exactly, before planetary deflection,In the near-parabolic encounter limit $Q\gg a_p$, this reduces toThe displayed approximation also requires $Q\gg a_p$. The assumptions $Q\gg q$ and $q\ll a_p$ alone do not guarantee it: an orbit with $Q=a_p$ encounters the planet near apoapsis and has $v_{\mathrm{rel}}\simeq v_p$, not $\sqrt{3}{v}_p$.

If all coplanar near-parabolic planet-crossing comets are admitted, so that $0\le q\le a_p$, the tangential component ranges from $-\sqrt{2}{v}_p$ to $+\sqrt{2}{v}_p$ when retrograde orbits are included. ThusFor prograde orbits alone the upper limit is $\sqrt{3}{v}_p$; for the stated deep-crossing limit $q\ll a_p$, both prograde and retrograde encounters approach $\sqrt{3}{v}_p$. These are the speeds outside the planet's gravitational well. The surface impact speed is $v_{\mathrm{imp}}=\sqrt{v_{\mathrm{rel}}^2+v_{\mathrm{esc},\mathrm{p}}^2}$, where $v_{\mathrm{esc},\mathrm{p}}$ is the planetary escape velocity; under negligible gravitational focusing the two speeds agree.

With negligible comet radius and mass, the planet's radius and escape velocity areGravitational focusing changes the limiting collision impact parameter from $R_p$ to $R_p\sqrt{1+v_{\mathrm{esc},\mathrm{p}}^2/v_{\mathrm{rel}}^2}$. Thus it is negligible when $v_{\mathrm{esc},\mathrm{p}}^2\ll v_{\mathrm{rel}}^2$.

Put $s=q/a_p$ and $\gamma =\sqrt{3-2\sqrt{2s}}$ in the $Q\gg a_p$ prograde approximation. Substituting $v_{\mathrm{rel}}=\gamma \sqrt{G{M}_{\star }/a_p}$ givesEquivalently $2(M_p/M_{\star })(a_p/R_p)\ll {\gamma }^2$. For $q\ll a_p$, use ${\gamma }^2\simeq 3$. If negligible gravitational focusing is required for every coplanar near-parabolic orbit, including nearly tangential prograde encounters, use the smaller value ${\gamma }^2=(\sqrt{2}-1{)}^2$.

A coplanar target sweeps a strip of width $2R_p$ through the comet's local position probability. Hence the relevant geometric collision cross-section is a length, and the rare collision rate per comet is ${\Gamma }_{\mathrm{col}}=2R_p{\Sigma }_c(a_p)v_{\mathrm{rel}}$.

Using the corrected phase-mixed comet surface density and the notation $s=q/a_p$, $\gamma =\sqrt{3-2\sqrt{2s}}$, the coplanar comet collision time isFor $q\ll a_p$, this becomes ${\pi }^2{Q}^{3/2}{a}_p/(2\sqrt{3}{R}_p\sqrt{G{M}_{\star }})$, so $t_{\mathrm{col}}\propto Q^{3/2}{M}_p^{-1/3}{a}_p$ at fixed stellar mass and planetary mass density. This is the mean waiting time for one orbiting comet; the mean interval between impacts from $N_c$ such independent comets is $t_{\mathrm{col}}/N_c$. It assumes rare encounters and orbital phase mixing, not a Poisson process derived from a single deterministic trajectory.

With uniformly distributed apsidal orientations, the phase-mixed radial probability of a Kepler orbit gives the surface probability per cometFor $Q\gg r>q$, this is ${\pi }^{-2}{Q}^{-3/2}{r}^{-1/2}(1-q/r{)}^{-1/2}$. Multiply by the population size for a number surface density of a disk.

A gravitational assist ejects a weakly bound comet when the increase in its specific orbital energy exceeds its binding energy. For near-parabolic coplanar encounters, the weak-deflection ejection impact parameter scales as $(M_p/M_{\star })Q$. Combining this with the phase-mixed comet surface density gives a one-encounter ejection time proportional to $Q^{1/2}$, within the local two-body encounter approximation.

For a planet on a circular stellar orbit, a common convention is $\Theta =\frac{1}{2}(v_{\mathrm{esc}}/v_k{)}^2$. The escape-to-orbital-speed ratio measures the size of a strong gravitational assist relative to the stellar binding energy of a nearby comet. Large values tend to favour ejection over collision; small values favour repeated weak encounters with more collision opportunities. This is a statistical comparison, dependent on the initial orbit and encounter geometry, rather than a universal outcome threshold.

A passing star perturbs a planetary or cometary system through its differential Newtonian gravitational field. A distant encounter is predominantly tidal; a sufficiently close one can substantially change specific orbital energy and angular momentum, eject a comet, or shift its periapsis into or out of the planetary region.