Pith. sign in

REVIEW 4 major objections 4 minor 2 references

Exploring the Small-scale Magnetic Fields in the Atmosphere of HD 49385 by Asteroseismic Analysis

T0 review · 4 major / 4 minor · reviewed 2026-08-01 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read Asteroseismic fits reveal ~80-gauss magnetic fields about 1,850 km above the surface of the evolved Sun-like star HD 49385, while pinning down its helium core and mass.

desk verdict Applies a solar-calibrated magnetic-atmosphere prescription to HD 49385 and fits well, but the 80 G/1850 km claim is not uniquely identified against standard surface effects. read the letter →

arxiv 2607.19724 v1 pith:52XE2XLB submitted 2026-07-22 astro-ph.SR

Yuetong Wang , Yaguang Li , Yan Li , Guifang Lin , Tao Wu This is my paper
classification astro-ph.SR
keywords asteroseismologysmall-scalemagneticfieldsstellaratmospheresp-modeoscillationsavoidedcrossingheliumcoresolar-likepulsatorsHD49385
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper tries to show that small-scale magnetic fields in the atmosphere of a remote Sun-like star can be detected through their imprint on stellar oscillation frequencies. The authors modify the standard Eddington temperature–optical-depth relation by adding a two-parameter exponential term that mimics the pressure of a magnetic field, then compute p-mode frequencies for the evolved star HD 49385. The best-fitting models require a field of roughly 80 gauss concentrated about 1,850 km above the photosphere, and they simultaneously determine the star's helium core (0.117 solar masses, 0.078 solar radii) and its mass (1.25 ± 0.02 solar masses). If correct, the approach gives astronomers a way to measure magnetic fields on stars where direct magnetic detection is impossible.

What carries the argument

The key machinery is the modified Eddington T–τ relation, T^4 = (3/4) T_eff^4 [τ + q(τ)] with q(τ) = 2/3 + a exp(−bτ). The second term creates a temperature bump in the upper atmosphere, raising the radiation pressure to mimic the pressure of a small-scale magnetic field; the parameters a and b control the field strength and height. The oscillation calculation imposes the boundary condition P' = 0 at a layer called the magnetic-arch splicing layer, meaning p-modes are completely reflected there and no acoustic wave travels higher. The avoided-crossing l = 1 mode at 748.48 μHz is the central diagnostic for the helium core, because its frequency is tightly connected to the core's mass and radi

What would settle it

Direct spectropolarimetric observations of HD 49385 that detect no magnetic field of roughly 80 G near 1,850 km above the photosphere would falsify the magnetic interpretation. Alternatively, demonstrating that an arbitrary two-parameter surface correction, without any magnetic meaning, fits the same frequencies equally well would show the frequency data alone do not require magnetic fields.

Watch

Extended reading notes

Core claim

The central claim is that the frequency spectrum of HD 49385, observed with space-based photometry, contains the signature of small-scale magnetic fields in its atmosphere. By inserting a modified Eddington T–τ relation, q(τ) = 2/3 + a exp(−bτ), into the stellar model and using the boundary condition P' = 0 at the magnetic-arch splicing layer, the calculated p-mode frequencies for l = 0, 1, and 2 match the observed ones within about 2 μHz for most modes, with no empirical surface correction. The two best-fit models, based on different standard chemical compositions, both require small-scale magnetic fields of about 80 G concentrated near 1,850 km above the photosphere. Requiring the avoided-

Load-bearing premise

The load-bearing assumption is that the extra temperature bump in the modified Eddington T–τ relation really represents magnetic pressure, and that p-modes are completely reflected at the magnetic-arch splicing layer; if the bump is just an arbitrary two-parameter surface correction, the frequency fit alone does not prove the presence or strength of a magnetic field.

Editorial extensions

If this is right

  • Asteroseismology becomes a practical way to estimate small-scale magnetic field strengths and heights on remote solar-like stars where spectropolarimetric measurement is impossible.
  • The helium core of HD 49385 is determined as 0.117 solar masses and 0.078 solar radii, identical for two independent chemical compositions, confirming that the avoided-crossing frequency is a precise core diagnostic.
  • The stellar mass is measured to 1.25 ± 0.02 solar masses, with ages of 4.1 Gyr for one composition and 4.5 Gyr for the other.
  • Fits without any empirical surface correction suggest that at least part of the conventional 'surface effect' in solar-like oscillators may be magnetic in origin.
  • Two stars with similar effective temperatures (the Sun and HD 49385) yield similar inferred field strengths (~80–90 G), while the field height scales with the more extended atmosphere of the evolved star.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the method holds up, the same two-parameter T–τ modification could be applied to large samples of solar-like oscillators observed by other space missions, producing a first statistical map of small-scale magnetic fields in stellar atmospheres across evolutionary states.
  • The interpretation of the temperature bump as magnetic pressure is not uniquely forced by the frequency data; an independent test would be to compare the method's predictions against spectropolarimetric measurements on stars with known fields, or against 3D magneto-convection simulations.
  • Because the two adopted chemical compositions yield nearly identical field parameters (a and b) and identical helium cores, the magnetic-field inference appears robust to composition, but the method may still be degenerate with other near-surface temperature perturbations, such as chromospheric heating.
  • A natural follow-up is to see whether adding convective-core overshooting to the models changes the helium-core mass or mass/age estimates; the paper notes this is a future step.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 4 minor

Summary. The paper reports an asteroseismic analysis of the subgiant HD 49385 using CoRoT frequencies. The authors modify the Eddington T–τ relation by adding a phenomenological term, a·b·exp(−bτ), intended to mimic the pressure effect of small-scale magnetic fields in the atmosphere. They compute adiabatic p-mode frequencies for ℓ=0,1,2 with MESA/adipls and fit four parameters (magnetic term parameters a, b; initial mass M; mixing-length α) in a two-step grid search for two chemical compositions (GS98, A09). They obtain best-fit models that match the observed frequencies within ~2 μHz and reproduce spectroscopic Teff and log g. From the fitted a,b they convert to a magnetic-field strength of ~80 G at a height of ~1850 km using the equipartition condition β≈1. They also use the avoided-crossing mode to infer a helium-core mass of 0.117 M☉ and radius 0.078 R☉, finally giving the stellar mass as 1.25 ± 0.02 M☉ and age 4.1–4.5 Gyr.

Significance. If the magnetic-field interpretation is correct, the paper offers a practical seismic probe of small-scale magnetic fields in resolved and unresolved stars, extending the solar work of Li et al. (2021) to an evolved star. The frequency extraction and standard grid-based chi-square fitting are methodologically sound, and the use of frequency ratios r_ij that are insensitive to outer layers gives independent core constraints. However, the central claim—the presence of ~80 G fields at ~1850 km—is not uniquely established: the magnetic field is a reinterpretation of the fitted T–τ bump parameters, no comparison is made against a conventional surface-effect correction, and the assumed rigid reflection boundary condition is not physically justified. The helium-core and mass results are more robust because they rely on the avoided-crossing mode and frequency ratios, but they inherit the model-selection uncertainties of the two-step fitting.

major comments (4)
  1. [§3.2, Eq. (2) and §5] The magnetic-field claim is load-bearing and currently circular. The modified Eddington relation in Eq. (2) introduces arbitrary parameters a and b; the temperature bump is stated to 'simulate' magnetic pressure, but no MHD or radiative-transfer derivation connects a,b to a physical field. The final conversion to B≈80 G and height≈1850 km is made in §5 using β=P_gas/P_mag≈1, i.e., the measured quantities are direct functions of the fitted parameters. Moreover, no comparison is made against a standard surface-effect correction (e.g., Ball & Gizon 2014) on the same grid. A generic temperature perturbation may produce similar frequency shifts without magnetism. I request: (i) a quantitative comparison with a two-parameter surface correction, and (ii) a demonstration that the frequency residuals uniquely favour the magnetic interpretation (e.g., via a model selection criterion).
  2. [§3.3, Eq. (5)] The boundary condition P'=0 is a crucial assumption. It is asserted that p-modes are completely reflected at the 'magnetic-arch splicing layer', but this layer is nowhere defined, and no physical argument shows that all modes are perfectly reflected. The standard isothermal boundary condition (Eq. 3) is replaced without a sensitivity study. Since the upper boundary condition directly affects the eigenfrequencies, the improvement in chi2 may be an artefact of a more reflective boundary rather than a magnetic field. Please test the sensitivity of the best-fit frequencies to the choice of boundary condition (e.g., interpolate between Eq. (3) and Eq. (5)) and justify the reflection assumption from the magnetic-arch model.
  3. [§4] The two-step fitting procedure is not statistically justified. First, a and b are fixed using single values of M and α (M=1.286, α=1.85), with the assertion that a,b are insensitive to M and α; no supporting evidence is shown. Then M and α are searched with a,b fixed. This neglects possible degeneracies (e.g., a and M both shift large separation). The chi2 map in Fig. 2 shows multiple local minima, and the global minimum at a≈154,b≈49 is barely lower than the others; no confidence intervals on a,b are given. A simultaneous least-squares search over all four parameters (or at least a cross-check with several (M,α) pairs) is needed to claim a unique magnetic-field solution.
  4. [§4, Table 2; §5] The improvement from including the magnetic-term parameters is not quantified. The text states in the Conclusions that the models 'show a much better fit', but no chi2 or Δchi2 is given for models without the magnetic term on the same grid. With two extra parameters (a,b), one expects a lower chi2; a model-comparison statistic (e.g., AIC/BIC or F-test) is required to demonstrate that the improvement is significant. Also, the helium-core mass (0.117 M☉) and radius (0.078 R☉) are quoted without uncertainties; the selection based on ν_AC is not propagated into error bars. Please provide a formal uncertainty estimate.
minor comments (4)
  1. [Table 3] The DM11 effective temperature is listed as '66080 ± 60' rather than '6080 ± 60'; please correct this typo.
  2. [§1] There is a typographical error: 'B y matching' should read 'By matching'. Some equations in §3.3 also appear poorly typeset in the displayed version; please ensure they are legible.
  3. [Figure 4] The caption does not explain the meaning of the error bars on the observed triangles, and the ordinate label is missing in the displayed version. Please add a clear axis label and a note on the error bars.
  4. [§3.2] The statement that a and b 'specify how strong the magnetic fields are and where they are located' is only made explicit in §5. It would help the reader to define the mapping earlier and to state explicitly that this mapping is a model assumption.

Circularity Check

2 steps flagged · score 6.0 of 10

Magnetic field strength/height are a re-parametrization of the fitted T–τ bump; the frequency fit alone does not identify the field, and the magnetic interpretation rests on a self-cited phenomenological ansatz.

  1. fitted input called prediction [Sec. 3.2 (Eq. 2); Sec. 4 (a,b scan); Sec. 5 (B and height derivation)]
    "we introduce two physical parameters, a and b, to specify how strong the magnetic fields are and where they are located in the atmosphere. ... Consequently, we have determined the parameters of magnetic fields as a = 154 and b = 49 ... Model E with GS98 composition has small-scale magnetic fields of approximately 80 G and spreading concentratively at a height of about 1850 km in the photosphere."

    The reported 80 G and 1850 km are not independently measured quantities: they are computed from the same fitted parameters a and b via the assumed T–τ bump and the β≈1 equipartition condition. The frequency fit only constrains a and b; the paper then presents B and height as a detection, but this is a re-statement of the fit under a chosen physical interpretation. No separate observable constrains B or height, and no nonmagnetic surface-correction model is compared on the same grid, so the magnetic-field measurement does not stand independently of the model calibration.

  2. ansatz smuggled in via citation [Abstract; Sec. 3.2 (Eq. 2)]
    "We incorporate a modified Eddington T-τ equation that phenomenologically mimics the effect of the magnetic fields in the atmosphere ... Y. Li et al. (2021) introduced an additional term in the Hopf-like function ... This rise in temperature leads to an increase in the radiation pressure, which can be used to simulate magnetic pressure, and hence the presence of magnetic fields in the upper atmosphere."

    The representation of magnetic fields as a temperature bump is adopted, not derived, from the authors' prior work, and the paper itself calls it 'phenomenologically mimics.' Therefore the improved frequency fit does not by itself establish that the bump is magnetic; any surface temperature correction with similar shape could produce comparable frequency shifts. The solar calibration in Li et al. (2021) provides external plausibility, but for HD 49385 no independent magnetic measurement anchors the interpretation, so the argument reduces to accepting the self-cited ansatz.

full rationale

The asteroseismic fitting pipeline is internally consistent: observed frequencies are extracted, stellar models are built with MESA, and a χ² fit selects two best models with GS98 and A09 compositions. The helium-core mass and radius (0.117 M☉, 0.078 R☉) and the stellar mass (1.25±0.02 M☉) are derived in a standard way from fitting the avoided-crossing mode and the overall frequency set; those parts are not circular. The circularity concerns focus on the headline magnetic-field claim. The field strength and height are not fitted observables but are constructed from the fitted T–τ parameters a and b, which the paper itself defines as representing field strength and location. The conversion to gauss and kilometers uses an assumed β≈1 reflection layer and imports the P′=0 boundary condition from the authors' prior work. Because the modified Eddington term is explicitly phenomenological and no alternative surface-correction model is tested, the claim that HD 49385 hosts ~80 G fields at ~1850 km is only as strong as the self-cited ansatz, not independently verified by the frequency data. This is partial circularity rather than a complete derivation-from-inputs, so the score is 6 rather than higher.

Assumptions & free parameters 4 free parameters · 6 assumptions · 3 invented entities

The central new content is an inference from fitted parameters: a,b are fit to frequencies, and from them the paper derives B and height under an assumed reflection condition. The stellar mass and helium core are also inferred from the same frequency fit using standard stellar physics. The main input the reader 'pays for upstream' is the phenomenological equivalence between a temperature bump and magnetic pressure.

free parameters (4)
  • a = 154 (GS98), 156 (A09)
    Controls the amplitude of the added temperature/magnetic-pressure term in the Hopf-like function q(tau); fitted by grid search minimizing chi2 to observed frequencies.
  • b = 49 (GS98), 59 (A09)
    Controls the height/optical-depth location of the temperature bump; fitted jointly with a by chi2 grid search.
  • Initial mass M = 1.252 Msun (Model E), 1.250 Msun (Model J)
    Grid-searched over 1.225-1.315 Msun to minimize frequency chi2.
  • Mixing-length parameter alpha = 1.73 (GS98), 1.60 (A09)
    Grid-searched over 1.51-1.99; standard free parameter in stellar modeling.
assumptions (6)
  • domain assumption Eddington gray atmosphere T-tau relation with Hopf-like function q(tau) = 2/3 + a*b*exp(-b*tau)
    Used as the atmospheric boundary condition in MESA; the standard Eddington relation is accepted, but the added exponential term is ad hoc.
  • ad hoc to paper Increased radiation pressure from the temperature bump simulates magnetic pressure
    Section 3.2: the temperature rise 'can be used to simulate magnetic pressure, and hence the presence of magnetic fields.' This is a phenomenological identification, not derived from MHD.
  • ad hoc to paper Complete reflection of p-modes at the magnetic-arch splicing layer, giving P' = 0
    Section 3.3, Eq. 5: replaces the standard isothermal boundary condition, following Li et al. (2021). If this boundary condition is wrong, the fitted frequencies and inferred field are affected.
  • domain assumption Magnetic field strength is inferred by setting P_gas = P_mag (beta ~ 1) at the reflection layer
    Section 5, citing Rosenthal et al. (2002) and Cally (2007), assumes reflection/transmission at gas-pressure/magnetic-pressure equality. The exact field strength depends on this equality.
  • domain assumption Neglect of core overshooting, element diffusion, and a traditional surface-effect correction does not bias the inferred parameters
    The models omit overshooting and diffusion, and no comparison is made with standard surface corrections; these choices could shift the best-fit parameters.
  • domain assumption Initial helium abundance Y = 0.248 + 2Z
    Section 3.4, using the Dotter et al. (2008) enrichment relation with adopted constants; variation in this relation is not explored.
invented entities (3)
  • Temperature bump term a*b*exp(-b*tau) in the T-tau relation
    purpose: Phenomenologically mimics magnetic pressure in the upper atmosphere and changes p-mode frequencies
    The term is fitted to the observed frequencies; there is no independent measurement on HD 49385 that confirms it represents a magnetic field rather than another surface effect.
  • Magnetic-arch splicing layer with complete p-mode reflection
    purpose: Provides the mechanical boundary condition P' = 0 for oscillation calculations
    Adopted from Li et al. (2021); its existence and location on HD 49385 are not independently established.
  • ~80 G small-scale magnetic field at ~1850 km in HD 49385
    purpose: Reported as the physical interpretation of the fitted temperature-bump parameters
    The field strength and height are derived from fitted parameters and the beta ~ 1 assumption; no direct spectropolarimetric or other magnetic measurement of HD 49385 is presented.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Exploring the Small-scale Magnetic Fields in the Atmosphere of HD 49385 by Asteroseismic Analysis." pith.science (2026). https://pith.science/paper/52XE2XLB

@misc{pith2026260719724,
  author       = {Pith},
  title        = {Pith review of: Exploring the Small-scale Magnetic Fields in the Atmosphere of HD 49385 by Asteroseismic Analysis},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/52XE2XLB}},
  note         = {Machine review of arXiv:2607.19724}
}
abstract

Recent asteroseismic studies have shown convincing evidences that magnetic fields may exist in the interior of some pulsating red giants. Inspired by this breakthrough, we explored the effect of small-scale magnetic fields on the p-mode oscillations in an evolved star, HD 49385. {\bf We incorporate a modified Eddington $T$-$\tau$ equation that phenomenologically mimics the effect of the magnetic fields in the atmosphere of HD 49385,} and calculate the frequencies of p-modes with $l=0$, 1, and 2. By comparing the calculated frequencies with the observed ones, we select two best-fit models with either GS98 or A09 chemical composition. Our best-fit models not only fit satisfactorily the observed frequencies, but also well reproduce some spectroscopically observed stellar parameters such as effective temperature and log\,$g$. Based on the two best-fit models, we have estimated that the small-scale magnetic fields possess a strength of approximately 80\,G and spread concentratively at approximately a height of 1850 km in the atmosphere. By selecting the best-fit models with special requirement on the avoided-crossing mode, we have confirmed that the frequency of the avoided-crossing mode is tightly related to the helium core of the star, and determined the size of the helium core as 0.117${\rm M}_\odot$ in mass and 0.078${\rm R}_\odot$ in radius. Based on the improvements of previous two sides, we can accurately determine the mass of HD 49385 to be $1.25\pm 0.02\,{\rm M}_\odot$ with an age of 4.1\,Gyr for GS98 composition and 4.5\,Gyr for A09 composition.

Figures

Figures reproduced from arXiv: 2607.19724 by the authors.

Figure 1
Figure 1. Differences between frequencies extracted in this work (black symbols) and in S. Deheuvels et al. (2010; green symbols) for HD 49385. The gray horizontal short lines represent the error bar of the observed frequencies. 2 The Astrophysical Journal, 973:44 (8pp), 2024 September 20 Wang et al [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Variations of call 2 with different combinations of a and b values in Equation (2) for the best-fitting asteroseismic models under the assumption of M = 1.286Me and α = 1.85 with GS98 composition [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Comparisons between observed and theoretical frequencies for HD 49385. Red symbols represent the results of Model E with GS98, and blue symbols represent those of Model J with A09. The black ones with gray error bars are the observed frequencies [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figures from the paper (2 more)
Figure 5
Figure 5. Figure 5: The variations of gas pressure and magnetic pressure with height. Panel (a) is for Model E with GS98, and panel (b) for Model J with A09. 6 The Astrophysical Journal, 973:44 (8pp), 2024 September 20 Wang et al [PITH_FULL_IMAGE:figures/full_fig_p006_5.png]
Figure 6
Figure 6. Figure 6: The propagation of some mixed modes in the interior of HD 49385. The black line represents the square of buoyancy frequency N2 , the blue dotted line corresponds to the distribution of He abundance, and the remaining three solid lines represent the normalized kinetic e…

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

2 extracted references

  1. [1]

    Aerts, C., Christensen-Dalsgaard, J., & Kurtz, D. W. 2010, Asteroseismology (Berlin: Springer ) Aizenman, M., Smeyers, P., & Weigert, A. 1977, A&A, 58, 41 Asplund, M., Grevesse, N., & Sauval, A. J. 2005, in ASP Conf. Ser. 336, Cosmic Abundances as Records of Stellar Evolution and Nucleosynthesis, ed. I. Barnes, G. Thomas, & F. N. Bash (San Francisco, CA: ...

    2010

  2. [6]

    The propagation of some mixed modes in the interior of HD 49385. The black line represents the square of buoyancy frequency N2, the blue dotted line corresponds to the distribution of He abundance, and the remaining three solid lines represent the normalized kinetic energy for l = 1 modes with frequencies of 748.48, 777.91, and 828.22 μHz. 7 The Astrophys...

    2024

Pith tools

Reviewed August 1, 2026 · model on record in the stance chip above.