pith. sign in

arxiv: 2606.02413 · v1 · pith:HRE5KIRLnew · submitted 2026-06-01 · 🌌 astro-ph.SR · astro-ph.IM

Modeling the Thermal Low-Frequency Radio Sun with Ray Tracing

Peijin Zhang , Bin Chen , Gregory Fleishman , Alexey Kuznetsov , Cooper Downs , Surajit Mondal , Sijie Yu This is my paper

Pith reviewed 2026-06-28 12:29 UTC · model grok-4.3

classification 🌌 astro-ph.SR astro-ph.IM
keywords solar radio emissionray tracingcoronal refractionfree-free emissionquiet Sunsynthetic radio mapspropagation effectsCarrington rotation
0
0 comments X

The pith

Ray tracing through 3D coronal models reproduces the quiet-Sun radio spectrum once refraction is included.

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

The authors build a forward-modeling method that traces radio rays through a global three-dimensional model of the corona while calculating how the plasma refracts the paths and changes the ray-tube area to keep flux conserved. They then solve the radiative transfer equation for thermal free-free emission and absorption along those rays to create synthetic brightness maps between 40 and 800 MHz. When the method is run on a full Carrington rotation, the quiet-Sun background spectrum matches observations only after the propagation effects are turned on. The same calculation under-predicts active-region brightness, indicating that other physics will be needed there.

Core claim

We develop a forward-modeling framework that combines refractive ray tracing through a global 3D coronal model with radiative transfer along each ray. The method tracks the ray-tube cross-sectional area S(s) using a step-wise perturbation retracing approach and incorporates a geometric magnification term proportional to d ln S / ds to enforce flux conservation under focusing/defocusing. Thermal free-free emission and absorption are then computed with the GRFF radiative transfer code to produce synthetic radio maps over 40--800 MHz. Applying the framework to Carrington rotation 2298, we find that including propagation effects allows the quiet-Sun background spectrum to be well reproduced.

What carries the argument

Ray-tube cross-sectional area S(s) tracked by step-wise perturbation retracing, together with the geometric magnification term d ln S / ds that enforces flux conservation during refraction.

If this is right

  • Synthetic radio maps generated this way become suitable for direct quantitative comparison with interferometric observations.
  • Propagation effects must be included to reproduce the observed quiet-Sun background at meter and decimeter wavelengths.
  • Active-region brightness remains under-predicted, so additional emission or absorption mechanisms are required for those features.
  • The framework can be used to test how refraction and focusing alter the apparent size and brightness distribution of the radio Sun.

Where Pith is reading between the lines

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

  • The same ray-tracing machinery could be applied to time-evolving coronal models to predict how radio morphology changes during solar rotation or flares.
  • Discrepancies remaining in active regions may point to the need for non-thermal electron populations that are not present in the current thermal-only calculation.
  • The method offers a route to forward-model radio images at frequencies where future arrays will provide high-resolution data, allowing direct tests of coronal density structure.

Load-bearing premise

The chosen global 3D coronal model supplies density and magnetic-field values accurate enough for the calculated refraction and free-free opacity to match the real quiet Sun.

What would settle it

A large mismatch between the modeled and observed quiet-Sun flux density spectrum across 40-800 MHz when the same framework is applied to an independent coronal density model would show that the refraction treatment does not capture reality.

Figures

Figures reproduced from arXiv: 2606.02413 by Alexey Kuznetsov, Bin Chen, Cooper Downs, Gregory Fleishman, Peijin Zhang, Sijie Yu, Surajit Mondal.

Figure 1
Figure 1. Figure 1: Ray-tube geometry and transverse basis; perturbed rays probe the cross-sectional area change. To form a stable transverse basis, we choose a reference axis a that is least aligned with ˆt i (e.g., a = zˆ unless |tˆz| ≈ 1, in which case a = yˆ), and set eˆ i 1 = a × ˆt i [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Plane of sky 2D slices through the PSI MAS model at z = 0, showing the electron density Ne, electron temperature Te, and magnetic-field strength |B| used for the ray-tracing calculations. (a) AIA 304 (b) AIA 171 (c) HMI magnetogram SDO context 2025-06-08T20:00 [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Context solar images from the Solar Dynamics Observatory (SDO). Here we use the MAS (Magnetohydrodynamic Algorithm outside a Sphere) model developed by Predictive Science Inc. (PSI) (Z. Miki´c & J. A. Linker 1996). PSI MAS is a global 3D, time-dependent MHD model of the solar corona [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Ray-tracing trajectories overlaid on an electron-density map from the PSI MAS coronal model. The 32 rays are launched from x = −3.5, y = 0, z = (−1.2, 1.2), colors indicate the variation of the magnification factor, positive indicates defocusing, negative indicates focusing. to focusing (d ln S/ds < 0) and positive (red) segments correspond to defocusing (d ln S/ds > 0) [PITH_FULL_IMAGE:figures/full_fig_p… view at source ↗
Figure 5
Figure 5. Figure 5: Concept figure for the ray-coordinate convention. The ray coordinate increases along the traced ray toward the Sun (observer to source), which is opposite to the radiative-transfer integration direction (source to observer). 0 1 2 3 4 Distance along ray (R_sun) 0.95 0.96 0.97 0.98 0.99 1.00 1.01 Magnification Factor S ray z0: -1.12 ray z0: -0.75 ray z0: -0.38 ray z0: 0.00 ray z0: 0.38 ray z0: 0.75 ray z0: … view at source ↗
Figure 6
Figure 6. Figure 6: Magnification factor (S) along the ray preceding toward the Sun. Upper panel shows the source area S along the ray path (equavlent to accumulated local magnification factor); values > 1 indicate defocusing and values < 1 indicate focusing along the ray, lower panel shows the variation of the magnification factor of rays along the ray path, the rays are lauched from x = −3.5, y = 0, z = (−1.2, 1.2) [PITH_F… view at source ↗
Figure 7
Figure 7. Figure 7: Synthetic quiet-Sun radio images at 40, 80, and 150 MHz computed from the PSI MAS coronal model with refractive ray tracing and thermal radiative transfer. The model-resolution images are convolved with a representative interferometric beam corresponding to a 15 km maximum baseline. The top row panels present the result with ray tracing, the bottom row panels present the result without ray tracing. Figures… view at source ↗
Figure 8
Figure 8. Figure 8: Synthetic quiet-Sun radio images at 280, 550, and 800 MHz computed from the PSI MAS coronal model using refractive ray tracing and free–free radiative transfer. The intrinsic images are convolved with a representative interferometric beam corresponding to a 5 km maximum baseline. The dotted circle marks the solar optical limb. paths sample a different coronal volume than the physically refracted trajectori… view at source ↗
Figure 9
Figure 9. Figure 9: Effect of the ray-tube magnification term in the radiative-transfer calculation. Shown are synthetic quiet-Sun brightness maps at 60 MHz computed with refractive ray tracing from the same coronal model, comparing solutions that include and omit the geometric magnification factor proportional to d ln S/ds. 3.4. Comparison with observations 3.4.1. Brightness temperature spectrum The near-center (< 0.5R⊙) qui… view at source ↗
Figure 10
Figure 10. Figure 10: Quiet-Sun brightness temperature spectrum. The solid curve shows the modeled < 0.5R⊙ averaged Tb derived from our synthetic images over 40–800 MHz. Symbols show representative measurements from previous quiet-Sun observations in the literature (e.g., K. R. Subramanian 2004; R. Ramesh et al. 2006; V. N. Melnik et al. 2018; P. Zhang et al. 2022) [PITH_FULL_IMAGE:figures/full_fig_p009_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Comparison between observed and modeled Sun brightness at decameter wavelengths. Panels (a1–a3) show OVRO–LWA quiet-Sun images from 2025-06-08T20:07:03 UT at 38.7, 52.4, and 84.6 MHz. Panels (b1–b3) show the corre￾sponding synthetic images computed with refractive ray tracing and GRFF radiative transfer using the PSI/MAS coronal model for Carrington rotation 2298, rotated to align the coronal configuratio… view at source ↗
Figure 12
Figure 12. Figure 12: Example image-level comparison between OVRO–LWA observations and the ray-tracing MAS–GRFF model at ∼66 MHz. (a) OVRO–LWA quiet-Sun brightness temperature map at 66.0 MHz. (b) Corresponding synthetic brightness temperature map from the model at 66.3 MHz. (c) Normalized difference map, (Iobs − Imodel)/Imax obs , shown with a diverging color scale. The black circle marks the solar disk, and the dashed lines … view at source ↗
read the original abstract

Incoherent radio emission at meter--decimeter wavelengths provides a key diagnostic of the coronal thermal plasma, but at frequencies below $\sim$\,1\,GHz coronal refraction can substantially bend ray paths and modify the apparent source size and brightness distribution. We develop a forward-modeling framework that combines refractive ray tracing through a global 3D coronal model with radiative transfer along each ray. The method tracks the ray-tube cross-sectional area $S(s)$ using a step-wise perturbation retracing approach and incorporates a geometric magnification term proportional to $d\ln S/ds$ to enforce flux conservation under focusing/defocusing. Thermal free--free emission and absorption are then computed with the \texttt{GRFF} radiative transfer code to produce synthetic radio maps over 40--800\,MHz. Applying the framework to Carrington rotation 2298, we find that including propagation effects allows the quiet-Sun background spectrum to be well reproduced. However, active region brightness is less accurately modeled, suggesting that additional physical factors should be considered in future work. These results establish a physics-based method for generating low-frequency quiet-Sun synthetic images suitable for quantitative comparison with interferometric observations and for assessing how propagation effects shape the observed morphology.

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

2 major / 2 minor

Summary. The paper develops a forward-modeling framework combining refractive ray tracing through an external global 3D coronal model with radiative transfer via the GRFF code. Ray-tube cross-sectional area S(s) is tracked via step-wise perturbation retracing, with a geometric magnification term d ln S/ds enforcing flux conservation. Thermal free-free emission and absorption are computed to generate synthetic maps at 40-800 MHz. Applied to Carrington rotation 2298, the framework is reported to reproduce the quiet-Sun background spectrum once propagation effects are included, while active-region brightness is less accurately modeled. The work positions itself as establishing a physics-based method for synthetic quiet-Sun images rather than claiming universal accuracy.

Significance. If the reproduction result holds with quantitative support, the framework supplies a reproducible, parameter-free (within the adopted 3D model) approach to generating low-frequency radio maps that incorporate refraction and free-free effects. This is useful for direct comparison with interferometric data and for isolating how propagation shapes observed morphology. The method builds directly on existing tools (global coronal model + GRFF) without introducing new free parameters or ad-hoc tuning, which strengthens its utility as a forward-modeling tool.

major comments (2)
  1. [Results (CR 2298 application)] Results section (application to CR 2298): the central claim that 'including propagation effects allows the quiet-Sun background spectrum to be well reproduced' lacks accompanying quantitative metrics (e.g., reduced chi-squared, RMS residuals, or direct spectral comparison with error bars). Without these, it is not possible to evaluate whether the match is within observational uncertainties or merely qualitative.
  2. [Method (ray tracing)] Method section on ray-tube tracking: the step-wise perturbation retracing for S(s) and the d ln S/ds magnification term are described at a high level, but no explicit test (e.g., conservation check on a known analytic case or comparison against full ray-bundle integration) is provided to confirm numerical accuracy of the flux-conservation implementation.
minor comments (2)
  1. [Abstract] Abstract: the phrasing 'well reproduced' is qualitative; replace with a brief statement of the level of agreement once quantitative metrics are added in the main text.
  2. [Notation] Notation: ensure consistent use of S(s) for ray-tube area throughout; a short table summarizing symbols would aid readability.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive review and the recommendation of minor revision. We agree that both major comments identify areas where the manuscript can be strengthened with additional quantitative support and validation. We address each point below and will incorporate the suggested changes in the revised version.

read point-by-point responses

  1. Referee: Results section (application to CR 2298): the central claim that 'including propagation effects allows the quiet-Sun background spectrum to be well reproduced' lacks accompanying quantitative metrics (e.g., reduced chi-squared, RMS residuals, or direct spectral comparison with error bars). Without these, it is not possible to evaluate whether the match is within observational uncertainties or merely qualitative.

    Authors: We agree that the current presentation of the CR 2298 results relies primarily on visual comparison of spectra. In the revised manuscript we will add quantitative metrics, including RMS residuals between the modeled and observed quiet-Sun spectra and a direct spectral plot that includes observational error bars, to demonstrate that the reproduction lies within the reported uncertainties. revision: yes

  2. Referee: Method section on ray-tube tracking: the step-wise perturbation retracing for S(s) and the d ln S/ds magnification term are described at a high level, but no explicit test (e.g., conservation check on a known analytic case or comparison against full ray-bundle integration) is provided to confirm numerical accuracy of the flux-conservation implementation.

    Authors: We acknowledge that an explicit numerical validation test is absent. In the revised methods section (or a new appendix) we will include a conservation check performed on a spherically symmetric analytic corona model, comparing the step-wise perturbation results against the expected analytic flux conservation to confirm the accuracy of the ray-tube area tracking and magnification term. revision: yes

Circularity Check

0 steps flagged

No significant circularity identified

full rationale

The paper's central result is a forward-modeling application of ray tracing plus the external GRFF code to an independent 3D coronal model for CR 2298. The reproduction of the quiet-Sun spectrum follows directly from the physics-based propagation and radiative-transfer calculation; no equation or parameter is fitted inside the paper and then relabeled as a prediction, no self-citation supplies a load-bearing uniqueness theorem, and no ansatz is smuggled via prior work by the same authors. The derivation chain is therefore self-contained against external inputs.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract-only review yields no explicit free parameters, axioms, or invented entities; the central claim rests on the unstated accuracy of the input 3D coronal model and the validity of the ray-tube perturbation method.

pith-pipeline@v0.9.1-grok · 5764 in / 1117 out tokens · 19626 ms · 2026-06-28T12:29:12.062961+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

20 extracted references · 19 canonical work pages

  1. [1]

    E., Bastian, T

    Alissandrakis, C. E., Bastian, T. S., & Nindos, A. 2023, A&A, 670, C5, doi: 10.1051/0004-6361/202243774e

    work page doi:10.1051/0004-6361/202243774e 2023

  2. [2]

    Bastian, T. S. 1994, ApJ, 426, 774, doi: 10.1086/174114

    work page doi:10.1086/174114 1994

  3. [3]

    Dulk, G. A. 1985, ARA&A, 23, 169, doi: 10.1146/annurev.aa.23.090185.001125

    work page doi:10.1146/annurev.aa.23.090185.001125 1985

  4. [4]

    D., Kuznetsov, A

    Fleishman, G. D., Kuznetsov, A. A., & Landi, E. 2021, ApJ, 914, 52, doi: 10.3847/1538-4357/abf92c

    work page doi:10.3847/1538-4357/abf92c 2021

  5. [5]

    P., Chen , X., Chrysaphi , N., et al

    Kontar, E. P., Chen, X., Chrysaphi, N., et al. 2019, ApJ, 884, 122, doi: 10.3847/1538-4357/ab40bb

    work page doi:10.3847/1538-4357/ab40bb 2019

  6. [6]

    A., Fleishman, G

    Kuznetsov, A. A., Fleishman, G. D., & Landi, E. 2026, Codes for computing the solar gyroresonance and free-free radio emissions, v1.0.1, Zenodo, doi: 10.5281/zenodo.19368163

    work page doi:10.5281/zenodo.19368163 2026

  7. [7]

    A., & Miki´ c, Z

    Lionello, R., Linker, J. A., & Miki´ c, Z. 2009, ApJ, 690, 902, doi: 10.1088/0004-637X/690/1/902

    work page doi:10.1088/0004-637x/690/1/902 2009

  8. [8]

    N., Shepelev, V

    Melnik, V. N., Shepelev, V. A., Poedts, S., et al. 2018, Solar Physics, 293, doi: 10.1007/s11207-018-1316-3

    work page doi:10.1007/s11207-018-1316-3 2018

  9. [9]

    2015, Astronomy & Astrophysics, 583, doi: 10.1051/0004-6361/201425540 Miki´ c, Z., & Linker, J

    Mercier, C., & Chambe, G. 2015, Astronomy & Astrophysics, 583, doi: 10.1051/0004-6361/201425540 Miki´ c, Z., & Linker, J. A. 1996, in American Institute of Physics Conference Series, Vol. 382, American Institute of Physics Conference Series, 104–107, doi: 10.1063/1.51370

    work page doi:10.1051/0004-6361/201425540 2015

  10. [10]

    B., Howard, R

    Mondal, S., Shaik, S. B., Howard, R. A., et al. 2026, ApJ, 999, 237, doi: 10.3847/1538-4357/ae4353

    work page doi:10.3847/1538-4357/ae4353 2026

  11. [11]

    S., Kathiravan, C., & Sastry, C

    Ramesh, R., Nataraj, H. S., Kathiravan, C., & Sastry, C. V. 2006, The Astrophysical Journal, 648, 707, doi: 10.1086/505677

    work page doi:10.1086/505677 2006

  12. [12]

    2020, The Astrophysical Journal, 903, doi: 10.3847/1538-4357/abb949

    Sharma, R., & Oberoi, D. 2020, The Astrophysical Journal, 903, doi: 10.3847/1538-4357/abb949

    work page doi:10.3847/1538-4357/abb949 2020

  13. [13]

    E., & Pohjolainen, S

    Shibasaki, K., Alissandrakis, C. E., & Pohjolainen, S. 2011, SoPh, 273, 309, doi: 10.1007/s11207-011-9788-4

    work page doi:10.1007/s11207-011-9788-4 2011

  14. [14]

    1971, A&A, 10, 362

    Boischot, A. 1971, A&A, 10, 362

    1971

  15. [15]

    Subramanian, K. R. 2004, Astronomy & Astrophysics, 426, 329, doi: 10.1051/0004-6361:20047120

    work page doi:10.1051/0004-6361:20047120 2004

  16. [16]

    J., & Gopalswamy, N

    Thejappa, G., MacDowall, R. J., & Gopalswamy, N. 2011, ApJ, 734, 16, doi: 10.1088/0004-637X/734/1/16

    work page doi:10.1088/0004-637x/734/1/16 2011

  17. [17]

    2018, A&A, 614, A54, doi: 10.1051/0004-6361/201630067

    Vocks, C., Mann, G., Breitling, F., et al. 2018, A&A, 614, A54, doi: 10.1051/0004-6361/201630067

    work page doi:10.1051/0004-6361/201630067 2018

  18. [18]

    2024, Nature Astronomy, 8, 50, doi: 10.1038/s41550-023-02122-6

    Yu, S., Chen, B., Sharma, R., et al. 2024, Nature Astronomy, 8, 50, doi: 10.1038/s41550-023-02122-6

    work page doi:10.1038/s41550-023-02122-6 2024

  19. [19]

    2022, The Astrophysical Journal, 932, doi: 10.3847/1538-4357/ac6b37

    Zhang, P., Zucca, P., Kozarev, K., et al. 2022, The Astrophysical Journal, 932, doi: 10.3847/1538-4357/ac6b37

    work page doi:10.3847/1538-4357/ac6b37 2022

  20. [20]

    M., & Hurford, G

    Zirin, H., Baumert, B. M., & Hurford, G. J. 1991, The Astrophysical Journal, 370, 779, doi: 10.1086/169861

    work page doi:10.1086/169861 1991