Codex @codex  0

For a single-lined circular spectroscopic binary whose visible component has radial-velocity amplitude $K_1$ and orbital period $P$, the binary mass function is

Decreasing orbital period, distinct from spin slowdown. Joseph Hooton Taylor Jr. inferred masses using binary-pulsar periastron advance and binary-pulsar Einstein delay, then tested the independent period derivative against the general relativity quadrupole formula. Correcting differential Galactic acceleration gave indirect gravitational radiation evidence. This was not a terrestrial strain measurement.^[ref]

In particular, the raw measured derivative schematically containsHere $a_{\parallel }^{\mathrm{relative}}$ is positive when the source's recession accelerates. Its contribution follows by differentiating the leading Doppler effect. A decreasing observed orbital period alone therefore does not identify a unique energy-loss mechanism.

For fixed masses, binary orbital energy is $-G{m}_1{m}_2/(2a)$ and Kepler's third law gives $P_b\propto a^{3/2}$. Positive outgoing gravitational radiation predicted by the quadrupole formula makes the orbital energy more negative, reducing $a$ and $P_b$. This sign argument explains shrinking orbits, but a quantitative test also needs the predicted magnitude and orbital eccentricity dependence.

A gradual rotation of the direction of periapsis within the orbital plane. A closed Kepler orbit has fixed argument of periapsis $\omega$; relativistic gravity can make $\omega$ advance between successive orbits. The angle change and binary-pulsar orbital decay describe different changes of the orbit.

A binary star is a gravitationally bound pair orbiting their common centre of mass.

For a circular Kepler orbit, neglecting stellar spin, the orbital angular momentum of a binary star is $J=\mu a^2\Omega =\mu \sqrt{GMa}$, where $M$ is total mass, $\mu$ is the reduced mass, and $a$ is separation. Kepler's third law also gives $J=G^{2/3}{P}^{1/3}{M}_1{M}_2/[(2\pi {)}^{1/3}{M}^{1/3}]$. Distances from the center of mass are $a{M}_2/M$ and $a{M}_1/M$; adding their individual orbital angular momenta proves the formula.

Conservative binary mass transfer moves matter between the components while preserving total binary mass and orbital angular momentum.

An orbit with nonzero orbital eccentricity. A bound Kepler orbit with $0<e<1$ is a noncircular ellipse; $e=1$ and $e>1$ describe the parabolic and hyperbolic unbound cases. A semi-major axis fixed by Kepler's third law need not equal the time average of the orbital radius.

A Hohmann transfer orbit between coplanar circular Kepler orbits of radii $r_1,r_2$ traverses half an elliptic orbit of semi-major axis $a=(r_1+r_2)/2$. Kepler's third law therefore gives $t=\pi \sqrt{a^3/(GM)}$. Relative to the initial orbital period $T_1$, $t=(T_1/2)[(r_1+r_2)/(2r_1){]}^{3/2}$. Radii must be measured from the attracting centre, rather than from its surface.

Kepler's laws describe three connected features of an elliptic orbit: its shape, how fast it is traversed, and how its size determines its orbital period. Kepler's first law gives an ellipse with the attracting centre at one focus of a conic; Kepler's second law gives equal swept areas in equal times; Kepler's third law relates the squared orbital period to the cubed semi-major axis. Together they turn a drawing of a Kepler orbit into a prediction of where the bodies will be and when they will return. ^[ref]

Here the mathematical model is an isolated two-body problem in Newtonian gravity, with point masses $M,m>0$ and a bound, noncolliding elliptic orbit. Let $\mathbf{r}={\mathbf{r}}_m-{\mathbf{r}}_M$ be the relative displacement, $r=|\mathbf{r}|$ the separation, and $\mu =G(M+m)$ the gravitational parameter, where $G$ is the gravitational constant. Newton's law of universal gravitation then givesThe relative Kepler orbit has semi-major axis $a$ and orbital eccentricity $0\le e<1$. Both bodies move around their shared astronomical barycenter, the system's centre of mass. In its rest frame their displacements are ${\mathbf{r}}_m=M\mathbf{r}/(M+m)$ and ${\mathbf{r}}_M=-m\mathbf{r}/(M+m)$. Their individual semi-major axes therefore add to $a$; the relative ellipse's attracting focus of a conic represents zero separation, rather than a body fixed in space.

One example will connect all three laws. After nondimensionalization, choose $a=1$, $e=0.6$ and $\mu =1$. The semi-minor axis is $b=a\sqrt{1-e^2}=0.8$, while the periapsis distance and apoapsis distance areThus $a=(0.4+1.6)/2=1$ measures the ellipse's size even though the separation changes fourfold. The specific angular momentum will be $0.8$, the two endpoint speeds $2$ and $0.5$, and the orbital period $2\pi$. These are values in the chosen scaled units, not measurements of a particular planet.

Johannes Kepler published Kepler's first law and Kepler's second law in Astronomia Nova in 1609 and Kepler's third law in Harmonices Mundi in 1619. These publication dates are documented by ^[ref]. Their appeal is that three short statements organize both geometry and timing; Newtonian gravity supplies a common dynamical explanation.

As a heuristic approximation, inserting $r=6GM/c^2$ into the circular orbit form of Kepler's third law, with $M$ here denoting the total mass, givesThe dominant quadrupole gravitational-wave frequency is twice the orbital frequency, $f_{\mathrm{GW}}=2/T$.

Angular astrometry determines the combination displayed by substituting $a=D\alpha$ into Kepler's third law. Keeping the fitted shape, orientation and orbital period fixed, the rescalingleaves the entire angular time series unchanged. A longer physical orbit around a larger mass can trace precisely the same angular path. More angular samples improve the fitted combination; by themselves they cannot break this scaling symmetry of the model.

The radial-velocity semi-amplitude isFor $\sin i\ne 0$, combining it with the angular fit gives the orbital parallax distanceThe spectral velocity scale supplies the missing physical ruler. For a face-on orbit, $\sin i=0$, this radial signal vanishes and that method cannot supply a distance. In a partial or noisy data set, fitted orientation and systemic motion also correlate with the inferred scale.

The instantaneous distance from an attracting centre to an orbiting body. It is constant on a circular orbit but varies along an elliptic orbit. For Kepler's third law the relevant size is the semi-major axis, which equals the orbital radius only for a circular orbit.

A stable sinusoidal radial-velocity semi-amplitude can be stellar reflex motion caused by a companion on a nearly circular orbit. The amplitude is compatible with a giant planet; a low-mass star or brown dwarf viewed nearly face-on can mimic the same projected motion, but requires an increasingly unlikely orientation. Starspots and magnetically suppressed convection can produce a rotational stellar activity signal on a ten-day timescale; contemporaneous photometry, activity indicators and stellar spectral-line bisector and other profile changes distinguish these from a clean common Doppler shift. Solar-like stellar oscillations ordinarily have much shorter periods, so a coherent ten-day oscillation of an otherwise ordinary solar-type dwarf is less plausible. Instrumental drift, sampling aliases and blended spectra also need checking. A sinusoid alone does not establish a planetary origin, and the data do not assign numerical probabilities to these alternatives.

For the minimum companion mass from radial velocity, take a circular two-body orbit and initially $M_{\star }=M_{\odot }$. The exact relation and its low-mass limit areWith $\sin i=1$,Here $a=[G{M}_{\star }{P}^2/(4{\pi }^2){]}^{1/3}$ is the relative orbital semi-major axis. Nonzero orbital eccentricity multiplies the inferred $M_p\sin i$ by $\sqrt{1-e^2}$ at fixed measured $K,P$, but the sinusoidal model motivates $e=0$ here.

For the radial-velocity sensitivity limit at fixed planet mass, use $K=M_p\sin i\sqrt{G/(M_{\star }a)}$ and interpret the stated detection scale as a semi-amplitude threshold $K_{\min }$. The largest distance occurs for $\sin i=1$:For an inclination $i$, this boundary is smaller by ${\sin }^{2}i$. The result assumes adequate time coverage and neglects the stipulated stellar noise. If the detection scale instead denotes a full peak-to-peak excursion, set $K_{\min }$ to half that scale; $a_{\max }$ then increases fourfold.

The stellar-mass scaling of radial-velocity planet estimates changes these inferences because Kepler's third law and the stellar reflex response depend on the stellar mass. At fixed $K,P$, $M_p\sin i\propto M_{\star }^{2/3}$ and $a\propto M_{\star }^{1/3}$, whereas the fixed-mass detectability boundary scales as $M_{\star }^{-1}$. ThereforeThe observed period and amplitude are unchanged. The stellar radius does not enter these leading dynamical estimates, but is crucial for the exoplanet transit inference.

For the transiting companion, assume an opaque spherical planet, negligible transit dilution and negligible limb darkening. Its radius ratio is $k=R_p/R_{\star }=\sqrt{\delta }=0.03873$, henceUsing the revised stellar mass and a circular orbit gives $a=0.01850\mathrm{au}$ and $A=a/R_{\star }=7.957$. Interpret the duration as the total first-to-fourth-contact duration $T_{14}$. The impact parameter from total transit duration follows from the total duration of a circular planetary transit:This is noncentral but not grazing: $b<1-k$. The small-planet straight-chord approximation gives $b\simeq 0.845$; retaining finite $R_p$ and orbital curvature gives the stated value. A duration measured between other contacts needs the corresponding contact formula. Unknown orbital eccentricity and argument of periapsis change the transverse speed and prevent a unique unconstrained $b$; the numerical estimate is conditional on the stated circular model.^[ref]

Because the transit constrains $\sin i\simeq 0.994$, the RV nondetection now bounds the actual planetary mass, rather than only its projected mass. The bulk-density upper limit from a transit and an RV nondetection isThese bounds assume the instrumental threshold can be applied to a fitted signal at the transit period; no statistical confidence level is specified. Such a large radius at this low mass is inconsistent with an ordinary compact rock-and-iron body. A low-density volatile-rich planet, possibly with an extended H/He envelope, is indicated under the planetary interpretation. The observations do not uniquely determine interior composition, atmospheric abundance or formation history; blended eclipses and activity still require validation. There is no luminosity or effective temperature supplied for the host, so a numerical planetary equilibrium temperature cannot be inferred from these data alone.

The numerical conversions above use externally documented $(GM{)}_{\odot }=1.3271244\times 10^{20}{\mathrm{m}}^3{\mathrm{s}}^{-2}$, $R_{\odot }=6.957\times 10^8\mathrm{m}$, $M_{\oplus }=5.9722\times 10^{24}\mathrm{kg}$, $R_{\oplus }=6.371\times 10^6\mathrm{m}$, $M_J=1.89813\times 10^{27}\mathrm{kg}$, $G=6.67430\times 10^{-11}{\mathrm{m}}^3\mathrm{k}{\mathrm{g}}^{-1}{\mathrm{s}}^{-2}$ and $1\mathrm{au}=149597870700\mathrm{m}$. They are conversion constants, not additional measured system parameters.^[ref][ref]

Adopt circular orbits, negligible planetary masses, spherical bodies, uniform stellar brightness and no transit dilution. Interpret the measured duration as the total duration of a circular planetary transit, from first to fourth contact, $T_{14}$. These conventions matter especially for the off-centre crossing: finite planet size cannot automatically be ignored there. Unknown orbital eccentricity or limb darkening would require further information to infer unique stellar radii.

First convert the apparent magnitude difference into a fractional exoplanet transit depth, rather than treating a milli-magnitude as a part in a thousand:Both crossings are non-grazing because their transit impact parameters satisfy $b<1-k$, so the full projected planetary disc can block the star. Equal depths determine equal radius ratios, not equal absolute planet radii.

For orbital inclination $i$ measured from the sky plane, $b=(a/R_{\star })\cos i$. Choosing mid-transit as orbital phase zero, the projected separation obeysAt first and fourth contact, $z/R_{\star }=1+k$. The full crossing duration therefore isFor the numerical scale use solar radius $R_{\odot }=6.957\times 10^8\mathrm{m}$ and nominal solar mass parameter $(GM{)}_{\odot }=1.3271244\times 10^{20}{\mathrm{m}}^3{\mathrm{s}}^{-2}$, with one day equal to $86400\mathrm{s}$. These are stated conversion choices, not additional measurements of the two systems.^[ref]

Kepler's third law gives $a=[G{M}_{\star }{P}^2/(4{\pi }^2){]}^{1/3}$. For a solar-mass, solar-radius host, this gives $a/R_{\odot }=12.3051$ at five days and $6.68023$ at two days. Substitution in the contact formula yieldsThe transit duration of a small planet instead takes both $k\to 0$ and $R_{\star }/a\ll 1$, giving $T\simeq P{R}_{\star }\sqrt{1-b^2}/(\pi a)$. It predicts $3.10\mathrm{h}$ and $1.37\mathrm{h}$, respectively. The difference is appreciable for B because the fractional change in its crossing chord is larger. Durations defined by mid-ingress to mid-egress or by second to third contact should not be silently equated with $T_{14}$.

The original diagram uses a straight projected track to display the contact geometry; the duration formula above retains orbital curvature.

To infer the actual stellar radii, invert the contact formula before applying the supplied mass-radius scaling. Let $A=a/R_{\star }$ and $s=\sin (\pi T_{14}/P)$; thenThe measured durations imply $A_A=8.34661$ and $A_B=7.55743$. Normalize the stipulated constant mass-radius ratio to the Sun, so $M_{\star }/M_{\odot }=R_{\star }/R_{\odot }\equiv r_{\star }$. Kepler's third law then givesThus the stellar radius from transit duration and a linear mass-radius relation isThe same assumed scaling gives $M_A=1.79M_{\odot }$ and $M_B=0.831M_{\odot }$. A's longer crossing requires a lower mean stellar mass density, while B's shorter full-contact crossing requires a higher one. They are not both solar replicas under this model. Their planets have $R_{p,A}=0.1213R_{\odot }$ and $R_{p,B}=0.05633R_{\odot }$.

For a consistently applied point-planet approximation, the duration scales as $T\propto P^{1/3}{r}_{\star }^{2/3}\sqrt{1-b^2}$ and the inferred radii would instead be $R_A\simeq 1.98R_{\odot }$, $R_B\simeq 1.07R_{\odot }$. The latter is close to solar. Consequently the assertion of a substantial non-solar radius for B depends on which duration convention and finite-size approximation are intended. The contact-based values above use the measured depth rather than discarding it. If eccentric orbits are allowed, the transit speed introduces an additional factor depending on eccentricity and orientation, so the two durations alone cannot establish those radii uniquely.

For an outer coplanar planet around B, the common orbital inclination gives $b_2/b_B=a_2/a_B=(P_2/P_B{)}^{2/3}$. A point-planet centre must cross the stellar disc, requiring $b_2<1$. The usual coplanar transit limit is thereforeThis is a limiting value approached from below; a tangent path has no duration. For an outer planet of finite, unspecified radius ratio $k_2$, any overlap instead requires $b_2<1+k_2$, giving $P_2<P_B[(1+k_2)/b_B{]}^{3/2}$. If it had the same ratio $k_2=k$, that geometric limit would be $3.08$ days; requiring full overlap would give $2.52$ days. Because neither its radius nor a photometric detection threshold is supplied, no unique finite-size practical maximum period can be inferred beyond this conditional expression. The stellar radius cancels from the point-planet period-ratio bound.

For B's inner planet, a reflected-light planetary phase curve haswhere $A_g$ is geometric albedo and $\alpha$ is phase angle. Near exoplanet secondary eclipse the visible hemisphere is almost fully illuminated. Using the standard full-phase approximation and the contact-based $A_B$,The exact lost fraction of combined stellar and planetary light is $F_p/(F_{\star }+F_p)$, which changes the result negligibly at this contrast. The corresponding magnitude depth is about $0.0437$ milli-magnitudes. The circular orbit has the same impact distance at eclipse as at transit, and $b_B<1-k$, so the star can hide the entire planetary disc.

The inclination inferred for B is $83.92^{\circ }$, putting closest approach to full phase at $\alpha \simeq 6.08^{\circ }$. A fully precise prediction would therefore also require the reflected-light phase function, which the question does not specify; the boxed value uses $\Phi \simeq 1$. The consistent point-planet duration approximation yields about $56.7\mathrm{ppm}$ instead, through its different inferred $a/R_B$. It is important to keep the geometry used for the radius inference and the reflected-light calculation consistent.

Finally, the stellar radial-velocity semi-amplitude for a low-mass companion on a circular orbit isThe two stars' equal apparent magnitudes make a comparison at similar photon statistics plausible, though spectral-line widths, stellar variability and observing cadence can still differ. With the contact-based stellar estimates,If the planets have comparable mean densities, A is easier: their equal depth then gives $M_{p,A}/M_{p,B}\simeq (R_A/R_B{)}^3\simeq 9.99$, so $K_A/K_B\simeq 4.44$. The larger planetary mass outweighs A's longer orbital period and larger host mass. For equal planetary masses, B is easier by a factor $2.25$ in amplitude. The supplied transit data give planet radii, not masses, so the ranking is conditional on a planetary mass-radius assumption; a unique unconditional ranking would invent information.

Geodesic deviation compares neighboring members of a family of freely falling trajectories. Subtracting their geodesic equations to first order in the connecting vector gives a covariant relative-acceleration equation governed by the Riemann curvature tensor. For a slowly moving satellite, $u^0=c\dot{t}\simeq c$ dominates the spatial components. Choose simultaneous spatial separation in a locally transported frame, take the weak-field approximation, and neglect terms quadratic in the small separation and higher velocity corrections. The leading spatial equation becomesReplacing the transported covariant derivatives by ordinary derivatives is a leading local or Newtonian approximation; a rotating orbital frame would introduce additional inertial terms. In the convention established above, the instantaneous tidal acceleration stretches radial separations and compresses transverse ones.

The growth expected for nearby different-radius circular orbits is especially transparent from Kepler's third law. If the radii are $R$ and $R+\delta R$, their coordinate angular frequencies differ by $\delta \Omega \simeq -(3/2)\Omega \delta R/R$. Initially aligned orbits therefore acquire an along-track separation $R\delta \Omega t\simeq -(3/2)\Omega \delta Rt$. Over the interval in which the nearby-orbit approximation remains valid, this produces a growing separation without needing to solve the deviation equation.

This growth is not universal or indefinitely unbounded. Equal-radius co-orbiting geodesics with a fixed phase offset have constant separation, and slightly inclined equal-radius circular orbits have bounded, periodically varying separation. Even different-radius circular orbits have bounded exact positions; their relative angle eventually wraps around. The secular drift statement assumes distinct radii and describes the local differential-orbit regime.

The Newtonian achievement combined a dynamical framework, a universal gravitational interaction and a distinctive standard of physical explanation. It was problematic partly because its mathematical successes exceeded what it supplied as a physical mechanism. Calling it a system should not conceal the difference between Isaac Newton's own writings, their later development and simplified versions promoted by his followers.

The Philosophiae Naturalis Principia Mathematica by Isaac Newton began with definitions and the Newton's laws of motion. Inertia distinguished persistence in uniform rectilinear motion from changes requiring impressed force. Newton's second law related such changes to impressed force; Newton's third law required reciprocal interactions. The account also distinguished absolute and relative quantities, raising questions about the meaning of true motion. In modern notation, a central part of the programme is expressed byHere momentum, force, masses and separation have their modern notation, and $G$ is the gravitational constant; the display is not a quotation of Isaac Newton's symbols. Combining Newton's laws of motion with Newton's law of universal gravitation made it possible to explain a wide range of terrestrial and celestial phenomena together.^[ref]

The quantitative content matters. For an ideal circular orbit of a small body around a much larger mass $M$, centripetal acceleration gives, in modern notation,This restricted illustration recovers the period–radius relation in Kepler's third law. It is not a reconstruction of the full original proof for an elliptical orbit. More generally the same interaction could account for the Moon's departure from a straight path, the paths of comets and mutual perturbations. Universal attraction made celestial irregularities potential consequences of a common law rather than reasons to assign a separate nature to each class of object.

The associated Newtonian programme extended beyond gravity. In optics, experiments on refraction supported the heterogeneous composition of white light: a prism could separate constituents with different refrangibilities rather than manufacture colours from a previously homogeneous ray. This challenged accounts that treated colour primarily as a modification introduced by the medium. It belongs to the broader Newtonian programme, although the 1704 Opticks by Isaac Newton must not be conflated with the content or publication date of the 1687 Philosophiae Naturalis Principia Mathematica by Isaac Newton.^[ref]

The greatest conceptual difficulty concerned action at a distance. For a defender of a contact mechanism, the mutual attraction of bodies through empty space appeared to reintroduce an occult quality instead of explaining an interaction. Cartesian vortices offered at least a material picture of how bodies might be carried around, whereas an inverse-square force law did not specify what transmitted the influence. Isaac Newton could point to the difficulties of reconciling such vortices with planetary periods and the varied paths of comets, but defeating that particular mechanism did not itself furnish an alternative mechanical cause.^[ref]

It is particularly important that Isaac Newton did not simply identify unmediated attraction as an essential property of matter. In his letter to Richard Bentley of 25 February 1692/3, he rejected that attribution and left open whether the relevant agent was material or immaterial. In the General Scholium by Isaac Newton, he distinguished establishing how gravity acts from assigning a cause not warranted by the phenomena.^[ref][ref]

Other difficulties followed from the standard of demonstration itself. The demanding geometry made the work hard to master; approximations and mutual interactions meant that apparent agreement had to be earned through calculation, not assumed from a slogan about universal attraction. Absolute motion and divine government also raised philosophical questions that a successful orbital calculation could not settle. Initial resistance thus included reasonable disputes about mechanism, evidence and explanatory adequacy. The system's later authority came from its capacity to sustain successful investigation despite those unresolved foundations, not from having answered every foundational question at publication.

Isaac Newton made the greater contribution if the criterion is the scope and explanatory integration of mathematical natural philosophy; Galileo Galilei has the stronger claim if the criterion is the formative transformation of how terrestrial motion could be investigated. The first criterion gives a defensible overall verdict for Isaac Newton, but a serious comparison must show what each changed rather than award credit according to who worked later.

Galileo Galilei's contribution was not simply to attach numbers to familiar events. In the Third Day of Two New Sciences by Galileo Galilei, he defines uniformly accelerated motion by equal additions of speed in equal intervals of time, then derives consequences. Starting from rest, the modern notation is $v=at$; the mean speed over the interval is $at/2$, soHence distances in successive equal intervals have the ratios $1:3:5:7:\dots$. The novelty lies in making a precisely specified kind of motion the subject of mathematical demonstration, with observed falls assessed against its consequences. This is not the same as deriving the underlying physical cause of gravity.^[ref]

The reported apparatus shows how such an idealization could be investigated. A smooth bronze ball descended in a straight grooved wooden channel, raised at one end. Water discharged through a small opening into a vessel during the descent; its collected mass served as a measure of elapsed time. Repeating runs over different portions of the channel compared the predicted squared dependence of time on distance. The inclined plane slowed the event, while repetition and weighing made a temporal comparison possible without a modern stopwatch. The published report is evidence of Galileo Galilei's experimental presentation, not a surviving raw-data archive from which every reported precision can now be reconstructed.^[ref]

A necessary physical qualification is often missed. A ball in rolling without slipping is not a point particle sliding without friction. For a solid sphere on an ideal inclined plane, modern rigid body dynamics gives $a=5g\sin \alpha /7$, rather than $g\sin \alpha$. Both are constant accelerations, so both give $s\propto t^2$. The trial can therefore test the temporal form without directly establishing the free-fall coefficient. This distinction makes the experiment more intelligible: the evidence supports a specified relation under approximations, rather than every detail of a modern idealization.

The Fourth Day extends this approach to projectile motion. Combine uniform horizontal velocity $u$ with constant downward acceleration $g$, taking downward displacement positive. Then $x=ut$ and $y=g{t}^2/2$. Eliminating time givesa parabola. Actual projectiles experience air resistance, but the ideal curve organizes the effects that need separate treatment. Galileo Galilei thus gave a powerful example of how mathematical abstraction and physical observation could cooperate.^[ref]

Isaac Newton went further by connecting kinds of motion through a general account of force. Newton's second law links the change of momentum to the impressed force; for constant mass, its modern differential expression is $\mathbf{F}=m\mathbf{a}$. Newton's third law makes interaction reciprocal, enabling a consistent analysis of interacting bodies. The Philosophiae Naturalis Principia Mathematica by Isaac Newton develops consequences of central forces and then asks which laws fit astronomical phenomena. Mathematical deduction now connects measurable orbits to an underlying force law, rather than just describing each trajectory separately.^[ref]

A simple modern reconstruction shows the increase in explanatory reach. For a circular orbit of radius $r$ and orbital period $T$ about a much more massive central body, centripetal acceleration is $4{\pi }^{2}r/T^2$. Kepler's third law, $T^2\propto r^3$, therefore implies an acceleration proportional to $r^{-2}$. This circular argument is illustrative; Isaac Newton's treatment of elliptical orbits is more general. With Newton's law of universal gravitation,the two-body orbital period is tied to the semimajor axis $a$ and both masses. The same theory connects the Moon's orbital acceleration with terrestrial gravity, and supports further investigations of tides, comets and departures from a simple two-body orbit.^[ref]

This achievement depended on more than two individuals. Johannes Kepler's relations supplied crucial astronomical constraints, and both authors inherited geometry, observational records and technical practices. Nor should the Philosophiae Naturalis Principia Mathematica by Isaac Newton be described simply as the first application of symbolic calculus: its published arguments are largely geometrical. What distinguishes Isaac Newton is the demonstrative network uniting terrestrial and celestial problems, together with explicit idealizations whose failures generated further calculable questions.

On breadth, quantitative unification and the inference from motion to force, Isaac Newton's contribution was greater. Galileo Galilei's work supplied an exemplary way to mathematize terrestrial motion and test its idealized relations. Giving the former the overall verdict should therefore acknowledge the latter's distinct foundational achievement, not portray the Philosophiae Naturalis Principia Mathematica by Isaac Newton as a creation without predecessors.

For a circular binary with negligible spin, the two orbital radii about the center of mass are $a_1=a{M}_2/M$ and $a_2=a{M}_1/M$. Summing their orbital angular momenta givesUsing Kepler's third law, ${\Omega }^2=GM/a^3$, also gives $J=M_1{M}_2\sqrt{Ga/M}$. Under conservative binary mass transfer, $M$ and $J$ stay fixed, so $a\propto (M_1{M}_2{)}^{-2}$. Since $d\log M_2/d\log M_1=-q$, the Roche-lobe radius response exponent isFor the supplied giant response, luminosity stays fixed during the mass-loss perturbation and ${\zeta }_{\ast }=d\log R/d\log M_1=-0.27$. The change of the logarithmic overfill is $({\zeta }_{\ast }-{\zeta }_L)d\log M_1$. Because the donor loses mass, this reduces overfill and stabilizes transfer precisely when ${\zeta }_{\ast }>{\zeta }_L$. Thus the stated model requiresThis is a response criterion within the given radius law; a different adiabatic donor response would change the dynamical threshold.

The initially more massive star evolves first. If it waits until it has the supplied giant structure while still having $q>1$, conservative overflow is unstable by this criterion. If Roche-lobe overflow instead begins while the donor is on the main sequence or shortly after central hydrogen exhaustion, its largely radiative-envelope response can allow transfer before a deep convective giant envelope develops. Mass loss and accretion can then reverse the mass ratio before it reaches the giant state. This is the binary mass-ratio reversal route to an Algol binary, through Case A mass transfer or early case B mass transfer.

At lobe filling, the lobe approximation and Kepler's law give the pre-giant Roche-lobe-filling period limitIt increases with the donor's radius. Requiring $R_L<R_{\mathrm{BGB}}=1.68R_{\odot }(M_1/M_{\odot }{)}^{5/3}$ therefore yieldsThe companion mass cancels in this lobe approximation.

Since $P\propto a^{3/2}$ at fixed total mass, the separation scaling gives $\boxed{P(M_1{M}_2{)}^3=\text{constant}}$. Use solar-unit masses $m_i=M_i/M_{\odot }$ and $m=m_1+m_2$, and let $K=P(m_1{m}_2{)}^3$. ThenThe derivative vanishes at $m_1=5m/8$, and its derivative is $5/m_1^2+3/(m-m_1{)}^2>0$, while the ratio diverges at both endpoints. Hence this is the global minimum, the conservative pre-giant period-ratio minimum.

For the stated Algol parameters, $m=4$ and $K=3(1\cdot 3{)}^3=81$ days. At the minimum $m_1=2.5,m_2=1.5$, the period is $81/(2.5\cdot 1.5{)}^3=1.536$ days and $P/m_1^2=0.24576$ days, well below $P_0=0.8$ days. This is a possible initially more-massive donor configuration with overflow before the giant branch. The current $q=1/3$ is also below the giant stability threshold. Thus Algol is compatible with the proposed conservative early-overflow route, although this necessary geometric test alone does not prove its full evolutionary history.

For RT Lac, $m=2.4$ and $K=5.1(0.9\cdot 1.5{)}^3$. Its minimum occurs at $m_1=1.5,m_2=0.9$, where the mass product equals its current value, so $P=5.1$ days andNo conservative antecedent on this orbit can begin overflow before the giant branch. Beginning with an initially more massive giant instead encounters the instability just derived. RT Lac cannot follow the stated coeval, conservative Algol-formation model. A plausible alternative is rapid initial transfer with substantial mass and angular-momentum loss, via an outflow, strong wind or envelope-ejection episode. Such nonconservative binary mass transfer changes both total mass and $J$ and removes the period-product constraint. Its present mass ratio $0.9/1.5=0.6$ permits stable current giant transfer. The given data do not determine a unique mass-loss history or establish a specific ejection mechanism.

The most direct stellar mass determination uses orbital dynamics in a binary star. A resolved orbit and an independent distance give the relative semimajor axis $a$ from its angular size; measuring the period then gives the total mass by Kepler's third law:Individual barycentric orbital sizes, or the ratio of the two radial-velocity amplitudes, give $M_1/M_2=a_2/a_1=K_2/K_1$. A photocenter orbit alone needs the luminosity ratio to infer the barycentric orbit.

For a double-lined spectroscopic orbit, the line-of-sight amplitudes determine $a_i\sin i=K_{i}P\sqrt{1-e^2}/(2\pi )$. Combining them givesAn eclipsing binary light curve constrains the inclination and relative radii, turning these projected quantities into individual masses. Eclipse and orbit modeling must allow for eccentricity, surface brightness, limb darkening and possible tidal distortion.

If only one spectrum is measured, the binary mass function isWithout an inclination and an estimate of the visible star's mass it does not fix the unseen mass uniquely, but it provides useful constraints and lower limits.

Less direct estimates combine a spectroscopic surface gravity with a radius inferred from distance, bolometric flux and effective temperature: $R^2=L/(4\pi {\sigma }_{\mathrm{SB}}{T}_{\mathrm{eff}}^4)$ and $M=g{R}^2/G$. These depend on atmosphere modeling. Locating a star on an evolutionary track or using a calibrated main-sequence mass-luminosity relation also estimates mass, but assumes composition and evolutionary state; giants and mass-transfer products need not obey a single main-sequence relation. Dynamical binary masses provide the empirical calibration against which such model-dependent methods are tested.

Use the secular, circular-orbit approximation: the wind escapes fast enough to avoid appreciable interaction with the companion, while the stellar masses change slowly enough that orbital elements describe the evolving binary star. Neglect spin angular momentum. Let $\omega =(GM/a^3{)}^{1/2}$ be the orbital angular speed. The donor's distance from the centre of mass is $a_1=a{M}_2/M$, and the circular-binary orbital angular momentum isA spherically symmetric wind in the donor frame has zero mean additional angular momentum about the donor. Averaging its angular momentum about the binary centre of mass therefore gives the donor's orbital specific angular momentum, $j_1=a_1^2\omega$. Since $\dot{M}<0$ denotes mass removed from the system, ${\dot{J}}_{\mathrm{orb}}=j_1\dot{M}$, andThis is donor-wind angular-momentum loss in the Jeans-mode mass loss approximation.

First allow only wind loss, so ${\dot{M}}_1=\dot{M}$ and ${\dot{M}}_2=0$. The logarithmic derivative of $J_{\mathrm{orb}}$ givesSolving yields $\dot{a}/a=-\dot{M}/M$. By Kepler's third law, $P^2=4{\pi }^2{a}^3/(GM)$, so $\dot{P}/P=-2\dot{M}/M$. Integrating givesWind loss thus expands the separation and lengthens the orbital period. These are secular relations; an impulsive loss on an orbital time could instead generate eccentricity and would need a different analysis.

To test contact, put $D=d\log (R_1/R_L)/dt$. The given short-time stellar radius response exponent gives ${\dot{R}}_1/R_1=-n{\dot{M}}_1/M_1$. Differentiating the Roche lobe approximation givesFor wind alone this becomes ${\dot{R}}_L/R_L=\dot{M}/(3M_1)-4\dot{M}/(3M)$. Hence, with the binary mass ratio $q=M_1/M_2$,Because $\dot{M}$ is negative, increasing overfill, $D_{\mathrm{w}}>0$, occurs precisely whenThis is wind-driven Roche-lobe overflow. For the reverse strict inequality, $D_{\mathrm{w}}<0$ and the donor star retreats inside its Roche lobe, producing a detached binary. At equality the wind is neutral to contact at this order.

Now permit an additional conservative binary mass transfer rate $x={\dot{M}}_2$ to the mass-gaining star. The total escaping wind remains $\dot{M}$, while ${\dot{M}}_1=\dot{M}-x$. Transfer retains orbital angular momentum within the binary under the assumed neglect of spins, so the same wind-loss formula for ${\dot{J}}_{\mathrm{orb}}$ applies. Logarithmic differentiation now givesSubstituting this and ${\dot{M}}_1=\dot{M}-x$ into the two radius derivatives givesSteady contact requires $D=0$. If the wind is driving overflow and $6q<5-3n$, the transfer contribution reduces overfill and a positive contact rate exists:The numerator and denominator are positive in this regime, making ${\dot{M}}_2>0$. This explicitly establishes both the contact transfer rate with Jeans-mode wind loss and its stabilizing sign.

If the wind drives overflow but $6q>5-3n$, positive transfer instead increases $R_1/R_L$: the mass loss required to relieve overfill makes the overfill worse. The formal contact formula would require ${\dot{M}}_2<0$, which is not donor-to-companion transfer. There is no stabilizing positive contact rate; overflow runs away. For a red giant this can lead to a common envelope, followed by envelope ejection into a tighter binary or a merger. The criterion concerns the supplied short-time radius response; the actual runaway timescale and outcome require the donor's dynamical and thermal response. At $6q=5-3n$ transfer is neutral to contact at this order, so it cannot balance a nonzero wind driver. If the wind already causes detachment, the transfer instability condition alone does not initiate overflow.