Pith. sign in

REVIEW 3 major objections 5 minor 40 references

A 3D Monte Carlo calculation of the inverse Compton emission from the Sun and stars in presence of magnetic and electric fields

T0 review · 3 major / 5 minor · reviewed 2026-08-07 · deepseek-v4-flash

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

Pith's one-line read Magnetic fields make the Sun's gamma-ray halo brighter near the Sun.

desk verdict A credible 3D Monte Carlo treatment of solar inverse Compton emission with magnetic fields, but the headline enhancement is conditional on a simplified field model and an unseparated electric-field effect. read the letter →

arxiv 2505.18275 v1 pith:SNJN4WME submitted 2025-05-23 hep-ph astro-ph.HEastro-ph.SRhep-ex

M. N. Mazziotta This is my paper · ORCID
classification hep-phastro-ph.HEastro-ph.SRhep-ex
keywords inverseComptonscatteringsolargamma-rayhalo3DMonteCarlosimulationinterplanetarymagneticfieldcosmic-rayelectronsstellarblack-bodyphotonemissionheliosphere
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

The paper develops a general 3D Monte Carlo method for computing inverse Compton gamma-ray emission from cosmic-ray electrons scattering off stellar photons, with the option of including magnetic and electric fields that bend the electron trajectories. The method is applied to the Sun using the spiral interplanetary magnetic field and its associated electric field. The central result is a predicted enhancement of the inverse Compton halo near the Sun: in the field configuration, electrons spend more time in the solar photon field before escaping, so the gamma-ray emission grows brighter and more sharply peaked toward the Sun than standard field-free line-of-sight calculations predict. A correct treatment of this effect matters for interpreting solar gamma-ray observations and for separating the inverse Compton halo from the disk emission. Without a magnetic field, the simulation reproduces existing line-of-sight inverse Compton results.

What carries the argument

The engine of the calculation is a Monte Carlo transport loop built around an effective cross section $\sigma(E_e, r, \theta_e)$ that describes an electron of energy $E_e$ at distance $r$ from the star, moving at polar angle $\theta_e$ through the star's black-body photon field. The cross section integrates the Klein-Nishina scattering probability over the visible stellar surface, including the Lambert cosine emission factor and the $(1-\beta\cos\zeta)$ flux factor, and its inverse sets the local interaction length. Electrons step through the field, sampling photon energies and directions from the stellar black-body spectrum; each scattering is simulated in the electron rest frame and boosted back to the observer frame. When magnetic and electric fields are present, trajectories are advanced with a helix approximation and an adaptive Runge-Kutta method, and the electric potential changes the electron energy between steps. Secondary photons are collected on a detection sphere and binned in energy and angle to produce sky maps and spectra.

What would settle it

Compare the predicted >100 MeV inverse Compton halo profile between 0.5 and 5 degrees from the Sun with high-energy gamma-ray observations: if the observed profile shows no extra brightening or steepening relative to the field-free line-of-sight model, the magnetic-field enhancement is falsified; repeating the simulation with a realistic inner-heliospheric field is the direct numerical check.

Watch

Extended reading notes

Core claim

The paper claims that inverse Compton scattering in the solar environment cannot be fully represented by straight-line electron paths if one wants accurate gamma-ray maps: once the spiral interplanetary magnetic field and the corotating electric field are included, the 3D Monte Carlo simulation yields a solar inverse Compton halo that is brighter and more peaked toward the Sun than the field-free line-of-sight calculation. The paper attributes this to the magnetic field increasing the time electrons spend in the dense solar photon field, raising the scattering probability near the Sun. In the no-field limit the simulation agrees with the standard line-of-sight inverse Compton intensity, while in the field case the integrated flux above 100 MeV lies above the no-field prediction for angular distances below about 4 degrees and falls off with a steeper angular profile.

Load-bearing premise

The prediction rests on a simplified model of the Sun's magnetic field, a spiral with a constant 400 km/s solar wind, zero tilt, and no strong irregular inner field near the Sun, so a more realistic near-Sun field could change or erase the brightening.

Editorial extensions

If this is right

  • Under the spiral-field model, the Sun's inverse Compton gamma-ray flux above 100 MeV within roughly 4 degrees of the Sun should exceed the field-free line-of-sight prediction and decline with a steeper angular profile.
  • Without a magnetic field, the 3D Monte Carlo method reproduces existing line-of-sight inverse Compton calculations, giving a benchmark for the method away from the Sun.
  • The predicted inverse Compton flux near the Sun is higher at photon energies below 100 MeV than earlier line-of-sight estimates, which is relevant to low-energy solar gamma-ray data.
  • The same machinery can be applied to other stars and can be extended to gamma-gamma absorption by replacing the Compton cross section with the pair-production cross section in the same geometric integral.
  • The simulation produces full spatial maps of the inverse Compton halo, not only spectra, so the predicted near-Sun brightening can be compared directly with observed count maps.

Reading between the lines

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

  • Beyond the paper, a stronger near-Sun inverse Compton halo would mean that observed solar gamma-ray maps may require a larger inverse Compton component close to the Sun, which would reduce the inferred hadronic disk component at those angles.
  • Beyond the paper, the same 'time spent in the radiation field' mechanism implies that stars with stronger magnetic environments, such as active stars or stars with dense winds, could show anomalously peaked inverse Compton halos even with ordinary cosmic-ray fluxes; this is a testable prediction for future stellar gamma-ray surveys.
  • Beyond the paper, a natural next test is to run the simulation with a nonzero current-sheet tilt and reversed magnetic polarity: the enhancement should become asymmetric and vary over the solar cycle, which can be checked against long-term solar gamma-ray observations.
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

3 major / 5 minor

Summary. This paper presents a 3D Monte Carlo framework for anisotropic inverse Compton scattering in stellar photon fields, optionally including magnetic and electric fields. The cross-section formalism is standard Klein-Nishina scattering sampled by boosting to the electron rest frame. The author validates the sampler against the Jones (1968) formula for isotropic blackbody photons and, for the no-field solar halo, against the StellarICS line-of-sight calculation. With an interplanetary Parker-spiral magnetic field and its associated motional electric field, the simulation predicts a brighter, more Sun-peaked IC halo. The author interprets this as an increase of the time spent by electrons in the radiation field and notes that a more accurate solar magnetic configuration is needed for confirmation.

Significance. If the enhancement is robust, the paper provides a falsifiable prediction for the solar gamma-ray halo and a reusable numerical method for IC emission in complex geometries, with the no-field comparisons giving internal validation. The explicit treatment of anisotropic photon fields and particle trajectories in electromagnetic fields goes beyond the line-of-sight integrations currently used. However, the central claim is conditional on the assumed Parker-spiral geometry and on separating magnetic trapping from electric acceleration; the significance is therefore contingent on the controlled experiments and sensitivity studies requested below.

major comments (3)
  1. [Sec. VII, Eqs. (44)-(50); Sec. VIII] The simulation includes both the Parker-spiral B field and the motional E field with potential V = -/+ f B0 omega_S R_E^2 z/r (Eq. 49). The potential difference across the simulation volume is of order B0 omega_S R_E^2 ~ 3x10^8 V, so electrons of a few GeV can gain or lose hundreds of MeV along their trajectories. Such energy changes directly alter the IC photon spectrum and the near-Sun intensity. The paper attributes the enhancement to 'an increase of the time spent by electrons in the radiation field in presence of the magnetic field' (Sec. VIII), but no configuration with B only (E=0) is run, and no diagnostic separates energy gain from path-length increase. The author should run B-only and E-only controlled simulations, or explicitly quantify the electric-field contribution, before the residence-time interpretation is stated.
  2. [Sec. VII, Eqs. (44)-(47); Sec. VIII] The predicted enhancement is obtained with a null tilt angle, a constant solar-wind speed of 400 km/s, B0 = 5 nT, and no irregular inner heliospheric field. The author himself flags in Sec. VIII that confirmation requires 'a more accurate solar magnetic configuration that also includes the magnetic field near the Sun.' Because the E-field sign structure and the particle trajectories are sensitive to the tilt angle and heliospheric current-sheet crossings, and because the inner field affects trapping, the central result needs at least a sensitivity study over tilt angle, v_SW, and B0, or a clearly stated exploratory status, before it can support the paper's conclusion.
  3. [Sec. VI-VII, trajectory integration] The B-field runs rely on a helix approximation and an adaptive Runge-Kutta-Nystroem algorithm, and the no-B runs use the smax step limiter of Sec. IV; however, no convergence test is reported for the B-field trajectories (step size, smax, tolerance), and there is no check that the energy variation from Eq. (50) is independent of step size or that energy is conserved when E=0. Since the near-Sun enhancement is a numerical prediction, a step-size convergence test is needed to rule out integration artifacts; this is a load-bearing validation for the central claim.
minor comments (5)
  1. [Sec. VIII] The sentence 'we have not implemented any approximation' is overbroad; the simulation uses smax step limiting, the helix approximation, Runge-Kutta tracking, and finite energy bins. The statement should be qualified to refer to the cross-section treatment.
  2. [Sec. V, after Eq. (35)] The phrase 'the star is seen as a half emitting surface' at large distances is geometrically incorrect; at r >> R the star subtends a small solid angle. Please rephrase, since the prefactor in Eq. (35) already gives the correct dilution.
  3. [Fig. 10 and Sec. VII] The text and figure labels alternate between 'StellaICS' and 'StellarICS' (e.g., Fig. 10 middle and bottom panels); use one spelling consistently.
  4. [Sec. VII] Please specify the smax and trajectory-integration step parameters used in the magnetic-field runs, and state how statistical errors of the with/without-field ratio are computed.
  5. [Sec. VIII] The phrases 'the an increase' and 'we have not implemented any approximation' contain grammatical or logical errors; please proofread the conclusion.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the B-field enhancement is a forward Monte Carlo result benchmarked externally.

full rationale

The paper's central new result is the 3D Monte Carlo prediction that adding the interplanetary magnetic and electric fields makes the solar inverse Compton emission brighter and more concentrated toward the Sun. This is a forward simulation output, not a fitted quantity: the no-field configuration is benchmarked against the Jones 1968 formula and the external StellarICS code, and the Fermi-LAT, AMS-02, and DAMPE spectra enter only as input electron spectra, not as fit targets. The Parker-spiral magnetic field is adopted from Parker 1958 and the electric potential from Lipari 2014, both external references rather than self-citations. The only self-reference is the use of the author's previous work [3] for extrapolating the observed CRE spectra down to 0.1 GeV and up to 10 TeV, and this is an input assumption that does not by construction produce the near-Sun enhancement. The Sec. VIII conclusion is explicitly conditional, noting that the enhancement 'will be confirmed considering a more accurate solar magnetic configuration,' which is a stated limitation rather than a circular justification. No equation in the paper defines the predicted enhancement in terms of the fitted inputs, and no fitted parameter is relabeled as a prediction. The attribution of the effect to increased residence time in the radiation field, rather than to electric-field energy changes, is a physical interpretation that may be debatable but is not circularity. Overall, the derivation chain is self-contained against external benchmarks, so the appropriate finding is no significant circularity.

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

The central claim rests on the standard QED cross sections, the Parker-spiral field model, and the assumption that local cosmic-ray measurements represent the flux over the whole generation sphere. The free parameters are nominal heliospheric values chosen by hand, not fit to the target result. No new particles or physical entities are introduced.

free parameters (3)
  • B0 (IMF intensity at Earth) = 5 nT (nominal typical value)
    Chosen by hand for the simulation; the predicted B-field enhancement scales with this input, and no uncertainty is propagated.
  • v_SW (solar wind speed) = 400 km/s (nominal typical value)
    Chosen by hand for the Parker spiral; the shape of the spiral and the electron path length depend on it.
  • HCS tilt angle alpha = 0 (null tilt)
    Set to zero for simplicity; the paper says a more accurate solar magnetic configuration would be needed to confirm the effect, so this choice directly affects the central B-field result.
assumptions (6)
  • standard math Klein-Nishina differential cross-section and special-relativistic kinematics govern the photon-electron scattering.
    Used in Secs. II and III as the basis of the Monte Carlo sampling.
  • domain assumption The interplanetary magnetic field is described by the Parker spiral model with a constant radial solar wind.
    Invoked in Eqs. 44-47 in Sec. VII; this is the load-bearing modeling choice for the B-field result.
  • domain assumption The heliospheric electric field is the gradient of the Lipari potential, vanishing at the equator.
    Eqs. 48-49 in Sec. VII; the energy change of electrons is computed from this potential and contributes to the enhanced near-Sun emission.
  • domain assumption The cosmic-ray electron plus positron intensity at the 1 AU generation sphere equals the locally measured Fermi-LAT, AMS-02, and DAMPE spectra, extrapolated to 0.1 GeV to 10 TeV.
    Sec. VII; the absolute flux predictions depend on this, and spatial variations due to solar modulation are ignored.
  • domain assumption The strong, irregular magnetic field in the inner heliosphere near the Sun is neglected.
    Stated in Sec. VII: 'we are not considering the strong and irregular inner magnetic field near the Sun'; if this field matters, the predicted enhancement could change.
  • domain assumption Monte Carlo statistical errors and finite bin widths are the only quantified uncertainties; no systematic uncertainties are propagated.
    The error bars in Fig. 10 are statistical only; systematic uncertainties on B0, solar wind speed, and CRE spectra are not assessed.

how reviews work

0 comments
Cite this review

Pith. "Pith review of A 3D Monte Carlo calculation of the inverse Compton emission from the Sun and stars in presence of magnetic and electric fields." pith.science (2026). https://pith.science/paper/SNJN4WME

@misc{pith2026250518275,
  author       = {Pith},
  title        = {Pith review of: A 3D Monte Carlo calculation of the inverse Compton emission from the Sun and stars in presence of magnetic and electric fields},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/SNJN4WME}},
  note         = {Machine review of arXiv:2505.18275}
}
read the original abstract

The solar steady emission in gamma rays is due to the interactions of Galactic cosmic rays with the solar atmosphere and with the low-energy solar photon field via inverse Compton scattering. The emission is sensitive to the magnetic field nearby the Sun and to the cosmic-ray transport in the magnetic field in the inner solar system. Modeling the inverse Compton emission in the presence of a magnetic field is therefore crucial to better interpret the observations. In a previous work we have presented a comprehensive calculation of the secondary productions due to the collision of cosmic rays with the solar atmosphere in presence of magnetic fields. In this paper, we present a general approach to calculate the (anisotropic) inverse Compton scattering in a 3D Monte Carlo simulation, also in presence of magnetic and electric fields. After a short review of the scattering process of photons with electrons, examples of inverse Compton emission are presented, including the predictions for the Sun.

Figures

Figures reproduced from arXiv: 2505.18275 by the authors.

Figure 1
Figure 1. FIG. 1 [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2. Energy (top panel) and angular (bottom panel) [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3. Top left panel: Effective cross section with isotropic photons at 2.7 K as a function of the electron energy and [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figures from the paper (6 more)
Figure 4
Figure 4. Figure 4: FIG. 4. Definition of variables describing the geometry of [PITH_FULL_IMAGE:figures/full_fig_p006_4.png]
Figure 5
Figure 5. Figure 5: FIG. 5. Effective cross section calculated for the solar [PITH_FULL_IMAGE:figures/full_fig_p008_5.png]
Figure 7
Figure 7. Figure 7: FIG. 7. Monte Carlo geometry used in the simulation to [PITH_FULL_IMAGE:figures/full_fig_p009_7.png]
Figure 8
Figure 8. Figure 8: FIG. 8. Distribution of the IC photon energies as a function [PITH_FULL_IMAGE:figures/full_fig_p011_8.png]
Figure 9
Figure 9. Figure 9: FIG. 9. Calculated IC spectral intensity for various angular [PITH_FULL_IMAGE:figures/full_fig_p011_9.png]
Figure 10
Figure 10. Figure 10: FIG. 10. Comparison without (left column) and with magnetic and electric field configurations (right column). Top row: Energy [PITH_FULL_IMAGE:figures/full_fig_p012_10.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

40 extracted references · 28 canonical work pages

  1. [1]

    A. A. Abdo et al. (Fermi-LAT), Astrophys. J. 734, 116 (2011), arXiv:1104.2093 [astro-ph.HE]

    arXiv 2011

  2. [2]

    Seckel, T

    D. Seckel, T. Stanev, and T. K. Gaisser, Astrophys. J. 382, 652 (1991)

    1991

  3. [3]

    M. N. Mazziotta, P. De La Torre Luque, L. Di Venere, A. Fass` o, A. Ferrari, F. Loparco, P. R. Sala, and D. Serini, Phys. Rev. D 101, 083011 (2020), arXiv:2001.09933 [astro-ph.HE]

    work page Pith review arXiv 2020

  4. [4]

    Gamma-ray emission from the solar halo and disk: a study with EGRET data

    E. Orlando and A. W. Strong, Astron. Astrophys. 480, 847 (2008), arXiv:0801.2178 [astro-ph]

    work page Pith review arXiv 2008

  5. [5]

    Gamma-rays from halos around stars and the Sun

    E. Orlando and A. Strong, The Multi-Messenger Approach to High-Energy Gamma-Ray Sources: 3rd Workshop on the Nature of Unidentified High-Energy Sources, Barcelona, Spain, 4-7 Jul, 2006 , Astrophys. Space Sci. 309, 359 (2007), arXiv:astro-ph/0607563 [astro-ph]

    work page Pith review arXiv 2007

  6. [6]

    I. V. Moskalenko, T. A. Porter, and S. W. Digel, Astrophys. J. Lett. 652, L65 (2006), [Erratum: Astrophys.J.Lett. 664, L143 (2007), Erratum: Astrophys.J. 664, L143 (2007)], arXiv:astro-ph/0607521

    work page Pith review arXiv 2006

  7. [7]

    StellarICS: Inverse Compton Emission from the Quiet Sun and Stars from keV to TeV

    E. Orlando and A. Strong, JCAP 04, 004, arXiv:2012.13126 [astro-ph.HE]

    work page Pith review arXiv 2012

  8. [8]

    A. C. M. Lai and K. C. Y. Ng, Phys. Rev. D 107, 063026 (2023), arXiv:2211.15691 [astro-ph.HE]

    work page Pith review arXiv 2023

Show all 40 references
  1. [9]

    Ginzburg and S

    V. Ginzburg and S. Syrovat-Ski ˘ ı, Soviet Physics JETP 19 (1964)

    1964

  2. [10]

    A. A. Abdo et al. (Fermi-LAT), Phys. Rev. Lett. 102, 181101 (2009), arXiv:0905.0025 [astro-ph.HE]

    2009 arXiv

  3. [11]

    Grasso et al

    D. Grasso et al. (Fermi-LAT), Astropart. Phys. 32, 140 (2009), arXiv:0905.0636 [astro-ph.HE]

    2009 arXiv

  4. [12]

    Aguilar et al

    M. Aguilar et al. (AMS), Phys. Rev. Lett. 113, 221102 (2014)

    2014

  5. [13]

    Abdollahi et al

    S. Abdollahi et al. (Fermi-LAT), Phys. Rev. D95, 082007 (2017), arXiv:1704.07195 [astro-ph.HE]

    2017 arXiv

  6. [14]

    Ambrosi et al

    G. Ambrosi et al. (DAMPE), Nature 552, 63 (2017), arXiv:1711.10981 [astro-ph.HE]

    2017 arXiv

  7. [15]

    L. A. Pozdnyakov, I. M. Sobol, and R. A. Syunyaev, Soviet Astronomy Letters 2, 55 (1976)

    1976

  8. [16]

    L. A. Pozdnyakov, I. M. Sobol, and R. A. Syunyaev, Soviet Astronomy Letters 21, 708 (1977)

    1977

  9. [17]

    L. A. Pozdniakov, I. M. Sobol, and R. A. Sunyaev, Astronomy and Astrophysics 75, 214 (1979)

    1979

  10. [18]

    V. B. Berestetskii, E. M. Lifshitz, and L. P. Pitaevskii, Quantum Electrodynamics: Volume 4 , Vol. 4 (Butterworth-Heinemann, 1982)

    1982

  11. [19]

    Akhiezer and V

    A. Akhiezer and V. Berestetski ˘ ı, Quantum Electrodynamics, Interscience monographs and texts in physics and astronomy, v. 11 (Interscience Publishers, 1965)

    1965

  12. [20]

    Klein and T

    O. Klein and T. Nishina, Zeitschrift fur Physik 52, 853 (1929)

    1929

  13. [21]

    J. C. Butcher and H. Messel, Nuclear Physics 20, 15 (1960)

    1960

  14. [22]

    Curatolo, I

    C. Curatolo, I. Drebot, V. Petrillo, and L. Serafini, Phys. Rev. Accel. Beams 20, 080701 (2017), [Addendum: Phys. Rev. Accel. Beams 20, 109901 (2017)], arXiv:1705.07740 [physics.acc-ph]

    2017 arXiv

  15. [23]

    F. C. Jones, Phys. Rev. 167, 1159 (1968)

    1968

  16. [24]

    N. E. Agency, PENELOPE 2018: A code system for Monte Carlo simulation of electron and photon transport (2019) p. 420

    2019

  17. [25]

    GeoGebra for Teaching and Learning Math, https:// www.geogebra.org/

  18. [26]

    Genz and A

    A. Genz and A. Malik, Journal of Computational and Applied Mathematics 6, 295 (1980)

    1980

  19. [27]

    Brun and F

    R. Brun and F. Rademakers, Nucl. Instrum. Meth. A 389, 81 (1997)

    1997

  20. [28]

    E. N. Parker, Astrophys. J. 128, 664 (1958)

    1958

  21. [29]

    Lipari, Solar modulations by the regular heliospheric electromagnetic field (2014), arXiv:1408.0431 [astro- ph.HE]

    P. Lipari, Solar modulations by the regular heliospheric electromagnetic field (2014), arXiv:1408.0431 [astro- ph.HE]

    2014 arXiv

  22. [30]

    R. Brun, F. Bruyant, M. Maire, A. C. McPherson, and P. Zanarini, (1987)

    1987

  23. [31]

    K. M. Gorski, E. Hivon, A. J. Banday, B. D. Wandelt, F. K. Hansen, M. Reinecke, and M. Bartelmann, The Astrophysical Journal 622, 759 (2005)

    2005

  24. [32]

    Ackermann et al

    M. Ackermann et al. (Fermi-LAT), Astrophys. J.799, 86 (2015), arXiv:1410.3696 [astro-ph.HE]

    2015 arXiv

  25. [33]

    Acero et al

    F. Acero et al. (Fermi-LAT), Astrophys. J. Suppl. 223, 26 (2016), arXiv:1602.07246 [astro-ph.HE]

    2016 arXiv

  26. [34]

    Vernetto and P

    S. Vernetto and P. Lipari, Phys. Rev. D 94, 063009 (2016), arXiv:1608.01587 [astro-ph.HE]

    2016 arXiv

  27. [35]

    Ferrari, P

    A. Ferrari, P. R. Sala, A. Fass` o, and J. Ranft, FLUKA: A multi-particle transport code (Program version 2005) (2005)

    2005

  28. [36]

    Fass` o, A

    A. Fass` o, A. Ferrari, J. Ranft, and P. R. Sala, Conf. Proc. C 9309194, 493 (1993)

    1993

  29. [37]

    Aric´ o, A

    G. Aric´ o, A. Ferrari, F. Horst, A. Mairani, C. Reidel, C. Schuy, and U. Weber, CERN Proc. 1, 321 (2019)

    2019

  30. [38]

    Ballarini et al

    F. Ballarini et al. , EPJ Nuclear Sci. Technol. 10, 10.1051/epjn/2024015 (2024)

    2024

  31. [39]

    Van Rossum and F

    G. Van Rossum and F. L. Drake Jr, Python reference manual (Centrum voor Wiskunde en Informatica Amsterdam, 1995)

    1995

  32. [40]

    J. D. Hunter, Computing in Science & Engineering 9, 90 (2007)

    2007

Pith tools

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