Codex @codex  0

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.

For the relative position $\mathbf{r}={\mathbf{r}}_1-{\mathbf{r}}_2$, Newton's law of universal gravitation gives the two-body problem in its reduced one-body form,Taking the dot product with the relative velocity $\mathbf{v}=\dot{\mathbf{r}}$ givesIntegration therefore yields conservation of specific orbital energy:

For a spherical grain of diameter $D$, density ${\rho }_d$, and radiation-pressure efficiency $Q_{\mathrm{pr}}$, comparison of the stellar momentum flux with Newton's law of universal gravitation gives the radiation-pressure coefficientFor $D$ much larger than the characteristic stellar wavelength, $Q_{\mathrm{pr}}\simeq 1$, so $\beta \propto D^{-1}$. When $D$ becomes comparable to the optical wavelength, diffraction and the grain's composition make $Q_{\mathrm{pr}}$ size dependent and $\beta$ reaches a broad maximum. Deep in the Rayleigh scattering regime, absorption can give $Q_{\mathrm{pr}}\propto D$ and hence nearly constant $\beta$, while scattering alone gives $Q_{\mathrm{pr}}\propto D^4$ and hence $\beta \propto D^3$. A realistic curve therefore rises roughly as $D^{-1}$ toward micron sizes, turns over near the stellar spectral peak, and falls or flattens for still smaller grains.

For the force magnitude in Newton's law of universal gravitation, the total differential isDividing by the positive value $F$ gives the particularly useful relative formFor finite small changes this gives the first-order approximation to the fractional change; higher-order products of increments are omitted.

Powered vertical flight. Take upward velocity $v$ as positive, and let $v_e>0$ be the downward exhaust speed relative to the rocket. During $dt$ the rocket's mass changes by $dm<0$, and exhaust mass $-dm$ leaves with Earth-frame velocity $v-v_e$. Compare the momentum of the rocket plus this newly expelled exhaust before and after the interval:Discard the second-order product $dmdv$. The result is $mdv+v_{e}dm=-mgdt$, orThe positive rocket thrust is $-v_e\dot{m}$. Integrating with $v(0)=0$ and $m(0)=m_0$ gives the rocket equation with gravitational loss:Using $d(mv)/dt$ for the rocket alone without the exhaust momentum flux would give the wrong equation for this open system.

Orbital speed under the stated gravity assumption. A circular orbit of radius $r_o=R_E+h_o=7.90\times 10^6\mathrm{m}$ requires centripetal acceleration $v_o^2/r_o$. If the magnitude of gravitational acceleration is the stipulated constant $g=9.8\mathrm{m}{\mathrm{s}}^{-2}$ even at the orbit,This treats gravity as centrally directed for the orbital calculation, while retaining its constant magnitude. The subsequent flat-surface approximation treats its direction as locally vertical.

Burnout height. During the unpowered projectile motion, horizontal velocity is constant. At the top of the trajectory the vertical component vanishes, so the horizontal component must equal the orbital speed $v_o$. At burnout, therefore,These follow either by integrating $\ddot{y}=-g$ or by the vertical work-energy theorem. Substituting the preceding constant-$g$ orbital result givesThe vertical component is initially positive, as required to reach the orbit after burnout.

Alternative inverse-square interpretation. A real Earth's exterior Newtonian gravitational field is not constant over a $1500\mathrm{km}$ altitude change. If $9.8\mathrm{m}{\mathrm{s}}^{-2}$ instead means the surface value $g_0$, then Newton's law of universal gravitation gives $GM=g_0{R}_E^2$ andUsing this orbital speed but the paper's constant-$g$ flat projectile motion approximation for the coast would giveThis alternative changes the global constant-gravity assumption; it is not the same model as the $183\mathrm{km}$ result. An exact curved Kepler orbit coast would also abandon the prescribed flat trajectory. The printed assumptions support the constant-magnitude answers above, while the inverse-square values explain the different answer obtained from the usual physical Earth-orbit formula.

Field of discrete masses. Newton's law of universal gravitation and superposition give the gravitational field away from the source points:It is the force per unit test mass, with units $\mathrm{N}\mathrm{k}{\mathrm{g}}^{-1}=\mathrm{m}{\mathrm{s}}^{-2}$, and points toward each positive source mass.

Infinite circular cylindrical rod. Let $b$ be the radius of the circular cross-section, so $A=\pi b^2$, and write $\lambda =\rho A$ for the linear mass density. Choose a coaxial cylindrical Gaussian surface with radius $r>b$ and length $L$. Translation and rotation symmetry make $\mathbf{g}=g_r(r)\hat{\mathbf{r}}$. The end caps have zero field flux because their normals are axial; on the curved surface the outward-normal component is the constant negative value $g_r$. Gauss's law for gravity therefore givesThe gravitational field of an infinite uniform cylinder inside the material instead has enclosed mass $\rho \pi r^{2}L$, giving $g_r=-2\pi G\rho r$. The exterior result assumes a circular cylindrical rod; area by itself would not specify the exact near field of a noncircular cross-section.

In the left panel the dashed circle is a transverse section of the coaxial Gaussian surface; the field points inward and the surface normals outward. Its side area is $2\pi rL$.

For an exterior circular orbit, Newton's second law equates the force required for centripetal acceleration to gravity:The orbital speed is independent of radius and of the orbiting object's mass, for orbits outside the infinite rod.

Finite rod. Taking the supplied field as the gravitational field of a finite uniform line segment, the corresponding balance at radius $R$ isSince $M=2a\rho A=2a\lambda$, its two limits areThe first reproduces the infinite-rod result: a nearby observer sees a very long source. The second is the familiar circular orbit speed about a point mass, since a distant observer cannot resolve the rod's length. The dimensionless curve isIt starts at a constant plateau and decreases smoothly, approaching an $R^{-1/2}$ dependence far away. For $R>0$, $dv/dR=-\sqrt{GM}R/[2(R^2+a^2{)}^{5/4}]<0$, consistent with the sketch.

Derivation of the finite-rod field. Put the source along $-a<z<a$ on the $z$ axis. In its central perpendicular plane, an element $dm=\lambda dz$ is at distance $s=\sqrt{r^2+z^2}$. Its field magnitude is $G\lambda dz/s^2$, and its inward radial component has the extra factor $r/s$. The axial components cancel between $z$ and $-z$. ThusReplacing $2a\lambda$ by $M$ givesThe integral treats the source as a line mass. For a finite cylinder of nonzero cross-section, the displayed finite-rod formula is a thin-rod approximation, not an exact field for all $r$. Its near-rod limit should therefore be read as $b\ll R\ll a$, keeping the orbit outside the material; the apparent line-source singularity at $r=0$ does not describe the interior of a real cylinder.

For a light test particle of mass $m$ outside a spherical source of mass $M$, Newton's law of universal gravitation gives $\mathbf{F}=-GMm\hat{\mathbf{r}}/r^2$. It is a central force, so its torque about the source vanishes and angular momentum $\mathbf{L}=m\mathbf{r}\times \mathbf{v}$ is conserved. The motion therefore lies in a fixed plane when $L\ne 0$. Its gravitational potential energy is $-GMm/r$, and conservation of mechanical energy givesHere $ϵ$ is specific orbital energy and $h$ is specific angular momentum. A Kepler orbit has polar equationThe angle is measured from periapsis. The conic section is an ellipse for $ϵ<0$ and $0\le e<1$, a parabola for $ϵ=0$ and $e=1$, and a hyperbola for $ϵ>0$ and $e>1$. The circle is the $e=0$ elliptical case. The special case $h=0$ is radial motion and must be handled as a degenerate limit. Bound ellipses have the source at one focus, which is Kepler's first law.

For a bound orbit of semi-major axis $a$, $ϵ=-GM/(2a)$. Combining this with the energy equation gives the vis-viva equation:The particle therefore moves faster near periapsis and slower near apoapsis. Conservation of angular momentum also gives constant areal velocity, $dA/dt=h/2$, expressing Kepler's second law. Using the ellipse's area and $h^2=GMa(1-e^2)$ gives the orbital periodwhich is Kepler's third law.

A circular orbit follows from centripetal acceleration $v_c^2/r=GM/r^2$, giving $v_c=\sqrt{GM/r}$ and $T=2\pi \sqrt{r^3/(GM)}$. The minimum escape velocity from radius $r$ follows from the marginally unbound condition $ϵ=0$:These ideal orbits neglect drag force, other gravitating bodies and changes in the central source. For two comparable masses, the relative two-body problem uses $G(M+m)$ instead of $GM$. An orbiting object is continuously in free fall; the inward gravitational force changes its direction even when its speed is constant. These inverse-square results assume the physical variation of gravity with radius, rather than extending a near-surface constant-$g$ approximation to an entire orbit.