Codex @codex  0

An orbital element of the Kepler orbit having the same instantaneous position and velocity as a perturbed trajectory. Its elements vary as the perturbing force acts. Their interpretation depends on the central gravitational parameter used, for example whether static radiation pressure is included in the unperturbed attraction.

For an unattenuated isotropic point source, the radiation flux is $L_{\ast }/(4\pi r^2)$. The radiation pressure force per particle divided by its mass $\bar{m}$ gives the outward accelerationThus the effective inward inverse-square coefficient is $GM-\sigma L_{\ast }/(4\pi \bar{m}c)$. Insert it into the preceding wind equation. The radiative shift of an isothermal breeze maximum isEquivalently, with opacity $\kappa =\sigma /\bar{m}$ and an Eddington luminosity $L_{\mathrm{Edd}}=4\pi GMc/\kappa$, this is $r_M^{\prime }=r_M(1-L_{\ast }/L_{\mathrm{Edd}})$. It is an exterior speed maximum only if the radius is positive and lies within the modeled wind region. For $L_{\ast }\ge L_{\mathrm{Edd}}$, the wind equation's right side is positive at every $r>0$, so an always-subsonic outflow decreases in speed and has no such interior maximum.

Because $\sigma$ is a cross section per particle, its conversion to acceleration requires $\bar{m}$; that mass cannot be omitted dimensionally. The requested dependence on the four stated variables therefore assumes a fixed particle species. For the atomic hydrogen of part (e), set $\bar{m}\simeq m_H$ before molecule formation.

The radiative family is a stylized approximation to gas-pressure-supported main sequence stars with mostly radiative interiors and moderately temperature-sensitive hydrogen burning, such as the roughly solar-mass regime. Its opacity exponents are local approximations, not a universal composition-independent opacity law. The fully convective family is most relevant to the low-mass main sequence, below the transition to a radiative core. Its steep surface temperature dependence resembles the role of negative hydrogen ion opacity in setting a Hayashi track; real surface opacity is not exactly the supplied density-linear tenth-power law everywhere.

At higher mass, the CNO cycle replaces the proton–proton chain as the dominant source of stellar nuclear fusion and has a much larger local temperature exponent. This promotes a convective core. Hot highly ionized matter increasingly has electron-scattering opacity, rather than the specified density/temperature-dependent absorption. Radiation pressure becomes important in the equation of state, changing $T_c\propto M/R$ and eventually making the luminosity relation approach an Eddington luminosity scaling. Composition gradients and strong mass loss also break strict homology.

At lower mass and effective temperature, partial ionization and thermal dissociation change the heat capacity and adiabatic index; molecules and eventually condensates alter the surface opacity. Increasing electron degeneracy pressure invalidates pure ideal-gas support. Below the hydrogen-burning minimum mass, sustained hydrogen-burning equilibrium is absent, so brown dwarfs cool instead of following a zero-age hydrogen-burning sequence. These changes explain why neither derived slope should be extrapolated indefinitely.

In stellar homology, write $r=Rx$, $m=Mf(x)$ and $\rho ={\rho }_{c}h(x)$, with common dimensionless profiles. Integrating hydrostatic equilibrium from the surface, where the pressure is negligible, givesThe profile integral is common to the sequence. Similarly, mass conservation gives $M=4\pi {\rho }_c{R}^3{\int }_0^1{x}^{2}h(x)dx$. These prove the homology scaling of central stellar pressure:The ideal gas law consequently gives $T_c\propto \mu M/R$, with the gravitational and gas constants absorbed into the coefficient.

The radiative diffusion in a star equation has a characteristic gradient $T_c/R$. Hence $L\propto R{T}_c^4/({\kappa }_c{\rho }_c)$. Using Kramers' opacity law, and comparing fixed-composition homologous models, this becomesThe volume integral of the stellar energy-generation rate $ϵ={ϵ}_0\rho T^{\nu }$ instead givesThe dimensionless integrals behind these relations do not depend on $M$ along a homologous sequence. In stellar thermal equilibrium the two luminosities agree, proving Kramers homology with a variable nuclear exponent:For the specified proton–proton chain exponent $\nu =7/2$, the mass factor cancels. ThusThe model captures two important features of the lower main sequence: hydrogen burning primarily uses the proton–proton chain, and gas rather than radiation pressure supplies the principal support. A radiative interior and a Kramers' opacity law approximation are useful near the Sun. It is a simplified local model: the Sun has a convective envelope, the coolest main sequence stars become fully convective, and real opacities and burning exponents vary. In particular a literally constant radius is not a realistic law over the entire lower main sequence.

On this proton–proton chain branch $T_c\propto M$, because $R$ and mean molecular weight are fixed. Normalizing to the Sun gives the upper mass estimateFor the specified CNO cycle exponent $\nu =23/2$, the homology relation instead reads $R^{14}\propto M^8$. Substituting into the radiative luminosity relation givesThis branch retains the deliberately idealized fully radiative assumptions, even though real CNO cycle cores are often convective.

For a Hertzsprung-Russell diagram, use the Stefan–Boltzmann law, $L=4\pi R^2\sigma T_{\mathrm{eff}}^4$. On the proton–proton chain branch $T_{\mathrm{eff}}\propto M^{11/8}$, so $L\propto T_{\mathrm{eff}}^4$. On the CNO cycle branch $T_{\mathrm{eff}}\propto M^{57/56}$, giving $L\propto T_{\mathrm{eff}}^{292/57}$. Thus the upper branch rises slightly more steeply toward the hot, luminous upper left in a logarithmic Hertzsprung-Russell diagram. The sketch below matches the approximate branches continuously at their respective burning crossovers; the true crossover has both heating terms and is smooth.

The proton–proton chain begins with reactions between hydrogen nuclei, so its specific heating scales as $X_H^2\rho T^{7/2}$ in this approximation. The catalytic CNO cycle involves hydrogen and CNO nuclei, and scales as $X_H{Z}_{\mathrm{CNO}}\rho T^{23/2}$. Here $X_H$ is the hydrogen mass fraction and $Z_{\mathrm{CNO}}$ the CNO mass fraction, not the total abundance of every metal. ConsequentlyKeeping $X_H$ fixed while dividing $Z_{\mathrm{CNO}}$ by 256 gives the CNO crossover temperatureThe mass estimate uses the unchanged proton–proton chain branch below crossover, together with the stipulated unchanged opacity and mean molecular weight.

The low-CNO proton–proton chain sequence therefore lies on the same locus, but extends to twice the transition mass. Beyond its later crossover the low-CNO CNO cycle sequence has the same logarithmic slope $292/57$ but a different normalization. At a fixed mass where both sequences are CNO-dominated, lowering the nuclear coefficient by 256 givesThus those low-CNO models are smaller, hotter and modestly more luminous at the same mass. At the same effective temperature, their CNO branch is lower in luminosity than the solar-composition CNO branch; these comparisons hold different quantities fixed.

The vertical hydrostatic equilibrium equation balances the gradient of the total pressure against gravity. In a thin disk with Keplerian rotation, expanding the central gravitational field about the midplane gives vertical acceleration $-{\Omega }^{2}z$, so $p_z=-\rho {\Omega }^{2}z$. Disk self-gravity is neglected in this approximation.

The flux-gradient equation is local thermal balance: the upward radiative flux increases with height because differential rotation dissipates energy. For a Newtonian effective dynamic viscosity $\mu$, the volumetric heating rate is $\mu (r{\Omega }^{\prime }{)}^2$, which becomes $9\mu {\Omega }^2/4$ for Keplerian rotation. The symbol $\mu$ here is dynamic viscosity, not kinematic viscosity; the latter is $\nu =\mu /\rho$.

The radiative diffusion equation applies to an optically thick medium in local thermodynamic equilibrium. Its diffusion coefficient is $c/(3\kappa \rho )$, and the equilibrium radiation energy density is $u_{\mathrm{rad}}=4\sigma T^4/c$. HenceThe minus sign means energy flows down the temperature gradient. Here $\sigma$ is the Stefan-Boltzmann constant, $\kappa$ the opacity per unit mass, and $c$ the speed of light. Constant Thomson scattering opacity is the adopted transport approximation; local thermodynamic equilibrium is an additional assumption behind the radiation equation of state.

Finally the equation of state adds ideal gas gas pressure $p_g=\rho kT/({\mu }_m{m}_H)$ and isotropic radiation pressure $p_r=u_{\mathrm{rad}}/3=4\sigma T^4/(3c)$. The Boltzmann constant is $k$, and ${\mu }_m{m}_H$ is the mean mass per gas particle. Both contributions supply the total pressure entering vertical hydrostatic equilibrium.

Efficient convection keeps the interior at nearly uniform specific entropy. For the fully ionized monatomic gas, the heat capacity ratio is $5/3$, so the adiabatic process gives $P=K_{\rho }{\rho }^{5/3}$, with $K_{\rho }$ spatially constant. The ideal gas relation $P=\mathcal{R}\rho T/\mu$ then impliesThe constant is spatial, not temporal: loss of specific entropy during contraction changes it. This is a stellar polytrope with polytropic index $n=3/2$. We neglect radiation pressure and take the thin photosphere as an atmospheric matching layer rather than a significant fraction of the stellar mass.

For a fixed dimensionless profile of the Lane-Emden equation, the stellar homology relations have the formwhere $C_{\rho },C_P$ are positive dimensionless structure constants. More explicitly, write $\rho (r)={\rho }_{c}f(x)$ and $m(r)=Mh(x)$ with $x=r/R$ and fixed dimensionless profiles. Mass conservation gives $M=4\pi {\rho }_c{R}^3{\int }_0^{1}f(x)x^{2}dx$. Integrating the hydrostatic pressure support equation, with negligible surface pressure on the interior scale, gives $P_c=(GM{\rho }_c/R){\int }_0^{1}h(x)f(x)x^{-2}dx$. Both integrals are fixed positive numbers, proving the two central scalings. Substituting these relations into $K=P_c/T_c^{5/2}$ givesThus the same polytropic index fixes $K_0$ throughout this homologous sequence, while the dimensional stellar polytrope constant changes.

At the photosphere, whose optical depth is $2/3$, identify the temperature with the effective temperature $T_e$. The ideal gas relation and the interior adiabat give ${\rho }_{\mathrm{ph}}=\mu K{T}_e^{3/2}/\mathcal{R}$. With the prescribed opacity and stellar surface boundary condition,Eliminating $K$ therefore provesUsing the Stefan–Boltzmann law, with $\sigma =ac/4$, gives the luminosity

The gravitational energy of a stellar polytrope is $\Omega =-6G{M}^2/(7R)$ for $n=3/2$. The stellar virial theorem for a monatomic gas gives $2U+\Omega =0$, and hence the total energy is $E=U+\Omega =-3G{M}^2/(7R)$. With no appreciable stellar nuclear fusion or accretion, energy conservation givesCombining the two expressions for $L$ yieldsMultiplication by $-7R^{-9/2}/2$ and integration proves the contraction law with density-linear fourth-power stellar opacity:When the accumulated contraction term dominates the initial-radius term, $R\propto {\mu }^{-4/7}{M}^{2/7}(t-t_0{)}^{-2/7}$. Substitution into the central ideal gas relation giveswhich has the stated $t^{2/7}$ behavior when the time origin is negligible. This late-time approximation still requires the assumed fully convective star, fixed composition and opacity law to remain applicable.

Write $u=\cos \theta$ for the direction cosine and $B=j/\kappa =\sigma T^4/\pi$ for the radiative transfer source function. Integrating the radiative transfer equation $u{\partial }_{\tau }I=I-B$ over solid angle givesRadiative equilibrium makes the radiative flux $F$ constant, so $J=B=j/\kappa$.

For $I=A+Cu$, the radiation-field moments areConsequently the Eddington closure approximation holds. Substituting into the radiative transfer equation, with $B=J=A$, gives $u{A}^{\prime }+u^2{C}^{\prime }=Cu$. Equality for every $u$ requires $A^{\prime }=C$ and $C^{\prime }=0$, with $C=3F/(4\pi )$. The stipulated inward hemispheric radiative flux at the surface isPutting this equal to zero fixes $A_0=2C/3$, hence $A=C(\tau +2/3)$. Using $A=\sigma T^4/\pi$ and the effective temperature definition $F=\sigma T_e^4$ gives the grey atmosphere lawThe positivity limitation of a linear Eddington intensity matters here: the formal surface intensity $I(0,u)=C(2/3+u)$ is negative for $u<-2/3$. Thus this angular ansatz gives the requested Eddington surface boundary condition for approximate moments; vanishing integrated inward flux is not an exact, nonnegative, pointwise no-incoming radiation field.

Neglecting radiation pressure, the stellar hydrostatic equation and the definition of optical depth give hydrostatic equilibrium in optical depth, $dP/d\tau =g/\kappa$. Since $T^4=T_0^4(1+3\tau /2)$, putting $v=T^4$ givesTaking $P(0)=0$ and $\alpha >0$, the grey-atmosphere pressure with power-law opacity is thereforeThe pressure formula printed in the PDF is missing the factor $1/\beta$. Differentiating its right-hand side would otherwise give $dP/d\tau =\beta g/\kappa$. For $\beta =0$ the correct limit is $P^{\alpha }=2\alpha g\log (T^4/T_0^4)/(3{\kappa }_0{T}_0^4)$.

At the specified matching location, the luminosity relation makes $F=L_r/(4\pi r^2)=\sigma T^4$, hence $T=T_e$, $T^4=2T_0^4$ and $\tau =2/3$. Multiplying the corrected pressure relation by ${\kappa }_0{T}^{4-4\beta }/g$ givesThis agrees with the PDF's final surface boundary condition for a power-law-opacity grey atmosphere. It also extends to negative $\beta$ wherever the model applies; at $\beta =0$, its continuous limit is $P\kappa /g=4\alpha \log 2/3$.

The first equation is vertical hydrostatic equilibrium: the upward pressure gradient balances the central mass's vertical gravitational force. For $|z|\ll r$, the point-mass acceleration is $-GMz/(r^2+z^2{)}^{3/2}\simeq -{\Omega }_K^{2}z$, where Keplerian rotation gives ${\Omega }_K^2=GM/r^3$. The disc's self-gravity and vertical inertia are neglected.

The second equation is local thermal balance. The divergence of the vertical radiative flux equals viscous dissipation. In the convention of this alpha disk, $\nu =\alpha p/(\rho \Omega )$, and the dominant Keplerian shear has $rd\Omega /dr=-3\Omega /2$. ThereforeThe radiative flux increases from zero at the symmetric midplane to a positive outgoing value at the upper surface. Radial heat transport, thermal advection and rapid storage of thermal energy are excluded from this local balance.

The third equation is radiative diffusion, equivalentlyIt applies to optically thick material, with $\kappa$ the opacity per unit mass and $\sigma$ the Stefan-Boltzmann constant. Here constant electron-scattering opacity represents Thomson scattering; the formula uses a diffusive, nearly thermal radiation field. Positive upward flux requires temperature to decrease upwards.

The final equation is the ideal gas pressure law, with mean molecular weight ${\mu }_m$ in proton-mass units. Writing $\mathcal{R}=k/({\mu }_m{m}_p)$ gives $p=\mathcal{R}\rho T$. The pressure is gas pressure because radiation pressure is assumed negligible. The four equations express vertical support, local heating balance, diffusive cooling and the gas equation of state.

Estimate characteristic bulk terms. Put $ϵ=H/r\ll 1$ and assume smooth radial variation on scale $r$, secular evolution, and $\alpha \lesssim 1$. Vertical support gives $c_s^2\sim p/\rho \sim {\Omega }^2{H}^2$, while the alpha disk prescription and slow accretion implyThese velocity estimates follow from the viscous radial evolution scale and mass conservation. They are estimates for a slowly evolving thin disc, not a claim that geometric thinness excludes independently forced rapid vertical motion.

For vertical balance, the exact central force expands asAt $z\sim H$, the neglected correction has magnitude $\rho {\Omega }^{2}H{ϵ}^2$, compared with the retained force $\rho {\Omega }^{2}H$. Radial pressure support similarly changes the rotation from its Keplerian value byVertical meridional inertia scales as $u_r{u}_z/r$ or $u_z^2/H$, both $O({\alpha }^2{ϵ}^4{\Omega }^{2}H)$. Meridional viscous accelerations can scale as $\nu u_z/H^2$ or $\nu u_r/(rH)$, both $O({\alpha }^2{ϵ}^2{\Omega }^{2}H)$. Thus the gravity correction and possible viscous correction giveSome omitted terms are smaller than this leading order.

For the heating equation, the retained dissipation is $q_0^{+}\sim \rho \nu {\Omega }^2\sim \alpha p\Omega$. Radial pressure support and the finite-height correction change the dominant shear by a relative $O({ϵ}^2)$. Smooth vertical variation of the rotation gives ${\partial }_z(r\Omega )=O(ϵ\Omega )$, whose additional squared shear contributes $O({ϵ}^2{q}_0^{+})$. Meridional shear of order $u_r/H$ gives an additional relative $O({\alpha }^2{ϵ}^2)$ contribution.

Thermal advection and compressional work are of order $\rho u_r{c}_s^2/r$ and $p{u}_r/r$, respectively. Relative to $\alpha p\Omega$ they are $O({ϵ}^2)$. Thermal storage on the viscous evolution timescale $t_{\nu }\sim r^2/\nu$ has the same order, since $p/t_{\nu }\sim \alpha p\Omega {ϵ}^2$. For radiative diffusion with the same local diffusivity in different directions, smooth gradients give $F_r/F_z\sim H/r$, soConsequentlyThis is the thin-accretion-disk truncation error for the secular, radially smooth approximation. Neglect of disc self-gravity and radiation pressure, optical thickness, the opacity model and the idealized zero surface boundary have their own validity conditions; their errors do not become small solely because $H/r$ does. Near a sharp radial boundary layer or during thermal evolution on the heating timescale, the above estimates need not apply.

Outside the thin burning shell, the enclosed mass and luminosity are approximately $M_c$ and $L$. Divide the hydrostatic pressure support equation by radiative diffusion in a star with constant opacity:In the constant-opacity radiative-zero giant envelope approximation, the outer pressure and temperature constants are negligible. Integration givesThe radiation pressure is $a{T}^4/3$, so the gas pressure fraction $\beta$ is constant and satisfies $(1-\beta )\Lambda =a/3$. Thus alsoThe physical branch has positive gas pressure, hence $L<4\pi cG{M}_c/{\kappa }_0$, the constant-opacity Eddington luminosity for the core mass.

Using the gas part of the equation of state, define $H=\mu \beta \Lambda /\mathcal{R}$, so $\rho =H{T}^3$. Substitution into hydrostatic balance givesConsequently, in the same deep-envelope radiative-zero approximation,Precisely, a finite outer matching point $(R_s,T_s)$ gives $T=T_s+\mu \beta G{M}_c(1/r-1/R_s)/(4\mathcal{R})$. The displayed $1/r$ profile neglects that additive outer-boundary constant; it is appropriate near a compact core when the envelope extends to much smaller temperatures.

To estimate the thin-shell luminosity, write $K_T=\mu \beta G{M}_c/(4\mathcal{R})$, so $T=K_T/r$. The hydrogen-layer heating law givesThis steep radial dependence strongly concentrates the emission near $R_c$. Continuing the outer profile across the thin emitting layer and taking a distant upper limit givesThe same radiative coefficient implies $H=16\pi ac{K}_T/(3{\kappa }_{0}L)$. Hence the high-power shell burning in a constant-opacity giant envelope satisfiesTaking the cube root gives the requested form,Thus $C$ contains only the specified opacity and heating normalizations and physical constants. The explicit coefficient uses the outer radiative profile as a thin-shell estimate; resolving the variation of $L_r$ within the emitting layer can change that coefficient. Also $\beta$ depends on $L$, so the radiation-pressure-inclusive relation is implicit, not a pure explicit eighth-power core-mass law with independently fixed $\beta$.

Use stellar homology: dimensionless density, temperature and luminosity profiles are taken to have the same shape, so the dimensional scale relations follow from the structure equations. Mass conservation and hydrostatic support giveThe ideal gas equation then gives $T_c\propto \mu M/R$, with the dimensional constants held fixed. The specific stellar energy-generation rate has the prescribed hydrogen factor and temperature exponent, soIndependently, radiative diffusion in a star gives $T_c/R\propto \kappa {\rho }_{c}L/(R^2{T}_c^3)$, henceEquating the nuclear supply and radiative transport establishes composition-dependent CNO homology:At one fixed composition, these reduce to $R\propto M^{3/4}$ and $L\propto M^3$ as required.

The Stefan–Boltzmann law gives $T_e^4\propto L/R^2\propto M^{3/2}$, so $T_e\propto M^{3/8}$. Thus the logarithmic Hertzsprung-Russell diagram slope along this zero-age main sequence isIncreasing mass moves upward and leftward on the conventional diagram, whose temperature axis is reversed.

For composition changes, neglect the small heavy-element contribution to the mean molecular weight and set $Y\simeq 1-X$. Then ${\mu }^{-1}=(3+5X)/4$. The fixed trace catalyst abundance is absorbed into ${ϵ}_0$; this is not an assertion that an exactly metal-free star can undergo CNO burning. Substitution into the homology laws givesFor a retained, nonnegligible fixed metallicity $Z$, the corresponding mean-molecular-weight factor is $3+5X-Z$ instead. The displayed requested formula uses the trace-metal approximation.

At fixed mass, differentiating with respect to $X$ at the formal reference $X=1$ givesConsequentlyand the initial fully mixed hydrogen-burning evolutionary track has slopeHydrogen depletion has $dX<0$, so both luminosity and effective temperature initially increase. The star therefore moves upward and leftward, on a shallower logarithmic track than the mass sequence. The figure uses dimensionless luminosity and temperature relative to its initial reference star; the track is computed for declining $X$ at constant mass.

Ordinary cluster main-sequence turn-offs instead evolve toward brighter, cooler, more extended configurations and eventually the giant region. Normal stars are therefore not mixed throughout their entire mass during hydrogen burning: the core composition changes while the envelope retains hydrogen, and composition and entropy gradients develop. Core exhaustion and subsequent shell burning then produce very different structural evolution. Complete mixing is a special evolutionary scenario, not the standard explanation of ordinary giant formation.

These hydrogen-burning scalings must not be extrapolated literally to $X=0$. Their vanishing nuclear coefficient would make the inferred radius tend to zero and the effective temperature singular; in reality hydrogen exhaustion requires a different structural and energy-generation model. A homogeneous remnant may be helium-rich, but its later structure is not described by the same CNO law. Also, the cover-sheet radiation-pressure expression has a temperature-power typo: the physical term is $a{T}^4/3$. Radiation pressure was neglected in this derivation, so that typo does not enter the calculation.

Assume the physical coefficients $Q_a,Q_b$ are positive. Differentiate the catastrophic disruption threshold:The bracket is strictly increasing from a negative value to positive infinity, so there is one critical point and it is the global minimum. In fact $Q_D^{\ast }$ tends to infinity at both endpoints $D\to 0$ and $D\to \infty$. The strength-gravity disruption transition is thereforeAt this diameter, $b{Q}_b{D}_w^b=a{Q}_a{D}_w^{-a}$, so $Q_w^{\ast }=[(a+b)/b]Q_a{D}_w^{-a}$. This minimizes the full catastrophic disruption threshold; it need not be where its two contributions are equal.

For equal-density spherical planetesimals, the specific impact energy per unit target mass is $\frac{1}{2}(D_{\mathrm{im}}/D_w{)}^3{v}_{\mathrm{rel}}^2$. Setting this to $Q_w^{\ast }$ gives the catastrophic projectile diameterHere $v_{\mathrm{rel}}$ is the impact speed, the threshold is defined per target mass, and energy-coupling changes are not included. If projectile and target densities differ, multiply the diameter by $({\rho }_{\mathrm{target}}/{\rho }_{\mathrm{projectile}}{)}^{1/3}$. For $a=1/2$, $b=3/2$,so collecting the powers givesThe result presumes that this projectile diameter is available in the population.

For the dust, $n(D)dD$ is the total number in a diameter interval, not a number per unit volume. The mass normalization of a power-law size distribution with dust mass $M/2$ givesThe approximation uses $(D_{\min }/D_{\max }{)}^{4-\alpha }\ll 1$. Merely having a small diameter ratio does not guarantee this accuracy if $\alpha$ approaches four, so retain the exact denominator in that case.

With $\gamma >0$, the stated Poynting–Robertson drag law integrates directly:Thus dust arrives in order of increasing diameter. This is burst transport under inverse-size radial drag, assuming an instantaneous release at $r_0$ and no intervening destruction or changes of diameter.

Dust that reaches the planet's orbital radius has several possible fates. It can pass without a close encounter and continue spiralling towards the star; undergo resonant trapping of dust in an exterior mean-motion resonance, where resonant torque temporarily balances the drag; be scattered onto a different stellar Kepler orbit or ejected; or collide with the planet and be accreted. Resonant trapping of dust can increase orbital eccentricity until escape or an encounter becomes possible. Temporary planet-bound capture can become permanent if drag or another dissipative process removes sufficient energy. Arrival at a radial circle alone does not guarantee a collision with the planet at the correct longitude.

A more massive planet has a larger Hill sphere, stronger resonant perturbations and stronger planetary scattering. Capture into a mean-motion resonance is favoured when drift across it is slow compared with its libration timescale. Since $|\dot{r}|\propto D^{-1}$ here, larger grains drift more slowly and are generally easier to trap; small fast-drifting grains more readily cross a resonance. Grain radiation pressure can also move resonance locations, though the prescribed nearly circular drift model omits that additional change. For direct impact, a larger physical radius increases the geometric cross-section, and gravitational focusing gives, in the two-body encounter approximation,A compact massive planet can nevertheless scatter or eject grains rather than hit them: its encounter region can be much larger than its physical radius. Slower dust migration permits more repeated encounter opportunities. These competing effects mean neither planet mass nor grain diameter alone specifies an accretion probability.

For the idealized complete-accretion calculation, put $T=t_p(D_{\max })$ and $\eta =D_{\min }/D_{\max }$. At time $t$, the arriving diameter is $D_t=D_{\max }t/T$ and $d{D}_t/dt=D_{\max }/T$. Hence the dust-burst accretion pulse follows from the arriving mass per diameter:Dropping the small lower-cutoff normalization correction gives the requested expressionThe rate is zero before the smallest grains arrive and after the largest arrive. Integrating the exact expression over its support returns $M/2$, as required by mass conservation. The approximate expression has that integral to the same lower-cutoff accuracy. This pulse assumes no resonant delay or scattering loss before accretion.

For the surface density of a size-sorted dust burst, take $0<r<r_0$ and writeAt fixed $t>0$, the size-radius map is monotone and $d{D}_r/dr=2r{D}_r/H_r$. Define $\Sigma (r,t)$ as the azimuthally averaged surface density, so the mass in an annulus is $2\pi r\Sigma dr$. Transforming the power-law size distribution into radial mass givesIt follows thatThe requested leading approximation omits $1-{\eta }^{4-\alpha }$ from the denominator. The surface density is zero where the inferred $D_r$ lies outside the original size interval. Equivalently, one can obtain it from the instantaneous inward mass flux divided by $2\pi r|\dot{r}(D_r)|$. At $t=0$ the idealized initial distribution is a delta-function ring at $r_0$; the formula above describes the subsequent radial spread. No uniform initial azimuthal distribution is required for this azimuthal average, but the local density of a clump would require its angular evolution. If the complete-absorption planet assumption is retained, this expression holds exterior to $r_p$ and $\Sigma =0$ interior to the planet.

The first relation is vertical hydrostatic equilibrium. In a thin disc, the central object's gravity near the midplane is $g_z=-{\Omega }^{2}z$; the pressure gradient balances it. The second is an equation of state combining ideal gas pressure with isotropic radiation pressure. Write $\mathcal{R}=k/({\mu }_m{m}_p)$, so $p_g=\mathcal{R}\rho T$ and $p_r=4\sigma T^4/(3c)$; $\sigma$ is the Stefan-Boltzmann constant and ${\mu }_m$ the mean molecular weight in proton-mass units.

The third is local energy balance: the increase of upward radiative flux with height equals viscous dissipation. Here $\mu$ is dynamic viscosity, not mean molecular weight. A Keplerian disc has shear $rd\Omega /dr=-3\Omega /2$, so its heating rate is $\mu (r{\Omega }^{\prime }{)}^2=9\mu {\Omega }^2/4$. The fourth is optically thick radiative diffusion, with opacity per unit mass $\kappa$ and diffusion directed down the temperature gradient. Midplane symmetry gives $F(0)=0$ and oppositely directed fluxes from the two faces. The four relations express vertical force balance, gas–radiation thermodynamics, local heating, and diffusive cooling.

Let $\Psi$ be the gravitational potential, with ${\nabla }^2\Psi =4\pi G\rho$, and let ${\varpi }^2=x^2+y^2$ measure distance from the rotation axis. Uniform solid-body rotation contributes the centrifugal potential $-{\Omega }^2{\varpi }^2/2$. In the corotating frame, static force balance therefore gives the uniformly rotating stellar hydrostatic equilibrium equationsThe final coefficient follows because ${\nabla }^2(x^2+y^2)=4$.