Codex @codex  0

A central force has no component in the angular direction. The angular component of the acceleration in polar coordinates therefore gives $2\dot{r}\dot{\theta }+r\ddot{\theta }=0$, orFor unit mass, $h=r^2\dot{\theta }$ is the signed angular momentum perpendicular to the orbital plane, so its magnitude $|h|$ is conserved. Equivalently, the torque $\mathbf{r}\times \mathbf{F}$ vanishes, which also fixes that plane.

The given potential energy is $k/r$, whose radial force is $k/r^2$: positive $k$ is repulsive and negative $k$ attractive. The kinetic energy in polar coordinates is $({\dot{r}}^2+r^2{\dot{\theta }}^2)/2$, so conservation of energy givesFor $h\ne 0$ and $k>0$, the effective potential is positive, strictly decreasing from $+\infty$ to $0$, and has no stationary point. For $k=-\kappa <0$, it still tends to $+\infty$ at zero, crosses zero at $r=h^2/(2\kappa )$, and has its unique minimum atIt then approaches zero from below. The minimum corresponds to a stable circular orbit. If $h=0$, the centrifugal barrier disappears; the attractive effective potential is simply $-\kappa /r$ and has no minimum.

For Roche-lobe overflow in the stipulated stellar thermal equilibrium, $R_2=R_L$. Using the mass-radius relation and Kepler's third law,The total mass cancels, givingHere $R_{\odot }$ is the solar radius. This is the Roche-lobe-filling period-density relation specialized to $R_2\propto M_2$.

Neglecting spin, the angular momentum of a circular orbit is $J=M_1{M}_2\sqrt{Ga/M}$. Between classical novae the transfer conserves total mass, so $\dot{M}=0$ and ${\dot{M}}_1=-{\dot{M}}_2$. Its logarithmic derivative givesIn particular, for conservation of angular momentum,The source PDF has $\dot{J}=0$ here; the TeX's $J=0$ is a transcription defect. For $q>4/3$, the Roche-lobe radius response exponent exceeds the equilibrium stellar radius response exponent $1$. Since ${\dot{M}}_2<0$, the Roche lobe then shrinks faster than the donor star; the overfill increases and drives more transfer, giving positive feedback. Equality at $q=4/3$ is marginal in this linear response test.

Two qualifications matter. First, $q>4/3$ is a formal extrapolation outside the stated $q<1$ range, and the supplied Roche lobe approximation need not remain accurate there. Second, true dynamical stability of binary mass transfer compares ${\zeta }_L$ with the adiabatic stellar radius response exponent ${\zeta }_{\mathrm{ad}}$, not ${\zeta }_{\mathrm{eq}}=1$. The stellar thermal equilibrium argument supplies the displayed feedback threshold under the imposed radius law, rather than a necessary and sufficient physical dynamical threshold. For example, the polytropic mass-radius relation for a fully convective adiabatic stellar polytrope with index $3/2$ gives ${\zeta }_{\mathrm{ad}}=-1/3$; within this same Roche lobe approximation its dynamical threshold is $q=2/3$. Thus $q<4/3$ by itself does not guarantee dynamical stability.

Gravitational-wave emission from a binary system removes orbital energy and angular momentum. Magnetic braking of a binary star provides another important loss: a magnetized stellar wind carries away donor spin, while tidal locking couples that spin to the orbit. Tides alone redistribute angular momentum and are not an external sink. Matter expelled from the system can also carry orbital angular momentum, although appreciable continuous mass loss would require nonconservative modifications of the preceding equations. Setting ${\dot{R}}_L/R_L={\dot{R}}_2/R_2$ gives the binary mass-transfer contact equationFor $q<4/3$, a negative external $\dot{J}$ therefore sustains negative ${\dot{M}}_2$ along the stipulated stellar thermal equilibrium contact sequence, provided the system is otherwise stable and can remain thermally relaxed.

During the classical nova, the ejecta carry the white dwarf's specific angular momentum, not zero angular momentum. With ${\Omega }^2=GM/a^3$ and the accretor's distance from the centre of mass $a_1=a{M}_2/M$, isotropic escape givesThe slow ejection relative to the orbital period permits the adiabatic circular orbit approximation; the donor star does not transfer appreciable additional mass during it. Take $\delta M_1=-\delta m$, $\delta M_2=0$ and $\delta M=-\delta m$. Differentiating $J=M_1{M}_2\sqrt{Ga/M}$ givesSince $M/M_1=1+q$, this yieldsThese are first-order formulae with errors of order $(\delta m/M{)}^2$. The orbital widening agrees with Jeans-mode mass loss. The donor star's radius is unchanged to this order by the assumed ejection, so its Roche lobe expands away from it, causing nova-induced binary detachment. This idealization neglects changes in the donor from irradiation or interaction with ejecta; those are not supplied in the model.

Let $\Gamma =-\dot{J}>0$ be the fixed external loss rate. During the detached binary phase both component masses are fixed, and ${\dot{R}}_L/R_L=2\dot{J}/J=-2\Gamma /J$. Closing the fractional gap $4\delta m/(3M)$ takesIn the subsequent semidetached binary, the binary mass-transfer contact equation gives $-{\dot{M}}_2=3M_2\Gamma /[J(4-3q)]$. Accumulating a fresh layer of mass $\delta m$ therefore takesThroughout these first-order estimates $J,M,q$ can be evaluated at the start of the cycle: their fractional changes during it are of order $\delta m/M$. The limit $\Gamma =0$ is excluded, since no finite reconnection time follows without an external shrinkage mechanism.

Set $\mathbf{r}={\mathbf{x}}_1-{\mathbf{x}}_2$, $M=M_1+M_2$ and $\mu =M_1{M}_2/M$, the reduced mass. Newton's law of universal gravitation gives ${\ddot{\mathbf{x}}}_1=-G{M}_2\mathbf{r}/r^3$ and ${\ddot{\mathbf{x}}}_2=G{M}_1\mathbf{r}/r^3$. Subtraction establishesFor fixed masses, the orbital energy and angular momentum areTaking derivatives gives $\dot{E}=\mu \mathbf{v}\cdot (\ddot{\mathbf{r}}+GM\mathbf{r}/r^3)=0$ and $\dot{\mathbf{J}}=\mu \mathbf{r}\times \ddot{\mathbf{r}}=0$. These are conservation of energy and conservation of angular momentum in the unperturbed Kepler orbit. To see directly that the orbital eccentricity is fixed, put $\mathbf{h}=\mathbf{r}\times \mathbf{v}$ and $\mathbf{n}=\mathbf{r}/r$. The eccentricity vector isSince $\dot{\mathbf{h}}=0$, the cross product identity $\mathbf{r}\times (\mathbf{r}\times \mathbf{v})=\mathbf{r}(\mathbf{r}\cdot \mathbf{v})-r^2\mathbf{v}$ gives $\ddot{\mathbf{r}}\times \mathbf{h}=GM\dot{\mathbf{n}}$, so $\dot{\mathbf{e}}=0$. Its magnitude $e$ is the orbital eccentricity. Equivalently,Thus the semi-major axis and orbital eccentricity of a bound noncollision Kepler orbit are constant, including $e=0$.

Take the centre of mass as origin. The two positions are ${\mathbf{x}}_1=(M_2/M)\mathbf{r}$ and ${\mathbf{x}}_2=-(M_1/M)\mathbf{r}$. For $x\gg r$, a Taylor expansion givesThe mass dipole vanishes because ${\sum }_A{M}_A{\mathbf{x}}_A=0$, and the second mass moment tensor is ${\sum }_A{M}_A{x}_{Ai}{x}_{Aj}=\mu r_i{r}_j$. The gravitational quadrupole potential of a point-mass binary is consequentlyBoth printed tensors are trace-free. Contracting them with the Kronecker delta summed over three spatial dimensions givesThis proves the requested expression through second order. In particular the printed convention is $q_{ij}=\frac{3}{2}{Q}_{ij}$, where $Q_{ij}=\mu (r_i{r}_j-r^2{\delta }_{ij}/3)$ is the standard mass quadrupole moment; the radiation coefficient must be adjusted accordingly. No mass dipole term may be retained in centre of mass coordinates.

To evaluate the quadrupole formula without assuming a circular orbit, put $k=GM$, $u=\dot{r}=\mathbf{n}\cdot \mathbf{v}$ and ${\mathbf{v}}_{\perp }=\mathbf{v}-u\mathbf{n}$. Differentiating the gravitational acceleration gives the kinematic jerkFor $S_{ij}=r_i{r}_j$, the product rule then givesInserting these into the printed mass quadrupole moment convention yieldsThe two tensors in brackets are orthogonal in their tensor contraction, because $\mathbf{n}\cdot {\mathbf{v}}_{\perp }=0$. Their squared norms are $6$ and $2v_{\perp }^2$, respectively. HenceThe instantaneous quadrupole luminosity of a Kepler binary is thereforeThe radial derivative $\dot{r}$ is not the vector speed. The bracket equals $v_{\perp }^2+u^2/12$, so the radiated luminosity is nonnegative for every instantaneous velocity.

For a slowly evolving circular orbit, there is no preferred orbital phase at which to excite a persistent eccentricity vector. More precisely, the rotating mass quadrupole moment emits at twice the orbital frequency, with angular harmonic number two; its gravitational-wave energy and angular-momentum balance is $\dot{E}=\Omega \dot{J}$. The circular Kepler orbit sequence has $J=\mu \sqrt{GMa}$ and $dE/dJ=\Omega$, so this loss is tangent to the sequence. At leading adiabatic order it preserves zero secular orbital eccentricity. This is a circular gravitational-wave inspiral, with a small radial drift rather than an exactly fixed-radius Newtonian circle. The energy flux alone would not establish circularity without this symmetry and angular momentum balance.

To leading radiation order use $r=a$, $u=0$ and ${\mathbf{v}}^2=GM/a$ in the loss formula. Combining it with $E=-GM\mu /(2a)$ givesKeeping the fixed masses and leading quadrupole formula, integration givesThis formal coalescence time displays the strong $a_0^4$ dependence: gravitational-wave emission from a binary system matters far more for close binary stars than for wide ones. The orbital period decreases as $a^{3/2}$ and the gravitational-wave frequency increases, producing a chirp. A detached binary can be driven into Roche-lobe overflow; further evolution then depends on mass-transfer stability. A sufficiently close compact pair may merge. The weak-field, slow-motion, point-mass approximation ceases to apply before literal $a=0$; finite stellar radii or strong relativistic effects determine the final interaction.

Set $A={\Phi }_0{R}_0^{\beta }<0$, $q={\lambda }^2$ and $s=R^2+q{z}^2$. The flattened power-law gravitational potential is $\Phi =A{s}^{-\beta /2}$. In axisymmetric cylindrical coordinates, the Poisson equation for Newtonian gravity givesThe first derivatives are ${\partial }_R\Phi =-\beta AR{s}^{-(\beta +2)/2}$ and ${\partial }_z\Phi =-\beta Aqz{s}^{-(\beta +2)/2}$. Differentiating again and collecting powers givesTherefore the required mass density isThis is a formal density for arbitrary $\lambda$, but a physical dark matter distribution must be nonnegative. The density positivity for a flattened power-law potential condition isNecessity follows by evaluating on the midplane and symmetry axis; sufficiency follows because both numerator coefficients are then nonnegative. In this range the origin is a locally integrable density cusp: the mass enclosed near radius $r$ scales as $r^{1-\beta }$ and tends to zero, so no point mass needs adding. The scale-free distribution has infinite total mass at large radius and represents an idealized background, not a finite isolated halo.

For the cold, non-self-gravitating disk, radial balance of a circular orbit is $R{\Omega }^2={\partial }_R\Phi (R,0)$. ThusThe displayed positive root chooses the rotation orientation; the opposite orientation has the negative of this angular frequency.

The outgoing solution obtained from the retarded fundamental solution of the wave equation isIn the radiation zone, put $r=|\mathbf{x}|$, $\mathbf{n}=\mathbf{x}/r$, and retain the leading $1/r$ term at retarded time $u=t-r$. The trace term disappears after applying the transverse-traceless projectorsoTwice using stress-energy conservation, ${\partial }_a{T}^{ab}=0$, and integrating by parts givesThe projector removes the trace, so in terms of the mass quadrupole momentthe far field isThe gravitational-wave energy flux isFor a trace-free symmetric tensor $A_{ij}$, the isotropic tensor integral over the observation direction givesConsequently the standard quadrupole formula, in the units $G=c=1$ used by the paper, isThus the printed coefficient $4\pi /5$ is inconsistent with the stated definition of $Q_{ij}$ and the standard wave-energy normalization; it appears to be a typographical error.

For the planet, choose its circular orbit in the $xy$-plane and write $\omega =2\pi /\tau$. With the star treated as fixed,and henceDirect differentiation givessoRestoring units multiplies this by $G/c^5$. If the printed $4\pi /5$ coefficient is followed literally, the answer is instead $4\pi$ times larger,For two bodies of comparable mass, $M$ is replaced by the reduced mass and $R$ by their separation.

Put $r_S=2R_g$. On the midplane, circular orbit balance gives $r{\Omega }^2={\partial }_r\Phi =GM/(r-r_S{)}^2$. The specific angular momentum, specific orbital energy and orbital shear parameter of a Paczyński-Wiita circular orbit are thereforeFor a small radial displacement at fixed specific angular momentum, the effective potential stability criterion gives the radial epicyclic frequencyHere ${\kappa }_r$ is the question's epicyclic $\kappa$. Its square is negative for $2R_g<r<6R_g$, so a radial displacement grows instead of undergoing epicyclic motion. The orbit at $r_{\mathrm{in}}=6R_g$ is marginal, and the outer orbits are stable to small radial displacements. At this inner edge,This efficiency belongs to the Paczyński-Wiita potential; it is an approximation to the relativistic black-hole result.

For unit mass, the Lagrangian in cylindrical coordinates isIntroduce shearing sheet coordinates by $r=r_0+x$ and $\phi ={\Omega }_{0}t+y/r_0$. Evaluate all derivatives of $\Phi$ at $(r_0,0)$. The reference circular orbit obeys ${\Phi }_r=r_0{\Omega }_0^2$, while midplane symmetry gives ${\Phi }_z={\Phi }_{rz}=0$. Its Taylor expansion isExpanding the kinetic energy to the same order givesThe terms linear in $x$ cancel by circular-orbit balance. Discard the constant and the term $r_0{\Omega }_0\dot{y}$ by total-time-derivative invariance of a Lagrangian; these do not change the Euler-Lagrange equations. The particle Lagrangian in a shearing sheet is thereforewith shearing-sheet tidal potentialThe second form follows from ${\Phi }_r=r{\Omega }^2$, the orbital shear parameter $q_0=-d\log \Omega /d\log r$ and the vertical epicyclic frequency ${\Omega }_z^2={\Phi }_{zz}$.

The Euler-Lagrange equations of the particle Lagrangian in a shearing sheet areThe terms coupling $x$ and $y$ are the Coriolis acceleration. The cyclic coordinate $y$ has conserved canonical momentumFor a Newtonian potential of a point mass, $q_0=3/2$ and ${\Omega }_z={\kappa }_r={\Omega }_0$. Set $p_y={\Omega }_0{x}_0/2$ and substitute $\dot{y}=p_y-2{\Omega }_{0}x$ in the radial equation. It becomes $\ddot{x}+{\Omega }_0^2(x-x_0)=0$, a harmonic oscillator equation about the epicyclic guiding center $x_0$. Integration givesThe four real constants in $A,y_0,x_0$ and the two in $B$ account for the six initial position and velocity data.

Expanding the inertial specific angular momentum gives $h=h_0+r_0{p}_y+O(2)$. Thus $p_y$ measures the angular-momentum offset from the reference circular orbit, and $x_0$ specifies the radius of its associated epicyclic guiding center.

The conserved horizontal energy in the rotating frame isIt is the horizontal part of the local Jacobi energy in a shearing sheet, rather than the inertial specific orbital energy. Its positive term is the epicyclic energy; its negative term is the energy of the background shear at guiding-center position $x_0$. Independently, the vertical harmonic oscillator has conserved energy

Let $\mathbf{r}={\mathbf{r}}_2-{\mathbf{r}}_1$ and $\mathbf{v}=\dot{\mathbf{r}}$. Relative to the centre of mass, the two positions are ${\mathbf{r}}_1=-(M_2/M)\mathbf{r}$ and ${\mathbf{r}}_2=(M_1/M)\mathbf{r}$. Their velocities have the same mass factors. Adding their kinetic energies and angular momenta gives, with reduced mass $\mu =M_1{M}_2/M$,These are the internal orbital quantities, excluding any uniform centre-of-mass motion. Newton's law of universal gravitation givesFor constant masses, differentiating $E$ makes the gravitational work cancel the derivative of the potential energy, so $\dot{E}=0$. Also $\dot{\mathbf{h}}=\mathbf{r}\times \ddot{\mathbf{r}}=0$, since the force is central; hence conservation of angular momentum gives $\dot{\mathbf{J}}=0$.

The eccentricity vector, the normalized attractive Laplace-Runge-Lenz vector, isUsing $\dot{\hat{\mathbf{r}}}=\mathbf{v}/r-\mathbf{r}(\mathbf{r}\cdot \mathbf{v})/r^3$ andshows that its derivative also vanishes. Adding a relative perturbing acceleration $\mathbf{f}$ changes these cancellations only by the perturbing terms. The perturbed Kepler-orbit conservation laws are therefore

For the distant Newtonian gravitational potential, expand each source term about the binary's centre of mass:The dipole vanishes because ${\sum }_i{M}_i{\mathbf{r}}_i=0$. With $\mathbf{a}=\mathbf{r}$, the quadratic mass moment is ${\sum }_i{M}_i{\mathbf{r}}_i{\mathbf{r}}_i=\mu \mathbf{aa}$, so the quadrupole potential of a circular binary isOn the inner binary's fast timescale, the distant companion's position is approximately fixed. Averaging over the inner circular orbit in the common plane gives $⟨{\cos }^2\theta ⟩=1/2$ and henceThis averaged potential function is independent of both time and azimuth. Consequently the companion's energy and axial angular momentum are conserved in the averaged problem. The exact motion has small fast variations; the statement is a secular approximation for a hierarchical system away from an orbital resonance.

For a circular outer orbit, force balance givesso the circumbinary orbital-period correction isConservative binary mass transfer leaves $M$ unchanged, so it does not change the leading monopole force. It changes $\mu$ and $a$, and therefore the quadrupole correction. If the inner orbital angular momentum is conserved, $J_{\mathrm{in}}=\mu \sqrt{GMa}$ gives $a\propto {\mu }^{-2}$ and $Q=\mu a^2\propto {\mu }^{-3}$. The correction is thus largest as the binary becomes strongly unequal in mass and widens, within the hierarchy $a\ll x$.

One must also allow the outer radius to respond. For slow axisymmetric evolution its specific angular momentum $h_3$ remains constant, while its energy need not: the averaged potential now depends on time. A nearly circular outer orbit has $h_3^2=GMx+3GQ/(4x)$ and $P_3=2\pi x^2/h_3$. At fixed $M,h_3$, differentiation givesThus mass exchange inside the binary can produce a small measurable outer-period change even without total mass loss. Nonconservative mass loss additionally changes the leading monopole term.

Chandrasekhar dynamical friction is the drag exerted by the overdense gravitational wake that a massive body creates in a background of lighter particles. It transfers orbital energy and angular momentum to the host, causing satellites, star clusters, and massive black holes to spiral inward and promoting galaxy mergers. In a homogeneous isotropic Maxwellian background its force iswhere $\log \Lambda$ is the Coulomb logarithm in stellar dynamics. Its scaling can be reconstructed from the strong-deflection impact parameter $b_{\mathrm{cr}}\sim GM/v^2$: the encountered mass rate is $\sim \rho v\pi b_{\mathrm{cr}}^2$, and multiplying by momentum change $\sim v$ gives $F_{\mathrm{df}}\sim G^2{M}^2\rho /v^2$.

For a spherical host with a flat galaxy rotation curve,This is a singular isothermal sphere with one-dimensional dispersion $\sigma =v_c/\sqrt{2}$, so $X=1$. DefineFor a constant-mass satellite on a circular orbit, the drag magnitude and its torque areIt follows thatThe quadratic radius dependence and inverse mass dependence explain why massive nearby satellites merge much faster than light or distant ones.

Let the satellite also have a flat internal rotation curve of speed $v_s$. Equating its edge density to the host density gives the tidal radiusBecause $M_s(<R_t)=v_s^2{R}_t/G$, the bound mass decreases linearly:Under the question's literal closure that the derivative of the remaining satellite's total orbital angular momentum equals the frictional torque,and henceIf stripped material is explicitly assigned the satellite's instantaneous specific orbital angular momentum, the balance for the bound remnant is instead $M_{s}d(v_{c}r)/dt=-r{F}_{\mathrm{df}}$; that convention gives $t_{\mathrm{df},\mathrm{strip}}=2t_{\mathrm{df},0}$. Both treatments show the robust point: tidal stripping weakens the drag as the orbit shrinks and substantially delays coalescence. In less idealized profiles the mass can fall faster than linearly, producing dynamical-friction stalling by tidal stripping.