{"id":"cbfbfc62-3ea1-4dc0-898b-e6b0a6294ba0","arxiv_id":"2608.05667","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"BORAY-3D is a new ray tracing code that combines broad frequency coverage, three-dimensional magnetic geometry, and relativistic electron cyclotron absorption in a single framework.","lead":"The authors present BORAY-3D, a ray tracing code that simulates radio frequency wave propagation in magnetized plasmas, covering frequencies from ion cyclotron to electron cyclotron in both tokamak and stellarator geometries. It unifies previously separate capabilities for broad frequency coverage, three-dimensional magnetic configurations, and relativistic electron cyclotron absorption in a single tool.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 'fully relativistic Maxwellian EC absorption' claim is internally inconsistent: the polarization entering the absorption coefficient Eq. (20) is computed from the weakly relativistic Krivenski–Orefice tensor (Sec. 2.3, Appendix A), so the absorption is a hybrid model, not fully relativistic.","rationale":"The central claim is that BORAY-3D is the first code to combine broad frequency coverage, arbitrary 3D geometry, and fully relativistic EC absorption. The cold-ray/hot-absorption split identified by the reader is a real limitation, but it is explicitly acknowledged in Sec. 4 and is shared by other ray-tracing codes (e.g., GENRAY), so it does not distinguish BORAY-3D from the state of the art. The relativistic polarization inconsistency is more load-bearing because it directly contradicts the third headline feature: the absorption coefficient in Eq. (20) is built from polarization components obtained from a weakly relativistic tensor. If 'fully relativistic absorption' is a key selling point, the implementation as described cannot support it. The concern is not merely semantic: α_ω is used for ECE and heating, and Fig. 7 shows the X3 band is sensitive to the absorption model. A concrete check with a fully relativistic polarization would quantify the error. The paper is otherwise carefully written, and the trajectory benchmarks are strong; the verdict stays CONDITIONAL, but the condition should require either a quantitative demonstration that the weakly relativistic polarization does not affect the results or a revised claim. The code is not released, but that is a reproducibility issue, not a physics flaw. Overall, the most load-bearing concern is the unsubstantiated 'fully relativistic' qualifier.","tokens_in":13562,"tokens_out":9046,"duration_ms":94515,"concrete_test":"Recompute the W7-X 140 GHz O-mode absorption and the X2/X3 ECE spectra of Fig. 7 using a polarization vector obtained from a fully relativistic dielectric tensor (e.g., the exact relativistic Maxwellian response of Bornatici et al. [1,22]) instead of the weakly relativistic Krivenski–Orefice tensor, keeping Eq. (20) unchanged. If the absorbed fractions or peak radiative temperatures shift by more than ~5%, the weakly relativistic polarization is quantitatively important and the 'fully relativistic' claim fails. If the results are indistinguishable, the authors should still rewrite the claim as 'relativistic resonance integral with weakly relativistic polarization' to avoid overstatement.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's third headline capability is 'fully relativistic Maxwellian EC absorption' (abstract; Sec. 2.3). However, Sec. 2.3 states that the wave polarization is calculated from the weakly relativistic Krivenski–Orefice tensor [23], and Appendix A provides the corresponding weakly relativistic Shkarofsky-function tensor (Eqs. 27–29). The absorption coefficient α_ω in Eq. (20) depends on the polarization components e_x, e_y, e_z through the harmonic coupling |C_n|^2 in Eq. (21). Consequently, the computed α_ω is not fully relativistic; it combines a fully relativistic resonance condition and Maxwell–Jüttner integral with a weakly relativistic polarization. This undercuts the novelty claim of the third feature, since many existing EC codes use weakly relativistic polarization with relativistic resonance. The practical risk is quantitative: the weakly relativistic polarization can deviate near cyclotron resonances, especially for higher harmonics (X3) or narrow resonances. The W7-X O-mode absorbed fractions in Sec. 3.2 spread from 73.3% to 87.1% across codes with different kinetic models, and Fig. 7 shows the X3 optical depth is sensitive to the absorption model; part of that spread may originate from the polarization approximation. At minimum, the phrase 'fully relativistic' is not supported by the described implementation.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":13760,"tokens_out":11636,"duration_ms":115010,"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":[{"comment":"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.","section":"Sec. 2.3, abstract"},{"comment":"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.","section":"Sec. 3.5, Table 1"},{"comment":"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.","section":"Sec. 3.2, Fig. 4"},{"comment":"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.","section":"Sec. 3.2, Fig. 2"}],"minor_comments":[{"comment":"The phrase 'an unified treatment' should be 'a unified treatment'.","section":"Abstract"},{"comment":"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.","section":"Eq. (20)"},{"comment":"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.","section":"Sec. 3.4, Fig. 6"},{"comment":"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.","section":"Sec. 3.2, Fig. 2"},{"comment":"Ref. [18] is listed as 'Yu et al. 2026'; please confirm the publication status and provide a DOI if available.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The main issue is the overstated 'fully relativistic' claim, which is a central novelty point and needs to be fixed either by implementation or by a careful revision of the wording. The ECE benchmark in Sec. 3.5 is also misrepresented as a validation when no comparison with the cited public data is shown. These are fixable within the manuscript's scope, so I recommend major revision rather than rejection. I do not see evidence of parameter tuning or circular benchmarking; the comparisons with GENRAY and Raytrax are credible for trajectories."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a solid code paper that genuinely extends BORAY to 3D equilibria and adds relativistic EC absorption, but the words 'fully relativistic' are doing more work than the implementation supports. The polarization tensor is weakly relativistic, so the absorption model is hybrid. That matters for the novelty claim and for near-resonance accuracy.\n\nWhat's new and good: the ray equations in cylindrical coordinates with varying n_phi (Sec. 2.1) are a real extension of axisymmetric BORAY, and the code handles arbitrary (r,phi,z) profiles with closed or open field lines. The benchmark suite is well chosen: GENRAY for LHW with ripple, Raytrax and TRAVIS for W7-X O-mode, plus HSX heating and LHD ICW. The trajectory agreements are very good — within 1–7 mm in the tokamak ripple case and 5 mm in W7-X — so the propagation part is convincing. The physical trends in the LHD and MUSE cases also look sensible.\n\nSoft spots, in proportion. First, the 'fully relativistic Maxwellian EC absorption' headline claim is not supported by Sec. 2.3. The absorption coefficient in Eq. (20) uses the polarization from the weakly relativistic Krivenski–Orefice tensor, derived in Appendix A. The resonance condition and Maxwell–Jüttner average are fully relativistic, but the damping depends on the polarization through the harmonic coupling. The model is hybrid. Many production codes do the same, but it should be labeled as such; the abstract and summary should say 'relativistic resonance with weakly relativistic polarization,' or the code should adopt a fully relativistic tensor.\n\nSecond, the split between cold ray and hot absorption is only indirectly acknowledged. The paper notes the geometric optics limit, but not that using a cold ray for the trajectory can misplace deposition near cutoffs or resonances, especially for low-frequency waves. The HSX first-pass absorption (0.29 vs 0.40) and the W7-X spread (73–87%) are consistent with that caveat, but the paper attributes all differences to 'different kinetic models' without testing the cold-ray sensitivity.\n\nThird, reproducibility is limited: no code release, and some benchmark data are unavailable. The HSX and LHD comparisons are trend-level. That is common for code papers, but it raises the bar for validation claims.\n\nBottom line: the propagation capability is well validated and the 3D extension is real. The absorption modeling is less clean, and the 'fully relativistic' wording needs correction. I would send this to a serious referee; it is a useful contribution to RF modeling. I would ask for the code or a test-case repository and a rewriting of the absorption claims.","headline":"A useful code paper that extends BORAY to 3D, but the 'fully relativistic' absorption claim is overstated because the polarization stays weakly relativistic.","tokens_in":14396,"tokens_out":3434,"would_cite":true,"duration_ms":32205,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["52.65.-y","52.50.Qt","52.25.Os"],"model":"deepseek-v4-flash","headline":"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…","keywords":["ray tracing","radio frequency waves","three-dimensional equilibrium","broad-frequency-range modeling","relativistic absorption","electron cyclotron emission","stellarator","tokamak"],"falsifier":"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.","tokens_in":13258,"feed_emoji":"⚡","tokens_out":4924,"duration_ms":48924,"temperature":0.7,"pith_summary":"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.","feed_headline":"One ray tracer now spans RF heating to ECE in 3D plasmas","feed_subtitle":"It reproduces established codes' trajectories within millimeters while covering ion cyclotron to EC emission.","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"The original BORAY code supplies the axisymmetric ray-tracing formulation and the non-relativistic hot-plasma absorption model that BORAY-3D extends to three dimensions.","marker":"[4]"},{"why":"Provides the circular-tokamak parameters and the analytic 18-coil ripple model used in the GENRAY comparison for lower-hybrid waves.","marker":"[9]"},{"why":"TRAVIS supplies the stellarator ECW/ECE modeling formalism and the benchmark trajectories against which BORAY-3D's W7-X results are compared.","marker":"[10]"},{"why":"Raytrax is the open-source stellarator ray tracing code used as the second independent benchmark for W7-X electron-cyclotron trajectories and absorption.","marker":"[12]"},{"why":"Supplies the fully relativistic Maxwellian absorption integral that underlies the EC absorption model and the ECE radiative transfer calculation.","marker":"[22]"},{"why":"Krivenski–Orefice tensor is used to compute the wave polarization for the relativistic absorption model.","marker":"[23]"},{"why":"HSX 28 GHz X2 heating results from Likin et al. provide the experimental density-dependence and deposition data used for the HSX comparison.","marker":"[25]"},{"why":"W7-X ECE measurements and receiver geometry from Marushchenko et al. are used to define the ECE simulation setup and to compare the computed spectra.","marker":"[26]"}],"fun_headline_variants":["BORAY-3D unifies RF heating and ECE in 3D plasmas","From ion cyclotron to ECE: one ray tracer for 3D","Relativistic ECE and broad RF now in a single 3D ray code","Ray tracing code spans 13 MHz to 220 GHz in 3D"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["BORAY-3D unifies RF heating and ECE in 3D plasmas","From ion cyclotron to ECE: one ray tracer for 3D","Relativistic ECE and broad RF now in a single 3D ray code","Ray tracing code spans 13 MHz to 220 GHz in 3D"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000689,"raw_usage":{"total_tokens":3166,"prompt_tokens":1031,"completion_tokens":2135,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":647,"completion_tokens_details":{"reasoning_tokens":2058}},"tokens_in":647,"tokens_out":2135,"duration_ms":15428,"temperature":1.0,"reasoning_tokens":2058,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-08T04:30:02.612111+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The original BORAY code supplies the axisymmetric ray-tracing formulation and the non-relativistic hot-plasma absorption model that BORAY-3D extends to three dimensions."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the circular-tokamak parameters and the analytic 18-coil ripple model used in the GENRAY comparison for lower-hybrid waves."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"TRAVIS supplies the stellarator ECW/ECE modeling formalism and the benchmark trajectories against which BORAY-3D's W7-X results are compared."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Raytrax is the open-source stellarator ray tracing code used as the second independent benchmark for W7-X electron-cyclotron trajectories and absorption."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the fully relativistic Maxwellian absorption integral that underlies the EC absorption model and the ECE radiative transfer calculation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Krivenski–Orefice tensor is used to compute the wave polarization for the relativistic absorption model."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"HSX 28 GHz X2 heating results from Likin et al. provide the experimental density-dependence and deposition data used for the HSX comparison."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"W7-X ECE measurements and receiver geometry from Marushchenko et al. are used to define the ECE simulation setup and to compare the computed spectra."}],"review_version":1}