Codex @codex  0

The central equation must be read as $\ddot{\mathbf{r}}=-\mu \mathbf{r}/r^3$; its printed left-hand expression lacks $=0$. Initially $\mu$ is constant. Dotting this equation with the velocity givesTaking its cross product with $\mathbf{r}$ gives $d(\mathbf{r}\times \mathbf{v})/dt=0$. In the fixed orbital plane the magnitude is $h=r^2\dot{\theta }$. These are conservation of specific orbital energy and specific angular momentum.

Let $f=\theta -\varpi$. Differentiating the polar equation of the Kepler orbit, with its elements initially constant, givesSubstitution into the specific orbital energy yields

For isotropic stellar mass loss with no recoil or drag, the force is still central, so $\mathbf{h}$ remains exactly constant, although $C$ does not. NowThe same algebraic relation between $C,h,e$ holds for the instantaneous osculating orbital elements. Differentiating it and using $h^2/(\mu r)=1+e\cos f$ givesand, for $e>0$,A particularly useful nonsingular description uses the eccentricity vectorThe vector equation remains meaningful at a circular orbit, where a longitude of pericentre is undefined. Its components along and perpendicular to $\mathbf{e}$ give the scalar formula above and

The osculating pericentre is $q=h^2/[\mu (1+e)]$. ThereforeFor mass loss, $\dot{\mu }<0$, this is nonnegative: pericentre does not decrease. It is positive except at pericentre itself, where $\cos f=1$ and $\dot{q}=0$. Thus a strictly opposite sign at every instant is not literally true; a pericentre passage with nonzero mass-loss rate is a counterexample to the strict wording. The intended monotonicity follows exactly, without requiring slow mass loss.

For adiabatic orbital expansion under isotropic mass loss, average over the unperturbed Kepler orbit and hold $\dot{\mu }/\mu$ constant to leading order during that orbit. If $E$ is the eccentric anomaly and $n=\sqrt{\mu /a^3}$, thenIt follows immediately thatAveraging with uniform true anomaly instead of uniform time would give the wrong result. Since $\mu a(1-e^2)=h^2$, conserved $h$ and constant secular orbital eccentricity giveThe pericentre and apocentre expand in the same secular proportion. These are leading adiabatic invariants, rather than exact invariants for arbitrary time-dependent $\mu$.

For the instantaneous precession fraction write $\theta =f+\varpi$ and $\dot{\theta }=h/r^2$. The adiabatic apsidal condition is controlled byIt is a signed fraction: apsidal advance contributes positively on the outward leg during mass loss, and negatively on the inward leg. Let ${\tau }_{\mu }=|\mu /\dot{\mu }|$. Ordinary orbit averaging requires ${\tau }_{\mu }\gg P$, together with negligible variation of the mass-loss law over one orbit. If the pericentre direction is also to vary negligibly relative to the instantaneous azimuthal motion, the maximum absolute fraction must be small. A precise condition for $0<e<1$ isThe maximizing cosine is $c=(1-\sqrt{1+8e^2})/(2e)$: differentiating the last factor gives $c+2e-e{c}^2=0$. Thus one can evaluate $F(e)$ explicitly by putting $\sqrt{1-c^2}/(1+ec{)}^2$ into the formula. Its useful limits areFor orbital eccentricity of order unity, mass loss much slower than an orbital period suffices; near a circular orbit, keeping the pericentre direction slowly varying additionally requires $n{\tau }_{\mu }e\gg 1$. The latter divergence is a singularity of the pericentre coordinate, not a divergence of the eccentricity vector dynamics. Indeed, at frozen $\dot{\mu }/\mu$, integrating the scalar evolution to first order gives $\delta e=-(\dot{\mu }/\mu )(1-e^2)\sin E/n$ up to an integration constant. Its absolute amplitude stays small for slow mass loss even when the scalar longitude of pericentre becomes unsuitable.

At each fixed semimajor axis, let $B=b_{3/2}^2(\alpha )/b_{3/2}^1(\alpha )$ and $Z_f=B{e}_{\mathrm{pl}}$. The planet's longitude of pericentre is the zero of longitude. The secular perturbation equation becomes $\dot{z}=iA(z-Z_f)$. Integrating it with the initial condition $z(0)=0$ givesThe forced eccentricity vector is the constant real $Z_f$; the free eccentricity vector is $-Z_f{e}^{iAt}$. Consequently the complex eccentricity follows a circle in the Argand plane centered at $Z_f$, with radius $Z_f$, starting at the origin and moving counterclockwise. Its initial velocity is downward, $\dot{z}(0)=-iA{Z}_f$. The secular period is $2\pi /A$, andThe longitude of pericentre is undefined at the origin, but the complex eccentricity passes smoothly through it. At $At=\pi$, $z=2Z_f$ and the orbit is maximally eccentric and aligned with the planet.

Using the small-$\alpha$ Laplace coefficient expansion givesAlso $n=n_{\mathrm{pl}}{\alpha }^{-3/2}$ to leading order in the planet-to-star mass ratio. The secular forcing by a distant outer planet is thereforeA relative $O({\alpha }^2)$ correction to $A$ can accumulate a phase error at very long times; the leading-frequency formula is not a uniform-in-time expansion.

If $a_2=a_1(1+\delta )$ with $\delta =\delta a/a_1\ll 1$, then $Z_{f,2}/Z_{f,1}=1+\delta +O({\delta }^2)$ and $A_2/A_1=1+3\delta /2+O({\delta }^2)$. The outer planetesimal therefore has a slightly larger forced eccentricity circle and a shorter secular period. Its accumulated secular phase lead is $\Delta (At)\simeq 3A_{1}t\delta /2$. Initially both circles start at the origin, but their eccentricity vectors progressively lose alignment. Linearizing the exponential in $\delta$ additionally requires $|A_{1}t\delta |\ll 1$.

The orbit crossing is caused by this accumulating differential secular phase, even though individual orbital eccentricities remain small. To first order in orbital eccentricity, the radius of an orbit at fixed inertial longitude $\theta$ isFor an infinitesimally neighboring orbit,This derivative is evaluated within the first-order orbital-shape model; derivatives of the neglected quadratic terms give relative $O(Z_f)$ corrections near crossing. The least separation over longitude is consequently $\delta a[1-|z+a{\partial }_{a}z|]$. The differential secular orbit-crossing criterion is $|z+a{\partial }_{a}z|=1$. Since $Z_f\propto a$ and $A\propto a^{3/2}$,The first term stays bounded by $4Z_f\ll 1$, whereas the second grows in proportion to time. To leading order in $Z_f$, crossing begins when $(3/2)Z_{f}At=1$. ThusMore precisely, in this first-order orbital-shape approximation the exact first crossing lies within a fractional $O(Z_f)$ window of this value, since the magnitude differs from $(3/2)Z_{f}At$ by at most $4Z_f$. Higher-order Laplace coefficients give additional corrections. The local derivation takes the neighboring-orbit limit before the long-time limit; for a finite separation it requires $\delta a/a_1\ll Z_f$ near crossing. A broad continuous disk contains such neighboring orbits. The crossings first appear in directions selected by the phase of $z+a{\partial }_{a}z$, not simultaneously at every longitude.

For the collision energy estimate adopt the phase-averaged mean orbital eccentricity $\overline{e}=4Z_f/\pi \simeq 5\alpha e_{\mathrm{pl}}/\pi$. The specified velocity estimate givesUsing the projectile kinetic energy per target mass convention for specific impact energy, equal-size and equal-density bodies have equal masses and $Q=\frac{1}{2}{v}_{\mathrm{col}}^2$. A catastrophic planetesimal collision requires $Q\gtrsim Q_D^{\ast }$, soThe numerical coefficient is only illustrative because the collision speed was specified only to order of magnitude. If the disruption threshold is instead defined using the centre of mass kinetic energy per combined mass, equal masses give $Q_R=v_{\mathrm{col}}^2/8$ and the coefficient becomes $8{\pi }^2/25$ for a threshold expressed in that convention. The robust result is a lower bound of order $Q_D^{\ast }{a}_{\mathrm{pl}}/(G{M}_{\star }{e}_{\mathrm{pl}}^2)$.

Within the inner-disk approximation, catastrophic collisions require sufficiently large $a_1/a_{\mathrm{pl}}$, not arbitrarily small radii. Although the Keplerian speed rises inward, the secularly induced orbital eccentricity falls faster: $v_{\mathrm{col}}^2\propto \alpha$. The lower bound must be much smaller than one to leave a domain compatible with $\alpha \ll 1$. A small planet orbital eccentricity, a large disruption threshold, or a distant planet can eliminate that domain. The planet mass controls the time to crossing but cancels from the eventual forced eccentricity amplitude and this collision-energy criterion. The threshold alone does not ensure that crossing has occurred within the disk's available lifetime.

With zero seed, both Bäcklund transformation equations concern ${\varphi }_a/2$. On a nonconstant branch, separation of variables uses $d\log |\tan ({\varphi }_a/4)|=d{\varphi }_a/[2\sin ({\varphi }_a/2)]$, givingFor $C>0$, absorb its magnitude into an additive constant $c=\log C$ in the exponent. A convenient smooth representative isThe constant-phase condition for this traveling profile determines its velocity:The profile is a traveling soliton with velocity $v_a=(1-a^2)/(1+a^2)$, strictly between $-1$ and $1$. The scalar-field vacua on the two sides differ by $2\pi$. With the topological charge convention $Q=[\varphi (+\infty )-\varphi (-\infty )]/(2\pi )$, this branch has $Q=\mathrm{sgn}a$: it is a Sine-Gordon kink for $a>0$ and an antikink for $a<0$. Its width is proportional to ${\gamma }_a^{-1}$, and its derivative decays exponentially away from its center, giving a finite-energy field configuration.

For completeness, the classical rest mass in this paper's coupling convention is $8m/\beta$, as derived in Question 2; a Lorentz boost gives energy $8m{\gamma }_a/\beta$. This localized, topologically protected traveling field is the required classical field-theory soliton. A negative $C$ reverses the field and hence the topological charge; $C=0$ gives the vacuum rather than a soliton. Vacuum shifts by $4\pi$ can be made without changing the displayed Bäcklund transformation equations.

To keep the paper's coupling convention explicit, write $g=\sqrt{\beta }$ and use physical coordinates $X=x/m$, $T=t/m$. Assume $\beta >0$ and weak coupling $g\ll 1$. The canonical real scalar field is $\phi =\varphi /g$, and the physical potential isThe symbol $g$ is the coupling usually appearing in canonical Sine-Gordon theory; the paper's $\beta$ is its square. The vacua have $g\phi =2\pi n$. A static Sine-Gordon kink has $g{\phi }_K=4\arctan e^{m(X-X_0)}$ and satisfies $\frac{1}{2}{\phi }_K^{\prime 2}=U({\phi }_K)$. Since $1-\cos (g{\phi }_K)=2{\mathrm{sech}}^2[m(X-X_0)]$, its classical mass isThe topological sector is important: this energy is measured relative to a vacuum, and the kink joins distinct vacua at the two spatial ends. Expanding around a spatially constant vacuum cannot construct this state by any finite-order perturbation in $g$.

For the general one-loop soliton mass correction, start with a canonical scalar field potential $U$ and a stable static soliton ${\phi }_K$. Write $\phi ={\phi }_K+\eta$. The term linear in $\eta$ vanishes by the classical Euler-Lagrange field equation. The quadratic fluctuation Hamiltonian isChoose a common large box and a common finite-mode regularization in quantum field theory. Expand the nonzero normal modes as $\eta ={\sum }_n{q}_n{f}_n$, with ${\mathcal{H}}_K{f}_n={\omega }_n^2{f}_n$ and normalized eigenfunctions. Each pair $(q_n,p_n)$ is a quantum harmonic oscillator contributing ground-state energy $\hslash {\omega }_n/2$. In the vacuum, replace ${\mathcal{H}}_K$ by ${\mathcal{H}}_0=-{\partial }_X^2+U^{\prime \prime }({\phi }_{\mathrm{vac}})$. Subtract the two ground-state energies and add the local counterterms evaluated on the soliton relative to the vacuum. This derives the general formulaThe prime excludes exact zero modes in field theory from oscillator quantization; their zero frequencies contribute no zero-point energy, but their role in mode counting must not be forgotten. The two sums mean a paired finite regulator, not separate subtractions of divergent answers. Equivalently the nonzero-mode term is the regulated difference of the square-root traces of the two fluctuation operators. Discrete bound states and continuum modes both contribute. This is a vacuum-subtracted soliton mass and the first term in the semiclassical soliton mass expansion.

Translation gives a zero mode in field theory because differentiating the static equation yields ${\mathcal{H}}_K{\phi }_K^{\prime }=0$. A Gaussian oscillator or an unprimed functional determinant is inappropriate along this flat direction. Replace its amplitude by the position collective coordinate $X_0(T)$ and require the residual fluctuation to obey $\int \eta {\phi }_K^{\prime }dX=0$, preventing double counting. The associated change-of-variables Jacobian supplies the zero-mode normalization. At low speed the collective-coordinate effective Lagrangian for a soliton is $-M_{\mathrm{cl}}+\frac{1}{2}{M}_{\mathrm{cl}}{\dot{X}}_0^2+\cdots$; quantizing the position gives the soliton momentum and its translational states, not an extra oscillator rest energy. More generally, every physical continuous modulus needs a collective coordinate; gauge directions require gauge fixing rather than additional physical states.

The Sine-Gordon kink fluctuation operator is particularly simple:while $D{D}^{†}=-{\partial }_X^2+m^2={\mathcal{H}}_0$. This supersymmetric factorization of the one-soliton potential shows stability and generates all nonzero eigenfunctions from vacuum plane waves. The sole normalizable bound state is the translational zero mode of a sine-Gordon kink, proportional to $\mathrm{sech}m(X-X_0)$; there is no positive-frequency internal bound oscillator. The continuum has $\omega (k)=\sqrt{k^2+m^2}$ and no reflection. Applying $D^{†}$ to $e^{ikX}$ gives a transmission amplitudeThis scattering phase shift changes the density of continuum modes. The high-frequency vacuum subtraction cancels the extensive vacuum contribution but still leaves a logarithmic ultraviolet divergence.

The finite part also requires consistent mode-number regularization of soliton masses. The periodic-box phase-shift quantization condition is $k_{n}L+\delta (k_n)=2\pi n$. Match $2N+1$ modes: the kink has its translation mode plus the two continuum modes at each $n=1,\dots ,N$, while the vacuum has the $k=0$ oscillator of frequency $m$ and the corresponding continuum pairs. With $\hslash =1$ and $\Lambda =2\pi N/L$, expanding $k_n-k_n^{(0)}=-\delta (k_n^{(0)})/L+o(L^{-1})$ givesIntegration by parts uses $\delta (0)=\pi$ and yieldsThe cutoff surface term for a Sine-Gordon kink tends to $-m/\pi$. It cannot be dropped merely because $\delta (\Lambda )\to 0$: $\omega (\Lambda )$ grows at the same time.

A renormalization condition must specify which mass and coupling are held fixed. For vacuum normal ordering, or cancellation of the vacuum tadpole diagram with the elementary mass fixed at $m$, the quartic interaction gives the mass countertermEvaluating this Sine-Gordon vacuum tadpole counterterm on the kink givesThe logarithmic ultraviolet divergence cancels, leaving the renormalized one-loop Sine-Gordon kink mass in the stated vacuum scheme:This example illustrates why a zero mode in field theory must be treated as a collective coordinate, why vacuum subtraction alone need not remove ultraviolet divergences, and why the finite relation between the two regulators matters. Different finite counterterms amount to different definitions of the renormalized parameters; an unregulated frequency difference without a renormalization condition is not a physical mass prediction.

Define $k_0=\omega /c_0$ and use the outgoing acoustic square-root branchFor positive real frequency reached from below, $\gamma =i\sqrt{k_0^2-k^2}$ on the propagating interval and is positive real for $|k|>k_0$. The outgoing field in the lower fluid has pressure amplitude $B{e}^{\gamma y}$, since $y<0$. If $\eta =\hat{\eta }{e}^{i\omega t-ikx}$, the shared normal velocity is $i\omega \hat{\eta }$. The linear homentropic acoustic equations therefore giveThis lower-fluid wave is outgoing; no additional incoming sound is included in defining the impedance seen by the upper fluid.

Put $K_s=T{k}^2-m{\omega }^2$. The sheet's force balance gives $K_s\hat{\eta }=-P+B$, hence $P=({\rho }_0{\omega }^2/\gamma -K_s)\hat{\eta }$. Its prescribed downward velocity amplitude is $V=-i\omega \hat{\eta }$. Thus the tensioned-sheet acoustic impedance isThe first term is the lower fluid's normal acoustic impedance, and the second is the sheet's inertial and elastic-sheet tension response. For a real propagating angle, $Z_f={\rho }_0{c}_0/\sin \theta$.

At fixed nonzero frequency and fixed wavenumber, $m\to \infty$ gives $|Z|\to \infty$, a zero-velocity, in-phase reflecting boundary. At fixed $k\ne 0$, $T\to \infty$ gives the same reflection limit, but there is an important exception: elastic-sheet tension does not resist the spatially uniform mode $k=0$. At normal incidence, the impedance remains ${\rho }_0{c}_0+im\omega$ however large the elastic-sheet tension is. With $m=0$ this mode is transparent. By contrast, arbitrarily large mass resists even a spatially uniform oscillation. These fixed-frequency limits exclude a simultaneously tuned structural resonance.

If $m=T=0$, $Z=Z_f$ and $R=0$. The identical fluids are effectively joined across a massless, untensioned interface: pressure and normal velocity continue without reflection. This is the matched case, rather than the pressure-release case.

The area element in plane polar coordinates is $rdrd\theta$. Thus the axial moment of inertia isFor the removed material, the angular width is $2\pi /13$, andIts moment of inertia is therefore $I_{\mathrm{cut}}=\pi {\rho }_0{a}^5/160$, leavingThe torque equation $I\dot{\Omega }=\tau$ gives constant angular acceleration, so the time to reach the prescribed angular speed from rest is

To find the centre of mass after removing material, let ${\mathbf{e}}_1$ point along the bisector of the missing sector in the disc. The original disc has zero first mass moment about its centre. The missing sector has zero transverse first moment by reflection symmetry, while its component along ${\mathbf{e}}_1$ isbecause the radial integral is $11a^4/192$. The remaining first moment is $-Q_{\mathrm{cut}}{\mathbf{e}}_1$. Dividing by the given mass $k{\rho }_0{a}^3$ gives the body-fixed centre of massIt lies away from the missing sector. In the inertial frame, let $t_{\ast }$ be the moment the applied torque stops and let $\psi ={\psi }_{\ast }+\Omega (t-t_{\ast })$ give the orientation of that sector's bisector. Then, with $d=11a\sin (\pi /13)/(96k)$,The centre of mass has inward centripetal acceleration of magnitude ${\Omega }^{2}d$. The required real force is the constraint reaction exerted by the fixed rod and its support. Its horizontal resultant is $k{\rho }_0{a}^3\ddot{\mathbf{R}}$; a zero applied axial torque does not imply a zero resultant force.

First apply momentum conservation to the rocket and the infinitesimal amount of expelled gas. During $dt$, the exhaust mass is $\alpha dt$ and its inertial velocity, to first order, is $v-u$. With no dust, the momentum balance isIt gives the rocket equation $\boxed{m\dot{v}=\alpha u}$, while $\dot{m}=-\alpha$.

For a dust-collecting rocket, the incoming dust mass is $\beta vdt$ with zero initial momentum. The same balance has rocket mass $m+(\beta v-\alpha )dt$ after the interval, soConsequentlyThe second term in the acceleration equation accounts for the momentum required to bring collected dust up to rocket speed.

Set $\lambda =\sqrt{\alpha /(\beta u)}$ and $V=\lambda u$, so $\alpha u=\beta V^2$. As long as the engines operate and $0\le v<V$, the velocity increases and time can be eliminated:Integrating from $v=0$, $m=m_0$, givesThus the required mass–velocity relation is

If $\lambda >1$, this formula would give $m\to 0$ as $v\to V^{-}$. In this interval, $\dot{m}=\beta (v-{\lambda }^{2}u)<0$. A real rocket retains positive body mass even after its fuel is gone, and collected dust only increases that residual mass. It therefore cannot follow this branch all the way to zero mass: its fuel must run out while $v<\lambda u$. Indeed the hypothetical time integral $dt/dv=m/[\beta (V^2-v^2)]$ is integrable at $V$ when $\lambda >1$, so even the formal zero-mass endpoint would occur at finite time.

If $0<\lambda <1$, the same relation instead gives $m\to \infty$ as $v\to V^{-}$. The mass first decreases, reaches its minimum at $v={\lambda }^{2}u$, and then increases because dust arrives faster than fuel is expelled. The time integral now diverges at $V$: under indefinitely maintained thrust, $v$ approaches $\lambda u$ asymptotically while stored dust grows without bound. At that limiting speed, thrust $\alpha u$ balances the dust-loading term $\beta v^2$, and $\dot{m}\to \beta u\lambda (1-\lambda )>0$. This is the physical reason the total rocket mass no longer signals fuel exhaustion. A rocket with a finite initial fuel supply still runs out at finite time $t=M_{\mathrm{fuel}}/\alpha$ and before reaching $V$; increasing total mass does not create new fuel.