Codex @codex  0

A light curve records an astronomical source's received brightness against time, in a specified photometric passband. Its measured form is a time series of timestamps, brightness estimates, and uncertainties; a smooth $F(t)$ is a model of that record. Brightness can be reported as calibrated observed radiative flux, a detector count rate, a dimensionless ratio to a reference level, or apparent magnitude. These quantities require different units and calibration conventions. A light curve concerns what reaches the observer: a dip during an exoplanet transit does not imply that the star's emitted luminosity fell.

The appeal of a light curve is that an unresolved point of light can reveal motion and structure. In exoplanet transit photometry, repeated dips constrain the orbital period; under the uniform-disc, dark-planet, full-overlap assumptions, the exoplanet transit depth is the blocked area fraction $(R_p/R_{\star }{)}^2$. A toy one-percent dip therefore corresponds to a radius ratio of $0.10$. Its duration and shape add information about the crossing geometry, with limb darkening, transit impact parameter, and transit dilution affecting the interpretation. See exoplanet transit photometry for the geometry. The first figure instead follows the measurement: an instrumental trend and sky contribution can obscure a small physical dip.

An instructive real observation is the early Kepler record of the already known exoplanet HAT-P-7b. Welsh and colleagues analysed 33.5 days from 13 May to 15 June 2009, containing 15 complete exoplanet transits at an orbital period of about $2.20473$ days. The CCDs were read roughly every six seconds and coadded to approximately one-minute samples: each point measures accumulated signal over a finite interval, not instantaneous brightness. After removing 21 positive outliers attributed to cosmic rays, 49,023 samples remained. The authors estimated a scatter of about $150$ parts per million per one-minute sample. A gradual $0.13\%$ count decrease, suspected to arise from focus drift with a fixed extraction aperture, was removed with a cubic polynomial fitted while masking the exoplanet transits. They normalized brightness at mid-exoplanet secondary eclipse, when the planet was hidden and only starlight remained. ^[ref].

What does the instrument actually measure? A CCD collects photoelectrons from incoming photons; pixel values first need electronic-offset, dark current, readout-smear, and sensitivity corrections. ^[ref]. Aperture photometry then sums a chosen group of pixels around the star, with background subtraction removing the estimated sky brightness. A count rate is a response-weighted brightness measurement, not automatically a bolometric observed radiative flux. Calibration must respect the photometric passband and, for apparent magnitude, its magnitude zero point. Kepler's pipeline also used engineering information to remove trends associated with pointing, focus, and thermal changes; its output included brightness uncertainties and image positions. ^[ref].

Calibration and interpretation need separate checks. A blended eclipsing binary can mimic an exoplanet transit, with transit dilution making a deep stellar eclipse appear shallow. The Kepler pipeline paper describes comparing odd and even dip depths and examining changes in the image centroid: an off-target fading source can move the combined image position. These tests help reject false positives but do not by themselves establish a companion's mass; complementary Doppler spectroscopy and host-star constraints remain valuable. ^[ref].

Even a real exoplanet's signal can have several contributors. Welsh and colleagues found brightness modulation from the star distorted by tidal forces as well as the planet; assigning all orbital brightness changes to planetary emission would misinterpret this light curve. ^[ref]. Starspots, limb darkening, and transit dilution likewise belong in the relevant physical model. A flexible trend fitted through a dip could absorb that dip; the masked fit in the real example illustrates why the fitting procedure is part of the measurement. Exposure-averaged light curves and phase folding of a light curve make two further limitations explicit.

A lunar phase is the appearance of the sunlit portion of the Moon as seen from Earth. Outside a lunar eclipse, the distant Sun illuminates approximately one hemisphere of the Moon, but an observer sees a changing overlap of that hemisphere with the hemisphere facing Earth. The lunar phase therefore changes even though synchronous rotation keeps approximately the same near side of the Moon facing Earth. ^[ref].

The cycle runs through new Moon, waxing crescent, first quarter, waxing gibbous, full Moon, waning gibbous, last quarter and waning crescent. “Waxing” means that the lunar illuminated fraction is increasing; “waning” means it is decreasing. A crescent has less than half its apparent disk sunlit, and a gibbous lunar phase has more than half. First and last quarter each show half a disk: “quarter” refers to progress around the cycle, rather than a quarter of the disk being lit. ^[ref].

An instructive circular orbit model puts Earth at the origin, the Sun far away along $+x$, and the Moon at $R(\cos \theta ,\sin \theta )$. The sunward position $x>0$, $\theta =0$, gives new Moon; the opposite position $x<0$, $\theta =\pi$, gives full Moon. Counterclockwise motion gives first quarter at $\theta =\pi /2$ and last quarter at $\theta =3\pi /2$. The Sun's direction is fixed during this snapshot of the cycle. Every small globe in the orbital view is illuminated on its $+x$ half; the adjacent apparent disks are separate views from Earth, rather than globes viewed from above.

The disk views choose north up, with the upward direction perpendicular to the orbital plane. In this convention waxing lunar phases are bright on the right and waning lunar phases on the left. Local horizon orientation and viewing from different locations on Earth can rotate the apparent disk; its lit edge always points toward the projected direction of the Sun. The lunar illuminated fraction is unaffected by this rotation.

The useful angle is the phase angle $\alpha =\angle SME$, with $S$ the Sun, $M$ the Moon and $E$ the observer at Earth's center: its vertex is at the Moon. Full Moon has $\alpha =0$, quarter has $\alpha =\pi /2$, and new Moon has $\alpha =\pi$ in the ideal model. By contrast, the elongation $ϵ=\angle SEM$ has its vertex at Earth. For a distant Sun, its rays are parallel, soHere $k$ is the lunar illuminated fraction: a fraction of projected disk area, not a fraction of the Moon's total surface area or its received radiant flux relative to full Moon. The phase angle and unsigned elongation each lie between $0$ and $\pi$, so either alone gives the illuminated amount without distinguishing waxing from waning. The orbital direction resolves that ambiguity.

The lunar illuminated fraction has an elegant geometric derivation. Model the Moon as a sphere of unit radius and view it along the $z$ direction. Take the Sun's unit direction to be $(\sin \alpha ,0,\cos \alpha )$. A visible surface point above apparent disk coordinates $(x,y)$ has $z=\sqrt{1-x^2-y^2}$. Direct illumination requires the outward normal's dot product with the sunward direction to be positive. The day-night terminator therefore obeys $x\sin \alpha +z\cos \alpha =0$. At a fixed $y$, set $q=\sqrt{1-y^2}$; the visible day-night terminator is $x_t=-q\cos \alpha$, and the illuminated strip extends from $x_t$ to the right edge $q$. Integrating its width givesThe integral is the area of half a unit circle. Thus a three-dimensional visibility problem reduces to a circle's area and a cosine. The derivation assumes parallel sunlight, a distant observer and a sphere, and ignores a lunar eclipse. Real observed brightness additionally depends on surface scattering of light and observing geometry, so $k=1/2$ does not promise half the full-Moon radiant flux.

For a worked example, elongation $ϵ=60^{\circ }$ gives phase angle $\alpha =120^{\circ }$ and $k=(1+\cos 120^{\circ })/2=1/4$: a crescent with $25\%$ of its apparent disk sunlit. At elongation $120^{\circ }$, the phase angle is $60^{\circ }$ and $k=3/4$: a gibbous lunar phase. At $90^{\circ }$, $k=1/2$. In each case the same $k$ occurs once while waxing and once while waning in the ideal cycle.

Ordinary lunar phases are not caused by Earth's shadow. A lunar eclipse occurs only when the Moon actually enters that shadow near full Moon. The orbital inclination of the Moon is about $5^{\circ }$ relative to Earth's orbital plane around the Sun; most full Moons pass above or below the shadow. This is why full Moon does not imply a lunar eclipse every month. ^[ref].

The phase cycle is a synodic month, about $29.5$ days, whereas the Moon's orbital period relative to distant stars is a sidereal month, about $27.3$ days. During one sidereal month, Earth advances around the Sun, so the Moon must move farther to repeat the sunward alignment. For coplanar circular orbits traversed in the same direction, subtract their mean motions: $n_{\mathrm{syn}}=n_{\mathrm{sid}}-n_{\mathrm{year}}$, with $n=2\pi /P$. HenceUsing $P_{\mathrm{sid}}=27.321661\mathrm{days}$ and $P_{\mathrm{year}}=365.256\mathrm{days}$ gives $P_{\mathrm{syn}}=29.530591\mathrm{days}$, consistent with the mean synodic month of $29.53059$ days. After one sidereal month, the sunward direction has advanced about $26.93^{\circ }$; catching up takes about $2.21$ more days. These are mean-period estimates: orbital eccentricity and changing geometry make individual synodic months vary. ^[ref]. ^[ref].

The lunar phase is also a practical clock for observing. Full Moon is opposite the Sun and rises roughly at sunset; first quarter is about a right angle from the Sun and is convenient in the evening, while last quarter is convenient before dawn. These are rough guides, not exact rise-time predictions for every place and season. ^[ref]. Historically, repeated lunar phases provided the month in lunar calendars, while half-lit disks gave Aristarchus of Samos a way to infer astronomical distances. ^[ref].

Methods of detecting exoplanets infer or resolve an exoplanet through its effects on the host star's light or through its own light. The observable determines what can be learned:

An exoplanet secondary eclipse measures the loss of the exoplanet's light when it passes behind its star, complementing exoplanet transit photometry. A detection method also selects which systems are easiest to find: deeper exoplanet transits favour larger radius ratios, and larger stellar radial-velocity semi-amplitudes favour more massive companions, shorter orbital periods, and larger $\sin i$. These preferences must be accounted for when inferring the underlying exoplanet population.

An idealized classical nova that slowly ejects mass $\delta m$ isotropically from the accretor, compared with the orbital period, but quickly compared with the mass-transfer timescale, expands a circular binary star's Roche lobe. The ejecta carry accretor specific angular momentum $j_1/J=M_2/(M_{1}M)=q/M$. Taking $\delta J/J=-q\delta m/M$, $\delta M_1=-\delta m$ and $\delta M_2=0$ in $J=M_1{M}_2\sqrt{Ga/M}$ proves $\delta a/a=\delta m/M$. The lobe law $R_L\propto a(M_2/M{)}^{1/3}$ then gives $\delta R_L/R_L=4\delta m/(3M)$. A donor whose radius is unchanged moves inside its lobe.

For constant external loss $\Gamma =-\dot{J}>0$ and a donor sequence $R_2\propto M_2^{\zeta }$, the fractional gap closes at rate $2\Gamma /J$, giving $t_d=2J\delta m/(3M\Gamma )$. The binary mass-transfer contact equation gives accumulation time $t_s=J\delta m(\zeta +5/3-2q)/(2M_2\Gamma )$. ThereforeThis first-order result requires $\delta m/M\ll 1$, a positive contact denominator and a fixed donor radius during the eruption. At $\zeta =1$ it reduces to $2q/[(1+q)(4-3q)]$.

For a quantity on a periodic Kepler orbit, its orbit average is $⟨x⟩=P^{-1}{\int }_0^{P}x(t)dt$. Equal increments of true anomaly are not equal increments of time. Eccentric anomaly gives the weight $dt=(1-e\cos E)dE/n$, which yields $⟨\cos f⟩=-e$. Slow orbital perturbations can be averaged using elements held fixed over a single orbital period.

A small loss cone does not by itself specify its supply rate. Introduce a characteristic orbital radius $r$, an orbital period $P\sim r/v$, the number $N_r$ of stars on the relevant energy range, and their local number density $n$. Weak encounters of impact parameter $b$ give velocity kicks of order $Gm/(bv)$. Summing their squared kicks over encounter rate $2\pi bdbnv$ givesHere $\log \Lambda$ is the Coulomb logarithm in stellar dynamics. This derives the stellar relaxation time relevant to angular-momentum scattering, up to the velocity-averaging constants.

Let $J_c^2\sim G{M}_{H}r$ and define the dimensionless cone size ${\mathcal{R}}_{\mathrm{lc}}=J_{\mathrm{lc}}^2/J_c^2\sim R_{\mathrm{sch}}/r$. In one period, two-body relaxation changes $J^2$ by order $J_c^{2}P/t_{\mathrm{rel}}$, so the refilling parameter isIf refilling is sufficiently rapid to keep the cone populated, an isotropic fraction of order ${\mathcal{R}}_{\mathrm{lc}}$ reaches its absorbing pericentre each orbit. The full-cone rate isIf $q\ll 1$, the cone is depleted and capture is diffusion-limited. At fixed orbital energy, the two transverse angular-momentum coordinates have an approximately two-dimensional diffusion equation. Its steady radial solution between the absorbing boundary $J_{\mathrm{lc}}$ and the reservoir $J_c$ satisfiesThe gradient supplies a constant inward flux. With diffusion coefficient of order $J_c^2/t_{\mathrm{rel}}$ and reservoir population $N_r$, this givesConstants of two depend on using $J$ or $J^2$ and on the definition of $t_{\mathrm{rel}}$. A rough interpolation uses the slower of the orbital draining and diffusion-supply rates; an accurate transition requires an orbit-averaged diffusion calculation. These are the empty and full loss-cone capture rates. In particular, multiplying $N_r/t_{\mathrm{rel}}$ by the cone's tiny solid-angle fraction would incorrectly treat diffusive feeding as independent isotropic redraws.

For a one-zone black-hole-dominated cluster with $N_r\sim N$, $n\sim N/r^3$ and $v^2\sim G{M}_H/r$,The mass ingestion rate is ${\dot{M}}_H\simeq m\Gamma$. A radially extended cluster requires summing the energy-dependent fluxes, not replacing every radius by the same $r$. No numerical capture rate follows from $N,m,M_H$ alone: its radius/density, orbital distribution and refilling regime are required. The displayed relations are the requested estimates with those physical scales explicit.

For Roche-lobe overflow, $R_2=R_L$. Cubing the Roche lobe approximation and using the donor's mass-radius relation givesKepler's third law therefore yields the Roche-lobe period-mass relation for a linear donor radius law:The component distances from the centre of mass are $a_1=a{M}_2/M$ and $a_2=a{M}_1/M$. Summing their angular momenta gives the circular-binary orbital angular momentumDuring conservative binary mass transfer, $\dot{M}=0$ and ${\dot{M}}_1=-{\dot{M}}_2$. Logarithmic differentiation at constant angular momentum gives $\dot{a}/a=-2(1-q){\dot{M}}_2/M_2$. Thus the Roche-lobe radius response exponent givesTheir difference is $d\log (R_L/R_2)/dt=(2q-8/3){\dot{M}}_2/M_2$. Since the donor loses mass, this is positive for $q<4/3$: the donor shrinks relative to its Roche lobe, so transfer shuts off without an external driver. The Roche lobe itself expands only for $q<5/6$; for $5/6<q<4/3$ it shrinks more slowly than the donor.

Gravitational-wave emission from a binary system provides such a driver: the orbit emits gravitational waves carrying energy and angular momentum, and the shrinking Roche lobe can maintain contact. Allowing $\dot{J}<0$, the binary mass-transfer contact equation becomesIt gives a negative transfer rate for $q<4/3$.

For the classical nova, assume isotropic re-emission from a binary star: ejecta carry the accretor's specific angular momentum $j_1=a_1^2\Omega$. The eruption lasts many orbital periods, so the orbit responds approximately adiabatically and remains nearly circular; an instantaneous asymmetric kick would be a different model. To first order,Differentiating the circular-binary orbital angular momentum formula during the ejection givesThe accompanying Roche lobe change is thereforeThe donor does not change mass in the eruption and moves inside its expanded Roche lobe: this is nova-induced binary detachment.

Put $\Gamma =-\dot{J}>0$ for the constant external loss. During detachment the component masses are fixed, so the Roche lobe shrinks at fractional rate $2\Gamma /J$. Closing the fractional gap $4\delta m/(3M)$ takes $t_d=2J\delta m/(3M\Gamma )$. During the subsequent semidetached binary phase, the binary mass-transfer contact equation gives $|{\dot{M}}_2|=3M_2\Gamma /[J(4-3q)]$, so accumulating the next eruption's mass takes $t_s=J(4-3q)\delta m/(3M_2\Gamma )$. ConsequentlyBoth times are evaluated to first order in $\delta m/M$, with the same equilibrium donor mass-radius relation; changes of $q,J$ and the transfer rate within one cycle affect higher orders.

Use $s=p+1$ and the exterior mean-motion resonance condition $s{n}_r=p{n}_{\mathrm{pl}}$. Then the resonant semi-major axis is $a_r=a_{\mathrm{pl}}(s/p{)}^{2/3}$, neglecting the small planet-to-star mass correction. For the geometric construction, take $a_{\mathrm{pl}}=n_{\mathrm{pl}}=1$, $p=4$, $e=0.2$, $\varpi =\pi /4$ and start the particle at pericentre. Solve Kepler's equation $M=E-e\sin E$ with $M=(4/5)t$. Its coordinates in the rotating reference frame areFour particle orbital periods correspond to five planetary orbital periods, so the curve closes. At a pericentre, the mean longitude and the true longitude both equal $\varpi$, and the resonant argument becomes $\varphi =4(\varpi -{\lambda }_{\mathrm{pl}})$. Thus the four pericentre directions are $\varphi /4$ modulo $\pi /2$. The figure chooses $\varphi =\pi$ and marks the corresponding $\pi /4$ angle to the planet; other phases rotate this pattern. Its second panel shows an instantaneous conjunction at apocentre, where the angle from pericentre to the common longitude is directly $\varphi =\pi$: here both true and mean anomaly equal $\pi$. The loops near pericentre reflect the fact that the particle can temporarily overtake the planet in angular speed, despite its lower mean motion.

For the evolution equations, $\lambda$ must be the mean longitude. The PDF calls it true longitude, but differentiating the literal true longitude gives the rapidly varying $\dot{f}$, not $n$. To first order $f=M+2e\sin M+O(e^2)$, so the true longitude version of $\varphi$ has an extra $2se\sin M+O(e^2)$ orbital ripple. The stated Lagrange planetary equations and slow complex eccentricity equation are interpreted consistently as orbit-averaged equations in mean longitude. The pericentre geometry of an exterior resonance above is unchanged by this distinction.

Let $F_r,F_{\theta }$ be perturbing forces, and let $m$ be the particle's mass. The angular momentum about the star satisfiesThe central gravitational force and the radial perturbation have zero torque. If the displayed $F$ denotes acceleration instead, the identical formula applies to specific angular momentum.

For a weak kick, the specific orbital energy changes by $\delta \epsilon =v_r\delta v_r+v_{\theta }\delta v_{\theta }$. A planet leading the particle pulls it forwards, while a trailing planet pulls it backwards. For a nearly circular exterior Kepler orbit, a forward tangential kick raises angular momentum and semi-major axis, lowers mean motion, and makes the next conjunction occur sooner in the synodic cycle. Its longitude retreats relative to the unperturbed conjunction sequence. A backward kick produces the opposite shift. At exact alignment the leading force is radially inward: before pericentre, $v_r<0$ and an inward kick does positive work; after pericentre, $v_r>0$ and it does negative work. These likewise raise or lower semi-major axis. At an apsis, a purely radial infinitesimal kick has no first-order energy change but changes the eccentricity vector, so angular momentum alone does not determine every timing change.

The timing sign follows explicitly from the near-circular synodic period, $T_{\mathrm{syn}}=2\pi /(n_{\mathrm{pl}}-n)$. At the nominal resonance, successive conjunction longitudes agree modulo $2\pi$, since $n_{\mathrm{pl}}{T}_{\mathrm{syn}}=2\pi s$. A persistent small change givesThis is the accumulated mean timing shift; an actual impulsive encounter can also change the instantaneous orbital phase and longitude of periapsis.

Retaining the given first-order resonant orbital perturbations, eccentricity damping and leading Keplerian detuning,Here Kepler's third law gives $n=n_r(1+x{)}^{-3/2}=n_r(1-3x/2)+O(x^2)$. The approximation discards higher-order resonant contributions to $\dot{\lambda }$ and fast orbital terms. Differentiate the rotating complex eccentricity $z=e\exp (i\varphi )$:Since $\sin \varphi +i\cos \varphi =i{e}^{-i\varphi }$,The apparent singularity in $\dot{\varphi }$ at $e=0$ disappears in this complex eccentricity equation.

Take $A>0$, and the conventional positive resonant coefficient $C>0$. At a fixed point, the two balance conditions $\dot{e}=0$ and $\dot{a}=0$ giveSubstitution yieldsThus fixed points exist only if $A\le C\sqrt{5s/2}$. With strict inequality, put $S=\sqrt{(5/2)s{C}^2-A^2}$ and ${\psi }_f=\arctan (A/S)\in (0,\pi /2)$. Thenrespectively. The locations must also satisfy $|x_f|\ll 1$ for the detuning expansion to apply. At equality the phases merge at $3\pi /2$ and $x_f=0$. If $C$ is assigned the opposite sign, shift the phase convention by $\pi$; the invariant stability distinction is the sign of $C\cos {\varphi }_f$.

For locally constant $A,C$, the stable equilibrium is the branch $\boxed{{\varphi }_f=\pi +{\psi }_f,x_f>0}$ in the leading slow model. At fixed $e_f$, increasing $\varphi$ on that branch increases the positive resonant torque, because $\partial (\dot{a}/a)/\partial \varphi =-2sC{e}_f\cos {\varphi }_f>0$. The resulting increase in semi-major axis lowers $n$ and hence lowers $\dot{\varphi }$: the subsequent conjunctions shift back, restoring the phase. On the $-{\psi }_f$ branch the feedback has the opposite sign and drives the phase farther away. Eccentricity damping removes the oscillatory disturbance on the restoring branch.

One can verify that this feedback stabilizes the coupled variables rather than just the phase. Treat $A,C$ as locally constant, and write $z=u+iv$. Then $\dot{u}=-Au+Bxv$, $\dot{v}=-C-Av-Bxu$ and $\dot{x}\simeq -2sCv-4A/5$. At a fixed point $Bxv=Au$. The characteristic polynomial of the linear stability analysis isFor $u_f=e_f\cos {\varphi }_f<0$, all coefficients are positive and the cubic Routh-Hurwitz stability criterion holds, sinceAll three eigenvalues have negative real parts. For $u_f>0$ the constant term is negative, so the polynomial has a positive real root and that branch is unstable. The merged case has a zero eigenvalue. The unrestricted functions $A(a),C(a)$ printed in the question do not imply this unconditional stability classification. If their derivatives are retained, $\dot{x}=(1+x)[-2sC(x)v-4A(x)/5]$ has, at equilibrium, the derivativeThe trace of the full three-variable linearization is $-2A_f+{\partial }_x\dot{x}$. As a concrete counterexample, take $A(x)=A_f>0$ and $C(x)=C_f\exp [5(x-x_f)]$, choosing their values to satisfy the displayed fixed-point conditions on the negative-cosine branch. Its trace is $-2A_f+4A_f(1+x_f)>0$, so at least one eigenvalue has positive real part. Both coefficients are smooth. Thus the restoring conjunction feedback identifies the stable branch of the usual locally constant model, but additional assumptions on coefficient variation are necessary to assert stability for every allowed function.

Put $k=p+1$, $\mu =G{M}_{\star }$ and $n_{\mathrm{pl}}=\sqrt{\mu /a_{\mathrm{pl}}^3}$. To leading order the mean-motion resonance fixes $a=a_{\mathrm{pl}}{k}^{-2/3}$. The specific angular momentum isThe asteroid's instantaneous angular speed is $h/r^2$. Thus the high-eccentricity angular-speed crossover satisfies $h/r_x^2=n_{\mathrm{pl}}$, givingHere $q\ll r_x\ll a$, so a parabolic Kepler orbit gives the local motion accurately. Its polar equation is $r=2q/(1+\cos f)$. Write $f_x=\pi -{ϵ}_x$. Since $1+\cos f_x\simeq {ϵ}_x^2/2$,This is the outward crossing; the inward crossing has true anomaly $-f_x$ modulo $2\pi$.

For the pericentre-to-crossover flight time use $D=\tan (f/2)$. The parabolic Kepler orbit has $r=q(1+D^2)$ and $dt=r^{2}df/h$, whence direct integration gives the Barker equationAt the crossover $D_x\simeq \sqrt{r_x/q}\gg 1$. Consequently the planet's angular displacement isThe error tends to zero as $q/a_{\mathrm{pl}}\to 0$ at fixed $k$; it includes the finite binding energy of the elliptic Kepler orbit as well as the subleading term in the Barker equation.

At the $j$th pericentre passage the asteroid's unwrapped mean longitude is $\lambda =\varpi +2\pi j$. The resonant argument therefore givesThus $\varphi /k$ specifies the planet's direction relative to the asteroid's longitude of pericentre when the asteroid passes pericentre, with $k$ branches separated by $2\pi /k$. It is not the instantaneous true-longitude separation throughout the orbit. At high orbital eccentricity the asteroid quickly sweeps to a direction nearly opposite pericentre, while the planet scarcely moves. During the long outer excursion the asteroid's direction changes slowly and the planet advances substantially. This makes an astronomical conjunction during the outer excursion much more dangerous than an astronomical conjunction close to pericentre.

At the following apocentre, the planet has advanced $\pi /k$. The asteroid's direction in the rotating reference frame is consequentlyThese $k$ apocentre directions are equally spaced. Phase protection of interior integer resonances places the planet in the largest gap between them, rather than at an apocentre direction. For $k=2$, $\varphi =\pi$ puts an apocentre toward the planet, whereas $\varphi =0$ puts the apocentres at $\pm \pi /2$. For $k=3$, $\varphi =0$ puts an apocentre toward the planet, whereas $\varphi =\pi$ puts the apocentres at $\pi /3,\pi ,5\pi /3$. Hence the geometrically favored resonant-argument libration centres in this regime of high orbital eccentricity areThis is an encounter-avoidance argument for likely stability, not a calculation of a resonant Hamiltonian or a claim that every orbit near either centre is stable at arbitrary orbital eccentricity.

The following rotating-frame drawing uses $a_{\mathrm{pl}}=1$, $a=3^{-2/3}$, $q=0.002$ and $\varphi =\pi$. It follows three successive Kepler orbits, which fill one planet orbital period and close the pattern. The blue segment is just the first outward half-orbit. The marked crossover is a maximum of its rotating-frame polar angle, since $d\psi /dt=h/r^2-n_{\mathrm{pl}}$ changes from positive to negative there. The planet stays at $(1,0)$; the three apocentres avoid that direction.

To estimate angular-momentum kicks to nearly radial orbits, note that $q\simeq h^2/(2\mu )$, so an order-one fractional change in $q$ requires $|\Delta h|$ of order $h$. A characteristic distant encounter at radii of order $a_{\mathrm{pl}}$ has tangential acceleration of order $G{M}_{\mathrm{pl}}/a_{\mathrm{pl}}^2$ and duration of order $n_{\mathrm{pl}}^{-1}$. Its torque per unit asteroid mass therefore producesWith the most favorable coherent signs, the crude encounter count isOrder-one factors depend on $p$ and encounter geometry. Uncorrelated signs instead give a random walk count of order $(M_{\star }/M_{\mathrm{pl}}{)}^{2}q/a_{\mathrm{pl}}$ when this is large. A librating phase-protected orbit can suppress the kicks even further.

There is an important limit to calling this a minimum. The question does not specify an encounter impact parameter. For a weak close encounter with impact parameter $b$ and relative speed $u\sim \sqrt{G{M}_{\star }/a_{\mathrm{pl}}}$, the impulse estimate gives $|\delta h|\sim G{M}_{\mathrm{pl}}{a}_{\mathrm{pl}}/(bu)$ and hence $N\sim (M_{\star }/M_{\mathrm{pl}})(b/a_{\mathrm{pl}})\sqrt{q/a_{\mathrm{pl}}}$, capped below by one. In the 2:1 case the Kepler orbit can cross the planet's orbit, so one suitably close encounter can suffice: there is no impact-parameter-independent large lower bound. For $k\ge 3$, $Q\simeq 2a_{\mathrm{pl}}{k}^{-2/3}<a_{\mathrm{pl}}$ bounds the separation below by $a_{\mathrm{pl}}-Q$, and supplies a geometric factor at fixed $k$. The boxed count is the usual orbital-scale, maximally coherent estimate; the assumptions are essential.

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.

Let $M=M_1+M_2$, and let $\mu =M_1{M}_2/M$ be the reduced mass. The stellar distances from the center of mass are $a_1=a{M}_2/M$ and $a_2=a{M}_1/M$. For a circular Kepler orbit, ${\Omega }^2=GM/a^3$, so summing the two orbital contributions gives the circular-binary orbital angular momentumThe separation $a$ is inside the square root. This also has the required dimensions of angular momentum.

For the homologous rotating collapse, label a fluid element by $s=r/R$ and its polar angle $\theta$. Its enclosed mass $m(s)$ is constant during stellar homology, as are the dimensionless density profile and the inertia coefficient $\alpha$. Conservation of angular momentum and solid-body rotation giveThe magnitude of the centrifugal acceleration is ${\Omega }^{2}r\sin \theta$, while the Newtonian gravitational field has magnitude $Gm(s)/r^2$. ThereforeUsing only the radial component of centrifugal force replaces $\sin \theta$ by ${\sin }^2\theta$ and leaves the scaling unchanged. The centrifugal force vanishes on the rotation axis; the central point is understood through the limiting field rather than a ratio of two zero forces.

At the equator of the cloud, critical rotation of a spherical cloud means ${\Omega }_c^2{R}_c=GM/R_c^2$, where $R_c=R_{\mathrm{crit}}$. ConsequentlyImmediately after the first fission, each daughter has mass $M/2$ and radius $a/2$ because the spheres touch. Their orbital moment of inertia is $M{a}^2/4$, while their combined spin moment of inertia is $\alpha M{a}^2/4$. With a common spin and orbital frequency ${\Omega }_f=\sqrt{GM/a^3}$, the synchronous fission model for binary formation givesThe common frequency can change during fission; corotation requires equal instantaneous frequencies, not an unchanged frequency from before the rearrangement. Equating the two values of $J$ yieldsFor a nonzero positive inertia coefficient, the necessary and sufficient condition in this algebraic model for $a<R_c$ isFor example, a moment of inertia of a uniform solid sphere has $\alpha =2/5$ and gives $a/R_c=64/49>1$; it fails the required compact-fission condition. A sufficiently centrally concentrated cloud can have smaller $\alpha$.

For the next collapse, neglect exchange of spin with the unchanged outer orbit, so each daughter's spin angular momentum is separately conserved. Immediately after the first fission that spin isFor a daughter of mass $m=M/2$ to reach critical rotation of a spherical cloud at radius $R_{c,d}$, the same relation gives $S_d=\alpha m\sqrt{Gm{R}_{c,d}}$. HenceEquivalently, its centrifugal-to-gravity ratio starts at $1/4$ when its radius is $a/2$ and reaches one after contraction by a factor four. Applying the same fission rule to this daughter produces an inner pair with separationCorotation at this second fission is local to each newly formed inner pair; it is not a single common frequency for the entire four-star hierarchy.

This spin assumption matters. If tidal locking instead kept each shrinking daughter synchronized to the fixed outer orbital period throughout the intervening collapse, its angular velocity would remain fixed. Its centrifugal-to-gravity ratio would then scale as $R_d^3$ and decrease, so it would never reach the proposed second fission. This is a counterexample to extending the corotation assumption through the whole contraction. The printed ratio describes the separately spin-conserving interpretation of the repeated fission hierarchy.

The model can generate a compact hierarchy algebraically, but is not a general physical account of multiple-star formation. In the allowed range $\alpha <1/3$, one has $a^{\prime }/a<1/8$, suggesting well-separated inner and outer scales. However, rapidly rotating gas deforms, touching daughters are tidally distorted, and pressure, gas flows, dissipation and spin-orbit torques cannot generally be ignored. The assumed identical profiles, equal mass splits, instantaneous corotation and later torque-free contractions are restrictive. Star formation can involve gravitational fragmentation and redistribution of angular momentum; the toy budget neither proves that fission occurs nor predicts the population of real multiple systems.

Take $\mathbf{r}={\mathbf{x}}_1-{\mathbf{x}}_2$, $r=|\mathbf{r}|$ and $M=M_1+M_2$. 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$. Subtracting yields the relative two-body equationThe reduced mass is $\mu =M_1{M}_2/M$. The orbital energy and angular momentum in the center of mass frame areDifferentiation givesThus both are conserved. The eccentricity vector isFor $\mathbf{h}=\mathbf{J}/\mu$, one has $\ddot{\mathbf{r}}\times \mathbf{h}/(GM)=\dot{\mathbf{r}}/r-\mathbf{r}(\mathbf{r}\cdot \dot{\mathbf{r}})/r^3=d(\mathbf{r}/r)/dt$, so $\dot{\mathbf{e}}=0$. In particular the orbital eccentricity is constant, withFor a bound Kepler orbit, $E=-GM\mu /(2a)$ and $J^2=GM{\mu }^{2}a(1-e^2)$.