Transitional evolution of the wave dispersion relation during helicon discharge ignition: an analytical theory from vacuum to steady operation

Lei Chang

REVIEW 3 major objections 3 minor 60 references

Transitional evolution of the wave dispersion relation during helicon discharge ignition: an analytical theory from vacuum to steady operation

T0 review · 3 major / 3 minor · reviewed 2026-08-14 · deepseek-v4-flash

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

Pith's one-line read The helicon and Trivelpiece–Gould waves are the two lobes of a single degenerate vacuum root, split continuously as density rises.

desk verdict 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. read the letter →

arxiv 2608.08649 v1 pith:3UATQ2YA submitted 2026-08-09 physics.plasm-ph

Lei Chang This is my paper · ORCID

classification physics.plasm-ph PACS 52.35.Hr52.50.Dg

keywords heliconwaveTrivelpiece-Goulddispersionrelationdischargeignitiondegeneratedoublerootcoalescencedensitycoldplasmadielectrictensortwo-timescaleadiabaticity

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

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.

What carries the argument

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.

What would settle it

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.

Watch

Extended reading notes

Core claim

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.

Load-bearing premise

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.

Editorial extensions

If this is right

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}$.

Reading between the lines

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

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.

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.

Desk Editor's Note

read the 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.

The 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'.

Table 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.

There 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.

Bottom 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.

Referee Report

3 major / 3 minor

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.

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 (3)

[Section III.D, abstract, Conclusions] 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.
[Table 4, Section VII vs Section III.F] 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 III.D and III.E] 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.

minor comments (3)

[Abstract and Table 2] 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 VII] 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.
[Figure 1] 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.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the dispersion-relation derivation is self-contained from the standard cold-plasma tensor; self-citations are contextual only.

full rationale

The paper's central derivation starts from the standard cold, collisional, two-species Stix tensor (Eqs. (4)-(7)), which is regular in the vacuum limit, and reduces the wave equation to the quartic in Eqs. (9)-(10). The double-root vacuum limit, the cancellation of the leading discriminant term, and the linear-in-density splitting are obtained by direct algebra (Eqs. (12)-(16)), not by assuming the result. The characteristic densities are defined by tensor-element conditions (n1: P=0; n2: |D|=1; n3: S=2), and the claims that the roots cross kz or remain near k0 at those densities are verified against numerical evaluation of the full quartic (Section 8.4), which is a consistency check on the same equations rather than an independent input. The thresholds ncut and nc follow explicitly from the roots and are cross-checked against the quartic. Self-citations ([27], [35]-[37], [43], [44], [49], [50]) appear only for context, comparison, or as examples of numerical/experimental verification; none is used as a premise in the derivation, and no uniqueness theorem or ansatz is imported from the author's prior work. The skeptical concern about Table 4 using the EMHD reduction at low density is an internal consistency/correctness issue about which reduction is valid where, not a circularity: those entries are not fitted inputs renamed as predictions. Therefore the derivation is self-contained and no circular step can be exhibited.

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

No free parameters are fitted to data; reference-case values are fixed inputs. No invented entities are introduced. The central results are explicit functions of B0, omega, kz, a, and the standard cold-plasma dielectric tensor.

assumptions (6)

domain assumption The cold, collisional, two-species Stix dielectric tensor with displacement current and ion inertia retained is valid from vacuum to steady state.
Invoked in Section II.D, Eqs. (6)-(7). All root trajectories follow from this tensor; kinetic effects are acknowledged as the principal limitation in Section VII.
domain assumption Fields are single Fourier modes exp[i(m*theta + kz*z - omega*t)] with a slowly varying envelope; partial_t goes to -i*omega at leading order.
Section II.A and Section V. Excludes broad-spectrum antenna coupling and mode-mode coupling beyond the two branches.
domain assumption Boundary conditions Jr(a)=0 with prescribed Bz(a), and radial eigenvalue T=3.832/a for m=1 with kz much less than T.
Section IV.B, Eqs. (34)-(36). Real antennas excite a kz spectrum, so the sharp resonance ladder is an idealization.
domain assumption The density and collision-frequency evolution are supplied by a global balance and diffusion transport model with tau_n approximately 22 microseconds.
Section VI, Eq. (55), and Figure 4. The non-adiabatic conclusions in Figure 5 depend on this modeled trajectory, not on measured data.
domain assumption Ionization source electrons are born with zero RF drift, so S_iz v_birth = 0 and no explicit dot_n_e/n_e term enters the current equation.
Section II.C, Eqs. (4)-(5). If newly born electrons carry RF momentum, the constitutive-law correction would differ.
standard math The Fresnel-integral solution of a linearly swept resonance and the Landau-Zener two-state formula are valid for the driven and conversion problems.
Section V.D-E, Eqs. (46)-(51). Standard results are applied without derivation.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Transitional evolution of the wave dispersion relation during helicon discharge ignition: an analytical theory from vacuum to steady operation." pith.science (2026). https://pith.science/paper/3UATQ2YA

@misc{pith2026260808649,
  author       = {Pith},
  title        = {Pith review of: Transitional evolution of the wave dispersion relation during helicon discharge ignition: an analytical theory from vacuum to steady operation},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/3UATQ2YA}},
  note         = {Machine review of arXiv:2608.08649}
}

abstract

The dispersion relations of the helicon and Trivelpiece--Gould (TG) waves in a steady, fully developed magnetised plasma column are well established, but the discharge does not begin in that state. During ignition the electron density rises through five to seven orders of magnitude and the electron--neutral collision frequency falls by two to three, so the dispersion relation is itself a time-dependent object whose very existence must be justified. This paper develops a closed analytical theory of that transition for a cylindrical column in a uniform axial field, using a cold, collisional, two-species description that remains regular in the vacuum limit. Three structural results emerge. First, the vacuum limit is a degenerate double root of the quartic governing the perpendicular wavenumber; the helicon and TG branches are not independent modes that appear at some threshold but the two lobes of a single vacuum root, split by the gyrotropic and anisotropic parts of the dielectric tensor. The leading-order contribution to the discriminant cancels identically, so the branches separate linearly in density rather than with the square-root behaviour generic to a perturbed repeated root; the separation is correspondingly suppressed, and the helicon root remains pinned at its vacuum value until the electron density reaches about 10$^{15}$ m$^{-3}$ while the TG root departs immediately.

Figures

Figures reproduced from arXiv: 2608.08649 by the authors.

**Figure 2.** FIG. 2. Regime map in the [PITH_FULL_IMAGE:figures/full_fig_p015_2.png] view at source ↗

**Figure 3.** FIG. 3. The ladder of bounded-mode resonant densities [PITH_FULL_IMAGE:figures/full_fig_p018_3.png] view at source ↗

Figures from the paper (4 more)

**Figure 4.** Figure 4: FIG. 4. Left: modelled density evolution [PITH_FULL_IMAGE:figures/full_fig_p019_4.png]

**Figure 5.** Figure 5: FIG. 5. The two adiabaticity parameters through the transient. The local parameter [PITH_FULL_IMAGE:figures/full_fig_p021_5.png]

**Figure 6.** Figure 6: FIG. 6. Exact solution of the swept-resonance problem for four values of the sweep parameter [PITH_FULL_IMAGE:figures/full_fig_p023_6.png]

**Figure 7.** Figure 7: FIG. 7. Two-branch driven eigenfunction (top) and power deposition profile (bottom) at four stages [PITH_FULL_IMAGE:figures/full_fig_p028_7.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

60 extracted references · 60 canonical work pages

[1]

The EMHD formulation therefore cannot be used to initialise a transient calculation; the full quartic, with displacement current and ion inertia retained, is required forne ≲n 3

The reason is that discarding the displacement current removes precisely the term that survives when the plasma current van- ishes. The EMHD formulation therefore cannot be used to initialise a transient calculation; the full quartic, with displacement current and ion inertia retained, is required forne ≲n 3. 13 G. Coalescence and the critical density The...

work page
[2]

Asn e →0the two computed roots converge to the single valuek 2 ⊥ =k 2 0 −k 2 z, and the splitting grows linearly inne, by a factor of ten for every tenfold rise in density. The linear law is asymptotic and should be tested over108 ≲n e ≲10 12 m−3, where it holds to better than 1%; byne = 1014 m−3 the decade ratio has fallen to 9.2 asn3 is approached and t...

work page
[3]

The upper root crosses|β|=k z atn 1 and the lower root atncut = (4/Λ)nc, the latter remaining within 1% ofk0 untiln e ≃n 3

work page
[4]

Ifν m is held aboveωce, no helicon–TG separation occurs at any density

work page
[5]

5.Rcrosses unity betweenp= 3andp= 4at 100 G, and no bounded mode resonates before coalescence ifB0 <19G

The coalescence point scales asB2 0k2 z; repeating a sweep at 300 G moves the branch exchange up by a factor of nine. 5.Rcrosses unity betweenp= 3andp= 4at 100 G, and no bounded mode resonates before coalescence ifB0 <19G

work page
[6]

A time-domain constitutive model and an instantaneous-tensor model differ by less than10 −3

work page
[7]

Quasi-static and transient calculations diverge exactly whereϵglob >1, that is during the first∼80µs

work page
[8]

The transient loading peak occurs at higher density than the static prediction and is followed by a decaying oscillation of rising frequency

work page

Show all 60 references

[9]

The fraction of power deposited outsider/a= 0.7peaks nearn e ∼3×10 17 m−3 at about 94% and falls thereafter
[10]

electron

A slow power ramp couples more energy into the helicon mode than a step of equal total energy, by approximately the overlap-integral factor of Section 5.6. 31 Predictions 3, 5, 8 and 10 are accessible to experiment through the loading resistanceRp(t) alone; the remainder requi...
[11]

R. W. Boswell, Plasma production using a standing helicon wave, Phys. Lett. A 33, 457 (1970)

1970
[12]

R. W. Boswell, Very efficient plasma generation by whistler waves near the lower hybrid fre- quency, Plasma Phys. Control. Fusion 26, 1147 (1984)

1984
[13]

Komori, T

A. Komori, T. Shoji, K. Miyamoto, J. Kawai, and Y. Kawai, Helicon waves and efficient plasma production, Phys. Fluids B 3, 893 (1991)

1991
[14]

R. W. Boswell and F. F. Chen, Helicons—the early years, IEEE Trans. Plasma Sci. 25, 1229 (1997)

1997
[15]

F. F. Chen and R. W. Boswell, Helicons—the past decade, IEEE Trans. Plasma Sci. 25, 1245 (1997)

1997
[16]

F. F. Chen, Helicon discharges and sources: A review, Plasma Sources Sci. Technol. 24, 014001 (2015)

2015
[17]

Isayama, S

S. Isayama, S. Shinohara, and T. Hada, Review of helicon high-density plasma: Production mechanism and plasma/wave characteristics, Plasma Fusion Res. 13, 1101014 (2018)

2018
[18]

Shinohara, Helicon high-density plasma sources: Physics and applications, Adv

S. Shinohara, Helicon high-density plasma sources: Physics and applications, Adv. Phys. X 3, 1420424 (2018)

2018
[19]

Chang, R

L. Chang, R. Boswell, E. Scime, S. Shinohara, K. Takahashi, S. Thakur, F. Filleul, A. Cal- darelli, S. Isayama, Y. Yu, M. Xu, H. B. Zhang, T. Y. Huang, M. Y. Wu, J. T. Wu, B. H. Xia, L. F. Lu, A. P. Sun, D. Du, Z. Y. Zhang, R. X. Yuan, A. D. Xu, X. Yang, D. Jing, Y. Xia, C. Wa...

2024
[20]

Chang, R

L. Chang, R. Boswell, and G. N. Luo, First Helicon Plasma Physics and Applications Work- shop, Front. Phys. 9, 808971 (2022)

2022
[21]

F. F. Chen, Plasma ionization by helicon waves, Plasma Phys. Control. Fusion 33, 339 (1991)

1991
[22]

K. P. Shamrai and V. B. Taranov, Volume and surface rf power absorption in a helicon plasma source, Plasma Sources Sci. Technol. 5, 474 (1996)

1996
[23]

A. W. Trivelpiece and R. W. Gould, Space charge waves in cylindrical plasma columns, J. Appl. Phys. 30, 1784 (1959)

1959
[24]

Li and G

W. Li and G. Wang, Damping characteristics of helicon and Trivelpiece–Gould waves in high density and low magnetic field helicon plasma, Phys. Plasmas 30, 023510 (2023)

2023
[25]

F. F. Chen and D. Arnush, Generalized theory of helicon waves. I. Normal modes, Phys. Plasmas 4, 3411 (1997)

1997
[26]

Arnush and F

D. Arnush and F. F. Chen, Generalized theory of helicon waves. II. Excitation and absorption, Phys. Plasmas 5, 1239 (1998)

1998
[27]

Arnush, The role of Trivelpiece–Gould waves in antenna coupling to helicon waves, Phys

D. Arnush, The role of Trivelpiece–Gould waves in antenna coupling to helicon waves, Phys. Plasmas 7, 3042 (2000)

2000
[28]

M. A. Lieberman and A. J. Lichtenberg, Principles of Plasma Discharges and Materials Pro- cessing, 2nd ed. (Wiley, New York, 2005)

2005
[29]

Charles, Plasmas for spacecraft propulsion, J

C. Charles, Plasmas for spacecraft propulsion, J. Phys. D: Appl. Phys. 42, 163001 (2009)

2009
[30]

Takahashi, Helicon-type radiofrequency plasma thrusters and magnetic plasma nozzles, Rev

K. Takahashi, Helicon-type radiofrequency plasma thrusters and magnetic plasma nozzles, Rev. Mod. Plasma Phys. 3, 3 (2019)

2019
[31]

Takahashi, Thirty percent conversion efficiency from radiofrequency power to thrust energy in a magnetic nozzle plasma thruster, Sci

K. Takahashi, Thirty percent conversion efficiency from radiofrequency power to thrust energy in a magnetic nozzle plasma thruster, Sci. Rep. 12, 18618 (2022)

2022
[32]

R. I. Pinsker, R. Prater, C. P. Moeller, J. S. deGrassie, C. C. Petty, M. Porkolab, J. P. Anderson, A. M. Garofalo, C. Lau, A. Nagy, D. C. Pace, H. Torreblanca, J. G. Watkins, and L. Zeng, Experiments on helicons in DIII-D: Investigation of the physics of a reactor-relevant no...

2018
[33]

Van Compernolle, M

B. Van Compernolle, M. W. Brookman, C. P. Moeller, R. I. Pinsker, A. M. Garofalo, R. O’Neill, D. Geng, A. Nagy, J. P. Squire, K. Schultz, C. Pawley, D. Ponce, A. C. Torrezan, J. Lohr, B. Coriton, E. Hinson, R. Kalling, A. Marinoni, E. H. Martin, R. Nguyen, C. C. Petty, M. Pork...

2021
[34]

A. W. Degeling, C. O. Jung, R. W. Boswell, and A. R. Ellingboe, Plasma production from helicon waves, Phys. Plasmas 3, 2788 (1996)

1996
[35]

C.M.Franck, O.Grulke, andT.Klinger, Modetransitionsinhelicondischarges, Phys.Plasmas 10, 323 (2003)

2003
[36]

B. N. Breizman and A. V. Arefiev, Radially localized helicon modes in nonuniform plasma, Phys. Rev. Lett. 84, 3863 (2000)

2000
[37]

J.J.Ma, L.Chang, M.Y.Wu, H.Zhou, Y.W.Zhang, I.I.Zadiriev, E.A.Kralkina, S.Isayama, and S. J. You, Revealing the transient ionization dynamics and mode-coupling mechanisms in helicon plasmas through a self-consistent multiphysics model, Sci. Rep. 16, 16529 (2026)

2026
[38]

Y. J. Chang, Z. Zhou, and G.-N. Luo, Numerical study on the temporal evolution of a helicon discharge, IEEE Trans. Plasma Sci. 49, 3733 (2021)

2021
[39]

J. T. Wu, C. F. Sang, C. J. Sun, B. Sun, S. S. Gao, Y. Bian, Y. Peng, Q. Wang, and D. Z. Wang, Experimental and simulation study of argon helicon discharge in multiple plasma simulation linear device (MPS-LD), Plasma Sources Sci. Technol. 33, 085007 (2024)

2024
[40]

Mouzouris and J

Y. Mouzouris and J. E. Scharer, Wave propagation and absorption simulations for helicon sources, Phys. Plasmas 5, 4253 (1998)

1998
[41]

T. H. Stix, Waves in Plasmas (AIP, New York, 1992)

1992
[42]

D. G. Swanson, Plasma Waves, 2nd ed. (Institute of Physics, Bristol, 2003)

2003
[43]

R. W. Boswell and D. Vender, An experimental study of breakdown in a pulsed helicon plasma, Plasma Sources Sci. Technol. 4, 534 (1995)

1995
[44]

G. D. Conway, A. J. Perry, and R. W. Boswell, Evolution of ion and electron energy distribu- tions in pulsed helicon plasma discharges, Plasma Sources Sci. Technol. 7, 337 (1998)

1998
[45]

Chang, B

L. Chang, B. N. Breizman, and M. J. Hole, Gap eigenmode of radially localized helicon waves in a periodic structure, Plasma Phys. Control. Fusion 55, 025003 (2013)

2013
[46]

Y. W. Zhang, L. Chang, and G. S. Xu, Wave penetration and power deposition of helicon current drive for magnetic fusion plasma with sharp edge gradient, Plasma Phys. Control. Fusion 67, 065022 (2025)

2025
[47]

Chang, M

L. Chang, M. J. Hole, J. F. Caneses, G. Chen, B. D. Blackwell, and C. S. Corr, Wave modeling in a cylindrical non-uniform helicon discharge, Phys. Plasmas 19, 083511 (2012)

2012
[48]

L. D. Landau, Zur Theorie der Energieübertragung II, Phys. Z. Sowjetunion 2, 46 (1932). 40

1932
[49]

Zener, Non-adiabatic crossing of energy levels, Proc

C. Zener, Non-adiabatic crossing of energy levels, Proc. R. Soc. London A 137, 696 (1932)

1932
[50]

G. B. Whitham, Linear and Nonlinear Waves (Wiley, New York, 1974)

1974
[51]

M. A. Lieberman and S. Ashida, Global models of pulse-power-modulated high-density, low- pressure discharges, Plasma Sources Sci. Technol. 5, 145 (1996)

1996
[52]

Ashida, C

S. Ashida, C. Lee, and M. A. Lieberman, Spatially averaged (global) model of time modulated high density argon plasmas, J. Vac. Sci. Technol. A 13, 2498 (1995)

1995
[53]

Chang, S

L. Chang, S. J. Zhang, J. T. Wu, Y. W. Zhang, C. Wang, Y. Peng, S. S. Gao, C. J. Sun, Q. Wang, C. F. Sang, S. C. Thakur, S. Isayama, and S. J. You, Transition of blue-core mode helicon discharge, Phys. Plasmas 32, 123504 (2025)

2025
[54]

C. Wang, L. Chang, L. F. Lu, S. Shinohara, Z. D. Zeng, I. Zadiriev, E. Kralkina, Z. Li, S. J. Zhang, Z. C. Kan, Y. Tao, and D. Z. Li, Spatial and temporal evolutions of blue-core helicon discharge driven by a planar antenna with concentric rings, Phys. Plasmas 32, 113509 (2025)

2025
[55]

Clarenbach, B

B. Clarenbach, B. Lorenz, M. Krämer, and N. Sadeghi, Time-dependent gas density and temperature measurements in pulsed helicon discharges in argon, Plasma Sources Sci. Technol. 12, 345 (2003)

2003
[56]

Cho, A self-consistent global model of neutral gas depletion in pulsed helicon plasmas, Phys

S. Cho, A self-consistent global model of neutral gas depletion in pulsed helicon plasmas, Phys. Plasmas 6, 359 (1999)

1999
[57]

Takahashi, Y

K. Takahashi, Y. Takao, and A. Ando, Neutral-depletion-induced axially asymmetric density in a helicon source and imparted thrust, Appl. Phys. Lett. 108, 074103 (2016)

2016
[58]

Isayama, S

S. Isayama, S. Shinohara, and T. Hada, Spatio-temporal behavior of density jumps and the effect of neutral depletion in high-density helicon plasma, Phys. Plasmas 26, 053513 (2019)

2019
[59]

Chang, J

L. Chang, J. F. Caneses, and S. C. Thakur, Wave propagation and power deposition in blue- core helicon plasma, Front. Phys. 10, 1009563 (2022)

2022
[60]

H. Zhou, L. Chang, G. S. Xu, Y. W. Zhang, M. Hole, D. Du, Z. S. Qu, and M. Q. Wu, Three- dimensional multiphysics modelling of helicon wave heating and antenna–plasma coupling for boundary density control in toroidal fusion plasmas, Nucl. Fusion 66, 076008 (2026). 41

2026

Pith tools

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

[1] [1]

The EMHD formulation therefore cannot be used to initialise a transient calculation; the full quartic, with displacement current and ion inertia retained, is required forne ≲n 3

The reason is that discarding the displacement current removes precisely the term that survives when the plasma current van- ishes. The EMHD formulation therefore cannot be used to initialise a transient calculation; the full quartic, with displacement current and ion inertia retained, is required forne ≲n 3. 13 G. Coalescence and the critical density The...

work page

[2] [2]

Asn e →0the two computed roots converge to the single valuek 2 ⊥ =k 2 0 −k 2 z, and the splitting grows linearly inne, by a factor of ten for every tenfold rise in density. The linear law is asymptotic and should be tested over108 ≲n e ≲10 12 m−3, where it holds to better than 1%; byne = 1014 m−3 the decade ratio has fallen to 9.2 asn3 is approached and t...

work page

[3] [3]

The upper root crosses|β|=k z atn 1 and the lower root atncut = (4/Λ)nc, the latter remaining within 1% ofk0 untiln e ≃n 3

work page

[4] [4]

Ifν m is held aboveωce, no helicon–TG separation occurs at any density

work page

[5] [5]

5.Rcrosses unity betweenp= 3andp= 4at 100 G, and no bounded mode resonates before coalescence ifB0 <19G

The coalescence point scales asB2 0k2 z; repeating a sweep at 300 G moves the branch exchange up by a factor of nine. 5.Rcrosses unity betweenp= 3andp= 4at 100 G, and no bounded mode resonates before coalescence ifB0 <19G

work page

[6] [6]

A time-domain constitutive model and an instantaneous-tensor model differ by less than10 −3

work page

[7] [7]

Quasi-static and transient calculations diverge exactly whereϵglob >1, that is during the first∼80µs

work page

[8] [8]

The transient loading peak occurs at higher density than the static prediction and is followed by a decaying oscillation of rising frequency

work page