REVIEW 3 major objections 4 minor 31 references
The effect of stellar magnetic activity on measurements of morning and evening asymmetry of planetary terminator
T0 review · 3 major / 4 minor · reviewed 2026-08-06 · deepseek-v4-flash
Pith's one-line read Magnetic activity on a host star can create up to 600 ppm of false asymmetry in a 10,000 ppm transit light curve.
desk verdict The effect is real and would matter if the assumed limb magnetisation contrast holds, but the 600 ppm headline is not yet backed by data. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is the limb-brightness contrast of magnetized small-scale surface concentrations, computed from 3D radiative magnetohydrodynamic simulations of stellar atmospheres with ray-by-ray spectral synthesis. Magnetic fields corrugate the stellar surface and expose bright walls near the limb, making magnetic regions brighter there than the quiet photosphere, which flattens the limb darkening; the paper then combines half-disk light curves with different magnetization levels to obtain asymmetric transits. The second element of the machinery is the wavelength dependence: because the Planck function's sensitivity to temperature changes decreases toward the infrared, the magnetic asymmetry shrinks at longer wavelengths, and this spectral signature is the paper's proposed tool for separating stellar from planetary asymmetry in observations.
What would settle it
Measure ingress-egress depth differences over several transits of a solar-like star while simultaneously imaging its surface magnetic field (for example with Zeeman Doppler imaging), and check whether the sign and amplitude of the asymmetry at about 0.6 µm track the inferred difference in magnetic flux between the eastern and western near-limb regions; a clear absence of correlation would falsify the mechanism.
Extended reading notes
Core claim
The paper's central claim is that stellar magnetic activity, through its effect on limb darkening, can generate an ingress-egress asymmetry in transit light curves that is independent of any planetary atmosphere. The asymmetry is driven by small-scale magnetic field concentrations, which appear as bright patches in intergranular lanes and whose brightness contrast grows toward the limb; a difference in magnetization between the two sides of the star near the transit chord therefore makes one half of the transit deeper than the other. In a 10,000 ppm transit at 0.6 µm, the computed ingress-egress depth difference reaches up to 600 ppm, grows nearly linearly with the east-west magnetization difference, and does not depend strongly on the absolute magnetization level. The signal peaks in the visible and weakens toward the infrared, following the temperature sensitivity of the Planck function, and the paper shows for the WASP-39b case that this wavelength dependence is distinct enough from the planetary morning-evening signal to allow disentangling the two. The authors also note that the relevant small-scale fields often form without spots and without detectable photometric variability, so a photometrically quiet star is not automatically safe from this contamination.
Load-bearing premise
The quantitative amplitude rests on assuming that the simulated magnetized patches, drawn from atmosphere boxes roughly nine million meters across with imposed field strengths of 100 to 300 gauss, faithfully represent real stellar limb regions, and that the two sides of the star can differ in near-limb magnetization by several hundred gauss.
Editorial extensions
If this is right
- An ingress-egress asymmetry in a transit light curve can no longer be treated as purely planetary; a stellar contribution of hundreds of ppm must be considered in the visible.
- Photometrically quiet host stars are not exempt, since the relevant small-scale magnetic features often appear without spots and without detectable broad-band variability.
- The stellar contamination can either amplify or mute the apparent planetary morning-evening temperature contrast, depending on which stellar limb is more magnetized during a given transit.
- For the WASP-39b system, the observed 4.4 µm asymmetry is unlikely to be explained by stellar magnetic activity alone, but the stellar signal can be comparable to the planetary signal at visible and near-infrared wavelengths.
- Averaging over many transits reduces the activity signal roughly as the inverse square root of the number of transits, except for tilted non-equatorial orbits, where the chord crosses different stellar latitudes at ingress and egress and the asymmetry persists.
Reading between the lines
- If the computed amplitudes hold, retrieval codes for high-precision transmission spectroscopy should fit a stellar limb-asymmetry component jointly with the planetary signal rather than subtracting a disk-averaged stellar baseline.
- The same mechanism should affect any measurement that compares stellar flux on either side of a transit chord, so combined analysis of ingress-only and egress-only spectra could constrain the east-west magnetization difference of the host star.
- A natural extension is to map the effect across spectral type and activity level, since the paper's models are solar-like; M dwarfs, with their different facular contrasts, could show a different wavelength dependence that either helps or complicates the disentangling.
- Observing the same planet at multiple wavelengths and fitting the two components simultaneously would test the proposed disentangling directly, and transit timing variations could provide an additional constraint on the limb magnetization asymmetry.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper proposes that small-scale stellar magnetic fields, by brightening the stellar limb relative to the non-magnetic case, can introduce an asymmetry between the ingress and egress of a planetary transit light curve. Using MURaM magnetoconvection simulations and MPS-ATLAS spectral synthesis, the authors model a star with different magnetization levels on the eastern and western hemispheres and find ingress-egress depth differences up to about 600 ppm for a 10,000 ppm transit at 0.6 µm. They compare this stellar signal with the morning-evening terminator asymmetry expected for WASP-39b and argue that the two have distinct wavelength dependencies, which could be used to disentangle them.
Significance. The paper identifies a genuinely new contribution to transit asymmetry that is not captured by standard 1D stellar atmosphere models and that is distinct from the usual unocculted-spot contamination. The modeling is based on established 3D MHD simulations and ray-by-ray spectral synthesis, and the forward-model nature of the calculation means there is no fitting to the target result. The proposed wavelength-dependence discriminant is a useful, falsifiable prediction. However, the quantitative amplitude of the effect is highly sensitive to the assumed magnetization contrast between the two stellar limbs, which is not observationally calibrated in this work; this limits the immediate certainty of the 'significant interference' claim.
major comments (3)
- [Section 2, Figure 3] The claim that near-limb east-west field differences 'can readily reach several hundred Gauss' lacks any supporting magnetogram statistics or reference. Since the computed asymmetry is nearly linear in the magnetization difference (as seen in Figure 2b), the headline 600 ppm value scales directly with this assumed contrast. If the typical contrast is only a few tens of gauss, the signal would be reduced to roughly the 50 ppm level, which is below the planetary signals of interest. Please provide quantitative constraints from solar or stellar magnetograms, or re-frame the result as a parametric scaling with the unknown field contrast.
- [Section 2] The two-hemisphere model assigns a uniform magnetization to each half of the stellar disk, but the transit asymmetry is set by the field distribution along the ingress/egress chord near the limb. The paper does not demonstrate that a hemispheric average correctly represents the near-limb region, where the magnetic brightening is strongest. A localized active region of finite size near one limb would produce a different asymmetry amplitude than a uniformly magnetized hemisphere. Please test the sensitivity of the amplitude to the spatial scale and location of the magnetized region, or provide a justification for why the hemispheric approximation is adequate.
- [Sections 2 and 3] The amplitudes reported (e.g., 600 ppm) are not accompanied by any uncertainty estimates. The manuscript does not state how many MURaM snapshots or time averages are used for each magnetization level, nor how the results vary across realizations. Given that the 9x9 Mm simulation boxes are small and the magnetic field configuration evolves, a discussion of snapshot-to-snapshot variability is needed to assess the robustness of the predicted amplitude.
minor comments (4)
- [Section 2, Figure 1] The references to 'Figure 1, right panel' and 'Figure 1, left panel' appear to be swapped relative to the caption, which identifies the left panel as quiet (SSD) and the right panel as active (300 G). Please correct the cross-references.
- [Figure 4 caption] The caption states 'The orange and blue lines show planetary asymmetry', but the stellar curve is also described as blue, making the color assignment ambiguous. Please specify distinct line styles and colors for each curve.
- [Section 2] The phrase 'To pinpoint the asymmetry of these transit light curves' might be more clearly rendered as 'To isolate the asymmetry' or 'To display the asymmetry'.
- [Section 3] The reference to Frazier (1971) for the typical field strength in solar faculae is dated; consider citing more recent measurements if available, though this does not affect the conclusions.
Circularity Check
No significant circularity: forward model with independently grounded inputs; only a minor self-citation for the limb-brightening input.
full rationale
The paper's chain is a forward model: (i) MURaM/MPS-ATLAS simulations yield limb intensity profiles for quiet (SSD) and magnetized (100/200/300 G) stellar atmospheres; (ii) these profiles are assigned to eastern and western hemispheres; (iii) transit light curves are computed by integrating the stellar disk. The resulting ingress–egress asymmetry is the sensitivity of the model to the assumed magnetization contrast, not a fitted quantity. No equation is defined in terms of the target asymmetry; the 600 ppm value is conditional ('can reach up to ... depending on the magnetization'), and Section 2 explicitly notes the asymmetry increases almost linearly with the imposed field difference. The claim that near-limb fields can differ by several hundred Gauss is asserted rather than calibrated, and Section 2 concedes hemispheric contrasts 'may be rare'; this is a realism/calibration concern, not circularity, because the derivation does not presuppose the conclusion. The main self-citation, Kostogryz et al. (2024), supplies the limb-darkening input; that prior work is an independent MHD simulation study with published validation against solar observations, so relying on it is standard use of prior evidence rather than a self-referential reduction. The planetary comparison uses published climate models, not fitted outputs. No pattern of self-definitional reasoning, fitted-input-called-prediction, imported uniqueness, ansatz-smuggling via citation, or renaming of known results is present. Score 1 reflects the minor load-bearing self-citation for the input limb-brightening effect, which is nevertheless externally grounded.
Assumptions & free parameters
free parameters (3)
- magnetization difference between limbs =
0-300 G (SSD vs 300 G)
- transit depth normalization =
10,000 ppm (rp/R*=0.1)
- G9 stellar model parameters =
Teff=5250 K, logg=4.4, M/H=0.0
assumptions (5)
- domain assumption MURaM 3D radiative MHD simulations with box-in-a-star setup reproduce the magnetic effect on stellar limb darkening.
- domain assumption MPS-ATLAS ray-by-ray spectral synthesis provides accurate limb spectra for the simulated atmospheres.
- domain assumption Small-scale dynamo (SSD) simulations represent the quiet stellar surface.
- domain assumption Each stellar hemisphere can be represented by a single MURaM snapshot with a homogeneous mean vertical field.
- domain assumption The WASP-39b planetary asymmetry models of Carone et al. 2023 (with CH4 omitted) adequately represent the true planetary signal for the comparison.
Cite this review
Pith. "Pith review of The effect of stellar magnetic activity on measurements of morning and evening asymmetry of planetary terminator." pith.science (2026). https://pith.science/paper/CEWZ6GST
@misc{pith2026250716451,
author = {Pith},
title = {Pith review of: The effect of stellar magnetic activity on measurements of morning and evening asymmetry of planetary terminator},
year = {2026},
howpublished = {\url{https://pith.science/paper/CEWZ6GST}},
note = {Machine review of arXiv:2507.16451}
}
read the original abstract
Differences in the ingress and egress shapes of transit light curves can indicate morning-evening temperature contrasts on transiting planets. Here, we pinpoint an alternative mechanism that can introduce asymmetries in transit light curves, potentially affecting the accurate determination of morning-evening differences. Small-scale magnetic field concentrations on the surfaces of the host star affect the visibility of stellar limb regions, making them brighter relative to the non-magnetic case. A difference in magnetization between the star's western and eastern limbs can thus create an asymmetry in limb brightness and, consequently, an asymmetry between transit ingress and egress. We model the limb darkening and stellar limb asymmetry in solar-like stars using the 3D radiative MHD code MURaM to simulate magnetized stellar atmospheres and the MPS-ATLAS code to synthesize spectra using ray-by-ray approach. Our results show that ingress-egress depth differences can reach up to 600 ppm for a 10000 ppm transit at 600 nm, depending on the magnetization of the stellar limbs--significantly interfering with planetary signals. Observations of the Sun show that such concentrations are often not accompanied by spots and do not manifest in photometric variability, indicating that even photometrically quiet stars can produce such asymmetries. However, planetary and stellar asymmetries exhibit distinct wavelength dependencies, which we propose to leverage for disentangling them.
Figures
Reference graph
Works this paper leans on
-
[1]
2021, An Introduction to Stellar Magnetic Activity, doi: 10.1088/2514-3433/ac2956
Basri, G. 2021, An Introduction to Stellar Magnetic Activity, doi: 10.1088/2514-3433/ac2956
-
[2]
Bhatia, T. S., Cameron, R. H., Solanki, S. K., et al. 2022, A&A, 663, A166, doi: 10.1051/0004-6361/202243607
-
[3]
Buehler, D., Lagg, A., & Solanki, S. K. 2013, A&A, 555, A33, doi: 10.1051/0004-6361/201321152
-
[4]
2023, arXiv e-prints, arXiv:2301.08492, doi: 10.48550/arXiv.2301.08492
Helling, C. 2023, arXiv e-prints, arXiv:2301.08492, doi: 10.48550/arXiv.2301.08492
-
[5]
Danilovic, S., Sch¨ ussler, M., & Solanki, S. K. 2010, A&A, 513, A1, doi: 10.1051/0004-6361/200913379 de Wit, J., Gillon, M., Demory, B. O., & Seager, S. 2012, A&A, 548, A128, doi: 10.1051/0004-6361/201219060
-
[6]
Espinoza, N., Steinrueck, M. E., Kirk, J., et al. 2024, Nature, 632, 1017, doi: 10.1038/s41586-024-07768-4
-
[7]
Faedi, F., Barros, S. C. C., Anderson, D. R., et al. 2011, A&A, 531, A40, doi: 10.1051/0004-6361/201116671
-
[8]
Frazier, E. N. 1971, SoPh, 21, 42, doi: 10.1007/BF00155772
Show all 31 references
-
[9]
P., Mather, J
Gardner, J. P., Mather, J. C., Clampin, M., et al. 2006, Space Science Reviews, 123, 485, doi: 10.1007/s11214-006-8315-7
2006 doi
-
[10]
2019, A&A, 631, A79, doi: 10.1051/0004-6361/201935771
Helling, C., Iro, N., Corrales, L., et al. 2019, A&A, 631, A79, doi: 10.1051/0004-6361/201935771
2019 doi
-
[11]
2023, A&A, 671, A122, doi: 10.1051/0004-6361/202243956
Helling, C., Samra, D., Lewis, D., et al. 2023, A&A, 671, A122, doi: 10.1051/0004-6361/202243956
2023 doi
-
[12]
E., Shapiro, A
Keers, R. E., Shapiro, A. I., Kostogryz, N. M., et al. 2024, ApJL, 977, L7, doi: 10.3847/2041-8213/ad8b51
2024 doi
-
[13]
A., et al
Kiefer, S., Samra, D., Lewis, D. A., et al. 2024, A&A, 690, A244, doi: 10.1051/0004-6361/202450526
2024 doi
-
[14]
Kipping, D. M. 2016, MNRAS, 455, 1680, doi: 10.1093/mnras/stv2379
2016 doi
-
[15]
M., Shapiro, A
Kostogryz, N. M., Shapiro, A. I., Witzke, V., et al. 2024, Nature Astronomy, 8, 929, doi: 10.1038/s41550-024-02252-5
2024 doi
-
[16]
G., Steffen, M., & Freytag, B
Ludwig, H. G., Steffen, M., & Freytag, B. 2023, A&A, 679, A65, doi: 10.1051/0004-6361/202346783
2023 doi
-
[17]
1995, Nature, 378, 355, doi: 10.1038/378355a0 Molli` ere, P., Wardenier, J
Mayor, M., & Queloz, D. 1995, Nature, 378, 355, doi: 10.1038/378355a0 Molli` ere, P., Wardenier, J. P., van Boekel, R., et al. 2019, A&A, 627, A67, doi: 10.1051/0004-6361/201935470
1995 doi
-
[18]
2017, ApJ, 834, 151, doi: 10.3847/1538-4357/aa4f6c
Rackham, B., Espinoza, N., Apai, D., et al. 2017, ApJ, 834, 151, doi: 10.3847/1538-4357/aa4f6c
2017 doi
-
[19]
V., Apai, D., & Giampapa, M
Rackham, B. V., Apai, D., & Giampapa, M. S. 2018, ApJ, 853, 122, doi: 10.3847/1538-4357/aaa08c
2018 doi
-
[20]
V., Espinoza, N., Berdyugina, S
Rackham, B. V., Espinoza, N., Berdyugina, S. V., et al. 2023, RAS Techniques and Instruments, 2, 148, doi: 10.1093/rasti/rzad009
2023 doi
-
[21]
2012, Living Reviews in Solar Physics, 9, 1, doi: 10.12942/lrsp-2012-1
Reiners, A. 2012, Living Reviews in Solar Physics, 9, 1, doi: 10.12942/lrsp-2012-1
2012 doi
-
[22]
2014, ApJ, 789, 132, doi: 10.1088/0004-637X/789/2/132 —
Rempel, M. 2014, ApJ, 789, 132, doi: 10.1088/0004-637X/789/2/132 —. 2018, ApJ, 859, 161, doi: 10.3847/1538-4357/aabba0
2014 doi
-
[23]
Seager, S., & Shapiro, A. I. 2024, ApJ, 970, 155, doi: 10.3847/1538-4357/ad509a
2024 doi
-
[24]
P., & Guillot, T
Showman, A. P., & Guillot, T. 2002, A&A, 385, 166, doi: 10.1051/0004-6361:20020101
2002 doi
-
[25]
K., D´ esert, J
Sing, D. K., D´ esert, J. M., Lecavelier Des Etangs, A., et al. 2009, A&A, 505, 891, doi: 10.1051/0004-6361/200912776
2009 doi
-
[26]
K., Inhester, B., & Sch¨ ussler, M
Solanki, S. K., Inhester, B., & Sch¨ ussler, M. 2006, Reports on Progress in Physics, 69, 563, doi: 10.1088/0034-4885/69/3/R02
2006 doi
-
[27]
K., Krivova, N
Solanki, S. K., Krivova, N. A., & Haigh, J. D. 2013, ARA&A, 51, 311, doi: 10.1146/annurev-astro-082812-141007 V¨ ogler, A., Shelyag, S., Sch¨ ussler, M., et al. 2005, A&A, 429, 335, doi: 10.1051/0004-6361:20041507
2013 doi
-
[28]
B., Shapiro, A
Witzke, V., Duehnen, H. B., Shapiro, A. I., et al. 2023, A&A, 669, A157, doi: 10.1051/0004-6361/202244771
2023 doi
-
[29]
I., Kostogryz, N
Witzke, V., Shapiro, A. I., Kostogryz, N. M., et al. 2022, ApJL, 941, L35, doi: 10.3847/2041-8213/aca671
2022 doi
-
[30]
I., Cernetic, M., et al
Witzke, V., Shapiro, A. I., Cernetic, M., et al. 2021, A&A, 653, A65, doi: 10.1051/0004-6361/202140275
2021 doi
-
[31]
I., Kostogryz, N
Witzke, V., Shapiro, A. I., Kostogryz, N. M., et al. 2024, A&A, 681, A81, doi: 10.1051/0004-6361/202346099
2024 doi
Reviewed August 6, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.