pith. sign in

arxiv: 2604.09901 · v1 · pith:YBB3HZI3new · submitted 2026-04-10 · 🌌 astro-ph.SR

Near-critical magnetic fields in Kepler red giants

S. Deheuvels , J. Ballot , F. Ligni\`eres , G. Li , M. Villate This is my paper

Pith reviewed 2026-05-10 16:40 UTC · model grok-4.3

classification 🌌 astro-ph.SR
keywords red giantsmagnetic fieldsasteroseismologyKeplerg-modesmixed modesstellar coresdynamo
0
0 comments X

The pith

Near-critical magnetic fields of 100 to 700 kG confined below the hydrogen-burning shell explain irregular g-mode periods in eight Kepler red giants.

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

The paper examines eight red giant stars from Kepler observations where l=1 modes appear as doublets showing deviations from the regular g-mode period spacing, with partial suppressions in three stars. Using an adapted non-perturbative formalism to treat strong magnetic effects on mixed modes, the authors obtain good fits only when the doublet components are identified as m=0 and m=1 rather than m=±1. This yields core fields ranging from 100 to 700 kG that are localized well below the hydrogen-burning shell. The best-fit models are low-mass stars of 1.1 to 1.2 solar masses whose field extents match the maximum size reached by their main-sequence convective cores, suggesting a dynamo origin. A sympathetic reader cares because these results provide direct constraints on internal magnetic fields and their possible role in angular momentum transport inside evolved stars.

Core claim

The paper claims that the observed deviations from regular g-mode period patterns and partial suppressions in eight Kepler red giants are produced by near-critical magnetic fields in their cores. Applying a non-perturbative formalism adapted from the traditional approximation for rotation, the authors show that identifying the l=1 doublet components as m=0 and m=1 produces very good fits to all observations. The corresponding fields have intensities of 100 to 700 kG confined well below the H-burning shell. The best-fit models indicate low masses of 1.1-1.2 solar masses, with the radial extent of the fields approximately matching the maximal size of the convective core during the mainsequence

What carries the argument

The adapted non-perturbative formalism for the effects of near-critical magnetic fields on mixed mode frequencies, which permits extraction of the radial profile of the radial magnetic field component Br.

If this is right

  • The detected fields are confined well below the hydrogen-burning shell with strengths between 100 and 700 kG.
  • The radial profile of the radial field component can be inferred for near-critical strengths, unlike weaker fields where only a weighted average of Br squared is measurable.
  • The best-fit models correspond to low-mass stars of 1.1-1.2 solar masses.
  • The radial extent of the fields matches the maximum size of the main-sequence convective core, consistent with a dynamo origin there.

Where Pith is reading between the lines

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

  • If these core fields persist, they could continue to influence angular momentum transport in later evolutionary phases.
  • The non-perturbative method could be applied to other stars showing mixed-mode irregularities to search for similar strong fields.
  • Stellar evolution models should incorporate core dynamo action on the main sequence to predict field strengths and extents matching these observations.

Load-bearing premise

The deviations from regular g-mode periods and the partial suppressions are produced by magnetic fields rather than other mechanisms, and the adapted non-perturbative formalism remains accurate for near-critical field strengths.

What would settle it

No magnetic field profile confined below the hydrogen-burning shell reproduces the observed mode frequencies, period deviations, and partial suppressions when the doublet components are assigned as m=0 and m=1.

Figures

Figures reproduced from arXiv: 2604.09901 by F. Ligni\`eres, G. Li, J. Ballot, M. Villate, S. Deheuvels.

Figure 1
Figure 1. Figure 1: Stretched period échelle diagram for KIC 7458743, obtained from 4-yr Kepler time series. Only peaks with a SNR above eight are shown, the size of the symbols being proportional to the peak ampli￾tude. For clarity, peaks in the vicinity of l = 0 and l = 2 modes are omitted. The échelle diagram is folded with an apparent period spacing of ∆Πmeas 1 = 76.5 s, which produces the best visual alignment of the rid… view at source ↗
Figure 2
Figure 2. Figure 2: Stretched period échelle diagrams of the eight Kepler red giants of our sample. The échelle diagrams were folded using apparent period spacings ∆Πmeas 1 that offer good visual alignment of the ridges. To guide the eye, the ridges were highlighted with different colors (red for the most distorted ridge, blue for the other). tiplet are observed suggests that the non-axisymmetry effects are indeed weak. We th… view at source ↗
Figure 3
Figure 3. Figure 3: Location of the stars under study in the (∆ν-∆Π1) plane. The ap￾parent period spacings of g modes ∆Πmeas 1 are shown as blue star sym￾bols and the red star symbols indicate the asymptotic period spacing ∆Π1 measured by accounting for magnetic fields in a non-perturbative manner (Sect. 6). Gray crosses correspond to all red giants studied in Li et al. (2024). Moreover, the identification as m = ±1 component… view at source ↗
Figure 4
Figure 4. Figure 4: Stretched échelle diagram of KIC 7458743. Filled blue circles show the observed mode frequencies. Open red circles show asymptotic mixed mode frequencies including rotation and magnetic field effects, assuming that the observed ridges correspond to m = ±1 components. The core rotation was determined to ensure a crossing of the ridges at νcross = 150 µHz and magnetic field properties were adjusted to repro￾… view at source ↗
Figure 6
Figure 6. Figure 6: Period échelle diagrams of l = 1 g modes for the example case described in Sect. 4.6. Left: Frequencies of m = 0 (circles) and m = 1 (triangles) modes computed for a field with a radial exten￾sion σB = 0.003 R⋆. Magnetic effects are treated either with a non￾perturbative approach (red symbols), or a perturbative approach (gray symbols). Right: Frequencies of m = 0 modes for values of σ = 0.0015 (red), 0.00… view at source ↗
Figure 5
Figure 5. Figure 5: Numerical solution of the horizontal problem for different values of the asymmetry parameter aasym. We show λ as a function of a for (l = 1, m = 0) modes (top panel) and for m = ±1 modes (bottom panel). The colored dots indicate that mode suppression is reached, and the black dashed line corresponds to b = 1. In [PITH_FULL_IMAGE:figures/full_fig_p008_5.png] view at source ↗
Figure 7
Figure 7. Figure 7: Impact of the radial profile of Br on the mode frequencies. We considered fields with radial extensions of σB = 0.0015 R⋆ (red lines) , 0.002 R⋆ (green lines), 0.003 R⋆ (blue lines) and a constant field (black lines). Left: Radial profile of Br . Middle: Profile of the parameter a(r) in the g-mode cavity for the mode closest to νmax (Eq. 12). Right: Profile of λ(r) deduced from the interpolation of the a(λ… view at source ↗
Figure 8
Figure 8. Figure 8: Stretched échelle diagrams of three stars of our sample. Blue circles and horizontal blue lines correspond to observed modes and their 1-σ uncertainties. Red empty circles correspond to the best asymptotic models including magnetic effects in a non-perturbative manner. perturbative approach, the non-rotating case holds an ambigu￾ity. The oscillation spectrum corresponding to δω0 and aasym is identical to t… view at source ↗
Figure 9
Figure 9. Figure 9: Corner plots of KIC 9227589 (left) and KIC 11408970 (right), restricted to the asymptotic period ∆Π1 of dipolar g modes, and the parameters characterizing the core magnetic field: the magnetic shift δν0 = δω0/(2π), σB, and aasym. The corner plots for other stars are given in Appendix D [PITH_FULL_IMAGE:figures/full_fig_p013_9.png] view at source ↗
read the original abstract

The recent seismic detection of magnetic fields in red giants cores has given the opportunity to characterize these fields, potentially giving information about their origin and their role in the internal transport of angular momentum. We detect strong deviations from the regular pattern of g-mode periods in eight Kepler red giants showing $l=1$ doublets. In three of these stars, the modes show partial suppression. We investigate the magnetic origin of these features and determine the characteristics of the core fields that can produce such signatures. We need to invoke strong, near-critical fields. Assessing the effects of such fields on the mixed mode frequencies requires a non-perturbative approach. We use and adapt a formalism that was recently proposed following a similar development as the traditional approximation for rotation (TAR). We then compute asymptotic expressions of mixed mode frequencies including magnetic effects and attempt to reproduce the observed oscillation spectra. We show that for near-critical fields, information can be obtained about the radial profile of the radial field $B_r$, as opposed to weaker fields for which only a weighted average of $B_r^2$ can be measured. For the eight targets, we find that the $l=1$ doublets cannot be identified as the $m=\pm1$ components. Instead, we show that very good fits to all the observations can be obtained by identifying the two components as $m=0$ and $m=1$. These solutions correspond to fields with intensities ranging from 100 to 700 kG that are confined well below the H-burning shell. Our best-fit models for the eight stars have low masses (1.1-1.2 $M_\odot$) and the maximal size of their convective core during the main sequence approximately corresponds to the radial extent of the measured magnetic fields. The detected fields could thus have been generated by dynamo action in the main-sequence convective core.

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 paper claims that deviations from regular g-mode periods and partial suppressions observed in l=1 doublets of eight Kepler red giants are produced by near-critical core magnetic fields. Adapting a non-perturbative TAR-like formalism, the authors show that identifying the doublet components as m=0 and m=1 (rather than m=±1) yields good fits with field strengths 100-700 kG confined well below the H-burning shell. The best-fit models have masses 1.1-1.2 M⊙, with the radial extent of the fields matching the maximum size of the main-sequence convective core, suggesting a dynamo origin.

Significance. If the magnetic interpretation and formalism hold, this would advance characterization of strong core fields in red giants, providing radial profile information on Br (beyond weighted averages of Br²) and linking fields to main-sequence dynamos. This has implications for internal angular momentum transport. The non-perturbative approach for near-critical regimes is a methodological step forward, though its unvalidated status limits current significance.

major comments (3)
  1. The adapted non-perturbative formalism (described in the methods section on mixed-mode frequencies) is not validated against direct numerical solutions of the full linear adiabatic MHD oscillation equations in the regime where magnetic and buoyancy frequencies are comparable. This is load-bearing for the central claim, as small inaccuracies could bias the reported field strengths, radial confinement, and preference for the m=0/m=1 identification over m=±1.
  2. In the results section on fits to the eight stars, the m=0/m=1 identification is adopted specifically because it permits good fits to the observed spectra; this introduces circularity between model choice and data, especially since no independent tests rule out m=±1 or non-magnetic mechanisms for the period deviations and suppressions. No quantitative fit metrics (e.g., reduced chi-squared) or error bars on derived parameters are supplied.
  3. The data-selection criteria for the eight targets and any validation tests of the overall procedure (e.g., recovery of injected signals or comparison to perturbative limits) are not detailed, which is load-bearing given that the solutions are fitted quantities for field strength and extent.
minor comments (2)
  1. The abstract states that 'very good fits' are obtained but does not define the quantitative threshold or metric used to assess fit quality to the oscillation spectra.
  2. Notation for 'near-critical' fields and the radial profile of Br should be clarified with explicit definitions or equations early in the text to aid readability.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for their detailed and constructive review of our manuscript. We address each major comment below, providing our responses and indicating planned revisions where appropriate to improve the paper's clarity and rigor.

read point-by-point responses

  1. Referee: The adapted non-perturbative formalism (described in the methods section on mixed-mode frequencies) is not validated against direct numerical solutions of the full linear adiabatic MHD oscillation equations in the regime where magnetic and buoyancy frequencies are comparable. This is load-bearing for the central claim, as small inaccuracies could bias the reported field strengths, radial confinement, and preference for the m=0/m=1 identification over m=±1.

    Authors: We acknowledge that direct numerical validation against the full linear adiabatic MHD equations would be the most rigorous test, particularly in the near-critical regime. Such computations are computationally intensive and were not feasible within the scope of this work. The formalism is an adaptation of the recently proposed magnetic TAR approach, and we have verified its reduction to the known perturbative limit for weaker fields. In the revised manuscript, we will add an expanded discussion section on the method's assumptions, validity range, and consistency checks with perturbative results to better contextualize potential limitations and their impact on the inferred parameters. revision: partial

  2. Referee: In the results section on fits to the eight stars, the m=0/m=1 identification is adopted specifically because it permits good fits to the observed spectra; this introduces circularity between model choice and data, especially since no independent tests rule out m=±1 or non-magnetic mechanisms for the period deviations and suppressions. No quantitative fit metrics (e.g., reduced chi-squared) or error bars on derived parameters are supplied.

    Authors: We agree that quantitative fit metrics and error estimates are needed. In the revised manuscript, we will report reduced chi-squared values for the adopted m=0/m=1 fits as well as for attempted m=±1 identifications, which produce significantly poorer matches to the observed period spacings and suppressions. Error bars on the derived field strengths and radial extents will also be included based on the fitting procedure. We will further discuss why non-magnetic mechanisms (e.g., rotation or structural variations) fail to reproduce the specific observed patterns, thereby reducing the appearance of circularity by presenting comparative evidence. revision: yes

  3. Referee: The data-selection criteria for the eight targets and any validation tests of the overall procedure (e.g., recovery of injected signals or comparison to perturbative limits) are not detailed, which is load-bearing given that the solutions are fitted quantities for field strength and extent.

    Authors: We will expand the methods section to detail the selection criteria used to identify the eight targets from the Kepler red-giant sample, specifically the requirements for clear l=1 doublets exhibiting period deviations and partial suppressions. We will also describe validation aspects of the procedure, including consistency checks against the perturbative limit for weaker fields and internal tests of the fitting robustness. These additions will provide the requested transparency on how the targets and solutions were obtained. revision: yes

Circularity Check

0 steps flagged

No significant circularity in derivation chain

full rationale

The paper applies an adapted non-perturbative formalism (modeled on TAR) to compute asymptotic mixed-mode frequencies under near-critical magnetic fields, then fits field strength, radial extent, and mode identification (m=0/m=1 rather than m=±1) to reproduce the observed g-mode period deviations and partial suppressions in eight stars. This constitutes standard inverse modeling and parameter estimation from data rather than a claimed first-principles derivation. No equations reduce to their inputs by construction, no fitted quantities are relabeled as independent predictions, and no load-bearing self-citations or ansatzes are quoted that would make the central results equivalent to the inputs. The reported field ranges (100-700 kG, confined below the H-burning shell) and mass estimates emerge from the fits but are presented as characterizations consistent with the data under the model, not as self-derived outputs. The derivation remains self-contained against external seismic observations.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The central claim rests on fitting magnetic-field parameters to seismic data under the assumptions that the observed period deviations are magnetic and that the non-perturbative formalism is valid near the critical field strength.

free parameters (2)
  • core magnetic field strength = 100-700 kG
    Adjusted per star to reproduce observed deviations in g-mode periods and partial suppressions; values range 100-700 kG.
  • radial extent of magnetic field
    Determined to lie well below the H-burning shell to match the data.
axioms (2)
  • domain assumption The observed deviations from regular g-mode period patterns and the partial mode suppressions are caused by magnetic fields.
    The paper states it investigates the magnetic origin of these features.
  • domain assumption The adapted non-perturbative formalism accurately computes mixed-mode frequencies for near-critical magnetic fields.
    Used to derive asymptotic expressions and perform the fits.

pith-pipeline@v0.9.0 · 5653 in / 1664 out tokens · 66853 ms · 2026-05-10T16:40:21.411169+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

4 extracted references · 4 canonical work pages

  1. [1]

    , " * write output.state after.block = add.period write newline

    ENTRY address archiveprefix author booktitle chapter edition editor howpublished institution eprint journal key month note number organization pages publisher school series title type volume year label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts #0 'before.all := #1 ...

    work page

  2. [2]

    write newline

    " write newline "" before.all 'output.state := FUNCTION n.dashify 't := "" t empty not t #1 #1 substring "-" = t #1 #2 substring "--" = not "--" * t #2 global.max substring 't := t #1 #1 substring "-" = "-" * t #2 global.max substring 't := while if t #1 #1 substring * t #2 global.max substring 't := if while FUNCTION word.in bbl.in " " * FUNCTION format....

    work page

  3. [3]

    , " * write output.state after.block = add.period write newline

    ENTRY address author booktitle chapter edition editor howpublished institution journal key month note number organization pages publisher school series title type volume year label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts #0 'before.all := #1 'mid.sentence := #2 '...

    work page

  4. [4]

    write newline

    " write newline "" before.all 'output.state := FUNCTION n.dashify 't := "" t empty not t #1 #1 substring "-" = t #1 #2 substring "--" = not "--" * t #2 global.max substring 't := t #1 #1 substring "-" = "-" * t #2 global.max substring 't := while if t #1 #1 substring * t #2 global.max substring 't := if while FUNCTION word.in bbl.in " " * FUNCTION format....

    work page