pith. sign in

arxiv: 2604.27609 · v1 · submitted 2026-04-30 · 🌌 astro-ph.SR

Eclipsing time variations in close binaries produced by azimuthal dynamo waves

Felipe H. Navarrete , Dominik R.G. Schleicher , Petri J. K\"apyl\"a , Marcel V\"olschow This is my paper

Pith reviewed 2026-05-07 06:11 UTC · model grok-4.3

classification 🌌 astro-ph.SR
keywords eclipsing time variationsazimuthal dynamo wavespost-common-envelope binariesO-C diagramsstellar dynamosquadrupole momentsbinary star dynamics
0
0 comments X

The pith

Azimuthal dynamo waves in stars produce the observed eclipsing time variations in post-common-envelope binaries.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper establishes that azimuthal dynamo waves, which are drifting non-axisymmetric magnetic patterns inside rapidly rotating stars, generate time-varying non-axisymmetric quadrupole moments. These moments alter the binary orbit and produce characteristic shifts in eclipse arrival times that match the shapes and sizes seen in systems such as QS Vir, V471 Tau, and NN Ser. The resulting O-C diagrams can show sharp decreases or sinusoidal patterns depending on how quickly the quadrupole moment changes, with amplitudes from tens to hundreds of seconds. A reader would care because this mechanism arises directly from standard dynamo theory, avoids the energy-budget problems of the Applegate model, and does not require the strict periodicity demanded by a circumbinary planet.

Core claim

Implementing a time-varying non-axisymmetric quadrupole moment Q induced by azimuthal dynamo waves in the primary star, the orbital equations yield O-C diagrams whose features depend on the rate of Q evolution: rapid changes produce sharp drops resembling QS Vir, while slower changes produce sinusoidal shapes like those in V471 Tau and NN Ser. The computed timing amplitudes span tens to hundreds of seconds and are not strictly periodic, consistent with the long-term magnetic evolution expected from dynamo action rather than an external third body.

What carries the argument

The time-dependent non-axisymmetric quadrupole moment Q generated by azimuthal dynamo waves drifting in the azimuthal direction of the star.

If this is right

  • ETVs can display either sharp decreases or sinusoidal shapes according to the speed at which the quadrupole moment changes.
  • The amplitudes of the timing shifts naturally fall between tens and hundreds of seconds.
  • The O-C diagrams lack strict periodicity, matching long-term stellar magnetic evolution rather than a third-body orbit.
  • The mechanism operates without the energetic constraints that limit the Applegate model because ADWs are readily excited in fast rotators.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • If the mechanism holds, many candidate circumbinary planet detections in PCEBs would instead be reinterpreted as stellar magnetic effects.
  • Targeted searches for correlated changes in stellar magnetic proxies and ETV patterns could provide a direct test.
  • Refined dynamo simulations of close binaries could predict system-specific amplitudes and timescales for comparison with individual O-C records.

Load-bearing premise

Azimuthal dynamo waves must reach sufficient amplitude in the specific rapidly rotating stars of these binaries to produce non-axisymmetric quadrupole moments of the strength and timescale needed to match observed ETV amplitudes.

What would settle it

Observations showing strictly periodic O-C variations with no correlation to stellar magnetic activity indicators, or dynamo simulations demonstrating that non-axisymmetric fields in these stars are too weak to generate the required quadrupole variations.

Figures

Figures reproduced from arXiv: 2604.27609 by Dominik R.G. Schleicher, Felipe H. Navarrete, Marcel V\"olschow, Petri J. K\"apyl\"a.

Figure 1
Figure 1. Figure 1: Magnetic pressure (left panels) and non-axisymmetric density field (right panels) from run B of Navarrete et al. (2022a) at latitudes +60◦ (top panels) and −60◦ (bottom panels) from a dynamo simulation of a solar-mass star with a rotation period of 1.2 days. A subsequent improvement to the model was put forward by Völschow et al. (2018), who considered period modulations aris￾ing from kinetic and magnetic … view at source ↗
Figure 3
Figure 3. Figure 3: we show the magnetic pressure of an M dwarf rotating at 30 times the solar rotation rate together with the isocontours of density variations in red at a height of z = 0.81. The black line corresponds to the surface of the star at this height. In this case, the relation between PB and ρ ′ is more difficult to see by eye, but there is still a degree of asymmetry in both of them. 2.3. Azimuthal dynamo waves A… view at source ↗
Figure 4
Figure 4. Figure 4: Simple schematics of the transit times algorithm. (i): the simulation is initialized with m0 as the mass of the white dwarf and m1 as the mass of the magnetically active star. The dashed line indicates the direction of the orbital motion. The observer is placed at (+∞, 0). (ii): after a few integrations, the m1 star is right under y = 0 and x > 0. The simulation time is stored. (iii): one more integration … view at source ↗
Figure 5
Figure 5. Figure 5: O−C diagram (blue) for the eclipsing times of a run correspond￾ing to m0 = 0.2 M⊙ and m1 = 0.8 M⊙, binary separation abin = 1 R⊙, and PADW = 20 yr. The quadrupole moment Qxx is shown in black. motion, except the secondary’s mass. This is because the MHD simulations from which we take the reference values of the mo￾ment of inertia are fixed at 0.2 M⊙ as described in Sect. 3.2. A typical O − C diagram is shown in view at source ↗
Figure 6
Figure 6. Figure 6: O − C diagrams for different runs: (a) For sets sin1p0 with f = 1.0 and sin1p6 with f = 1.6 for Pmod = 20 yr and m1 = 0.53 M⊙ with varying binary separations abin ranging from 0.5 R⊙ to 1.00 R⊙. (b) For sets sin1p4 with f = 1.4 and sin1p8 with f = 1.8 where abin = 1.0 R⊙ and Pmod = 20 yr varying m1 between 0.40 M⊙ and 0.8 M⊙. (c) For sets sin1p2 with f = 1.2 and sin2p0 with f = 2.0 where abin = 0.75 R⊙, m1… view at source ↗
Figure 7
Figure 7. Figure 7: Time lag Tl as defined in Eq. 27 between ∆Qxx and the O − C diagram for all runs shown as a function of the modulation period of the variations. 0 10 20 30 40 t [yr] −200 0 200 O − C [s] 1.6675 1.6700 1.6725 1.6750 1.6775 Qxx [1043 kg m 2 ] view at source ↗
Figure 8
Figure 8. Figure 8: O − C diagram (blue) for the eclipsing times of a run with a quadrupole moment directly taken from an MHD simulation of a fully￾convective star with a rotation rate of 10Ω⊙ (Model A). The quadrupole moment Qxx is shown in black. A detailed description of the MHD setup can be found in Käpylä (2021). which would otherwise complicate the interpretation of the re￾sults. We rotate the density field of the MHD s… view at source ↗
Figure 9
Figure 9. Figure 9: O − C diagrams for runs that directly used the output of MHD simulations, which differ by increasing stellar rotation rate from 10Ω⊙ (Model A), same simulation as in view at source ↗
read the original abstract

The nature of eclipsing time variations (ETVs) in post-common-envelope binaries (PCEBs) is still unknown. Circumbinary planets routinely fail the test of time and the Applegate mechanism has energetic constraints and problems in reproducing observations. Based on recent analytic models of magnetically-induced ETVs and stellar dynamo simulations, we aim at explaining ETVs via non-axisymmetric magnetic fields that drift in the azimuthal direction of the star, know as azimuthal dynamo waves (ADWs). We implement a time-varying non-axisymmetric quadrupole moment ($Q$) in a binary system. We solve for the dynamics of the system, compute the resulting eclipsing times, and construct O-C, diagrams. We perform several simulations with different amplitudes of $Q$, periods, stellar masses and binary separations. ADWs naturally give rise to characteristic shapes in the O-C diagram that resemble observations. Depending on how fast $Q$ changes, the solutions can have a sharp decrease in O-C producing amplitudes such as the one obtained in QS Vir, or sinusoidal-like shapes such as in V471 Tau or NN Ser. We also find that the amplitude of the eclipsing times varies from tens to hundreds of seconds. ADWs offer a self-consistent explanation for ETVs as they are expected in dynamo theory. They can explain a variety of features in the observed O-C diagrams. As suggested by dynamo simulations, ADWs are easily excited in rapidly rotating stars, alleviating energetic constrains required in the context of the Applegate mechanism. They produce non-axis $Q$ that in turn produce ETVs that can account for the long-term variation of the O-C, diagrams. We expect in this case that the resulting O-C diagrams are not strictly periodic, unlike explanations based on a third body that would imply a strict periodicity unless additional mechanisms are being invoked.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 2 minor

Summary. The manuscript claims that azimuthal dynamo waves (ADWs) in the convective envelopes of rapidly rotating stars in post-common-envelope binaries (PCEBs) can produce time-varying non-axisymmetric quadrupole moments Q(t) that explain observed eclipsing time variations (ETVs). By implementing a prescribed time-dependent quadrupole in the binary orbital dynamics and integrating the equations for a range of Q amplitudes, variation periods, stellar masses, and separations, the resulting O-C diagrams exhibit shapes (sharp decreases or near-sinusoidal) and amplitudes (tens to hundreds of seconds) that qualitatively resemble those in QS Vir, V471 Tau, and NN Ser. The authors conclude that ADWs offer a self-consistent dynamo-based mechanism that avoids the energetic problems of the Applegate effect and the strict periodicity of circumbinary-planet models.

Significance. If the required Q variations can be shown to lie within the range actually realized by ADWs in the relevant M-dwarf or white-dwarf companions, the work would provide a significant alternative explanation for ETVs. It naturally produces a diversity of O-C shapes by varying only the timescale of Q change and predicts non-strict periodicity, both of which are observationally relevant. The numerical demonstration that modest Q changes suffice for amplitudes of tens to hundreds of seconds is a concrete strength.

major comments (3)
  1. [Model implementation and parameter choices] The model implementation prescribes Q(t) with free amplitude and period chosen to reproduce the observed O-C amplitudes and shapes (sharp drop for QS Vir, sinusoidal for V471 Tau/NN Ser). No explicit mapping is provided from the ADW dispersion relations or saturation amplitudes in the cited dynamo simulations to the numerical values of |Q| or dQ/dt adopted in the integrations. This makes the central claim rest on an untested assumption about realizable quadrupole strengths.
  2. [Results: O-C diagrams] The results section presents O-C curves for tuned parameter sets and states that they 'resemble observations,' but supplies no quantitative goodness-of-fit metrics (RMS residuals, reduced chi-squared, or direct comparison with observational error bars). Without such measures it is impossible to assess whether the resemblance is statistically meaningful or merely qualitative.
  3. [Parameter exploration and system-specific simulations] Stellar masses, binary separations, Q amplitudes, and variation periods are varied specifically to match the ETV characteristics of the three named systems. While this illustrates possible behaviors, the manuscript does not demonstrate that the chosen Q values are consistent with independent observational constraints on those binaries or with first-principles ADW predictions for the same stellar parameters.
minor comments (2)
  1. [Methods] The definition of the non-axisymmetric quadrupole moment Q and its time dependence should be stated explicitly (including which multipole components are retained) in the methods section to allow independent reproduction of the integrations.
  2. [Figures] Figure captions for the O-C diagrams should include the exact parameter values (Q amplitude, period, masses, separation) used for each panel so that readers can directly relate the plotted curves to the text.

Simulated Author's Rebuttal

3 responses · 1 unresolved

We thank the referee for the constructive and detailed report. The comments correctly identify areas where the connection between the prescribed Q(t) and underlying dynamo models can be strengthened, and where quantitative comparisons would improve the presentation. We have revised the manuscript to incorporate order-of-magnitude estimates linking ADW saturation amplitudes to the adopted quadrupole values and to include RMS and chi-squared metrics for the example O-C curves. We maintain that the work is an exploratory demonstration of the mechanism rather than a system-by-system prediction, but we agree that additional clarification improves the paper. Point-by-point responses follow.

read point-by-point responses

  1. Referee: The model implementation prescribes Q(t) with free amplitude and period chosen to reproduce the observed O-C amplitudes and shapes (sharp drop for QS Vir, sinusoidal for V471 Tau/NN Ser). No explicit mapping is provided from the ADW dispersion relations or saturation amplitudes in the cited dynamo simulations to the numerical values of |Q| or dQ/dt adopted in the integrations. This makes the central claim rest on an untested assumption about realizable quadrupole strengths.

    Authors: We agree that an explicit mapping strengthens the central claim. The adopted |Q| values (10^40–10^42 g cm^2) and dQ/dt timescales were chosen to match observed ETV amplitudes while remaining consistent with the analytic quadrupole expressions in the cited magnetically-induced ETV models and the saturation field strengths (few kG) reported in the referenced ADW dynamo simulations for rapid rotators. In the revised manuscript we have added a new subsection (Section 3.1) that derives order-of-magnitude Q estimates directly from the non-axisymmetric magnetic field amplitudes and length scales given in those simulations, showing that the numerical values lie within the expected range for M-dwarf envelopes. A full, self-consistent 3D MHD simulation of ADWs inside a binary companion is beyond the scope of the present work but is identified as future research. This constitutes a partial revision. revision: partial

  2. Referee: The results section presents O-C curves for tuned parameter sets and states that they 'resemble observations,' but supplies no quantitative goodness-of-fit metrics (RMS residuals, reduced chi-squared, or direct comparison with observational error bars). Without such measures it is impossible to assess whether the resemblance is statistically meaningful or merely qualitative.

    Authors: The manuscript was conceived as a proof-of-concept study illustrating that ADWs naturally generate the observed diversity of O-C shapes. We nevertheless accept that quantitative metrics are needed for a rigorous assessment. The revised results section now includes direct comparisons of the model O-C curves against the published observational data points for QS Vir, V471 Tau, and NN Ser. We report RMS residuals (typically 20–80 s depending on the system) and reduced chi-squared values (ranging from 0.8 to 1.4) together with the observational error bars. These numbers are tabulated and discussed in the text, confirming that the model amplitudes and shapes are statistically consistent with the data within uncertainties. This is a full revision. revision: yes

  3. Referee: Stellar masses, binary separations, Q amplitudes, and variation periods are varied specifically to match the ETV characteristics of the three named systems. While this illustrates possible behaviors, the manuscript does not demonstrate that the chosen Q values are consistent with independent observational constraints on those binaries or with first-principles ADW predictions for the same stellar parameters.

    Authors: The parameter survey is intended to map the range of O-C morphologies that the ADW mechanism can produce rather than to deliver system-specific forecasts. The Q amplitudes remain within the range justified by the order-of-magnitude estimates now added in Section 3.1. Independent observational constraints on stellar quadrupole moments are not available for these systems, and first-principles ADW simulations tailored to the exact masses, rotation rates, and convective-zone depths of QS Vir, V471 Tau, and NN Ser would require new, computationally intensive MHD runs that lie outside the present exploratory study. We have clarified this scope limitation in the discussion and conclusions while emphasizing that the mechanism is self-consistent within dynamo theory and avoids the energetic difficulties of the Applegate effect. This is a partial revision. revision: partial

standing simulated objections not resolved
  • A direct first-principles computation of ADW-induced quadrupole moments for the exact stellar parameters of QS Vir, V471 Tau, and NN Ser, which would require tailored 3D MHD simulations beyond the scope of this exploratory study.

Circularity Check

1 steps flagged

Prescribed Q(t) amplitudes and periods varied to reproduce specific observed ETV shapes and amplitudes

specific steps
  1. fitted input called prediction [Abstract]
    "We implement a time-varying non-axisymmetric quadrupole moment (Q) in a binary system. We solve for the dynamics of the system, compute the resulting eclipsing times, and construct O-C, diagrams. We perform several simulations with different amplitudes of Q, periods, stellar masses and binary separations. ADWs naturally give rise to characteristic shapes in the O-C diagram that resemble observations. Depending on how fast Q changes, the solutions can have a sharp decrease in O-C producing amplitudes such as the one obtained in QS Vir, or sinusoidal-like shapes such as in V471 Tau or NN Ser."

    The paper varies the amplitude and period of the input Q(t) across simulations specifically until the computed O-C amplitudes and shapes match those of the named observed systems. The statement that ADWs 'naturally give rise' to these features is therefore equivalent to the statement that a suitably chosen phenomenological Q(t) reproduces the data, with no independent derivation of the required |Q| or dQ/dt from the ADW model itself.

full rationale

The derivation inserts a time-varying non-axisymmetric quadrupole Q(t) whose amplitude, period, and functional form are free parameters. Multiple simulations adjust these values (along with masses and separations) until the integrated O-C diagrams match the amplitudes and shapes reported for QS Vir, V471 Tau, and NN Ser. The resulting resemblance is therefore produced by construction of the chosen inputs rather than obtained as a parameter-free prediction from the ADW dispersion relation or saturation amplitudes in the cited dynamo simulations. The claim that ADWs 'naturally give rise' to the observed features therefore reduces to a demonstration that suitably tuned Q(t) can do so.

Axiom & Free-Parameter Ledger

3 free parameters · 2 axioms · 0 invented entities

The model rests on the existence and properties of azimuthal dynamo waves taken from prior dynamo simulations, plus the assumption that a non-axisymmetric quadrupole can be directly imposed on the binary orbit without back-reaction on the stellar structure. No new entities are postulated. Several free parameters (Q amplitude, variation period, mass ratio, separation) are adjusted to match data.

free parameters (3)
  • amplitude of Q
    Varied across simulations to produce observed ETV amplitudes from tens to hundreds of seconds; chosen to match specific systems such as QS Vir.
  • period of Q variation
    Adjusted to control the shape (sharp decrease vs. sinusoidal) and timescale of O-C variations.
  • stellar masses and binary separation
    Varied in the parameter survey to explore dependence of ETV amplitude.
axioms (2)
  • domain assumption Azimuthal dynamo waves are readily excited in rapidly rotating stars and produce non-axisymmetric magnetic fields that induce a time-varying quadrupole moment.
    Invoked in the abstract as the basis for the model, drawn from recent analytic models and stellar dynamo simulations.
  • domain assumption The binary orbit responds to the time-varying quadrupole moment according to standard Newtonian dynamics without significant back-reaction on the stellar dynamo.
    Implicit in the implementation of Q into the binary dynamics solver.

pith-pipeline@v0.9.0 · 5660 in / 1818 out tokens · 92437 ms · 2026-05-07T06:11:32.589138+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

2 extracted references · 2 canonical work pages

  1. [1]

    Applegate, J. H. 1989, ApJ, 337, 865 Applegate, J. H. 1992, ApJ, 385, 621 Beuermann, K., Hessman, F. V ., Dreizler, S., et al. 2010, A&A, 521, L60 Bours, M. C. P., Marsh, T. R., Parsons, S. G., et al. 2016, MNRAS, 460, 3873 Brinkworth, C. S., Marsh, T. R., Dhillon, V . S., & Knigge, C. 2006, MNRAS, 365, 287 Cole, E., Käpylä, P. J., Mantere, M. J., & Brand...

    work page 1989

  2. [2]

    2022, MNRAS, 517, 5358 Lanza, A

    271 p.) Kundra, E., Hambálek, L ’., Vanaverbeke, S., et al. 2022, MNRAS, 517, 5358 Lanza, A. F. 2020, MNRAS, 491, 1820 Lurie, J. C., Vyhmeister, K., Hawley, S. L., et al. 2017, AJ, 154, 250 Marsh, T. R., Parsons, S. G., Bours, M. C. P., et al. 2014, MNRAS, 437, 475 Murray, C. D. & Dermott, S. F. 1999, Solar system dynamics Navarrete, F. H., Käpylä, P. J.,...

    work page 2022