{"id":"3a76d58c-3c89-4893-b0eb-efd6ed82882b","arxiv_id":"2501.18193","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"For SST-1, refractive beam deflection limits the interferometer wavelength to below 1.2 mm, and at 432.6 μm the deflection stays within the beam diameter for vertical and lateral chords.","lead":"This paper calculates how much a far-infrared probing beam bends when it crosses the plasma in the SST-1 tokamak, for circular, elliptical, and D-shaped cross sections, and concludes the probing wavelength should stay below 1.2 mm. It is a diagnostic design check: if the bending is too large, the beam misses the detector and density measurements fail.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central 1.2 mm limit is not reproducible because Fig. 3 compares refraction shift with beam diameter without giving the beam waist, focusing, or port-to-window distances; a Gaussian beam diameter is itself wavelength-dependent, so the crossover can move substantially under different optics.","rationale":"The reader's weakest assumption is the density-profile model, and that is a real limitation: the deflection pattern and the 1.2 mm limit are computed for n0 = 3e19 m^-3 and the specified D-shape parametrization, not for measured SST-1 profiles. My concern is more immediate: even if the density model is taken as given, the numerical value of the wavelength cutoff is not derivable from the text because the comparison criterion in Fig. 3 requires an unstated beam-diameter model. A beam diameter at a window or retroreflector is a function of the beam waist, focusing arrangement, and wavelength, and different choices move the crossover. This directly controls the central recommendation that the 432.6 micron laser is safe. In addition, Eqs. (4)-(5), used to validate the ray tracer in Fig. 5, appear dimensionally wrong as printed, so the only internal check on the numerics is not reproducible. The paper makes a plausible device-specific design claim, and a conditional verdict is appropriate, but the missing optical-beam specification and the misprinted analytical formulas need to be corrected before the 1.2 mm limit can be accepted as stated.","tokens_in":4816,"tokens_out":13111,"duration_ms":137636,"concrete_test":"Reproduce Fig. 3 using the actual SST-1 port geometry and the Gaussian beam parameters of the 432.6 micron HCOOH interferometer from Ref. [6]: specify w0 and focal plane, compute d(lambda) = 2 w(z, lambda) at the exit window and at the retroreflector, and find the crossover with the ray-traced maximum shift. Then vary w0 by +/-30% and repeat; if the crossover drops below 432.6 microns, or shifts from 1.2 mm by more than 20%, the paper's design limit is not robust. Independently re-derive Eq. (4) from the parabolic profile; if it is not corrected to alpha = 2 (n_e/n_c) x sqrt(r0^2 - x^2)/r0^2, the analytical validation cannot be used.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The claim that vertical and lateral SST-1 interferometry must stay below lambda = 1.2 mm rests on Sec. 5/Fig. 3: beam displacement due to refraction is compared with the beam diameter d at the first optical component after the beam exits the plasma, and the wavelength where the two cross is quoted as 1.2 mm. The paper never states the values used for d, the waist w0 and focal position of the probing beam, or the distances from the plasma edge to the exit window and to the retroreflectors. For a real Gaussian beam, d is not fixed: at a given component it depends on w0, focusing geometry, and lambda through diffraction. Different plausible focusing choices can move the crossover by a large factor, and the claim that 432.6 microns is safe depends on that choice. A related weakness is that the analytical formulas used to validate the ray tracer, Eqs. (4)-(5), are printed with wrong dimensions (the factor sqrt(r0^2 - x^2) appears in the denominator instead of the numerator), so the validation in Fig. 5 cannot be checked as written. Both issues leave the quantitative central claim unsupported by the manuscript alone.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper computes refractive beam deflection for interferometry beams in circular, elliptical, and D-shaped plasma cross-sections, using analytical formulas for the first two geometries and ray tracing (eikonal approximation) for the D-shaped SST-1 plasma. It presents deflection results for a 432.6 μm HCOOH probing beam and concludes that the probing wavelength for vertical and lateral SST-1 interferometry should not exceed 1.2 mm, based on comparing the maximum refraction-induced beam displacement with the beam diameter at the first optical component. The manuscript also compares ray-tracing results with analytical formulas for validation.","tokens_in":5138,"tokens_out":5010,"duration_ms":45440,"significance":"If the central 1.2 mm limit were reproducible from the manuscript alone, this would be a useful design constraint for the SST-1 far-infrared interferometer and a reasonable demonstration of ray-tracing methodology for D-shaped plasmas. The work has the strength of using no fitted parameters: n0, γ, κ, and δ are stated as design inputs, and the ray tracing equations are standard. However, the load-bearing quantitative claim is currently not verifiable because the manuscript omits the optical layout and beam parameters used to define the beam diameter, and because the analytical formulas in Eqs. (4) and (5) are dimensionally inconsistent as printed, undermining the validation shown in Fig. 5. The general approach is sound, but the paper does not yet support its headline wavelength limit.","major_comments":[{"comment":"Equations (4) and (5) are dimensionally inconsistent as printed. In Eq. (4), the right-hand side has units of inverse length (x/(r0^2 sqrt(r0^2-x^2)) ≈ L/(L^2 L) = 1/L), whereas α must be dimensionless. The standard result for the parabolic profile in Eq. (3) is α = 2 (n_e/n_c) x sqrt(r0^2 - x^2)/r0^2, with sqrt(r0^2 - x^2) in the numerator. The same missing factor appears in Eq. (5). Because Fig. 5 validates the ray tracer against these analytical formulas, the validation cannot be checked as written. The authors should correct the formulas and either show the integration leading to them or cite a source.","section":"Section 3, Eqs. (4) and (5)"},{"comment":"The central claim that λ < 1.2 mm is derived by comparing the maximum refraction-induced beam displacement with the beam diameter d at the first optical component, but the manuscript never specifies d, the beam waist w0, the focusing geometry, or the distances from the plasma edge to the exit window (vertical viewing) and to the retroreflectors (lateral viewing). Since a Gaussian beam diameter is itself wavelength-dependent and the displacement scales with propagation distance, the crossover point at 1.2 mm cannot be reproduced from the paper alone. Please state the assumed optical parameters and geometry, or provide a sensitivity analysis showing how the limit changes with reasonable variations in these quantities.","section":"Section 5, Fig. 3"},{"comment":"The conclusion that 'the probing wavelength for vertical and lateral viewing of interferometer should not exceed 1.2 mm' is stated without any uncertainty or sensitivity range. Given that the design criterion depends on a somewhat arbitrary choice of 'beam diameter' and on the density profile model (γ between 0.5 and 2), the paper should either quantify how the limit varies across that range or explicitly state that the 1.2 mm value corresponds to a particular profile and optics choice. As written, the headline limit is presented as a universal property of SST-1, which is not justified by the information provided.","section":"Section 5, Fig. 3 and Section 6"}],"minor_comments":[{"comment":"The notation α(x) is used in Eqs. (4) and (5), but x is a scalar impact parameter while x is also used as a coordinate vector elsewhere; consider using b for the impact parameter for clarity.","section":"Section 3"},{"comment":"The captions of Figs. 4 and 5 are nearly identical and both refer to 'vertical paths'; the text should clarify that Fig. 4 is for the D-shaped SST-1 plasma only, while Fig. 5 compares circular, elliptical, and D-shaped geometries. The caption of Fig. 6 also appears truncated ('...at the pivot point y tracing calculations').","section":"Captions of Figs. 4 and 5"},{"comment":"The symbol n_e in Eqs. (2), (4), and (5) is used for both the local density and the central density; for the parabolic profile of Eq. (3), the central density is n_o. Please distinguish these to avoid confusion.","section":"Section 3, Eq. (2)"}],"recommendation":"major_revision","confidential_remarks":"The paper's core methodology is standard and the topic is suitable for a plasma diagnostics journal, but the lack of optical layout parameters and the dimensionally inconsistent analytical formulas make the main quantitative claim unsupported as submitted. The issues are fixable without changing the scope, so major revision seems appropriate rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: the genuinely new piece is the D-shaped, finite-κ, finite-δ ray-tracing calculation for SST-1, and the engineering takeaway is that 432.6 μm HCOOH is safe while the wavelength should stay below 1.2 mm. The ray-tracing equations and the validation strategy are sound—reducing the D-shape to circular and elliptical limits and checking against analytical points is a fair self-consistency test, and no parameter is fitted to the final deflection claim. The heavy self-citation to the 1995 design report is appropriate; that is the actual design context.\n\nThe soft spots are real but manageable. Equations (4) and (5), as printed, are dimensionally wrong: they give 1/length, not radians. The missing factor is a chord length in the numerator, so these are almost certainly transcription errors, but as published they cannot be checked. The derivation is not shown either; since these are standard small-deflection formulas, a citation or a two-line derivation would fix it. More importantly, the central 1.2 mm limit in Fig. 3 compares refraction shift with beam diameter d at the first optical component, but the paper never states d, the beam waist, focusing geometry, or the port-to-window distances. For a Gaussian beam d is itself wavelength-dependent, so the crossover can move substantially under different optics. The authors need to supply those numbers before the limit can be verified.\n\nSmaller points: no code or raw numerical data are provided, so the ray-tracing results in Figs. 4–6 cannot be reproduced independently; Fig. 6 has a caption typo; and the claim about sharper inboard gradients is reasonable but not backed by a quantitative comparison.\n\nOverall this is a legitimate engineering-validation paper, not new physics. It deserves a serious referee because the method is standard and the SST-1 result is useful to the plasma-diagnostics community, but the authors should be asked to correct Eqs. (4)–(5), provide the optics parameters behind Fig. 3, and release the numerical tables. As written, the central quantitative claim is plausible but under-supported.","headline":"Useful SST-1 design calculation with a standard ray-tracing method; the specific 1.2 mm cutoff is plausible but not yet reproducible from the manuscript as written.","tokens_in":5586,"tokens_out":2310,"would_cite":false,"duration_ms":22971,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Refraction by transverse density gradients in SST-1 sets the upper usable wavelength for interferometry at $\\lambda < 1.2$ mm, so the 432.6 $\\mu$m HCOOH laser is compatible with the machine geometry.","keywords":["plasma diagnostics","D-shaped plasma","transverse gradient","refraction","ray tracing","SST-1 tokamak","far-infrared interferometry","beam deflection"],"falsifier":"Probe SST-1 with a beam at $\\lambda$ just below 1.2 mm during a discharge reaching $n_0=3\\times10^{19}$ m$^{-3}$ with a peaked profile; if the beam displacement at the exit window (vertical) or inner-wall reflector (lateral) exceeds the beam diameter, the stated limit is wrong. Conversely, a deflection measurement at 432.6 $\\mu$m that matches the ray-tracing curves would support the limit.","tokens_in":4622,"feed_emoji":"⚛️","tokens_out":6389,"duration_ms":55421,"temperature":0.7,"pith_summary":"Far-infrared interferometry measures electron density through phase shift, but the plasma's transverse density gradients refract the probing beam. This paper asks whether the beam deflection stays small enough for the SST-1 tokamak geometry, and answers with a clear design rule: for both single-pass vertical and double-pass lateral viewing, the probing wavelength must stay below $\\lambda < 1.2$ mm. The authors derive the refraction limit analytically for circular and elliptical plasmas, then numerically ray-trace the D-shaped SST-1 plasma. Their deflection curves for a 432.6 $\\mu$m laser sit comfortably below the beam diameter at the exit window and at the retroreflectors, which is why that source is considered suitable. If the calculations are right, the same method gives a simple, reusable compatibility test for any interferometer wavelength and density profile.","feed_headline":"Refraction caps SST-1 probe wavelength at 1.2 mm","feed_subtitle":"Ray tracing and analytic calculations show the 432.6 μm HCOOH laser stays safe for density measurements.","key_machinery":"The quantitative engine is a set of four coupled ray-tracing ODEs derived from the Eikonal approximation of the ordinary-mode dispersion relation $D = k^2 c^2 + \\omega_p^2 - \\omega^2$. The electron density is modeled on flux surfaces as $n(R,z)=n_0(1-\\rho^2/a^2)^\\gamma$, with $\\rho$ defined by $R=R_0+\\rho\\cos(\\theta+\\delta\\sin\\theta)$ and $z=\\kappa\\rho\\sin\\theta$, taking $\\kappa=1.7$, $\\delta=0.6$, $n_0=3\\times10^{19}$ m$^{-3}$, and $\\gamma\\in[0.5,2]$. The deflection is the angle between the initial and final wave vectors, and the key design comparison is between the resulting beam displacement and the beam diameter at the first optical component. That comparison is what turns a refraction curve into the wavelength bound $\\lambda < 1.2$ mm.","core_discovery":"The paper's central claim is that the maximum usable probing wavelength for SST-1 interferometry is set by refractive beam displacement, not by any other instrument constraint, and that this limit is $\\lambda < 1.2$ mm. The displacement is computed by comparing the transverse shift of the beam at the first optical surface (exit window for vertical chords, inner-wall retroreflectors for lateral chords) with the beam diameter. For the 432.6 $\\mu$m and 337 $\\mu$m far-infrared lasers considered, the computed shift remains smaller than the beam diameter for all three density peaking cases ($\\gamma = 0.5$, $1$, $2$) at $n_0 = 3\\times10^{19}$ m$^{-3}$. The ray-tracing calculation is validated by reducing the D-shape to circular and elliptical cross-sections and matching the analytic deflection formulas.","pith_inferences":["The 1.2 mm limit is tied to the assumed $n_0 = 3\\times10^{19}$ m$^{-3}$; if SST-1 later operates at higher central densities, the usable wavelength would shrink, possibly below the 432.6 $\\mu$m choice.","The analytic circular and elliptical formulas could be inverted to infer line-averaged density gradients from measured beam displacement, turning the parasitic effect into a diagnostic.","The criterion 'displacement $\\le$ beam diameter' is conservative for Gaussian beams, since part of the beam edge can still reach the detector; a threshold based on the detector's field of view might allow somewhat longer wavelengths.","The validation covers only axisymmetric limits; extending the approach to measured, turbulence-distorted density profiles would test the safety margin more realistically."],"forward_implications":["The 432.6 $\\mu$m HCOOH laser and the 337 $\\mu$m HCN laser both remain below the 1.2 mm limit for SST-1's design parameters, so the multichannel interferometer can be built without changing optics.","The computed deflection pattern gives the maximum expected transverse shift at each vertical and lateral chord, guiding port and detector placement to avoid cross talk.","For peaked density profiles ($\\gamma = 2$) the deflection is largest, and in D-shaped plasmas it is asymmetric with stronger inboard-side bending, so channel spacing must account for this asymmetry.","The same beam-diameter-versus-displacement criterion can be applied to any future probing wavelength or density profile to test compatibility before installation.","Refraction imposes not only an upper wavelength limit but also a constraint on where the beam can enter and exit, since compensating for the shift during a discharge is impractical."],"supporting_citations":[{"why":"Establishes the phase-shift principle of multichannel far-infrared interferometry that motivates the whole calculation.","marker":"[1]"},{"why":"Provides the multichannel interferometer/polarimeter system context from another tokamak, supporting the diagnostic design.","marker":"[2]"},{"why":"Introduces the refractive bending constraint $\\nabla n_e L < 10^{14}$ cm$^{-3}$ that sets the upper-limit philosophy.","marker":"[3]"},{"why":"Gives the analytical line-integration method for refraction in cylindrical plasma that the paper extends to elliptical cross-sections.","marker":"[4]"},{"why":"Supplies a two-color interferometer ray-tracing comparison on another tokamak, the numerical approach the paper follows.","marker":"[5]"},{"why":"Sets the lower wavelength limit of 120 $\\mu$m and is the SST-1 interferometer conceptual design whose compatibility is being assessed.","marker":"[6]"},{"why":"Defines the SST-1 tokamak and its D-shaped cross-section, the geometry used in the ray tracing.","marker":"[7]"},{"why":"Provides the Eikonal ray-tracing equations from which the four coupled ODEs in Section 4 are derived.","marker":"[11]"}],"fun_headline_variants":["1.2 mm ceiling for SST-1 interferometry beams","Refractive deflection sets SST-1 probe wavelength cap","SST-1 THz beams safe up to 1.2 mm wavelength","Bending limits wavelength for SST-1 density probe"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The computed deflections assume a smooth analytic density profile $n(R,z)=n_0(1-\\rho^2/a^2)^\\gamma$ with $n_0=3\\times10^{19}$ m$^{-3}$, $\\gamma$ between 0.5 and 2, $\\kappa=1.7$, $\\delta=0.6$; these are model inputs, not measured SST-1 profiles, and a different peaking, shape, or density would change the deflection curves and the 1.2 mm limit.","fun_headline_variants_meta":{"raw":{"variants":["1.2 mm ceiling for SST-1 interferometry beams","Refractive deflection sets SST-1 probe wavelength cap","SST-1 THz beams safe up to 1.2 mm wavelength","Bending limits wavelength for SST-1 density probe"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000605,"raw_usage":{"total_tokens":2758,"prompt_tokens":819,"completion_tokens":1939,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":435,"completion_tokens_details":{"reasoning_tokens":1868}},"tokens_in":435,"tokens_out":1939,"duration_ms":13800,"temperature":1.0,"reasoning_tokens":1868,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T00:20:16.388856+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Probe SST-1 with a beam at $\\lambda$ just below 1.2 mm during a discharge reaching $n_0=3\\times10^{19}$ m$^{-3}$ with a peaked profile; if the beam displacement at the exit window (vertical) or inner-wall reflector (lateral) exceeds the beam diameter, the stated limit is wrong. Conversely, a deflection measurement at 432.6 $\\mu$m that matches the ray-tracing curves would support the limit.","supporting_citations":[{"cited_title":"M ultichannel far-infrared laser interferometer for electron density measurements on the tokamak fusion test reactor,","cited_arxiv_id":null,"evidence_quote":"Establishes the phase-shift principle of multichannel far-infrared interferometry that motivates the whole calculation."},{"cited_title":"Multichannel interferometer/polarimeter system for the RTP tokamak,","cited_arxiv_id":null,"evidence_quote":"Provides the multichannel interferometer/polarimeter system context from another tokamak, supporting the diagnostic design."},{"cited_title":"CO2 interferometer operation in Dou ble III,","cited_arxiv_id":null,"evidence_quote":"Introduces the refractive bending constraint $\\nabla n_e L < 10^{14}$ cm$^{-3}$ that sets the upper-limit philosophy."},{"cited_title":"Proposed Diagnostic Method for Cylindrical Plasmas,","cited_arxiv_id":null,"evidence_quote":"Gives the analytical line-integration method for refraction in cylindrical plasma that the paper extends to elliptical cross-sections."},{"cited_title":"Two-color interferometer system for Alcator C- MOD,","cited_arxiv_id":null,"evidence_quote":"Supplies a two-color interferometer ray-tracing comparison on another tokamak, the numerical approach the paper follows."},{"cited_title":"Conceptual design of a far infrared interferometer for SST-1","cited_arxiv_id":null,"evidence_quote":"Sets the lower wavelength limit of 120 $\\mu$m and is the SST-1 interferometer conceptual design whose compatibility is being assessed."},{"cited_title":"Present Status of SST-1 project,","cited_arxiv_id":null,"evidence_quote":"Defines the SST-1 tokamak and its D-shaped cross-section, the geometry used in the ray tracing."},{"cited_title":"Eikonal Method in Magnetohydrodynamics,","cited_arxiv_id":null,"evidence_quote":"Provides the Eikonal ray-tracing equations from which the four coupled ODEs in Section 4 are derived."}],"review_version":1}