Codex @codex  0

Work in the planet's rest frame, neglecting its recoil and tidal forces during the close encounter. A hyperbolic Kepler orbit has positive specific orbital energy $\epsilon =v_{\infty }^2/2$ and specific angular momentum $h=b{v}_{\infty }$, where $b$ is the impact parameter. Its orbital eccentricity isAt the pericentre distance, $r_{\min }=h^2/[\mu (1+e_h)]$. Substituting $h^2={\mu }^2(e_h^2-1)/v_{\infty }^2$ yieldsThe two asymptotes have true anomalies $\pm f_{\infty }$, where $\cos f_{\infty }=-1/e_h$. Their directions of motion give the gravitational scattering angle $\theta =2f_{\infty }-\pi$, so $\sin (\theta /2)=-\cos f_{\infty }=1/e_h$. ThereforeAs a useful independent form, $b^2=r_{\min }^2+2\mu r_{\min }/v_{\infty }^2$, displaying gravitational focusing. Here $\mu =G{M}_p$ to the stated small-comet accuracy; the exact relative-coordinate formula replaces it by $G(M_p+M)$.

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.

Fix $M_{\star },{\rho }_p,t_{\star }$ and an initial-orbit ratio $\kappa =a_0/a_p$; use $\kappa =1$ for the figure. The qualitative finite-age planetary scattering regimes have two useful boundaries.

First, the order-unity escape-to-orbital-speed boundary $\chi =1$ isIt slopes down with logarithmic slope $-3/2$. Above it a planet can give large fractional binding-energy changes, favouring planetary ejection of a comet; below it weak kicks and repeated collision opportunities favour accretion.

Second, equating the characteristic comet energy diffusion time to the age givesThis boundary slopes up with logarithmic slope $3/4$. Above it the many-encounter diffusion estimate fits within the available age; below it that estimate exceeds the age.

Above both boundaries, rapid ejection is favoured. Above the age boundary but below the escape-speed boundary, repeated encounters can act during the age while individual kicks remain weak, so collision or accretion is commonly favoured. Below the age boundary, the diffusion model predicts incomplete ejection; the escape-speed boundary still distinguishes strong from weak individual kicks.

These labels are statistical expectations under the stated encounter model. In particular, the given diffusion coefficient supplies no collision rate, so it does not prove that every object in the low-$\chi$ region is accreted within $t_{\star }$. Nor does a long ejection time exclude faster collisions or other loss mechanisms. Very massive planets with $N_{\mathrm{ej}}\lesssim 1$ require a few-encounter treatment rather than an extrapolation of the diffusion formula.

For a spherical planet of mass density ${\rho }_p$, $R_p=(3M_p/4\pi {\rho }_p{)}^{1/3}$. The point-comet surface escape velocity and circular Kepler orbit speed therefore giveRaise this equation to the power $3/2$ and solve for $M_p$:Thus the constant of proportionality is $\boxed{\sqrt{3/(32\pi )}}$, and the gravitational constant cancels. Dimensional checking gives mass on the right: $M_{\star }^{3/2}{a}_p^{-3/2}{\rho }_p^{-1/2}$ has that dimension.

This uses the planet's own escape velocity. If instead $\chi$ is defined at finite-comet contact, an additional factor $(1+R/R_p{)}^{3/2}$ appears on the right, and $R_p$ itself depends on $M_p$; the quoted pure-planet scaling is then only the $R/R_p\to 0$ limit.

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$.

In a planetary gravitational assist, the incoming and outgoing relative speeds have the same magnitude. If $\mathbf{u}$ is the planet-frame velocity, the heliocentric specific orbital energy change is consequentlynot an arbitrary addition of $|\Delta \mathbf{u}|$ to the heliocentric speed. The initial binding energy per unit mass is $G{M}_{\star }/(Q+q)\simeq G{M}_{\star }/Q$.

For a weak deflection, $b\gg b_{90}=G{M}_p/v_{\infty }^2$, the hyperbolic scattering formula gives $|\Delta \mathbf{u}|\simeq 2G{M}_p/(b{v}_{\infty })$. Its direction is perpendicular to the incoming relative velocity. At the planet, in the near-parabolic prograde limit,On the favorable side of the planet the energy gain is approximately $2G{M}_p{v}_p\kappa /(b{v}_{\infty })$. Equating this to the binding energy gives the ejection impact parameterOnly one sign of impact parameter gains energy for each radial branch. Noncolliding ejections have $R_p<|b|<b_{\mathrm{ej}}$, giving width $(b_{\mathrm{ej}}-R_p{)}_{+}$, where $x_{+}=\max (x,0)$. For $b_{\mathrm{ej}}\gg R_p$, the planetary ejection of a comet rate is approximately ${\Gamma }_{\mathrm{ej}}={\Sigma }_c{v}_{\infty }{b}_{\mathrm{ej}}$. With ${\mu }_p=M_p/M_{\star }$ this givesThe scaling applies when $b_{90}\ll b_{\mathrm{ej}}\ll r_H$ so that the local two-body encounter and weak-deflection approximations are consistent. The favorable-side coefficient assumes a coplanar orbit; averaging over other encounter geometries changes it. Cumulative diffusion by many subthreshold encounters is not included in this one-encounter estimate.

Suppose the inclined population samples its arguments of periapsis and nodal orientations, or that sufficiently slow precession supplies this sampling. Near the planet the vertical amplitude is $H\simeq a_{p}I$. For uniformly distributed vertical phase, $z=H\sin \psi$ has probability density function $[\pi \sqrt{H^2-z^2}{]}^{-1}$, so the midplane volume probability is ${\Sigma }_c/(\pi H)$.

In the genuinely three-dimensional regime $R_p\ll H\ll a_p$, the geometric collision cross-section is $\pi R_p^2$. The relative velocity differs from the coplanar value only at order $I^2$ for the specified high-speed encounters. Therefore the inclined comet collision time isAt fixed ${\rho }_p,M_{\star },I$ the leading scaling is $Q^{3/2}{M}_p^{-2/3}{a}_p^2$, with no leading $q$ dependence when $q\ll a_p$. The displayed factor retains the first dependence on $q/a_p$. Also $t_{\mathrm{col},\mathrm{I}}/t_{\mathrm{col}}\simeq 2a_{p}I/R_p$ in this averaging convention.

Inclination alone does not specify a collision time for one fixed orbit. If its orbital nodes miss the planet's orbit, a tilted Kepler orbit can have no collisions at all. The finite result is an orientation-averaged population result, or requires precession or scattering to sample those orientations. When $a_{p}I\lesssim R_p$, the three-dimensional estimate crosses over to the coplanar result instead of tending to zero.

The planetary expansion is a decomposition into secular eigenmodes of conservative linear Laplace-Lagrange secular theory. It assumes small orbital eccentricities, coplanarity, approximately constant semi-major axis values, and averaging over fast orbital phases away from important mean-motion resonances and close encounters. Dissipation, migration, and nonlinear secular effects are absent from this constant-coefficient representation.

The $g_i$ are real eigenfrequencies of the Laplace-Lagrange secular matrix, describing apsidal precession. For mode $i$, the signed real coefficients $e_{ji}$ give its amplitude in planet $j$; their relative magnitudes are fixed by an eigenvector, and their signs distinguish aligned and anti-aligned apsides. The common phase ${\beta }_i$ specifies that mode's orientation at the chosen time origin. A planet's actual orbital eccentricity is $|z_j(t)|$, which generally varies through interference of modes; it is not a sum of the positive mode amplitudes.

For one isolated perturbing planet, its Kepler orbit has fixed longitude of periapsis, so $z_1$ is constant and its secular forcing frequency is zero. A constant forced eccentricity satisfies $0=iA{z}_f+i{A}_1{z}_1$. Since $A=-A_{1}f({\alpha }_1)$,This assumes $A\ne 0$. If an external process made the planet precess at frequency $g$, the particular solution would instead be $z_f=A_1{z}_1/(g-A)$; the stated formula relies on the isolated single-planet case.

Write $ϵ={\mu }^{1/3}$ and ${\delta }_j=(a_j-a_0)/a_0$. The leading Keplerian shear isThe distinction between total mass $m$ and the individual planet's two-body central mass is order $\mu$, below the retained order $ϵ$.

At one common reference epoch, the local Kepler orbit expansion takes the formwhere ${\lambda }_{j,0}$ is the initial mean longitude relative to the rotating reference ray. Comparing its sine and cosine coefficients with the free solution of Hill equations givesThe required orbital elements from Hill coordinates are thereforeThese equalities have the first-order accuracy of the local approximation. If $e_j=0$, the longitude of periapsis is undefined. The combination $D_{4j}+2D_{2j}$ fixes the initial true longitude divided by $ϵ$. The special initial alignment imposed in part v would require $D_{4j}=-2D_{2j}$, but the general two-planet solution need not impose it.

Substitution also verifies ${\ddot{\xi }}_j-2n_0{\dot{\eta }}_j-3n_0^2{\xi }_j=0$ and ${\ddot{\eta }}_j+2n_0{\dot{\xi }}_j=0$: these are the unforced Hill equations. Close encounters add mutual-gravity forcing and can change the constants. Replacing $n_j$ by $n_0$ in the oscillation is consistent at leading order over local orbital times; accumulated phase differences must be retained when following much longer evolution.

For the single-line greenhouse deficit, use hydrostatic balance, fixed surface pressure $p_s$, constant gravity $g$, and the dry adiabatThe gas constant and specific heat refer to the background gas. Let ${\mu }_d=\cos \tilde{\theta }>0$ be the supplied effective ray cosine. The vertical optical depth and slant transmittance from level $p$ to space areThis permits pressure-dependent line widths. With local thermodynamic equilibrium, negligible infrared scattering, a black thermal surface, and no incident thermal radiation from space, the formal radiative transfer solution isThe weights add to one. Subtract the transparent-surface radiance and put $B(p)=B_{{\nu }_0}(T(p))$, using the narrow-line approximation only in the line deficit. Define its pressure-dependent equivalent widthApproximating each hemispheric flux by $\pi$ times the effective-ray radiance preserves the black-surface flux. The outgoing longwave flux is thereforeIntegration by parts gives the particularly transparent formBoth boundary terms vanish: $E(0,q)=0$ and $B(p_s)-B(p_s)=0$. The dry adiabat and monotonic Planck function give $dB/dp>0$. Atmospheric emission replaces part of the warm surface radiation by emission from cooler levels, reducing emission to space for a specified $T_s$. The inequality is strict when opacity occurs over a nonzero range with temperature contrast; an isothermal absorbing layer at the surface temperature would not produce such a contrast.

For a spherical planet the mean absorbed instellation is $F_{\ast }=(1-\alpha )L_{⊛}/4$. The required global energy balance isConsequently every equilibrium of this model has $T_s\ge T_{\mathrm{bare}}=[(1-\alpha )L_{⊛}/(4{\sigma }_{\mathrm{SB}}){]}^{1/4}$, with strict warming when $D>0$. At fixed temperature profile, increasing $q$ increases $E$ and hence $D$. This argument does not assume that the greenhouse gas changes the stellar-band albedo, surface pressure, or background thermodynamics; those are held fixed. For sufficiently broad or overlapping saturated wings the full frequency-dependent Planck function must replace the narrow-line approximation.

Interpret habitability as the ability to sustain life with broadly terrestrial biochemistry. Accessible liquid water, usable free energy, suitable chemical ingredients and sufficiently persistent conditions are the central requirements. Earth-like geography, a breathable atmosphere and actual habitation are not part of that definition.

A viable local physical environment. The relevant temperature and pressure must permit liquid water where the organisms live. For a surface habitat this normally requires an appropriate atmosphere and surface temperature; the roughly 6.1 millibar triple-point pressure of pure water is a minimum phase constraint, not a guarantee of a temperate surface. High salinity, extreme acidity and low water activity can make apparently wet environments unsuitable for known life. A subsurface ocean beneath ice can provide liquid conditions even when the surface is frozen, so the entire body need not lie in the conventional surface-water habitable zone.

A suitable orbit and rotation. For a surface-water planet, orbital distance must give an absorbed radiation budget compatible with liquid conditions after the body's Bond albedo, atmosphere and greenhouse effect are included. The simple planetary equilibrium temperature estimateuses orbital irradiance $S$, Bond albedo $A$ and the Stefan-Boltzmann constant $\sigma$, assumes unit thermal emissivity and averages over the whole sphere. It is not a prediction of surface temperature beneath a atmosphere with a greenhouse effect. Excessive absorbed energy risks a runaway greenhouse; inadequate surface heating can freeze surface water, although heat from the interior or tidal heating can maintain a subsurface habitat.

The orbit must persist without collision or ejection for a biologically useful interval. Large orbital eccentricity can produce extreme variations in heating, but neither exactly circular motion nor Earth's particular obliquity is essential. Rotation period and obliquity govern illumination and seasons; atmospheric and ocean heat flux determine how severe those variations become. Tidal locking does not automatically make a body uninhabitable, because sufficient heat transport can moderate the day–night contrast. A natural satellite also needs a stable orbit around its primary and a tolerable combination of irradiation and tidal heating: enough heating may maintain an ocean, whereas excessive heating can make the environment hostile. These are constraints on the body's location and path, without requiring a particular kind of surrounding planetary system.

An adequate chemical inventory and access to it. Known life uses carbon, hydrogen, nitrogen, oxygen, phosphorus and sulfur, together with catalytic trace elements. A rocky planet supplies mineral surfaces and reservoirs, but total inventory is not enough: chemical weathering, water-rock interaction or other transport must make nutrients accessible. A water-rich body with an ocean isolated from its rocky interior by a thick high-pressure ice layer could have difficulty recycling some nutrients; this is a possible limitation, not proof of sterility.

A usable energy supply. Photosynthesis can exploit light in an illuminated habitat. A dark subsurface habitat instead needs sustained chemical disequilibrium, for example oxidized and reduced species brought together by water-rock interaction, hydrothermal circulation or radiolysis. Warm water at complete chemical equilibrium is not, by itself, an energy supply for sustained metabolism.

Retention, protection and longevity. The body's mass, radius, surface gravity and thermal history affect retention of water and an atmosphere. Escape velocity must be considered together with gas temperature, molecular mass, atmospheric escape and replenishment: there is no universal Earth-mass threshold. An atmosphere, ice or rock overburden can protect life from damaging ionizing radiation and short-wavelength light. An intrinsic planetary magnetic field can modify charged-particle exposure and escape, but it is neither a universal requirement nor a guarantee of retention.

Finally, conditions must last long enough. Stored formation heat, radiogenic heating and possibly tidal heating can maintain circulation and chemical exchange. Plate tectonics can recycle materials and support a carbonate-silicate cycle, but it is not established as a necessary condition for every habitable body: conditional long-term carbon cycling is possible on a stagnant-lid planet as well.^[ref] The necessary end result is persistence of a suitable environment, rather than a checklist requiring every feature of the modern Earth.

Eighteenth-century support for the Newtonian system accumulated through successful predictions, the resolution of apparent discrepancies and increasingly precise terrestrial tests. The relevant system included laws of motion, universal gravitation and a programme of mathematical celestial mechanics; accepting one element did not automatically settle every optical or metaphysical dispute.

Edmond Halley identified earlier comet appearances with the same orbit and used Newtonian gravitation to predict a return. Its recovery in 1758, with the subsequent perihelion in 1759 and perturbative corrections for planetary attractions, made an especially persuasive test of a forecast made long before the event. It also brought comets into the same quantitative dynamics as planets rather than treating them as exceptional signs.^[ref]

Geodetic expeditions compared meridian arcs at different latitudes. The observed polar flattening agreed with the qualitative Newtonian expectation for a rotating self-gravitating Earth, against competing prolate interpretations. The inference required both accurate surveying and a theory of how rotation, gravity and internal mass distribution affect shape; exact flattening could not be predicted from rotation alone without an internal-structure model.^[ref]

Lunar and planetary perturbation theory provided a sustained, more demanding test than unperturbed Keplerian ellipses. Discrepancies in the Moon's apsidal motion initially appeared troublesome, but improved approximations within inverse-square gravitation recovered the missing contribution. Agreement achieved by a better controlled calculation supported the theory while demonstrating why a failure of an early approximation was not automatically a failure of the force law.

James Bradley's annual stellar aberration furnished evidence for Earth's orbital motion and finite light propagation; later nutation measurements also connected careful positional astronomy to the dynamics of Earth's axis. Aberration supports the moving-Earth framework, though it is not by itself a unique test of an inverse-square force.^[ref]

Finally, Henry Cavendish's torsion-balance experiment in 1798 measured attraction between laboratory masses and inferred Earth's mean density. It extended gravitation from celestial fitting to a quantitatively controlled terrestrial attraction; describing the original aim as determining Earth's density is historically more accurate than saying he set out to measure a modern separately defined $G$. Together these tests supplied converging support across scales. They did not prove Newtonian mechanics exactly or establish every part of Newton's natural philosophy, but they made its gravitational research programme substantially more credible.

Hydrostatic equilibrium means that the local pressure gradient balances gravity, with no large acceleration or strength-supported departure from that balance. In a spherical planet,Here $M(r)$ is the enclosed mass, not a constant total mass at every interior radius. If composition is uniform and the radial thermodynamic path is adiabatic, the relevant bulk modulus is $K_s=\rho (dP/d\rho {)}_s$. The chain rule then gives the Adams-Williamson relationA phase transition or composition jump must be treated separately; its density change is not obtained by integrating one smooth homogeneous-adiabatic equation through it.

The seismic parameter isThe shear modulus $\mu$ measures resistance to shear deformation. Solids with $\mu >0$ transmit shear waves; an ideal fluid has zero static shear modulus, so it cannot carry a bulk S wave, while compressional P waves can still propagate through its nonzero bulk modulus.

Limitations include uncertain seismic coverage, temperature and composition tradeoffs, nonadiabatic boundary layers, phase changes, lateral heterogeneity, rotation and finite material strength. Small asteroids and comet nuclei commonly retain irregular shapes because their self-gravity is too weak to overcome strength and impact history; a few larger bodies, such as Ceres, are much closer to global hydrostatic figures. A spherical shape alone is not a complete proof of present hydrostatic equilibrium.

Two complementary remote constraints are precise spacecraft tracking, which measures gravitational harmonics and constrains mass distribution when combined with shape and rotation; and tidal or rotational measurements, such as Love numbers, librations and spin precession, which constrain deformability, moment of inertia and possible fluid layers. Magnetic induction measurements in a time-varying external field can provide an additional conductivity constraint, but are not a unique density inversion.

For the June 2018 lunar picture, start with a regolith and fractured crust, a largely anorthositic crust tens of kilometres thick, a silicate mantle, a deep partially molten region and a small iron-rich core. GRAIL and Apollo-based models gave a global mean crustal thickness around $34$-$43$ km.^[ref] A representative seismic model has a partial-melt boundary near radius $480$ km, a liquid outer-core boundary near $330$ km and a solid inner-core radius near $240$ km; these are model-dependent radii from the centre, not depths below the surface.^[ref]

The long-running Apollo network had only four principal stations on the near side, with poor global coverage and strong scattering in the fractured near-surface material. Far-side structure and weak deep reflected phases were therefore difficult to constrain. The historical date matters: the answer does not import later Mars seismology or later lunar-core refinements into the 2018 premise.

In core accretion, solid material grows from dust and aggregates into planetesimals or rapidly accreted pebbles. A sufficiently large core binds an atmosphere; beyond the snow line, additional condensed volatiles aid solid growth. Cooling allows the envelope to contract and acquire more gas. Once gas self-gravity becomes important, rapid gas accretion can build a giant planet until disk supply, gap formation or disk dispersal limits growth. Migration and interactions with other bodies can alter its final orbit.

Heavy-element enrichment in giant planets, evidence for substantial interior heavy material, the prevalence of smaller solid planets and the correlation of giant-planet occurrence with host-star metallicity all support an important role for solid accumulation. They do not uniquely establish one growth pathway. A classic difficulty is forming cores and cooling envelopes before a gas disk disperses, especially far from the star; migration, envelope opacity, solid-delivery rates and fragmentation complicate the calculation. Pebble-assisted growth can change the timescale, and disk gravitational instability remains a distinct possible route in some environments.

Planetary luminosity is the total emitted power. For the Kelvin-Helmholtz budget one needs the intrinsic luminosity, excluding reradiated absorbed starlight. With orbital distance $r_p$, Bond albedo $A_B$, solar luminosity $L_{\odot }$ and present-like radius $R_J$,The given ratio of emitted to absorbed power is $1.67$, so $L_{\mathrm{int}}=0.67P_{\mathrm{abs}}$. Using the supplied radius, $r_p=5.2$ AU and a solar luminosity around $3.8\times 10^{26}$ W givesconsistent with the quoted $3.42\times 10^{17}$ W. Taking $1.67P_{\mathrm{abs}}$ as intrinsic power would incorrectly include reradiated sunlight.

The Hill radius is the scale at which stellar tidal acceleration competes with the planet's attraction. It bounds the approximate region in which gas can remain gravitationally associated with the planet and helps set feeding and envelope boundary scales. With a Jupiter mass,Using this as the stipulated extended proto-planet radius and the present-like intrinsic luminosity givesUsing the compact present radius instead would give an even longer value, of order $3\times 10^{11}$ years at that same luminosity. Neither number is a complete formation time: mass, radius, luminosity and envelope structure evolve together. The Hill radius is not automatically the radius holding most of the binding energy.

A constant present-like luminosity with this simple energy budget is too slow compared with a gas-disk lifetime of a few million years. At the Hill-radius energy scale, radiating in $3$ Myr would require an average intrinsic luminosity of order $5\times 10^{19}$ W, much brighter than the adopted value. This diagnoses a failure of the constant-current-Jupiter assumptions, not a general disproof of core accretion. Formation models must compute a much more luminous, evolving accretion and contraction history.

A planet's total emitted radiant power. Its intrinsic luminosity excludes reradiated absorbed starlight and is the appropriate cooling power in a simple contraction-energy timescale.

Planetary migration changes a planet's orbital semimajor axis through torques from a disc or exchanges of energy and angular momentum with other bodies, often followed by tidal dissipation. Short-period giant planets, eccentric progenitors, resonant architectures, and stellar obliquity provide complementary constraints rather than one unique migration diagnostic.