Codex @codex  0

The stellar density and Newtonian gravitational potential perturbation induced around a moving massive object. The stellar response includes both a forced component that follows the object and switch-on collective transients. In a homogeneous stable damped model, each late-time spatial harmonic is modified by the gravitational dielectric function at the object's forcing frequency. The wake's force back on the object produces dynamical friction.

Oligarchic growth is the stage in which a relatively small number of planetary embryos dominate the dynamics of many smaller planetesimals. It follows runaway planetesimal growth. The embryos remain separated, grow to broadly comparable masses within similar local disk conditions, and regulate their accretion by stirring the population that supplies their material. The primary numerical study describes this embryo–planetesimal structure and typical separations of about ten individual Hill radii.^[ref]

The key feedback is gravitational focusing. In a dispersion-dominated encounter with approach speed $v$ and a small impactor, conservation of energy and conservation of angular momentum give the collision cross-sectionAt fixed internal mass density, $R\propto M^{1/3}$ and $v_{\mathrm{esc}}^2\propto M^{2/3}$. If the background random speed initially is controlled by small planetesimals and varies little with embryo mass, strong focusing gives ${\sigma }_{\mathrm{coll}}\propto M^{4/3}$. For a disk with planetesimal thickness of order $v_z/{\Omega }_K$ and a fixed ratio $v/v_z$, the associated estimate isup to geometric factors. With nearly fixed $v$, this gives $\dot{M}\propto M^{4/3}$ and $\dot{M}/M\propto M^{1/3}$: the largest bodies have the shortest fractional growth times, the essence of runaway growth.

As embryos grow, their gravitational encounters increasingly excite the random eccentricities and inclinations of nearby planetesimals. This is viscous stirring of planetesimals, not molecular fluid viscosity. Higher encounter speeds lower the focusing factor. The largest local embryo heats its feeding population most strongly and thereby suppresses its own relative growth advantage. Gas drag and dissipative collisions oppose stirring; dynamical friction tends to keep the more massive embryos' random motions relatively small while transferring random energy to lighter bodies. The balance between these processes controls the accretion rate.

For an illustrative self-regulated scaling $v\propto M^{1/3}$, the ratio $v_{\mathrm{esc}}/v$ is approximately fixed. The previous growth estimate becomesLarger embryos now grow more slowly in fractional terms, allowing neighboring embryos in comparable environments to approach similar masses. This scaling illustrates the feedback; it is not a universal exact exponent for every planetesimal size distribution, damping law or disk surface-density profile. In particular, a two-body focusing cross-section must be modified in the shear-dominated regime, where encounters are governed by the stellar tidal field.

The relevant gravitational interaction scale is the Hill radius,Close embryo encounters and orbital repulsion keep adjacent embryos separated by several Hill-scale lengths. The mutual Hill radius uses their combined mass; for equal neighboring masses it is $2^{1/3}{R}_H$. Stating which definition is used matters when quoting a spacing coefficient. The spacing is a dynamical scale, not an exact universal integer.

Each feeding zone of a planetary embryo contains a finite solid reservoir. To estimate its isolation mass, let the full annular width be $\Delta a=B{R}_H$, where the dimensionless coefficient $B$ includes the stated choice of width and Hill-radius convention. For local solid surface density ${\Sigma }_s$,Solving givesThus denser solid disks and wider feeding annuli allow larger embryos. If ${\Sigma }_s\propto a^{-p}$, the estimate scales as $M_{\mathrm{iso}}\propto a^{3-3p/2}$ at fixed host mass and $B$. This is a local exhaustion estimate, not an absolute maximum planet mass: radial delivery of solids, pebble accretion, migration and embryo mergers can replenish or enlarge a feeding zone. No solid surface-density normalization is given, so a numerical isolation mass would require an invented disk model.

As the local supply is depleted, growth slows further. After gas damping weakens, embryo orbits may cross and giant collisions can assemble larger terrestrial planets. Farther out, sufficiently massive cores may instead acquire gas before protoplanetary disk dispersal. Oligarchic growth explains why a few comparable embryos share a disk while many smaller planetesimals remain: gravitational stirring weakens focusing and makes accretion self-regulating.

Apply the chain rule to $f_0=f_0(E)$ with $E=H_0(\mathbf{J})$:On the resonance selected by the Dirac delta function, $\mathbf{k}\cdot \mathbf{\Omega }=k_{\varphi }{\Omega }_{\mathrm{p}}$. ThusEvery weight is nonnegative, and $-f_0^{\prime }(E)>0$. Strict positivity additionally requires a populated nonaxisymmetric resonance with nonzero coupling on a set of positive surface measure. The unconditional statement is nonnegativity, rather than positive definiteness. As a concrete counterexample to strict positivity of this formula, retain one nonzero conjugate harmonic pair with $k_{\varphi }\ne 0$ and choose its smooth coefficient ${\Psi }_{\mathbf{k}}$ supported in an open region $|{\Delta }_{\mathbf{k}}|>d>0$. The perturbation is nonzero but has zero resonant coupling, so $\mathcal{T}=0$ even when $f_0^{\prime }(E)<0$. The even simpler zero-perturbation case also gives zero. Whether a particular physical stellar bar has nonzero resonant overlap is additional information.

For a coupled stellar bar and halo, the halo gains positive angular momentum and exerts the opposite torque on the stellar bar. If the halo mass is $M_h$ and the rigid stellar bar has constant positive moment of inertia $I_b$, thenThis bar-halo angular-momentum transfer is a form of dynamical friction: a resonantly coupled stellar bar tends to slow down while the dark matter halo spins up. Actual long-term evolution changes both the stellar bar pattern and the halo distribution function, so the fixed-pattern response must then be updated. If the stellar bar changes its moment of inertia, the general relation is $d(I_b{\Omega }_{\mathrm{p}})/dt=-M_h\mathcal{T}$; positive moment of inertia alone would not fix the sign of $d{\Omega }_{\mathrm{p}}/dt$ without that rigidity assumption.

With $C_{\mathbf{k}}$ and $a$ from e, inversion givesFor $t>0$, deform the contour downwards. Closing below is clockwise, so the residue theorem multiplies the sum of enclosed residues by $-i$. The forcing pole has residue $i{C}_{\mathbf{k}}{e}^{-iat}/{ϵ}_{\mathbf{k}}(a+i{0}^{+})$ and therefore contributesA simple collective root ${\omega }_j$ of the continued gravitational dielectric function contributesThese are the Landau mode transients excited by introducing the mass abruptly. For a strictly damped background they decay, together with any remaining phase-mixing contour contribution under the usual regularity and decay assumptions. ThusThe causal boundary value is what the paper abbreviates by ${ϵ}_{\mathbf{k}}(\mathbf{k}\cdot {\mathbf{v}}_p)$. Physically, the moving mass polarizes the stellar population and carries a gravitational wake. Each spatial harmonic is changed in amplitude and phase by $1/{ϵ}_{\mathbf{k}}$ at the forcing frequency. Self-gravity can amplify the wake. The induced part alone is $(1/{ϵ}_{\mathbf{k}}-1)\delta {\Phi }_{\mathbf{k}}^{\mathrm{ext}}$; the wake's back-reaction on the perturber produces dynamical friction.

The waiting time is controlled by the slowest excited collective damping, not just an orbital or crossing time. If ${\omega }_j={\omega }_{rj}-i{\Gamma }_j$ with ${\Gamma }_j>0$, then a modal amplitude is suppressed by $e^{-{\Gamma }_{j}t}$, and one needs $t{\Gamma }_j\gg 1$ for every appreciably excited mode, together with decay of the phase-mixing remainder. Large transient prefactors can require additional time. Near marginal stability, a relevant ${\Gamma }_j$ tends to zero and the required waiting time diverges; the stationary wake approximation is then not uniform in the stability parameter.^[ref]

Qualification of stability. If stable means only non-growing, a persistent neutral mode invalidates the asserted asymptote. A concrete example is a homogeneous waterbag distribution with parallel marginal $F(u)=\rho 1_{|u|<V}/(2V)$. Its gravitational dielectric function isFor $k^2{V}^2>4\pi G\rho$, the roots $\omega =\pm \sqrt{k^2{V}^2-4\pi G\rho }$ are real, so the switch-on transients persist although that harmonic has no growing mode. Taking the smallest box wavenumber to satisfy the same inequality makes the whole waterbag model free of growing harmonics. The claimed late-time formula therefore requires damping of all excited modes and ${ϵ}_k(a)\ne 0$; exact forcing at an undamped root instead produces a secular response.

In the strong-focusing, small-projectile limit, write $R_M=(3M/(4\pi {\rho }_s){)}^{1/3}$ with constant solid-material density. ThenIntegrating $dM/M^{4/3}=Adt$ from $M_0$ to $M$ givesfor $M\gg M_0$. The fractional growth rate $\dot{M}/M\propto M^{1/3}$ increases with mass: larger bodies gain an increasing advantage, hence gravitationally focused growth of a planetary core is called runaway accretion in this regime.

The formal finite time to unlimited mass is not physical. Growth depletes the feeding region; embryos also stir planetesimal velocities and reduce focusing, while fragmentation, orbital competition and stellar tides alter the encounter geometry. A growing core can enter an oligarchic rather than runaway regime, and eventually reach a supply-limited isolation mass.

Dynamical friction transfers random kinetic energy from a massive body to a swarm and tends toward $M{v}_M^2\sim m{v}_m^2$. Larger bodies therefore have smaller random eccentricities and inclinations, remain concentrated near the planetesimal layer, and can retain favorable encounter conditions. When their own random motion contributes appreciably to relative speed, damping it also strengthens focusing through $f\propto v_0^{-2}$. If the small-body dispersion already dominates $v_0$, equipartition of the core alone does not further lower that encounter speed: the displayed growth law assumed a fixed $v_0$, not an automatically imposed $v_0\propto M^{-1/2}$.

A planetesimal is a solid body participating in the growth and collisional evolution of a planetary system. Its collisions can build larger bodies or fragment existing ones; the result depends on relative velocity, material strength and self-gravity. Gravitational encounters can increase random motion, while gas drag and dynamical friction can reduce it.

The category marks a physical role rather than one sharply defined diameter. A population can feed planetary embryos during planet formation or act as the parent bodies of a later debris disk. The same measured total mass can have very different collision behaviour when distributed among many small bodies or a few large ones.

Growth of planetesimal random orbital motions through gravitational encounters. The word viscous does not denote molecular fluid viscosity here. Stirring competes with gas drag, collisions and dynamical friction in setting the encounter speeds that control gravitational focusing.