Modified-Gravity Calculation Reshapes Roche Lobes—But Not Where the Two Lobes Meet

In a close binary, each star governs a teardrop-shaped pocket of space. Inside it, gas and plasma stay bound to that star; outside it—across a narrow gravitational saddle called the inner Lagrange point, L1—material can slip toward the companion. That boundary is the Roche lobe, and its size decides whether a star stays intact or begins donating mass to its neighbor.

For four decades, the workhorse for that calculation has been a compact formula published by Peter Eggleton in 1983. It gives the lobe's volume-equivalent radius: the radius of a sphere that would hold the same volume as the teardrop. Binary-evolution and population-synthesis codes rely on such fits to avoid repeatedly solving the full equipotential geometry. It assumes, of course, that gravity follows the inverse-square law.

A new preprint by A. M. Velásquez-Toribio, of the Federal University of Espírito Santo in Brazil, asks what happens to that teardrop if gravity does not quite follow Newton. The paper, posted to arXiv, extends the Roche-lobe calculation to two families of alternative gravity and builds Eggleton-style fitting formulas for them.

Why doubt the inverse-square law?

Newtonian gravity has been tested to extraordinary precision across the Solar System. But that precision does not make the inverse-square law inevitable: unification schemes involving extra dimensions, Randall–Sundrum-type braneworld models, and scalar-tensor theories of gravity can predict small departures at some distance scale. The paper considers two phenomenological shapes for those departures.

The first is the Yukawa potential, which modifies the usual 1/r potential with a strength α and a range λ. If α is large and λ is comparable to the binary's separation, the correction matters; if λ is tiny compared with the binary's orbital separation, gravity snaps back to Newtonian. The second is a power-law correction, of the kind Randall–Sundrum models suggest, adding a term that grows toward small distances.

What the author actually computed

The setup is the classic restricted circular three-body potential: two point masses in a circular orbit, with a massless test particle feeling their gravity. Working in a frame that rotates with the binary, the author writes down an effective potential—the two stars' gravity plus the centrifugal term—and finds the point on the line between them where the effective force vanishes in the rotating frame. That is L1.

From there, the Roche lobe is the closed equipotential surface passing through L1. To trace it, the author uses a tangent-field method: at each point, step along a direction perpendicular to the gradient of the potential, which keeps the trajectory on the same level surface, and periodically apply a Newton projection step to correct any drift. The resulting two-dimensional contour is revolved to get the full three-dimensional lobe and its volume. The volume-equivalent radius follows directly.

Lobes that look almost Newtonian

The paper shows contours for mass ratio q—the ratio of the two masses—from 0.04, a very unequal pair, to 1.2. For the parameters chosen, the modified-gravity lobes sit close to the Newtonian ones. At extreme mass ratios the shifts are small but visible; as the two stars approach equal mass, the three potentials produce nearly indistinguishable lobes.

The three-dimensional reconstructions tell a similar story. The overall geometry remains close, but sizes, volumes, and some boundary details shift. In the Yukawa example (range parameter Λ = 0.3, strength α = 1.5), the low-mass lobe inflates slightly because the near-field attraction is stronger, while the companion's lobe adjusts to keep the L1 balance. In the q = 1.200 Yukawa case, the equivalent-volume radii deviate at the percent level. For the power-law case (index p = 1.0, amplitude β = 0.3), volumes shift systematically with the correction factor.

The biggest warp is not at the throat

Because the lobes look so similar, the author added a diagnostic: measure the boundary as a function of angle around each star, with zero degrees pointing toward L1, and compare the modified-gravity radius with the Newtonian one.

The largest deformations often appear not at the throat but in the lateral sectors, around 80° to 110° from that direction. Changing the potential shifts the transverse width and curvature of the critical surface, not just the location of the point where the two lobes touch. For the parameters shown, the power-law correction distorts the boundary more than the Yukawa model, and the lobe that deforms most depends on mass ratio: the lighter star's lobe at small ratios, the companion's at larger ones.

Fast formulas, small residuals

The second half of the paper builds practical replacements for Eggleton's formula. The strategy is multiplicative: take Eggleton's Newtonian expression and multiply it by a series in fractional powers of the mass ratio q, with coefficients fitted to the numerical solutions.

For the Yukawa case the residuals stay within about 1 percent across the mass ratios explored, comparable to the Newtonian deviations from Eggleton's formula used as a benchmark. The power-law fit shown in Figure 6 (p = 1.0, β = 0.05) tracks its numerical points to a few percent. Both keep the one-third exponent that governs how the radius scales for a very low-mass star, and both respect the symmetry that swapping the two stars should map one lobe onto the other.

That matters because these fits are meant to be dropped into stellar-evolution and population-synthesis codes, where fast evaluations matter across large sets of binary histories.

What this calculation does not include

The work is deliberately narrow. It treats the stars as point masses and isolates only the effect of changing the two-body potential. Real close binaries also contain tides, rotational distortions, finite stellar size, asynchronous rotation, eccentricity, and non-perturbative distortions of the stars themselves. In the tightest systems, the paper notes, those effects may be comparable to or larger than any modified-gravity contribution.

The power-law correction is also only meaningful as a small perturbation: it must remain small over the relevant portion of the Roche surface, otherwise it should be read as an effective description of a stronger short-range modification rather than a small nudge to Newton. And to use the fitting formulas in real binary evolution, researchers will need a physically motivated range for the new parameters—something the paper does not supply.

So this is not a detection of modified gravity, and it is not yet a full model of interacting binaries. It is a controlled baseline: a way to ask how the critical Roche surface responds when the law of gravity itself is changed, with everything else held fixed.

That is the quiet appeal of the result. The Roche lobe is a gravitational fingerprint, and the new work shows the fingerprint does change—not dramatically, and not where you might first look, but measurably in these calculations, in the outward shoulders of the teardrop where the star bulges sideways. If gravity ever departs from Newton, that is where the deviation may wait.