Pith. sign in

REVIEW 4 major objections 5 minor 34 references

BORAY-3D: A ray tracing code for three-dimensional magnetized plasma configurations

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

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

Pith's one-line read BORAY-3D is a ray tracing code that aims to model radio-frequency wave propagation and absorption in any magnetized plasma configuration—axisymmetric or fully three-dimensional, closed or open field lines—across the frequency range from…

desk verdict A useful code paper that extends BORAY to 3D, but the 'fully relativistic' absorption claim is overstated because the polarization stays weakly relativistic. read the letter →

arxiv 2608.05667 v1 pith:XOPDXLFM submitted 2026-08-06 physics.plasm-ph

Yuxuan Wang , Huasheng Xie This is my paper · ORCID
classification physics.plasm-ph PACS 52.65.-y52.50.Qt52.25.Os
keywords raytracingradiofrequencywavesthree-dimensionalequilibriumbroad-frequency-rangemodelingrelativisticabsorptionelectroncyclotronemissionstellaratortokamak
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

BORAY-3D is a ray tracing code that aims to model radio-frequency wave propagation and absorption in any magnetized plasma configuration—axisymmetric or fully three-dimensional, closed or open field lines—across the frequency range from ion cyclotron and helicon waves through lower hybrid to electron cyclotron waves and emission. The paper argues that combining these three capabilities in one code is feasible and useful, and validates the code against GENRAY, Raytrax, TRAVIS, and published HSX and W7-X results. If correct, it gives researchers a single tool for RF heating, current drive, and ECE diagnostics in devices like tokamaks and stellarators.

What carries the argument

The central object is the geometric-optics dispersion function $D(\omega, k_\parallel^2, k_\perp^2, r, \phi, z)$ that defines the ray Hamiltonian. BORAY-3D solves the Hamilton ray equations in cylindrical coordinates $(r, \phi, z)$ with wave-vector components $(k_r, n_\phi/r, k_z)$, treating $n_\phi$ as a dynamic variable because the dispersion relation depends explicitly on $\phi$. Propagation uses the cold-plasma dispersion relation; absorption is then computed separately along the ray, with the non-relativistic Maxwellian hot-plasma tensor for ion-cyclotron, helicon, and lower-hybrid waves and a fully relativistic Maxwellian absorption coefficient for electron-cyclotron waves. The same absorption coefficient feeds a reciprocal radiative-transfer equation, yielding the ECE radiative temperature spectrum via Kirchhoff's law.

What would settle it

Run a low-frequency helicon or ion-cyclotron case in a compact device where the wavelength is comparable to the machine size, and compare the predicted deposition profile against a full-wave simulation: if the cold-ray trajectory misses the resonance layer that the full-wave solution finds, the geometric-optics-plus-cold-ray split would be shown to fail in that regime.

Watch

Extended reading notes

Core claim

The paper's central claim is that three capabilities—broad frequency coverage, arbitrary three-dimensional geometry with closed and open field lines, and fully relativistic electron-cyclotron absorption—can be integrated into a single ray tracing code without sacrificing accuracy. It establishes this by deriving the three-dimensional ray equations in cylindrical coordinates with a non-conserved toroidal mode number, coupling cold-plasma ray propagation to species-resolved hot-plasma absorption, and benchmarking against GENRAY, Raytrax, TRAVIS, and public HSX and W7-X results. Trajectory agreement is at the level of millimeters for W7-X electron-cyclotron cases, and absorbed fractions are consistent with the spread among comparison codes, while the relativistic EC model is shown to be necessary for correct higher-harmonic optical depths and ECE spectra.

Load-bearing premise

The ray path is computed from the cold-plasma dispersion relation while absorption is evaluated afterward along that same path, so the argument holds only when thermal effects do not significantly bend the ray—valid when the wavelength is much shorter than the equilibrium scale length, and unreliable near cutoffs, resonances, or for low-frequency waves with long wavelengths.

Editorial extensions

If this is right

  • The same code can model heating, current drive, and ECE diagnostics in tokamaks with toroidal-field ripple, in stellarators like W7-X and HSX, and in permanent-magnet devices like MUSE, without switching codes or coordinate systems.
  • The relativistic EC absorption model correctly captures higher-harmonic X3 optical depths and ECE spectra, which non-relativistic models underestimate by roughly an order of magnitude.
  • Toroidal-field ripple can shift ray trajectories by up to about 0.7 m in lower-hybrid cases, and BORAY-3D reproduces this shift in agreement with GENRAY.
  • For W7-X electron-cyclotron cases, ray trajectories agree within 5 mm with TRAVIS and Raytrax, and differences in final absorbed fractions are attributed mainly to differing kinetic absorption models rather than to geometry handling.
  • Because density and temperature are not required to be flux functions, the code can be applied to open-field-line regions and numerical equilibria as readily as to analytic axisymmetric models.

Reading between the lines

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

  • A natural extension, which the authors list as future work, is to add collisional damping, wall reflection, current drive, and non-Maxwellian distributions; coupling to full-wave solvers could then address diffraction and edge reflection that geometric optics cannot.
  • The cold-ray/hot-absorption split could be tested more rigorously by comparing against full-wave simulations for low-frequency waves in compact devices where the wavelength is not small compared with the equilibrium scale length, a regime the paper itself flags as needing caution.
  • The cylindrical-coordinate, flux-free formulation should make it straightforward to couple BORAY-3D with numerical equilibria from VMEC, Biot–Savart coil models, or experimental reconstructions, consistent with the paper's use of public design equilibria for MUSE and HSX.
  • The demo of MUSE helicon waves is explicitly an application example rather than a benchmark; a quantitative study would require measured profiles and antenna spectra, which the paper notes are not yet published.
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 / 5 minor

Summary. The paper presents BORAY-3D, a geometric-optics ray tracing code for magnetized plasmas in three-dimensional cylindrical coordinates. The code uses the cold-plasma dispersion relation for ray propagation (Sec. 2.1), a non-relativistic hot-plasma absorption model for ion-cyclotron, helicon, and lower-hybrid waves (Sec. 2.2), and a relativistic Maxwellian absorption model for electron-cyclotron waves and ECE (Sec. 2.3), with benchmarks against GENRAY, Raytrax, TRAVIS, and public results from HSX and W7-X, plus application examples for LHD-like ICW and MUSE helicon (Sec. 3). The abstract claims three combined features: broad frequency coverage, unified treatment of 2D and 3D configurations without flux coordinates, and fully relativistic Maxwellian EC absorption.

Significance. If the claims hold, BORAY-3D occupies a useful niche: a single ray tracing framework spanning IC, helicon, LH, EC, and ECE in arbitrary three-dimensional equilibria, with cylindrical-coordinate input not restricted to flux functions. The trajectory benchmarks show millimeter-level agreement with established codes, the relativistic absorption module reproduces expected X2/X3 ECE features, and the paper is generally clear. The code implements standard and documented physics models without parameter tuning. However, the novelty claim of 'fully relativistic' absorption is overstated because the polarization is weakly relativistic, and some benchmark comparisons (HSX, W7-X ECE) are qualitative or not quantitatively compared with the cited public data. These issues need to be addressed before the paper can be recommended for acceptance.

major comments (4)
  1. [Sec. 2.3, abstract] The headline claim of 'fully relativistic Maxwellian EC absorption' is not supported by the implementation: the wave polarization entering the absorption coefficient in Eq. (20) via Eq. (21) is computed from the weakly relativistic Krivenski–Orefice tensor (Sec. 2.3 and Appendix A, Eqs. (27)–(29)). Thus the absorption model is hybrid, combining a fully relativistic resonance condition and Maxwell–Jüttner integral with a weakly relativistic polarization. Because the X3 optical depth is sensitive to the absorption model (Fig. 7), the authors should either implement fully relativistic polarization or explicitly revise the claim and discuss the expected quantitative error for higher harmonics. This is a load-bearing point for the third listed capability.
  2. [Sec. 3.5, Table 1] The W7-X ECE case is listed as a benchmark against Marushchenko et al. [26], but Fig. 7 shows only BORAY-3D spectra with relativistic and non-relativistic models; no published spectrum is overlaid or numerically compared. The reference is used only for receiver geometry and profile fitting. The authors should either add a quantitative comparison with the published data or reclassify this case as an application example, not a benchmark.
  3. [Sec. 3.2, Fig. 4] The HSX comparison reports a first-pass absorbed fraction of 0.29 against the published 0.40 at n_e = 1.5×10^18 m^-3, a 27% relative difference. The paper correctly notes that the machine-readable equilibrium and profiles are unavailable, but as presented the comparison supports only the density dependence and central deposition, not a quantitative benchmark. Please quantify the sensitivity of the absorbed fraction to the assumed equilibrium scaling and profile choices, or explicitly state that the HSX comparison is a qualitative cross-check rather than a quantitative benchmark.
  4. [Sec. 3.2, Fig. 2] In the W7-X O-mode benchmark, the final absorbed fractions are 85.7% (BORAY-3D), 87.1% (TRAVIS), and 73.3% (Raytrax). The text attributes the spread to non-identical kinetic absorption models, but no evidence is given for why Raytrax differs by about 14 percentage points. A comparison of the optical-depth profiles along the ray would help the reader judge whether the difference is understood and whether it affects the code's validation.
minor comments (5)
  1. [Abstract] The phrase 'an unified treatment' should be 'a unified treatment'.
  2. [Eq. (20)] The symbol F is used as both a vector (Eq. (33)) and a scalar in the denominator of Eq. (20); please clarify that |F| is the magnitude of the dimensionless power-flux vector.
  3. [Sec. 3.4, Fig. 6] The MUSE helicon absorption values (0.27% and 0.09%) are for an assumed profile set and may be dominated by geometric-optics limitations at this frequency; the paper's caveat in Sec. 4 should be reiterated in the figure discussion.
  4. [Sec. 3.2, Fig. 2] The label 'BORAY hot' in Fig. 2(c) refers to the non-relativistic model of Sec. 2.2; consider renaming to 'non-relativistic' to avoid confusion with the hot-plasma tensor used in the relativistic module.
  5. [References] Ref. [18] is listed as 'Yu et al. 2026'; please confirm the publication status and provide a DOI if available.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the paper implements explicit dispersion and absorption models and validates them against external codes and published results.

full rationale

The paper does not fit parameters to benchmarks or rename its inputs as predictions. The ray trajectories are defined by the explicit cold-plasma dispersion relation (Eq. 14) together with the Hamiltonian ray equations (Eqs. 1-13), and absorption is computed from independent kinetic models: the non-relativistic hot-plasma tensor of Sec. 2.2 (Eq. 16) and the relativistic Maxwellian integral of Sec. 2.3 (Eq. 20). The 3D capability is a direct extension of BORAY's cylindrical-coordinate equations in which phi-dependence is retained, and the derivation is presented in the paper rather than imported solely by citation. The benchmarks are external: GENRAY for ripple LHW, Raytrax and TRAVIS for W7-X ECW, and public HSX and W7-X data and equilibria. Where the HSX and LHD comparisons use only published trends because original data are unavailable, the paper states this limitation openly and does not tune the code to force agreement. The self-citations to BORAY, PDRF/PDRK, and Yu et al. supply the prior implementation heritage, but the new 3D and relativistic results are checked against independent codes and published values, so the self-citations are not load-bearing in a circular sense. One wording concern is that the 'fully relativistic Maxwellian EC absorption' label is not fully supported by Sec. 2.3, where polarization is taken from the weakly relativistic Krivenski-Orefice tensor; however this is an internal consistency and novelty-claim issue, not a circular reduction of a prediction to its inputs. The cold-ray/hot-absorption split is also an explicitly stated modeling assumption with a documented geometric-optics limitation, not a circular step.

Assumptions & free parameters 0 free parameters · 5 assumptions · 0 invented entities

The paper introduces no new physical entities and fits no free parameters. It relies on established dispersion relations and kinetic absorption models from the prior literature, and the validation is performed against external codes and experiments.

assumptions (5)
  • domain assumption Geometric optics approximation: the wavelength is much shorter than the equilibrium scale length.
    Invoked in Sec. 1 and Sec. 4; the entire ray tracing framework depends on this condition.
  • domain assumption Ray propagation is governed by the cold-plasma dispersion relation; thermal effects are only applied to absorption.
    Sec. 2.1 Eq. (14) and Sec. 2.2; this is the standard ray tracing split used in BORAY and GENRAY.
  • domain assumption The plasma distribution is Maxwellian (or Maxwell-Juttner for the relativistic EC model).
    Sec. 2.3 Eq. (19); the absorption and ECE models assume thermal equilibrium distributions.
  • domain assumption The weakly relativistic Krivenski-Orefice tensor provides sufficient accuracy for the EC wave polarization.
    Sec. 2.3 and Appendix A; the paper does not quantify the error from this approximation.
  • domain assumption The equilibrium data are smooth enough for numerical partial derivatives in the ray equations.
    Sec. 2.1 Eqs. (8)-(13) require derivatives of B and n_s0 with respect to r, phi, and z; the paper does not describe the interpolation or smoothing method.

how reviews work

0 comments
Cite this review

Pith. "Pith review of BORAY-3D: A ray tracing code for three-dimensional magnetized plasma configurations." pith.science (2026). https://pith.science/paper/XOPDXLFM

@misc{pith2026260805667,
  author       = {Pith},
  title        = {Pith review of: BORAY-3D: A ray tracing code for three-dimensional magnetized plasma configurations},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/XOPDXLFM}},
  note         = {Machine review of arXiv:2608.05667}
}
abstract

Ray tracing codes are useful tools for studying electromagnetic wave propagation and absorption using the geometrical-optics approximation. Existing codes commonly provide either broad radio-frequency coverage in axisymmetric equilibria or three-dimensional capability specialized for electron-cyclotron (EC) applications. BORAY-3D integrates three desirable features in a single version. First, it has a broad frequency range of validity regime from ion-cyclotron, helicon and lower-hybrid waves to EC waves and emission. Second, it provides a unified treatment of arbitrary two- and three-dimensional magnetic-plasma configurations, including both closed and open field-line regions. Third, it incorporates fully relativistic Maxwellian EC absorption. The code extends the axisymmetric BORAY formulation by solving the ray equations in cylindrical coordinates $(r,\phi,z)$ while allowing the toroidal mode number $n_\phi$ to vary. Magnetic-field, density and temperature data are directly described in $(r,\phi,z)$ coordinates without the restriction of flux functions, so that numerical equilibria and analytic field models can be handled in the same form. The non-relativistic hot-plasma model inherited from BORAY is used for lower-hybrid, ion-cyclotron and helicon absorption, whereas the relativistic model is coupled to reciprocal radiative transfer for electron cyclotron emission (ECE). Practical applications include 13.56 MHz helicon and 50 MHz fast waves, a 3.7 GHz lower-hybrid wave, and 115--220 GHz EC emission. BORAY-3D has been systematically benchmarked against GENRAY for tokamak toroidal-field ripple, Raytrax and TRAVIS for W7-X, as well as public HSX heating and W7-X ECE results.

Figures

Figures reproduced from arXiv: 2608.05667 by the authors.

Figure 1
Figure 1. Comparison of BORAY-3D and GENRAY for a 3.7 GHz lower-hybrid wave in a tokamak with strong toroidal-field ripple. Solid and dashed curves denote GENRAY and BORAY-3D, respectively. (a) Poloidal ray trajectories, (b) normalized radial position, (c) remaining power and (d) density and electron-temperature profiles. The ray trajectories agree well. The difference in power absorption is due mainly to the different hot-pl… view at source ↗
Figure 2
Figure 2. Comparison of BORAY-3D, Raytrax and TRAVIS for a W7-X 140 GHz O-mode ray. The trajectories agree within 5 mm. The differences in remaining power are mainly due to the different absorption models [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 3
Figure 3. Comparison of BORAY-3D and Raytrax for three W7-X Port-A D1 launches with poloidal aiming angles −20◦, −10◦ and +5◦. The ray trajectories and final absorbed fractions agree well for all three angles [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figures from the paper (4 more)
Figure 4
Figure 4. Figure 4: Comparison of BORAY-3D with the HSX 28 GHz X2 heating results of Ref. [25]. (a) The 18 Gaussian-beam rays at ¯ne = 2 × 1018 m−3 , (b) first-pass absorption versus line-averaged density, (c) averaged power from the plasma edge toward the magnetic axis and (d) density an…
Figure 5
Figure 5. Figure 5: BORAY-3D results for a 50 MHz fast wave in an LHD-like helical field. (a) Top views of the rays for ne0 = 1, 3 and 5 × 1020 m−3 , (b) parallel wave number, (c) first-pass tritium (solid) and electron (dashed) absorption and (d) density and temperature profiles. 7 [PIT…
Figure 6
Figure 6. Figure 6: BORAY-3D results for a 13.56 MHz helicon wave in the MUSE design equilibrium. (a) Top-view ray trajectories, (b) parallel wave number, (c) first-pass collisionless absorption and (d) assumed plasma profiles. The two rays launched at ϕ0 = 0 pass close to the magnetic ax…
Figure 7
Figure 7. Figure 7: W7-X ECE spectra calculated by BORAY-3D. (a) The ϕ = 5.88◦ equilibrium cross-section and 140 GHz central rays, (b) LFS spectrum, (c) HFS spectrum and (d) density and temperature profiles. Blue circles and orange triangles denote the relativistic and non-relativistic re…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

34 extracted references · 34 canonical work pages

  1. [26]

    Marushchenko N B, Erckmann V, Hartfuss H J, Hirsch M, Laqua H P, Maassberg H and Turkin Y 2007Plasma and Fusion Research2S1129

    work page

  2. [1]

    Bornatici M, Cano R, De Barbieri O and Engelmann F 1983Nuclear Fusion231153–1257

    work page

  3. [2]

    Hartfuss H J, Geist T and Hirsch M 1997Plasma Physics and Controlled Fusion391693–1769

    work page

  4. [3]

    Stix T H 1992Waves in Plasmas(New York: American Institute of Physics)

    work page

  5. [4]

    Xie H s, Banerjee D, Bai Y k, Zhao H y and Li J c 2022Computer Physics Communications 276108363

    work page

  6. [5]

    Xie H s 2014Computer Physics Communications185670–675

    work page

  7. [6]

    Xie H and Xiao Y 2016Plasma Science and Technology1897–107

    work page

  8. [7]

    Xie H s 2019Computer Physics Communications244343–371

    work page

Show all 34 references
  1. [8]

    Beidler C D, Smith H M, Alonso Aet al.2021Nature596221–226

  2. [9]

    Peysson Y, Decker J and Morini L 2012Plasma Physics and Controlled Fusion54045003

  3. [10]

    Marushchenko N B, Turkin Y and Maassberg H 2014Computer Physics Communications185 165–176

  4. [11]

    Hirshman S P and Whitson J C 1983Physics of Fluids263553–3568

  5. [12]

    Proxima Fusion Raytrax: a ray tracing code for stellaratorshttps://github.com/ proximafusion/raytrax

  6. [13]

    Ii Tsujimura T, Kubo S, Takahashi H, Makino R, Seki R, Yoshimura Y, Igami H, Shimozuma T, Ida K, Suzuki C, Emoto M, Yokoyama M, Kobayashi T, Moon C, Nagaoka K, Osakabe M, Kobayashi S, Ito S, Mizuno Y, Okada K, Ejiri A, Mutoh T and the LHD Experiment Group 2015Nuclear Fusion55123019

  7. [14]

    Yanai R, Tsujimura T, Kubo S, Yoneda R, Yoshimura Y, Nishiura M, Igami H, Takahashi H and Shimozuma T 2021Plasma and Fusion Research162402084

  8. [15]

    Tribaldos V, Jim´ enez J A, Guasp J and van Milligen B P 1998Plasma Physics and Controlled Fusion402113–2130

  9. [16]

    Castej´ on F, Tereshchenko M A, Sarksyan K A, Fern´ andez´A, Cappa ´A, Batanov G M, Sakharov A S and Mart ´ ın R 2004Fusion Science and Technology46327–334

  10. [17]

    Yanagihara K, Kubo S, Dodin I Y and the LHD Experiment Group 2021Nuclear Fusion61 106012

  11. [18]

    Yu W, Xie H, Mao A, Ma H, Liu T and Wang Z x 2026Physics of Plasmas33042511

  12. [19]

    R¨ onnmark K 1983Plasma Physics25699–701

  13. [20]

    Xie H 2024AIP Advances14075007

  14. [21]

    R¨ onnmark K 1982 WHAMP: Waves in homogeneous, anisotropic, multicomponent plasmas Tech. Rep. 179 Kiruna Geophysical Institute Kiruna, Sweden 11 IOP PublishingJournalvv(yyyy) aaaaaa Wang and Xie

    1982

  15. [22]

    Albajar F, Bertelli N, Bornatici M and Engelmann F 2007Plasma Physics and Controlled Fusion 4915–29

  16. [23]

    Krivenski V and Orefice A 1983Journal of Plasma Physics30125–131

  17. [24]

    Saito K, Seki T, Kumazawa R, Kasahara H, Shimpo F, Nomura G and Mutoh T 2010Plasma and Fusion Research5S1030

  18. [25]

    Likin K M, Abdou A, Almagri A F, Anderson D T, Anderson F S B, Brower D, Canik J, Deng C, Gerhardt S P, Guttenfelder W, Oh S, Radder J, Sakaguchi V, Schmitt J, Tabora J, Talmadge J N and Zhai K 2003Plasma Physics and Controlled Fusion45A133–A142

  19. [27]

    Landreman M, Medasani B, Wechsung F, Giuliani A, Jorge R and Zhu C 2021Journal of Open Source Software63525

  20. [28]

    LHD LINES OF FORCE.EXE

    Watanabe T, Yoshida M, Masuzaki S, Emoto M and Nagayama Y 2006 How to use the numerical system “LHD LINES OF FORCE.EXE” for the visualization of lines of force in the large helical device Tech. Rep. NIFS-TECH-014 National Institute for Fusion Science URLhttps://www. nifs.ac.jp...

    2006

  21. [29]

    Wang E, Morita S, Goto M, Kobayashi M and Dong C 2012Plasma and Fusion Research7 2402059

  22. [30]

    Qian T M, Chu X, Pagano C, Patch D, Zarnstorff M C, Berlinger B, Bishop D, Chambliss A, Haque M, Seidita D and Zhu C 2023Journal of Plasma Physics89955890502

  23. [31]

    Granetzny M and Schmitz O 2025Plasma Sources Science and Technology34115002

  24. [32]

    Marushchenko N B, Dinklage A, Hartfuss H J, Hirsch M, Maassberg H and Turkin Y 2006 Fusion Science and Technology50395–402

    2006

  25. [33]

    Zanini M, Laqua H P, Thomsen Het al.2020Nuclear Fusion60106021

  26. [34]

    Xie H, Ma H and Bai Y 2021 Plasma waves accessibility diagrams: A tutorial to include the fluid and kinetic thermal effects (Preprint2111.05669) URLhttps://arxiv.org/abs/2111.05669 12

    2021 arXiv

Pith tools

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