{"id":"8da7348d-0358-494e-9918-a07c5523bec2","arxiv_id":"2608.08649","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"The helicon and Trivelpiece-Gould wave branches are shown to be two lobes of a single degenerate vacuum root, and their ignition-time splitting, thresholds, and non-adiabatic response are derived analytically.","lead":"This paper gives an analytical theory for how the wave dispersion relation of a helicon plasma discharge evolves from vacuum to steady operation during ignition. It derives threshold densities, a magnetic-field lower bound for clean helicon coupling, and separates local from global adiabatic effects, with testable predictions.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 'pinned helicon root' claim is an internal inconsistency: the lower root of the full quartic stays near the vacuum value in k_perp^2, not in beta=(k_perp^2+kz^2)^{1/2}, and Table 4's low-density rows use the EMHD reduction that Section III.F says is invalid at ignition.","rationale":"The reader's verdict is CONDITIONAL, and that remains the right level: the discriminant cancellation in Section III.C is real, the vacuum double root is elementary but correctly identified, and the collisionally demagnetized regime is explicitly discussed in Section III.E. The reader's weakest-assumption concern about the cold collisional tensor and Landau damping is legitimate, but it is partly anticipated by the paper's own demagnetization and limitation statements. The more concrete and load-bearing problem is internal: the pinning claim is stated for beta=(k_perp^2+kz^2)^{1/2}, but the actual full-quartic behavior is that k_perp^2 stays near the vacuum value -kappa^2 while beta itself does not, because beta^2 is the small difference of two large numbers. Table 4's low-density rows appear to use the EMHD beta_- that Section III.F explicitly excludes during ignition, so the quantitative early-transient figures inherit that inconsistency. This does not overturn the structural double-root result, but it does mean the paper's quantitative claims about the pinned helicon root and the early eigenfunctions are not currently supported by a consistent set of definitions. The proposed check settles the discrepancy directly. Since the reader already demanded clarification and independent confirmation, no verdict change is needed.","tokens_in":24035,"tokens_out":37304,"duration_ms":380522,"concrete_test":"Recompute the two roots of Eq. (9) using Eq. (7) at ne = 8.3e14 m^-3 and at the first two rows of Table 4 (ne=1e14, 3e9 s^-1 and ne=1e15, 3e9 s^-1), using the same full-quartic code that produced Table 6. For each case report beta_Helicon = (k_perp^2 + kz^2)^{1/2} and compare with k0=0.284 m^-1 and with the tabulated beta_Helicon. If the full quartic does not give beta within 1% of k0 at ne=8.3e14, or if the tabulated entries reproduce only the EMHD beta_- of Eq. (24) rather than the full quartic roots, then the pinning claim and the early-transient quantitative results require correction or explicit relabeling of k_perp^2 as the pinned quantity.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section III.D and the abstract claim that the lower root remains pinned at beta=k0 to within 1% until ne~8.3e14 m^-3, with beta defined in Section III.D as (k_perp^2+kz^2)^{1/2}. Direct evaluation of Eq. (9) with Eq. (7) in the collisionless limit contradicts this. At ne=0.1 n1, the lower root has k_perp^2 ~ -988.7 m^-2, so beta^2 = k_perp^2 + kz^2 ~ -2.7 m^-2, i.e. the root is radially evanescent rather than pinned at beta=k0~0.284 m^-1. At ne=8.3e14 m^-3 the lower root gives beta ~ 0.46 m^-1, again not within 1% of k0. What actually remains close to its vacuum value is k_perp^2, not beta; beta^2 = k0^2 + delta with delta small only against kz^2, not against k0^2. Table 4 is harder to reconcile: the first rows list beta_Helicon = 0.0+0.00i at ne=1e14 and 0.1+0.00i at ne=1e15, while the full quartic at those densities gives a real part of beta not near k0. Numerically those table entries match the EMHD root beta_- = omega mu0 n_e e/(k_z B0), which at low density is small and is exactly the reduction that Section III.F (Eq. 27) shows fails for ne < n3 because it does not reproduce the vacuum root. The central double-root framework and the linear splitting are not questioned, but the pinning statement, Figure 1, and Table 4 rely on incompatible definitions or reductions unless a mislabeling is corrected.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops an analytical cold-collisional two-species fluid theory for the evolution of the helicon/TG dispersion relation during discharge ignition. Beginning from the standard Stix tensor, it shows that the vacuum limit of the perpendicular-wavenumber quartic is a degenerate double root, that the leading discriminant term cancels so the branches split linearly in density, and that the helicon root remains near its vacuum value until n_e~10^15 m^-3 while the TG root departs immediately. It introduces a ladder of characteristic densities, a resonance-coalescence ratio for bounded modes, a two-timescale adiabaticity analysis separating local dielectric accuracy from global cavity response, and a set of numerical benchmarks and experimental predictions.","tokens_in":26,"tokens_out":29776,"duration_ms":841590,"significance":"If the central derivation is accepted, the double-root viewpoint is a valuable unifying picture: the discharge has no mode-creation threshold, only continuous deformation of the dispersion relation with a suppressed, linear-in-density splitting. The derivation is self-contained from the Stix tensor and standard algebra, with no fitted parameters, and the paper supplies explicit falsifiable predictions (linear splitting law, n1 and ncut crossings, the ν_m>ω_ce separation condition, the resonance-coalescence crossing at p*=3.85, and the edge-absorption peak) together with verification residuals against the full quartic. These are genuine strengths. However, two quantitative claims that the paper elevates to headline status—the beta-pinning of the helicon root and the low-density rows of Table 4—do not survive direct evaluation of Eq. (9)-(10), so the paper needs substantial revision before its central narrative can be accepted.","major_comments":[{"comment":"Section III.D (Eqs. (19)-(20)), the abstract, and the Conclusions state that the lower root remains pinned at beta=k0 to within 1% until n_e≃n3. Direct evaluation of Eq. (9) with Eq. (7) in the collisionless reference case contradicts this. At n_e=8.3×10^14 m^-3 the quartic gives k_{⊥,-}^2≈-986.5 m^-2, so beta_-^2=k_{⊥,-}^2+k_z^2≈-0.5 m^-2 (vacuum value k0^2=0.0807 m^-2); the lower root is radially evanescent with imaginary beta≈0.7 m^-1, not beta≈0.284 m^-1. What remains within about 1% of its vacuum value is k_⊥^2 (equivalently the evanescence rate κ), not the total wavenumber beta defined in Section III.D. The 'pinned helicon root' claim should be reworded to refer to k_⊥^2/κ, and the corresponding statements in the abstract, Table 2, Figure 1, and the Conclusions should be corrected.","section":"Section III.D, abstract, Conclusions"},{"comment":"Table 4 lists beta_Helicon=0.0+0.00i, 0.1+0.00i, and 0.5+0.01i at n_e=10^14, 10^15, and 10^16 m^-3. These entries reproduce the EMHD lower root beta_-≈ω μ0 n_e e/(k_z B0) of Eqs. (24)-(25), not the roots of the full quartic Eq. (9). The EMHD root vanishes as n_e→0 (Eq. (27)), whereas the full quartic's lower root tends to beta=k0=0.284 m^-1 in the vacuum limit and is of order 1 m^-1 at n_e=10^14 m^-3 in the collisionless limit. Section III.F explicitly states that the EMHD reduction is invalid for n_e≲n3, so these rows cannot represent the ignition transient. Because Table 4 underpins the stage I-II field profiles in Figure 7, the edge-absorption percentages, and the resolution estimates of Section VIII.B, those low-density results should be recomputed from Eq. (9), or the table should be restricted to the density range where the EMHD reduction has been validated and the restriction stated in the caption.","section":"Table 4, Section VII vs Section III.F"},{"comment":"The statement that the TG root 'ceases to be radially evanescent precisely at n1' (Eq. (19)) and the abstract's 'Two thresholds acquire exact meaning' are collisionless statements. For finite ν_m, P=1-ω_pe^2/[ω(ω+iν_m)] is complex and never vanishes at any real density, so there is no exact real-density crossing; the quoted five-figure agreement must refer to a low-collision evaluation. Section III.E itself notes that ν_m can exceed ω_ce at the time n1 is crossed early in the transient. The threshold discussion should be explicitly restricted to the collisionless (or demagnetisation-free) limit, with a stated collision-dependent replacement for the early phase, for example the density where Re beta_+=k_z; otherwise the 'exact' framing in the abstract and Section III.D overstates the result.","section":"Section III.D and III.E"}],"minor_comments":[{"comment":"The abstract says the helicon root is pinned until about 10^15 m^-3, while Section III.D gives 8.3×10^14 m^-3; these should be aligned after the pinning statement is corrected to refer to k_⊥^2 rather than beta.","section":"Abstract and Table 2"},{"comment":"The final paragraph of Section VII correctly identifies Landau damping as the principal limitation, but the quantitative predictions in Table 4 and prediction 9 (94% edge absorption) are collisional-only results; the text should state explicitly that these are lower-bound estimates for late-transient edge absorption, where ω/(k_z v_te)=3.7 makes kinetic damping non-negligible for the TG branch.","section":"Section VII"},{"comment":"Figure 1 plots Re beta and Im beta, while the text of Section III.E describes the lower root as having |beta|≃k0 at low density; when beta is imaginary, Re beta is zero, so the figure and the text should make clear whether the plotted quantity is |beta|, Re beta, or Im beta for the evanescent lower root.","section":"Figure 1"}],"recommendation":"major_revision","confidential_remarks":"The double-root framework and the linear-splitting result appear sound and are the paper's main contribution. The pinning error and the EMHD-based low-density rows of Table 4 are substantive but local; they can be fixed by reworking the wording and recomputing the affected tables and figures. The 'exact threshold' claims also need explicit collisionless qualification. The manuscript would benefit from a careful pass over the abstract and conclusions to remove overstatements before resubmission."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Lei Chang's paper is worth a serious read, but the headline claim needs surgery. The core structural result is real: the vacuum limit of the cold-plasma quartic is a degenerate double root, and the helicon/TG pair is the splitting of that root by D and P-S, with the first-order discriminant term cancelling so the separation is linear in density. That part checks out, and the paper's careful walk through the discriminant is the best thing in it. The resonance-coalescence ratio R=(4/Lambda)sqrt(1+T^2/kz^2) and the omega_ce>4omega bound are compact and, as far as I know, new. The two-timescale separation of local vs global adiabaticity is well argued, and the swept-resonance Fresnel solution is a nice piece of work.\n\nThe soft spots are real, though. The abstract and Section III.D claim the helicon root stays pinned at beta=k0 to within 1% until ne about 10^15 m^-3, with beta defined as sqrt(k_perp^2+kz^2). That is not what the full quartic gives. At ne=0.1 n1, the lower root already has beta^2 negative, |beta| around 0.7 m^-1 rather than 0.284; by ne=1e15 the lower root has |beta| around 30 m^-1, nowhere near k0. What stays near its vacuum value is k_perp^2 relative to kz^2, not beta. The 'pinned at beta=k0' statement is wrong as written. The paper should either redefine the claim in terms of k_perp^2 or drop the word 'pinned'.\n\nTable 4 has a related problem: the low-density rows (1e14 and 1e15 m^-3) list beta_Helicon = 0.0 and 0.1 m^-1, which match the EMHD small root, not the full quartic. The full quartic at those densities gives an evanescent lower root with |beta| of order tens of m^-1. Since Section III.F correctly warns that EMHD is invalid below n3, the table contradicts the paper's own analysis. The table needs to state which model each row uses, or be recomputed with the full quartic.\n\nThere is also a moderate caveat, already acknowledged in the paper: the exact n1 crossing is collisionless, while the early transient has nu_m > omega_ce, which blurs the threshold. That is a limitation, not a fatal flaw.\n\nBottom line: the double-root framework and the new ratio are solid and deserve referee time. The pinning claim and Table 4 need fixing before publication. I would send it to review with a request for major revision, and ask the author to reconcile Figure 1, Table 4, and the abstract with the full quartic.","headline":"A genuinely new double-root framework for helicon ignition, but the 'pinned helicon root' headline claim and Table 4 are inconsistent with the paper's own full quartic.","tokens_in":24979,"tokens_out":24039,"would_cite":false,"duration_ms":199619,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["52.35.Hr","52.50.Dg"],"model":"deepseek-v4-flash","headline":"The helicon and Trivelpiece–Gould waves are the two lobes of a single degenerate vacuum root, split continuously as density rises.","keywords":["helicon wave","Trivelpiece-Gould wave","dispersion relation","discharge ignition","degenerate double root","coalescence density","cold plasma dielectric tensor","two-timescale adiabaticity"],"falsifier":"A direct test: in a 13.56 MHz argon discharge at about 100 G, measure the two perpendicular wavenumbers while ramping density from $10^8$ to $10^{12}$ m$^{-3}$. The theory predicts the two roots converge to a single value $k_\\perp^2=k_0^2-k_z^2$ as $n_e\\to0$ and that their separation grows linearly with density, a factor of ten per decade; if instead the separation grows as $\\sqrt{n_e}$, or if the roots do not converge to the vacuum evanescent value, the degeneracy-breaking picture fails. A second test: at $B_0<19$ G and 13.56 MHz, the theory says no bounded mode resonates before coalescence, so the loading resistance should show no clean resonance peak; observing one would falsify the $\\omega_{ce}>4\\omega$ bound.","tokens_in":23714,"feed_emoji":"⚡","tokens_out":6346,"duration_ms":59733,"temperature":0.7,"pith_summary":"This paper tries to show that the helicon and Trivelpiece–Gould (TG) waves of a magnetised plasma discharge are not two separate modes that switch on at some density threshold, but the two halves of a single wave solution that already exists in vacuum. The argument starts from the quartic that fixes the perpendicular wavenumber and observes that, as the electron density goes to zero, the quartic degenerates into one repeated root: the evanescent near field of the antenna. Plasma formation then breaks that degeneracy through the gyrotropic and anisotropic pieces of the cold-plasma dielectric tensor, so the two branches separate continuously from the vacuum state. The paper works out the exact densities at which each branch stops being evanescent and the condition under which a bounded mode can resonate before the branches coalesce, and it shows that a helicon source can be followed analytically through the whole ignition transient if it is treated with an instantaneous-tensor fluid model.","feed_headline":"One vacuum wave splits into helicon and TG as plasma forms","feed_subtitle":"A continuous dispersion relation from vacuum to steady discharge, with exact threshold densities and a sharp field bound.","key_machinery":"The central object is the quartic (quadratic in $k_\\perp^2$) obtained from $\\nabla\\times\\nabla\\times\\mathbf{E}=k_0^2 \\boldsymbol{\\varepsilon}\\cdot\\mathbf{E}$ after eliminating the polarization, with coefficients $A=S$, $B=(S+P)(k_z^2-S k_0^2)+D^2 k_0^2$, and $C=P[(k_z^2-S k_0^2)^2-D^2 k_0^4]$. Its vacuum limit is a degenerate double root because $S=P=1$, $D=0$; the discriminant's leading term cancels identically when density is switched on, forcing a linear-in-density branch separation. The mechanism doing the work is the cold, collisional Stix dielectric tensor with elements $S$, $D$, and $P$ that reduces continuously to the identity as $n_e\\to0$, together with the two invariants of the quartic, the sum and product of the $k_\\perp^2$ roots, which track the transient without branch-cut ambiguity.","core_discovery":"At the heart of the paper is a structural statement: in the vacuum limit the quartic $A k_\\perp^4 + B k_\\perp^2 + C=0$ governing the perpendicular wavenumber has a double root $k_\\perp^2 = k_0^2 - k_z^2$, so the discriminant vanishes identically. The helicon and TG branches are therefore the two lobes of one vacuum root, split by the gyrotropic element $D$ and the anisotropy $P-S$ of the dielectric tensor. The leading-order term in the discriminant cancels, so the splitting grows linearly with electron density $n_e$ rather than as $n_e^{1/2}$; the TG root leaves the vacuum value immediately and becomes radially propagating exactly at $n_1 = \\epsilon_0 m_e \\omega^2 / e^2$, while the helicon root remains pinned near its vacuum value until $n_e \\simeq 10^{15}$ m$^{-3}$ and propagates only at $n_{\\rm cut} = (4/\\Lambda) n_c$. The paper also derives the resonance-to-coalescence ratio $R = n_{\\rm res}/n_c = (4/\\Lambda)\\sqrt{1+T^2/k_z^2}$, the sharp bound $\\omega_{ce} > 4\\omega$ for a clean bounded helicon resonance, and a two-timescale result: the local dielectric response is adiabatic to better than $10^{-3}$ throughout, while the driven cavity response is strongly non-adiabatic for the first tens of microseconds.","pith_inferences":["The linear-in-density splitting law gives a direct experimental handle: measuring the two perpendicular wavenumbers between $10^8$ and $10^{12}$ m$^{-3}$ should show a factor-ten separation increase per decade, and a deviation from that law would expose where the cold-fluid description starts to fail.","The resonance-to-coalescence ratio $R$ could serve as an antenna-design criterion: choosing $k_z$ and $B_0$ so that the excited axial modes lie above $k_z^*$ should select clean, well-separated helicon resonances rather than broad coalescence-embedded ones.","Because the magnetisation condition $\\nu_m<\\omega_{ce}$ is crossed by neutral heating and depletion rather than by ionisation, the theory implies that gas temperature and neutral depletion timing control the E–H–W sequence; a testable extension would compare discharges in gases with very different momentum-transfer cross-sections."],"forward_implications":["The helicon and TG branches exist in a well-defined sense at every density, including zero; there is no threshold at which modes are created, only densities at which each branch stops being evanescent.","The TG branch becomes radially propagating at $n_1\\simeq2.28\\times10^{12}$ m$^{-3}$ in the reference case, about five orders of magnitude before the helicon branch leaves its vacuum value and more than two before it propagates at $n_{\\rm cut}=(4/\\Lambda)n_c$.","No bounded axial mode can resonate before helicon–TG coalescence unless $\\omega_{ce}>4\\omega$; at 13.56 MHz this puts clean helicon operation at $B_0$ above about 19 G, independent of density, power, and gas.","The instantaneous dielectric tensor is accurate to better than $10^{-3}$ throughout ignition, whereas the driven cavity response is non-adiabatic during roughly the first 80 $\\mu$s; the observable signatures are a delayed, reduced, ringing resonance peak.","Fractional edge absorption peaks near 94% at $n_e\\simeq3\\times10^{17}$ m$^{-3}$ and falls as collisions decrease, because the Landau–Zener conversion amplitude at coalescence scales as $(\\nu/\\omega)^{1/2}$."],"supporting_citations":[{"why":"Supplies the Stix-form cold collisional dielectric tensor whose vacuum limit is the identity, the constitutive foundation of the quartic.","marker":"[31]"},{"why":"Identifies the helicon–TG coalescence condition that the paper uses for $n_c$ and the resonance-coalescence ratio $R$.","marker":"[12]"},{"why":"Defines the Trivelpiece–Gould space-charge branch in cylindrical columns that the paper tracks through ignition.","marker":"[13]"},{"why":"Establishes the helicon branch relation $\\beta k_z = \\omega\\mu_0 n_e e/B_0$ used for the helicon cutoff density $n_{\\rm cut}$.","marker":"[11]"},{"why":"Provides the two-branch normal-mode formulation that underlies the bounded eigenvalue problem and the $T$ eigenvalue.","marker":"[15]"},{"why":"Supports the claim that TG waves mediate antenna coupling to helicons, motivating the edge-absorption analysis.","marker":"[17]"}],"fun_headline_variants":["Single vacuum root splits into helicon and TG waves at ignition","Helicon and TG branches are lobes of a single vacuum root","Dispersion relation evolves continuously from vacuum to plasma","One wave becomes two: helicon and TG as plasma forms"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The paper's load-bearing assumption is that the cold, collisional, two-species fluid dielectric tensor, evaluated with the instantaneous density and collision frequency and with newly born electrons carrying zero drift velocity, correctly describes the radio-frequency response at every stage from vacuum to steady state.","fun_headline_variants_meta":{"raw":{"variants":["Single vacuum root splits into helicon and TG waves at ignition","Helicon and TG branches are lobes of a single vacuum root","Dispersion relation evolves continuously from vacuum to plasma","One wave becomes two: helicon and TG as plasma forms"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000464,"raw_usage":{"total_tokens":2402,"prompt_tokens":1111,"completion_tokens":1291,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":727,"completion_tokens_details":{"reasoning_tokens":1222}},"tokens_in":727,"tokens_out":1291,"duration_ms":10157,"temperature":1.0,"reasoning_tokens":1222,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T04:30:05.071988+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct test: in a 13.56 MHz argon discharge at about 100 G, measure the two perpendicular wavenumbers while ramping density from $10^8$ to $10^{12}$ m$^{-3}$. The theory predicts the two roots converge to a single value $k_\\perp^2=k_0^2-k_z^2$ as $n_e\\to0$ and that their separation grows linearly with density, a factor of ten per decade; if instead the separation grows as $\\sqrt{n_e}$, or if the roots do not converge to the vacuum evanescent value, the degeneracy-breaking picture fails. A second test: at $B_0<19$ G and 13.56 MHz, the theory says no bounded mode resonates before coalescence, so the loading resistance should show no clean resonance peak; observing one would falsify the $\\omega_{ce}>4\\omega$ bound.","supporting_citations":[{"cited_title":"Takahashi, Thirty percent conversion efficiency from radiofrequency power to thrust energy in a magnetic nozzle plasma thruster, Sci","cited_arxiv_id":null,"evidence_quote":"Supplies the Stix-form cold collisional dielectric tensor whose vacuum limit is the identity, the constitutive foundation of the quartic."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Identifies the helicon–TG coalescence condition that the paper uses for $n_c$ and the resonance-coalescence ratio $R$."},{"cited_title":"Komori, T","cited_arxiv_id":null,"evidence_quote":"Defines the Trivelpiece–Gould space-charge branch in cylindrical columns that the paper tracks through ignition."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the helicon branch relation $\\beta k_z = \\omega\\mu_0 n_e e/B_0$ used for the helicon cutoff density $n_{\\rm cut}$."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the two-branch normal-mode formulation that underlies the bounded eigenvalue problem and the $T$ eigenvalue."},{"cited_title":"Isayama, S","cited_arxiv_id":null,"evidence_quote":"Supports the claim that TG waves mediate antenna coupling to helicons, motivating the edge-absorption analysis."}],"review_version":1}