Pith. sign in

REVIEW 2 major objections 2 minor 121 references

The Parker transport equation overestimates low-energy galactic cosmic ray intensities by about 30 percent at Earth compared to the focused transport equation.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

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

T0 review · grok-4.3

2026-06-27 14:09 UTC pith:BUNR6JSU

load-bearing objection This first numerical comparison finds the Parker TPE overestimates GCR intensities by 30-40% versus the focused version under normalized diffusion, but the result rests on that normalization choice. the 2 major comments →

arxiv 2606.09298 v1 pith:BUNR6JSU submitted 2026-06-08 physics.space-ph

Are the Parker and Focused Transport Equations Equivalent for Galactic Cosmic Ray Modulation?

J. P. van den Berg , N. E. Engelbrecht , R. D. Strauss This is my paper
classification physics.space-ph
keywords galactic cosmic raysParker transport equationfocused transport equationsolar modulationdiffusion coefficientsheliospherecosmic ray anisotropystochastic differential equations
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

The paper compares the Parker transport equation, used for decades to model galactic cosmic ray modulation, against the focused transport equation that tracks particle anisotropy instead of assuming isotropy. Under matched diffusion conditions obtained by scaling the pitch-angle dependent coefficients to the isotropic ones, the Parker equation produces intensities roughly 30 percent higher at Earth's orbit for low energies and 40 percent higher over the poles. The difference arises because the Parker equation averages particle streaming over the poles, where scattering is weaker and the field less wound, into an effective parallel diffusion term that proves too strong. The focused equation shows that this streaming allows particles easier inner-heliosphere access, and that results stay nearly the same across different pitch-angle forms of the diffusion coefficients.

Core claim

The Parker transport equation overestimates the galactic cosmic ray intensity at Earth's orbit for low energies by ~30%, and by ~40% over the poles. This stems from a small first-order anisotropy caused by particle fluxes over the poles. Particles gain easier access to the inner heliosphere by streaming in over the poles, where pitch-angle scattering is generally weaker, and the magnetic field is typically less wound. The focused transport equation also yields nearly identical results for different pitch-angle dependencies of the diffusion coefficients. The description of particle streaming and weak pitch-angle scattering as effective parallel diffusion in the Parker transport equation makes

What carries the argument

Normalization of the pitch-angle-dependent diffusion coefficients in the focused transport equation to the isotropic diffusion coefficients in the Parker transport equation, so that both solve the same modulation problem under identical scattering strength.

Load-bearing premise

The pitch-angle-dependent diffusion coefficients used in the focused transport equation can be normalized to the isotropic diffusion coefficients used in the Parker transport equation to produce a fair comparison under identical diffusion conditions.

What would settle it

Measurements of galactic cosmic ray intensity spectra at low energies both at Earth and at high heliographic latitudes during solar minimum that match one equation's prediction but deviate systematically from the other.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • The Parker transport equation treats particle streaming and weak pitch-angle scattering as effective parallel diffusion and is therefore overly diffusive.
  • Diffusion coefficients obtained by fitting the Parker transport equation to observations are likely underestimated.
  • Galactic cosmic ray spectral shape and anisotropy data alone cannot distinguish scattering theories that share similar mean free paths but differ in pitch-angle dependence.
  • The focused transport equation produces nearly identical intensities for different assumed pitch-angle forms of the diffusion coefficients.

Where Pith is reading between the lines

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

  • Modulation models that rely on the Parker equation may require lower diffusion coefficients or added anisotropy terms to match low-energy observations accurately.
  • Extending the comparison to include particle drifts could test whether the intensity difference grows or shrinks when drift effects are restored.
  • The result suggests that polar regions act as a preferred entry channel whose importance scales with the weakness of scattering there.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

2 major / 2 minor

Summary. The paper claims that the Parker and focused transport equations (TPEs) are not equivalent for galactic cosmic ray (GCR) modulation without drifts. Using a new proton model for solar minimum conditions solved via stochastic differential equations, with pitch-angle-dependent diffusion coefficients (DCs) in the focused TPE normalized to the isotropic DCs of the Parker TPE, the Parker TPE is found to overestimate GCR intensity at 1 AU by ~30% at low energies and by ~40% over the poles. The discrepancy arises from a small first-order anisotropy due to polar particle streaming (where scattering is weaker and the field less wound). The focused TPE yields nearly identical results across different D(μ) forms, implying that the Parker TPE's effective treatment of streaming as parallel diffusion is overly diffusive and that fitted DCs from Parker-based studies are likely underestimated. Spectral and anisotropy data alone cannot distinguish scattering theories with similar mean free paths but differing μ-dependencies.

Significance. If the result holds after addressing the normalization validation, the work would challenge the 60-year dominance of the Parker TPE as the standard for GCR modulation modeling by demonstrating that it misses first-order anisotropy effects from polar access. This has implications for interpreting fitted diffusion coefficients and for the information content of GCR observations. Strengths include the first direct numerical comparison under controlled (no-drift) conditions, the use of SDEs for both equations, and the explicit test of robustness to the functional form of D(μ).

major comments (2)
  1. [§3] §3 (normalization of DCs): The central claim of inequivalence under 'identical diffusion conditions' rests on normalizing the pitch-angle-dependent DCs of the focused TPE to the isotropic DCs of the Parker TPE. Because the focused TPE retains explicit μ-dependence, focusing, and first-order anisotropy (especially weaker polar scattering), an integral or mean-free-path matching may not produce equivalent effective transport; the reported 30-40% intensity differences could be an artifact of this mapping rather than a fundamental inequivalence. The intra-focused insensitivity test does not validate the cross-TPE normalization.
  2. [Results] Results (intensity profiles at 1 AU and poles): The quantitative overestimation figures (~30% at low energies at Earth's orbit, ~40% over poles) are presented without reported sensitivity tests to the normalization factor (a free parameter) or direct comparison of effective κ_|| or mean free paths between the two TPEs. This is load-bearing for the inequivalence conclusion.
minor comments (2)
  1. [Abstract] The abstract and introduction would benefit from a brief explicit statement of the normalization procedure and its limitations to aid readers in interpreting the equivalence test.
  2. [Figures] Figure captions and axis labels should clarify whether intensities are normalized or absolute, and whether the polar overestimation refers to a specific latitude or integrated over the polar cap.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful and constructive review. The comments raise valid points about the normalization procedure and the need for additional validation of the quantitative results. We address each major comment below and will incorporate revisions to strengthen the manuscript.

read point-by-point responses

  1. Referee: [§3] §3 (normalization of DCs): The central claim of inequivalence under 'identical diffusion conditions' rests on normalizing the pitch-angle-dependent DCs of the focused TPE to the isotropic DCs of the Parker TPE. Because the focused TPE retains explicit μ-dependence, focusing, and first-order anisotropy (especially weaker polar scattering), an integral or mean-free-path matching may not produce equivalent effective transport; the reported 30-40% intensity differences could be an artifact of this mapping rather than a fundamental inequivalence. The intra-focused insensitivity test does not validate the cross-TPE normalization.

    Authors: The normalization in §3 is performed by matching the pitch-angle-averaged parallel diffusion coefficient (equivalently the mean free path λ_||) between the focused TPE and the isotropic κ_|| of the Parker TPE, which is the standard approach for comparing the equations under comparable diffusion conditions. The insensitivity of focused-TPE results to the functional form of D(μ) indicates that the intensity differences are driven by the explicit treatment of first-order anisotropy and polar streaming rather than details of the μ-dependence. We nevertheless agree that cross-equation validation of effective transport coefficients is needed and will add this in revision. revision: yes

  2. Referee: [Results] Results (intensity profiles at 1 AU and poles): The quantitative overestimation figures (~30% at low energies at Earth's orbit, ~40% over poles) are presented without reported sensitivity tests to the normalization factor (a free parameter) or direct comparison of effective κ_|| or mean free paths between the two TPEs. This is load-bearing for the inequivalence conclusion.

    Authors: We agree that sensitivity tests to the normalization factor and explicit comparison of effective κ_|| (or λ_||) are required to confirm the robustness of the 30–40% differences. In the revised manuscript we will include: (i) a direct side-by-side comparison of the effective parallel diffusion coefficients realized in both models, and (ii) sensitivity plots showing how the reported intensity overestimations change when the normalization factor is varied by ±20%. revision: yes

Circularity Check

0 steps flagged

No significant circularity in the TPE equivalence test

full rationale

The paper performs a direct numerical comparison of the Parker and focused transport equations via stochastic differential equations under explicitly normalized diffusion conditions. The normalization of pitch-angle-dependent DCs to isotropic DCs is stated as a deliberate design choice to enable the test of equivalence, and the reported intensity differences are presented as numerical outcomes arising from the focused TPE's retention of first-order anisotropy and polar streaming. No derivation reduces a claimed result to its inputs by construction, no parameters are fitted and then relabeled as predictions, and no self-citation chain or uniqueness theorem is invoked as load-bearing. The finding that results are insensitive to the form of D(μ) further indicates an independent numerical exploration rather than a tautological mapping.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The central claim rests on the normalization procedure for diffusion coefficients and the domain assumptions of no drifts and solar minimum conditions. Full details on specific parameter values or implementation are unavailable from the abstract alone.

free parameters (1)
  • Normalization factor for pitch-angle-dependent DCs
    Applied to align the focused TPE diffusion with the isotropic DCs of the Parker TPE for consistent comparison.
axioms (2)
  • domain assumption Particle drifts are neglected in both models
    Stated to isolate diffusion and anisotropy effects in the comparison.
  • domain assumption The proton model represents solar minimum conditions
    Basis for the new numerical model developed to solve the TPEs.

pith-pipeline@v0.9.1-grok · 5867 in / 1561 out tokens · 34372 ms · 2026-06-27T14:09:19.686856+00:00 · methodology

0 comments
read the original abstract

The Parker transport equation (TPE) has been the equation of choice for the past 60 years in studies of galactic cosmic ray (GCR) modulation. Conversely, the focused TPE describes the same processes on a more fundamental level than the Parker TPE by modelling an anisotropic distribution rather than an isotropic one. It is usually assumed that the Parker TPE is valid for modelling GCRs, but the two TPEs have not been tested against each other in this context. We conduct a first-of-its-kind comparison of these TPEs without particle drifts to test whether they produce the same results under identical diffusion conditions. A new model for protons during solar minimum conditions is developed to numerically solve the TPEs using stochastic differential equations. The TPEs are designed to be as consistent as possible for diffusion by normalising the pitch-angle-dependent diffusion coefficients (DCs) used in the focused TPE to the isotropic DCs used in the Parker TPE. The Parker TPE overestimates the GCR intensity at Earth's orbit for low energies by ~30%, and by ~40% over the poles. This stems from a small first-order anisotropy caused by particle fluxes over the poles. Particles gain easier access to the inner heliosphere by streaming in over the poles, where pitch-angle scattering is generally weaker, and the magnetic field is typically less wound. The focused TPE also yields nearly identical results for different pitch-angle dependencies of the DCs. The description of particle streaming and weak pitch-angle scattering as effective parallel diffusion in the Parker TPE makes it overly diffusive. This suggests that DCs derived from fitting the Parker TPE to observations are likely underestimated. Furthermore, GCR spectral and anisotropy data alone cannot distinguish between scattering theories with similar mean free paths but different pitch-angle dependencies.

Figures

Figures reproduced from arXiv: 2606.09298 by J. P. van den Berg, N. E. Engelbrecht, R. D. Strauss.

Figure 1
Figure 1. Figure 1: Modified focusing length (Eq. 11; grey lines) and various MFPs for 1 GeV protons as a function of radius in the equatorial [PITH_FULL_IMAGE:figures/full_fig_p006_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Top: Proton spectrum at a radial distance of 1 au in the [PITH_FULL_IMAGE:figures/full_fig_p011_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Top: Percentage deviation between the results of the [PITH_FULL_IMAGE:figures/full_fig_p011_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Normalised PADs averaged over all latitudes, ra [PITH_FULL_IMAGE:figures/full_fig_p012_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: First- (top) and second-order (bottom) anisotropies calculated from fitting Eqs 47 (dashed orange) or C.6 (solid blue) to the [PITH_FULL_IMAGE:figures/full_fig_p013_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Radial advective (top left), radial diffusive (top right), total radial (bottom left), and polar diffusive (bottom right) fluxes throughout the heliosphere resulting from the Parker TPE for four different energies (indicated in the quadrants). Hatches in the right panels indicate flux directed in the opposite direction than expected, and dots in the bottom left panel indicate inward flux. Still, the anisot… view at source ↗

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

121 extracted references

  1. [1]

    U., Alfaro, R., Alvarez, C., et al

    Abeysekara, A. U., Alfaro, R., Alvarez, C., et al. 2019, , 871, 96

    2019

  2. [2]

    C., Bazilevskaya, G

    Adriani, O., Barbarino, G. C., Bazilevskaya, G. A., et al. 2013, , 765, 91

    2013

  3. [3]

    C., Bazilevskaya, G

    Adriani, O., Barbarino, G. C., Bazilevskaya, G. A., et al. 2015, , 811, 21

    2015

  4. [4]

    & Vainio, R

    Agueda, N. & Vainio, R. 2013, J. Space Weather Space Clim., 3, A10

    2013

  5. [5]

    2008, , 675, 1601

    Agueda, N., Vainio, R., Lario, D., & Sanahuja, B. 2008, , 675, 1601

    2008

  6. [6]

    2019, , 883, 33

    Ajello, M., Baldini, L., Barbiellini, G., et al. 2019, , 883, 33

    2019

  7. [7]

    J., et al

    Amenomori, M., Ayabe, S., Bi, X. J., et al. 2006, Sci., 314, 439

    2006

  8. [8]

    J., et al

    Bartoli, B., Bernardini, P., Bi, X. J., et al. 2018, , 861, 93

    2018

  9. [9]

    & Wibberenz, G

    Beeck, J. & Wibberenz, G. 1986, , 311, 437

    1986

  10. [10]

    H., Smith, C

    Bieber, J., Matthaeus, W. H., Smith, C. W., et al. 1994, , 420, 294

    1994

  11. [11]

    H., Minnie, J., et al

    Breech, B., Matthaeus, W. H., Minnie, J., et al. 2005, , 32, L06103

    2005

  12. [12]

    A., Kr \"u ger, T

    Burger, R. A., Kr \"u ger, T. P. J., Hitge, M., & Engelbrecht, N. E. 2008, , 674, 511

    2008

  13. [13]

    A., Moraal, H., & Webb, G

    Burger, R. A., Moraal, H., & Webb, G. M. 1985, , 116, 107

    1985

  14. [14]

    A., Potgieter, M

    Burger, R. A., Potgieter, M. S., & Heber, B. 2000, , 105, 27447

    2000

  15. [15]

    Caballero-Lopez , R. A. & Moraal, H. 2004, (Space Phys.), 109, A01101

    2004

  16. [16]

    2024, , 961, 87

    Chakraborty, M., Ahmad, S., Chandra, A., et al. 2024, , 961, 87

    2024

  17. [17]

    L., Conlon, T

    Chenette, D. L., Conlon, T. F., Pyle, K. R., & Simpson, J. A. 1977, , 215, L95

    1977

  18. [18]

    2010, (Space Phys.), 115, A10106

    Dunzlaff, P., Kopp, A., & Heber, B. 2010, (Space Phys.), 115, A10106

    2010

  19. [19]

    & Litvinenko, Y

    Effenberger, F. & Litvinenko, Y. E. 2014, , 783, 15

    2014

  20. [20]

    2025, , 221, 75

    Effenberger, F., Walter, D., Fichtner, H., et al. 2025, , 221, 75

    2025

  21. [21]

    L., Engelbrecht, N

    Els, P. L., Engelbrecht, N. E., Lang, J. T., & Strauss, R. D. 2024, , 975, 134

    2024

  22. [22]

    Engelbrecht, N. E. 2017, , 849, L15

    2017

  23. [23]

    Engelbrecht, N. E. 2019, , 880, 60

    2019

  24. [24]

    Engelbrecht, N. E. 2024, , 975, 227

    2024

  25. [25]

    Engelbrecht, N. E. & Burger, R. A. 2013 a , , 772, 46

    2013

  26. [26]

    Engelbrecht, N. E. & Burger, R. A. 2013 b , , 779, 158

    2013

  27. [27]

    Engelbrecht, N. E. & Di Felice , V. 2020, , 102, 103007

    2020

  28. [28]

    E., Effenberger, F., Florinski, V., et al

    Engelbrecht, N. E., Effenberger, F., Florinski, V., et al. 2022 a , , 218, 33

    2022

  29. [29]

    E., Herbst, K., Strauss, R

    Engelbrecht, N. E., Herbst, K., Strauss, R. D. T., et al. 2024, , 964, 89

    2024

  30. [30]

    E., Vogt, A., Herbst, K., Strauss, R

    Engelbrecht, N. E., Vogt, A., Herbst, K., Strauss, R. D. T., & Burger, R. A. 2022 b , , 929, 8

    2022

  31. [31]

    Ferreira, S. E. S., Potgieter, M. S., Burger, R. A., Heber, B., & Fichtner, H. 2001, , 106, 24979

    2001

  32. [32]

    E., Bosilca, G., et al

    Gabriel, E., Fagg, G. E., Bosilca, G., et al. 2004, in Proc. 11th European PVM/MPI Users' Group Meeting, Budapest, Hungary, 97

    2004

  33. [33]

    Gardiner, C. W. 1994, Handbook of stochastic methods for physics, chemistry and the natural sciences (Springer-Verlag)

    1994

  34. [34]

    & Jokipii, J

    Giacalone, J. & Jokipii, J. R. 1999, , 520, 204

    1999

  35. [35]

    Gleeson, L. J. & Axford, W. I. 1968, , 154, 1011

    1968

  36. [36]

    Gleeson, L. J. & Urch, I. H. 1971, , 11, 288

    1971

  37. [37]

    Gradshteyn, I. S. & Ryzhik, I. M. 2007, Table of integrals, series, and products , 7th edn., ed. A. Jeffrey & D. Zwillinger (Academic Press)

    2007

  38. [38]

    2007, , 75, 062003

    Guillian, G., Hosaka, J., Ishihara, K., et al. 2007, , 75, 062003

    2007

  39. [39]

    R., Millman, K

    Harris, C. R., Millman, K. J., van der Walt , S. J., et al. 2020, , 585, 357

    2020

  40. [40]

    & Wibberenz, G

    Hasselmann, K. & Wibberenz, G. 1970, , 162, 1049

    1970

  41. [41]

    & Schlickeiser, R

    He, H.-Q. & Schlickeiser, R. 2014, , 792, 85

    2014

  42. [42]

    R., Bykov, A., et al

    Herbst, K., Baalmann, L. R., Bykov, A., et al. 2022, , 218, 29

    2022

  43. [43]

    Hunter, J. D. 2007, Computing Sci. Eng., 9, 90

    2007

  44. [44]

    Isenberg, P. A. 1997, , 102, 4719

    1997

  45. [45]

    Isenberg, P. A. & Jokipii, J. R. 1979, , 234, 746

    1979

  46. [46]

    Jokipii, J. R. & Parker, E. N. 1970, , 160, 735

    1970

  47. [47]

    Kloeden, P. E. & Platen, E. 1995, Numerical solution of stochastic differential equations , 2nd edn. (Springer-Verlag)

    1995

  48. [48]

    D., & Potgieter, M

    Kopp, A., B \"u sching, I., Strauss, R. D., & Potgieter, M. S. 2012, Comp. Phys. Comm., 183, 530

    2012

  49. [49]

    2013, , 176, 391

    K \'o ta, J. 2013, , 176, 391

    2013

  50. [50]

    2011, PhD thesis, University of Osnabrück, Germany

    Lampa, F. 2011, PhD thesis, University of Osnabrück, Germany

    2011

  51. [51]

    T., Strauss, R

    Lang, J. T., Strauss, R. D., Engelbrecht, N. E., et al. 2024, , 971, 105

    2024

  52. [52]

    le Roux , J. A. & Webb, G. M. 2012, , 746, 104

    2012

  53. [53]

    A., Webb, G

    le Roux , J. A., Webb, G. M., Florinski, V., & Zank, G. P. 2007, , 662, 350

    2007

  54. [54]

    A., Webb, G

    le Roux , J. A., Webb, G. M., & Ye, J. 2014, in Astro. Society Pacific Conf. Series, Vol. 484, Outstanding problems in heliophysics: From coronal heating to the edge of the heliosphere, ed. Q. Hu & G. P. Zank, 110

    2014

  55. [55]

    Litvinenko, Y. E. & Noble, P. L. 2013, , 765, 31

    2013

  56. [56]

    Litvinenko, Y. E. & Schlickeiser, R. 2013, , 554, A59

    2013

  57. [57]

    Maalal, N. D. & Zhang, M. 2025, , 992, 46

    2025

  58. [58]

    Malkov, M. A. 2017, , 95, 023007

    2017

  59. [59]

    Malkov, M. A. 2018, Nuc. Part. Phys. Proc., 297-299, 152

    2018

  60. [60]

    2018, , 854, L2

    Martucci, M., Munini, R., Boezio, M., et al. 2018, , 854, L2

    2018

  61. [61]

    H., Qin, G., Bieber, J

    Matthaeus, W. H., Qin, G., Bieber, J. W., & Zank, G. P. 2003, , 590, L53

    2003

  62. [62]

    J., Barraclough, B

    McComas , D. J., Barraclough, B. L., Funsten, H. O., et al. 2000, , 105, 10419

    2000

  63. [63]

    B., Zhang, M., Heber, B., Kunow, H., & Sanderson, T

    McKibben , R. B., Zhang, M., Heber, B., Kunow, H., & Sanderson, T. R. 2007, , 55, 21

    2007

  64. [64]

    L., Rodgers-Lee , D., & Vidotto, A

    Mesquita, A. L., Rodgers-Lee , D., & Vidotto, A. A. 2021, , 505, 1817

    2021

  65. [65]

    W., Matthaeus, W

    Minnie, J., Bieber, J. W., Matthaeus, W. H., & Burger, R. A. 2007, , 670, 1149

    2007

  66. [66]

    D., Engelbrecht, N

    Moloto, K. D., Engelbrecht, N. E., & Burger, R. A. 2018, , 859, 107

    2018

  67. [67]

    2013, , 176, 299

    Moraal, H. 2013, , 176, 299

    2013

  68. [68]

    1987, , 313, 471

    Moses, D. 1987, , 313, 471

    1987

  69. [69]

    O'Neill, M. E. 2014, PCG: A family of simple fast space-efficient statistically good algorithms for random number generation , Tech. Rep. HMC-CS-2014-0905, Harvey Mudd College, Claremont, CA

    2014

  70. [70]

    Owens, M. J. & Forsyth, R. J. 2013, Living Rev. Solar Phys., 10, 5

    2013

  71. [71]

    Palmer, I. D. 1982, Rev. Geophys. Space Phys., 20, 335

    1982

  72. [72]

    Parker, E. N. 1958, , 128, 664

    1958

  73. [73]

    Parker, E. N. 1965, , 13, 9

    1965

  74. [74]

    Parker, E. N. 1967, , 15, 1723

    1967

  75. [75]

    W., Burger, R

    Pei, C., Bieber, J. W., Burger, R. A., & Clem, J. 2010, (Space Phys.), 115, A12107

    2010

  76. [76]

    Potgieter, M. S. 2013, Living Rev. Solar Phys., 10, 3

    2013

  77. [77]

    S., Vos, E

    Potgieter, M. S., Vos, E. E., Boezio, M., et al. 2014, , 289, 391

    2014

  78. [78]

    H., Teukolsky, S

    Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, B. P. 1992, Numerical recipes in C: The art of scientific computing (Cambridge University Press)

    1992

  79. [79]

    Pyle, K. R. & Simpson, J. A. 1977, , 215, L89

    1977

  80. [80]

    2007, , 656, 217

    Qin, G. 2007, , 656, 217

    2007

Showing first 80 references.