Magnetic flux ropes have been successfully created with plasma guns in the newly commissioned PHAse Space MApping (PHASMA) experiment. The flux ropes exhibit the expected m = 1 kink instability. The observed threshold current for the onset of this kink instability is half of the Kruskal–Shafranov current limit, consistent with predictions for the non-line tied boundary condition of PHASMA. The helicity, paramagnetism, and growth rate of the observed magnetic fluctuations are also consistent with kink instability predictions. The observed fluctuation frequency appears to be a superposition of a real frequency due to a Doppler shift of the kink mode arising from plasma flow ( kHz) and a contribution from a wave mode ( kHz). The dispersion of the wave mode is consistent with an Alfvén wave. Distinct from most previous laboratory studies of flux ropes, the working gas in PHASMA is argon. Thus, the ion cyclotron frequency in PHASMA is quite low and the frequency of the Alfvénic mode plateaus at of the ion gyro frequency with increasing background magnetic field strength.
I. INTRODUCTION
A magnetic flux rope is a current-carrying plasma with an embedded helical magnetic field. The field is composed of a background axial magnetic field Bz and an azimuthal magnetic field , generated by an axial plasma current.1 When the plasma current exceeds a threshold value, the m = 1 kink instability appears. The kink instability arises from the imbalance in magnetic pressure forces between inner bunched and outer spread magnetic field lines during transverse perturbations of the flux rope column. Sakurai2 first identified kinked structures on the Sun during remote imaging of solar flare events.3,4 The kink instability drives a variety of explosive solar phenomena, e.g., corona mass injections and magnetic reconnection.5,6 Kinked structures are also observed in magnetized jets emanating from astrophysical sources.7 The kink instability plays important roles in a variety of laboratory plasma devices. In tokamaks, kink instabilities distort the plasma torus, drive sawtooth oscillations, and also drive other instabilities.8 Kink instabilities also facilitate the conversion of toroidal magnetic flux into poloidal magnetic flux during the formation process of spheromak discharges.9
Because of widespread interest, a number of laboratory experiments have focused on the study of kink instabilities. A sampling of kink instability experiments in linear devices is listed in Table I. To excite the kink instability, the flux rope plasma current, Ibias, should exceed a threshold value, , where IKS is the Kruskal–Shafranov limit and L0 is the flux rope length.10 The coefficient α = 1 for line tied (LT) boundary conditions and for non-line tied (NLT) boundary conditions.11,12 At a (non) line tied boundary in a linear experiment, the end of the flux rope is (free) fixed to a (resistive) conducting surface. Kink instabilities with LT boundary conditions were investigated in the Resistive Wall Mode (ReMW) and Rotating Wall Machine (RWM) devices at the University of Wisconsin-Madison. Internal kink modes and the stabilizing effects of different vessel walls were identified.13,14 At Caltech, the growth rate of kink instabilities under LT boundary conditions and multiscale cascades arising from the kink instability were investigated.9,15,16 Kink instabilities with NLT boundary conditions were studied on the Reconnection Scaling eXperiment (RSX) device at Los Alamos National Laboratory17–19 and Large Plasma Device (LAPD) at the University of California, Los Angeles.20 In these two very different NLT kink studies, the threshold current was experimentally verified to be one-half of IKS. In both NLT studies, the observed frequency of the fluctuations was attributed to a Doppler frequency shift, fD, of a zero frequency kink mode (in the plasma frame) due to either axial flow, fVZ, or E × B drift, .
List of fundamental plasma parameters regarding kink instabilities on different devices.
. | PHASMA . | RSX21 . | ReWM13 . | RWM14 . | LAPD20 . | Caltech9 . |
---|---|---|---|---|---|---|
Axial boundary | NLT | NLT | LT | LT | NLT | LT |
Gas species | Argon | Hydrogen | Hydrogen | Hydrogen | Helium | Hydrogen |
Electron density ne ( ) | 5 | 0.9 | 4 | 25 | 0.1 | 10 |
Electron temperature Te (eV) | 5 | 12 | 2 | 10 | 4 | 12 |
Magnetic field Bz (Gauss) | 375 | 120 | 270 | 600 | 660 | 800 |
Length L0(m) | 1–1.7 | 0.92 | 1.2 | 1.2 | 11 | 0.25 |
Radius R0 (cm) | 2 | 2 | 6 | 9 | 2.5 | 5 |
Bias current Ibias (A) | 500 | 320 | 5000 | 5000 | 75 | 105 |
Kink frequency f (kHz) | 5–10 | 25 | 25 | 20 | 5 | >100 |
Doppler shift fD (kHz) | fVZ, 1–3 | fVZ, 20–40 | fVZ, | ∼ | ||
0.7 | 0.1 | 0.06 | 0.02 | 0.02 | 0.1 |
. | PHASMA . | RSX21 . | ReWM13 . | RWM14 . | LAPD20 . | Caltech9 . |
---|---|---|---|---|---|---|
Axial boundary | NLT | NLT | LT | LT | NLT | LT |
Gas species | Argon | Hydrogen | Hydrogen | Hydrogen | Helium | Hydrogen |
Electron density ne ( ) | 5 | 0.9 | 4 | 25 | 0.1 | 10 |
Electron temperature Te (eV) | 5 | 12 | 2 | 10 | 4 | 12 |
Magnetic field Bz (Gauss) | 375 | 120 | 270 | 600 | 660 | 800 |
Length L0(m) | 1–1.7 | 0.92 | 1.2 | 1.2 | 11 | 0.25 |
Radius R0 (cm) | 2 | 2 | 6 | 9 | 2.5 | 5 |
Bias current Ibias (A) | 500 | 320 | 5000 | 5000 | 75 | 105 |
Kink frequency f (kHz) | 5–10 | 25 | 25 | 20 | 5 | >100 |
Doppler shift fD (kHz) | fVZ, 1–3 | fVZ, 20–40 | fVZ, | ∼ | ||
0.7 | 0.1 | 0.06 | 0.02 | 0.02 | 0.1 |
Here, we report the appearance of a wave mode with Alfvénic features in a kinking flux rope in the new PHASMA (PHAse Space MApping) experiment at West Virginia University. In PHASMA, the observed frequency of the fluctuations associated with the kink instability is much larger than any plausible Doppler frequency shift and the characteristics of the differential frequency, , imply the presence of a distinct wave mode. Because the working gas in PHASMA is argon (to facilitate laser-induced fluorescence measurements of ions22,23), the observed fluctuation frequency is a significant fraction of the ion gyro frequency fci. Since the kink instability can be thought of as one part of the spectrum of surface Alfvén waves,21 coupling between low-frequency plasma modes and the kink instability is not entirely unexpected. Kink instabilities have been observed to provide free energy to excite other wave modes in plasmas. For example, the kink instability was recently found to excite whistler modes in the Caltech jet experiment.24
This paper is organized as follows: The PHASMA device, the experimental apparatus used to create flux ropes, and the fluctuation diagnostics are described in Sec. II. The observed fluctuations are identified as kink instabilities in Sec. III, including measurements of the instability threshold in terms of the safety factor , the right (left)-handed helicity for ( ), and the paramagnetic nature of the fluctuations. As noted previously, the frequency of the fluctuations is much larger ( ) than the Doppler frequency shift of a kink due to axial flow as described by Ryutov et al.12 In Sec. IV, we suggest that a plasma mode with a frequency equal to the differential frequency co-exists with kink instability. The Alfvénic nature of the plasma mode is verified by comparison with expectations for Alfvén wave dispersion. A brief summary of the results is given in Sec. V.
II. EXPERIMENTAL APPARATUS
A. PHASMA Device
Construction of the new PHASMA plasma experiment at West Virginia University was completed in 2019. PHASMA is designed to study space plasma-relevant phenomena, including particle heating and acceleration at kinetic scales during magnetic reconnection,25 electromagnetic instabilities driven by ion temperature anisotropy and plasma pressure,26 ion acceleration in expanding plasmas,27 and cross field ion flows near plasma-material interfaces immersed in a magnetized, high-density plasma.28 A key feature of PHASMA is the availability of volumetric, non-perturbative, laser diagnostics for ion and electron velocity distribution function measurements with spatial resolution at the kinetic scale (∼ mm for the electron inertial length). State-of-the-art laser-induced fluorescence (LIF) schemes are available to measure velocity distribution functions of ions22,23 and neutrals.29–31 An incoherent Thomson scattering system32 provides measurements of the electron velocity distribution function.
As shown in Fig. 1, PHASMA consists of two main sections, a helicon plasma source housed on a vacuum chamber 1.7 m in length and 15 cm in diameter and another vacuum vessel 2.7 m in length and 0.4 m in diameter that contains two plasma gun sources and a movable anode. The total three turbomolecular pumps located at the ends of the system provide a total pumping speed of 3800l/s that maintain a base pressure of Torr. A static axial magnetic field up to 2200 Gauss in the helicon plasma source (375 Gauss for plasma gun) is generated with 22 electromagnet coils. Various magnetic field configurations and uniform, flared, and magnetic mirror are created by varying the relative currents in the coils. The helicon source operates at an antenna frequency of 9–15 MHz and an RF power up to 2 kW to create steady-state and pulsed plasmas in argon, helium, and xenon.23,33,34 The plasma guns form pulsed high-density argon, helium, and hydrogen plasmas with a repetition rate of 1 shot per minute. The two plasma sources are fully independent and can operate simultaneously. The use of a background helicon generated plasma during plasma gun discharges provides a pre-ionized, target plasma (for manipulation of the neutral density), introduces flows in the interaction region between the plasma gun discharges, and controls density gradients across the interaction region. In addition, background helicon source plasmas provide both steady-state and high repetition rate (> 1 Hz) target plasmas for testing of laser diagnostics. PHASMA accesses a wide range of magnetized plasma operational regimes, e.g., plasma beta and Lundquist number , by independently controlling the plasma density.
(a) Solidworks rendering of the PHASMA facility highlighting the helicon source and plasma gun sections. Three turbomolecular pumps maintain a base pressure of Torr. Twenty-two solenoid coils are used to generate the background axial magnetic field. (b) Photo of PHASMA with the helicon plasma source operating. (c) Photo of an argon plasma gun discharge.
(a) Solidworks rendering of the PHASMA facility highlighting the helicon source and plasma gun sections. Three turbomolecular pumps maintain a base pressure of Torr. Twenty-two solenoid coils are used to generate the background axial magnetic field. (b) Photo of PHASMA with the helicon plasma source operating. (c) Photo of an argon plasma gun discharge.
B. Flux rope formation
Figure 2 shows the experimental geometry for the flux rope experiments. Horizontal and Vertical plasma guns are inserted along the x and y directions at the z = 0 plane to launch plasma columns along . In these single flux rope studies, only the vertical plasma gun along the y direction is used. The plasma column center is adjustable over 6 cm radially and ±6° azimuthally. A conical anode (bias plate) is placed at the end of PHASMA to drive axial current along the initial seed plasma columns, thereby forming flux ropes. The conical anode has a diameter of 250 mm and a half angle 60°. The magnetic field penetration time of the anode is approximately s. The ratio between the Alfvén transit time τA and inductive decay time is denoted by , confirming that non-line tied boundary condition models are appropriate for these experiments.12,21 The conical anode is designed to be translated axially from m to –1.73 m.
Schematic of the experimental geometry used for these experiments. The vertical plasma gun and bias plate (anode) were used in these single flux rope studies. Four magnetic probes A, B, C, and D are placed equidistant azimuthally, while a fifth probe, E, is placed on a different axial plane at the same azimuthal location as probe B. When both the background magnetic field, Bz, and the bias current density, Jz, are along , fluctuations propagate azimuthally and axially.
Schematic of the experimental geometry used for these experiments. The vertical plasma gun and bias plate (anode) were used in these single flux rope studies. Four magnetic probes A, B, C, and D are placed equidistant azimuthally, while a fifth probe, E, is placed on a different axial plane at the same azimuthal location as probe B. When both the background magnetic field, Bz, and the bias current density, Jz, are along , fluctuations propagate azimuthally and axially.
To characterize the flux rope, four magnetic probes A-D are placed at m, rad and a fifth magnetic probe E is placed at m, . The propagation features of the fluctuations are determined from the phase differences between each probe. All magnetic signals are integrated in real time with active integrator circuits that have a time constant of 10 s and sampled at 5 MSamples/s. The plasma density and electron temperature are measured with triple Langmuir probes located at the z = 0 plane.
Figure 3 presents the typical temporal evolution of a single flux rope. Neutral gas injection starts at ms and ends at ms when an arc plasma is generated with a discharge current of approximately 700 A. At t = 0 ms, a bias current Ibias of 500 A is fired during the plateau of the arc plasma. A flux rope with both axial and azimuthal magnetic fields lasts for more than 10 ms, much longer than the axial Alfvén time s. The azimuthal magnetic field at r = 2 cm is roughly 15 Gauss. Large fluctuations appear before the bias current reaches the peak value and persist for the remainder of the discharge. The plasma density and electron temperature at r = 0 cm also increase due to the additional Ohmic heating and ionization when the bias current is triggered.
Typical temporal evolution of a single argon plasma gun discharge in PHASMA. Solid blue traces are with bias current, and black dashed lines are without bias current. (a) Arc current; (b) bias current; (c) azimuthal magnetic field at r = 2 cm; (d) plasma density at r = 0 cm; (e) electron temperature at r = 0 cm.
Typical temporal evolution of a single argon plasma gun discharge in PHASMA. Solid blue traces are with bias current, and black dashed lines are without bias current. (a) Arc current; (b) bias current; (c) azimuthal magnetic field at r = 2 cm; (d) plasma density at r = 0 cm; (e) electron temperature at r = 0 cm.
III. IDENTIFICATION OF KINK INSTABILITY
A. Safety factor, q, profile
The measured magnetic field and pressure profiles are remarkably consistent with the classic Bennett profile for a z-pinch in magnetohydrodynamic (MHD) equilibrium; similar results were reported for flux ropes in the RSX device.21 The azimuthal magnetic field profile is plotted in Fig. 4(a) and is well fit with a Bennett profile of with a characteristic radius of 1.7 cm. The axial current density derived from the measured azimuthal magnetic field is [see Fig. 4(b)]. Figures 4(c) and 4(d) show the measured electron density and electron temperature profile for the flux rope. The resultant pressure profile is shown in Fig. 4(e), where and kB is the Boltzmann constant. The ion pressure is ignored assuming an ion temperature much smaller than the electron temperature. From fits of the Jz and P profiles, we determine the radius of the flux rope to be 2 cm; note that the black dotted lines in the Jz and P profile plots represent the background levels for both quantities. The safety factor profile, , is plotted with red squares in Fig. 4(f). Extending to the edge of the current channel, a radius of 2 cm, q < 2 is satisfied. Therefore, the entire non-line tied flux rope is vulnerable to the kink instability.12 Indeed, the magnetic fluctuations observed in PHASMA are distributed globally [see the open black circles in Fig. 4(f)]. Globally distributed fluctuations with a clear internal m = 1 mode structure were also observed in ReWM.14 If q < 2 were to be satisfied in a restricted region of the flux rope, e.g., , where R0 is the radius of flux rope, an internal kink mode would grow without significant perturbation of the plasma edge.13 Here, the entire flux rope is unstable. The kink mode is, therefore, an external kink mode, and significant perturbations at R0 are observed.
Radial profile of a flux rope for Bz = 375 Gauss and Ibias = 500 A. (a) Azimuthal magnetic field (black circles) and a fit assuming a Bennett profile (dashed line); (b) Axial current density, Jz; (c) plasma density, ne; (d) electron temperature, Te; (e) Plasma pressure ; (f) Black circles are the amplitude of fluctuations, and the red circles are the safety factor, q. Kink instability threshold q = 2 is marked by the red dashed line.
Radial profile of a flux rope for Bz = 375 Gauss and Ibias = 500 A. (a) Azimuthal magnetic field (black circles) and a fit assuming a Bennett profile (dashed line); (b) Axial current density, Jz; (c) plasma density, ne; (d) electron temperature, Te; (e) Plasma pressure ; (f) Black circles are the amplitude of fluctuations, and the red circles are the safety factor, q. Kink instability threshold q = 2 is marked by the red dashed line.
The global nature of this kink instability is also evident in high-time resolution visible light images of the flux rope. Images obtained using a fast camera at a sample rate of 75 000 frames/s are shown in Fig. 5. Figure 5(a) shows the temporal evolution of a slice of the plasma along the y axis at m. Black and blue dashed lines show the amplitude of the image fluctuations at the center and edge of the flux rope, respectively. The frequency of the fluctuations is consistent with those observed in . The plasma edge is clearly undergoing significant perturbation. Comparing the displacement of the flux rope edge at m [Fig. 5(b)] with that at , we observe a slight increase in the displacement closer to the bias plate. The increase in perturbation amplitude closer to the bias plate is consistent with a non-line tied boundary condition at the bias plate.
Contour plots of visible emission light intensity during a flux rope discharge recorded at m (a) and m (b) with a fast camera operating at 75 000 frames/s. The black lines show the amplitude of the fluctuations at the center Y = 0, while the blue lines show the amplitude of the fluctuations in the outer region of the flux rope.
Contour plots of visible emission light intensity during a flux rope discharge recorded at m (a) and m (b) with a fast camera operating at 75 000 frames/s. The black lines show the amplitude of the fluctuations at the center Y = 0, while the blue lines show the amplitude of the fluctuations in the outer region of the flux rope.
B. Helicity and paramagnetism
For a flux rope with both and Bz, the equilibrium magnetic field has a specific helical pattern, right (left) handedness for ( ). Fluctuations with a helical structure that matches the equilibrium magnetic field structure experience less twisting than those with the opposite helicity, i.e., less stabilizing magnetic tension force. Therefore, the helical perturbations that match the equilibrium structure are more vulnerable to the kink instability.10 In other words, kink instabilities with should dominate.9 For the experimental configuration shown in Fig. 2, , kink fluctuations with right-handed helicity, , should dominate.21
Figure 6(a) shows the temporal evolution of fluctuations obtained with the azimuthally distributed magnetic probes (A–D) at the same radial location. A coherent large-amplitude mode is clearly propagating in the azimuthal direction. From the measured phase difference, the azimuthal number m is determined. The frequency of the mode is kHz and propagates along with , the electron gyration direction. The phase difference between probes A and E located at similar θ but different axial positions is shown in Fig. 6(b), indicating that the mode propagates along , from the plasma gun tip to the bias plate. This mode has right-handed helicity, as expected for a kink instability. With the axial magnetic field reversed, , as is shown in Fig. 6(c), the mode helicity flips to left handed, with along the electron gyration direction and kz pointing from the gun tip toward the bias plate [Fig. 6(d)], .
Propagation of a coherent perturbation along the azimuthal and axial directions for two different experimental scenarios: background magnetic field Bz parallel and anti-parallel to Jz. Contour plots (a) and (c) of fluctuations as measured with the four magnetic probes, A to D, demonstrate that the fluctuation propagates along when and along when . (b,d) show time series of fluctuations obtained from axially spaced probes B and E, indicating for both scenarios. Thus, the fluctuations exhibit right (left)-handed helicity for ( ), consistent with a kink instability.
Propagation of a coherent perturbation along the azimuthal and axial directions for two different experimental scenarios: background magnetic field Bz parallel and anti-parallel to Jz. Contour plots (a) and (c) of fluctuations as measured with the four magnetic probes, A to D, demonstrate that the fluctuation propagates along when and along when . (b,d) show time series of fluctuations obtained from axially spaced probes B and E, indicating for both scenarios. Thus, the fluctuations exhibit right (left)-handed helicity for ( ), consistent with a kink instability.
One consequence of the right (left)-handed helicity of the kink instability is its paramagnetic effect.9,35 In laboratory experiments, proper accounting of contributions from other factors is also needed to reveal the paramagnetic effect of kink instability. For example, diamagnetic effects due to the plasma pressure gradient must be quantified and included in any analysis. For PHASMA, the expected decrease in the axial magnetic field due to diamagnetic effects is given by Gauss. As is shown in Fig. 7, the net change in the axial magnetic field during the flux rope evolution is - 10 Gauss, meaning that the paramagnetic contribution from the kink instability is roughly 4 Gauss, close to the estimated increase in the axial magnetic field that should result from the azimuthal component of the bias current Gauss. For context, we note that for high-current operation in the RWM, Pa led to a diamagnetic effect of Gauss, while twisting of their 5 kA bias current introduced Gauss. The overall paramagnetic effect was, therefore, Gauss, close to their reported result.14 Based on these measurements, the coherent fluctuation observed in PHASMA exhibits the expected paramagnetic feature of a kink instability.
The spatial-temporal evolution of the change in the axial magnetic field, . The change of Gauss results from the paramagnetic effect of the kink instability (5–6 Gauss) superimposed on the diamagnetic effect due to the radial pressure gradient (−14 Gauss).
The spatial-temporal evolution of the change in the axial magnetic field, . The change of Gauss results from the paramagnetic effect of the kink instability (5–6 Gauss) superimposed on the diamagnetic effect due to the radial pressure gradient (−14 Gauss).
C. Growth rate
A kink instability is expected to grow on an MHD timescale, in PHASMA.15,36 Since τA is smaller than one oscillation of the observed m = 1 structure, we have focused our attention on the growth of the amplitude of the very first oscillation of . There are two contributions to in the initial phase of the flux rope formation. One contribution arises from the increasing bias current in the flux rope, and the other comes from the increase (if any) in the amplitude of the m = 1 perturbation. Figure 8(a) shows the experimentally measured (blue line) as a function of time. The black dashed line shows the calculated increase in at the fixed radial location, r, of the probe due to the axial plasma current. The difference between these two signals results from the perturbation of the kink instability shifting the axis of the plasma current, dr/dt. , where RP is magnetic probe's initial radial location.
(a) Measured (blue solid line) and the calculated based on the bias current (black dashed line) for the location of the magnetic probe; (b) difference between measured and calculated (blue line) and the quadratic fit (red dashed line); the dotted vertical red line denotes the time at which the kink instability begins to grow.
(a) Measured (blue solid line) and the calculated based on the bias current (black dashed line) for the location of the magnetic probe; (b) difference between measured and calculated (blue line) and the quadratic fit (red dashed line); the dotted vertical red line denotes the time at which the kink instability begins to grow.
The initial phase of the perturbation, shown with a dotted red quadratic fit line in Fig. 8(b), demonstrates that the first quarter of first oscillation of , from t = 0.160 ms to t = 0.183 ms, is not a sinusoidal wave but a signal that is growing in time. Note that the bias current is 160 A when the kink instability starts to grow at t = 0.160 ms, completely consistent with a threshold value for q = 2. The bias current is growing linearly during this short period, at a rate of A/μs. Therefore, the quadratic nature of the growth implies that dr/dt is constant, i.e., the kink instability is growing linearly. Linear growth of the shift in r is independently verified with camera images, in which mm/μs is observed. Using measured values, the predicted amplitude of the quadratic term in is Gauss/ms2, comparable to the value of 300 Gauss/ms2 obtained from the quadratic fit. Linear growth of a kink instability has also been reported in other experimental systems.15,16
IV. EVIDENCE FOR AN ALFVÉNIC MODE
The real frequency of a kink instability in the plasma frame is zero.37 In the laboratory frame, observed finite frequencies are typically attributed to Doppler shifts introduced by either axial flow or E × B drift. The observed frequency of the kink instability reported here cannot be explained by E × B drift. The observed frequency vs axial magnetic field strength is shown in Fig. 9(a). As is shown in Fig. 9(b), the calculated E × B drift frequency is one order larger than the observed frequency and decreases monotonically with increasing background axial magnetic field strength Bz. The radial electric field is derived from calculations of the plasma potential38,39 based on the values of Vf and Te as measured with the triple Langmuir probe, . According to measurements in RSX,11,21 which operated with similar plasma parameters as used in this study, the Doppler shift from axial flow contributes to the observed frequency of the kink instability. Figure 9(c) shows the estimated Doppler shift due to axial flow,
where the last term includes the contribution from a bias current larger than the threshold value12,36 and axial flow speed is given by , where cs is the ion sound speed. While comparable in magnitude to the observed frequencies, the scaling of fVZ with the increasing magnetic field is completely inconsistent with the observations. Thus, the discrepancy between f and fVZ suggests the existence of another mode with finite real frequency.
(a) Observed frequency of coherent fluctuation as a function of applied axial magnetic field strength. (b) Estimated Doppler shifted frequency that would result from E × B drift. (c) Estimated Doppler shift frequency fvz that would result from axial flow Vz.
(a) Observed frequency of coherent fluctuation as a function of applied axial magnetic field strength. (b) Estimated Doppler shifted frequency that would result from E × B drift. (c) Estimated Doppler shift frequency fvz that would result from axial flow Vz.
A. Fluctuation frequency spectra
The evolution of the frequency spectrum of the fluctuations [Fig. 10(a)] during the relatively long 10 ms discharges in PHASMA (compared to ∼1 ms in RSX) provides an opportunity to identify the wave mode at later times after the amplitude of the global kink instability decreases. We use the inductance of the flux rope as one indicator of amplitude of global kink instability. The global twisting (deformation) of a flux rope introduces extra inductance into the external bias current circuit. The bias circuit inductance is related to the circuit impedance by
where Lcircuit and Rcircuit are the inductance and resistance in the bias circuit. Based on Eq. (2), Fig. 10(b) shows and in coordinates of along the x axis and along the y axis. Different colors identify different time intervals in the discharge. The slope of this plot is Lcircuit and is shown in Fig. 10(c) as a function of discharge time. The decrease in Lcircuit from ms to 2 ms is a clear indication that the amplitude of the global kink instability is decreasing. For comparison, the time evolution of Lcircuit for Ibias = 150, i.e., for a discharge current below the kink instability threshold, , is shown in Figs. 10(e) and 10(f). Below threshold, no additional inductance appears in the bias circuit from ms to 2 ms. Thus, the measured inductance provides a crude indicator of the amplitude of the global kink instability in our experiment.
(a) The frequency spectrum of for A; (b) temporal evolution of and in coordinates of along the x axis and along the y axis. The colors indicate time intervals during the discharge. The slope of these data equals the effective inductance Lcircuit in the bias circuit. (c) Temporal evolution of Lcircuit. (d)–(f) Similar to (a)–(c), but for A.
(a) The frequency spectrum of for A; (b) temporal evolution of and in coordinates of along the x axis and along the y axis. The colors indicate time intervals during the discharge. The slope of these data equals the effective inductance Lcircuit in the bias circuit. (c) Temporal evolution of Lcircuit. (d)–(f) Similar to (a)–(c), but for A.
The evolution of the frequency spectrum shown in Fig. 10(a) is obtained through a wavelet analysis of the fluctuations. As the amplitude of global kink instability decreases from ms to 2 ms, both the amplitude and the frequency of the fluctuations decrease. After t = 2 ms, the contribution from global kink instability to the observed fluctuation frequency is negligible since the global deformation of the flux rope, as indicated by the dramatically smaller inductance, has ceased. Note that there is no change in the fluctuation amplitude and frequency without excitation of the kink instability [Fig. 10(d)]. The presence of very weak fluctuation, still discernible in Fig. 10(d), implies that the observed wave modes are preexisting fundamental modes in the plasma and the kink instability simply provides more free energy to this mode. Also note that frequency and amplitude of the mode remain almost constant even after the bias current has dropped by roughly 30% from t = 2 ms to 8 ms—further evidence that this wave mode is not a global kink instability.
From t = [0.5 ms, 2 ms] to t = [2 ms, 8 ms], the frequency of the fluctuations decreases by kHz, consistent with the expected Doppler shift of a global kink instability arising from axial flow. The predicted value of kHz, close to observed frequency, is consistent with Ryutov's model for non-lined tied kink instabilities in a flowing plasma.12 While the kink instability only persists until t = 2 ms, 2 ms is still much longer than the Alfvén transit time, , and is comparable to the “long-lifetime” kink instabilities observed in RSX.21 These PHASMA discharges last long enough that other plasma modes are readily identifiable.
B. Coupling to the kink instability
To investigate the relation between the kink instability and this new wave mode, we investigated the dependence of the amplitude of the mode on bias current from t = 2 ms to 8 ms, shown as blue circles in Fig. 11(a). When Ibias < 180 A, the amplitude of the mode increases with increasing bias current, but the rate of increase is much smaller than that for Ibias > 180 A. The predicted threshold current for the kink is A. The amplitude of the global kink instability is also calculated by subtracting the amplitude of the additional mode from the total amplitude of fluctuations from t = 0.5 ms to 2 ms [red squares in Fig. 11(a)]. The resultant amplitude of the kink instability is shown as black circles in Fig. 11(b), verifying that Ibias = 180 A is the threshold bias current for the global kink instability. Therefore, this additional large amplitude mode draws energy from the kink instability once it is excited.
(a) The dependence of amplitude of fluctuations from t = 2 ms to 8 ms (blue circles) and for t = 0.5 ms to 2 ms (red squares) as a function of bias current. (b) The dependence of global kink instability amplitude, deduced from the difference in the two amplitudes in (a) as a function of current.
(a) The dependence of amplitude of fluctuations from t = 2 ms to 8 ms (blue circles) and for t = 0.5 ms to 2 ms (red squares) as a function of bias current. (b) The dependence of global kink instability amplitude, deduced from the difference in the two amplitudes in (a) as a function of current.
C. Alfvénic nature of the new mode
The kink instability is just one branch of the Alfvén surface wave spectrum.40 Thus, it seems reasonable that other Alfvénic wave modes could couple to and draw energy from the kink instability. Two features of this mode that point toward an Alfvénic nature are the mode's frequency being less than the ion gyrofrequency and the strong perpendicular magnetic field fluctuations, .
Very low frequency Alfvénic waves in an ideal MHD plasma propagate along the background magnetic field at the Alfvén speed, VA. In PHASMA, the wave frequency is comparable to the ion gyrofrequency, , and thus, finite frequency effects must be included in the wave dispersion relation,
where kz is the parallel wavevector and is the average Alfvén speed, .21,41 The phase angle difference between the axially separated magnetic probes B and E is obtained from a correlation analysis and used to calculate kz. Figure 12(a) shows kz as a function of bias current, black squares for Bz = 225 Gauss and blue squares for Bz = 375 Gauss. The coherence found in the correlation analysis is large, > 0.9, only for large bias currents, and so smaller bias current Ibias < 300 A is not included in these plots. The plasma density and electron temperatures measured with the triple Langmuir probe do not increase with bias current when Ibias > 300 A, consistent with the relatively constant values of kz seen in Fig. 12(a) if the mode is Alfvénic. The corresponding frequency of this mode fAW is shown in Fig. 12(b). Similar to the results for kz, the frequency remains almost constant with increasing bias current Ibias, but increases with increasing background magnetic field Bz. The theoretical predictions of fAW based on the Alfvén dispersion relation of Eq. (3) are shown in Fig. 12(c) as a function of Bz. The uncertainty in fAW shown by magenta shading is due to an assumed 10% uncertainty in both kz and ne. The experimental measurements of fAW are plotted as black and blue squares and are consistent with the Alfvén dispersion relation.
(a) Parallel wavevector kz and (b) wave frequency as a function of bias current for Bz = 225 Gauss (black squares) and 375 Gauss (blue squares). (c) Range of theoretical predictions for the Alfvén wave frequency fAW as a function of background magnetic field Bz (shaded region) assuming 10% uncertainties in ne and kz. Black and blue circles are the measured values of fAW for Bz = 225 Gauss and 375 Gauss.
(a) Parallel wavevector kz and (b) wave frequency as a function of bias current for Bz = 225 Gauss (black squares) and 375 Gauss (blue squares). (c) Range of theoretical predictions for the Alfvén wave frequency fAW as a function of background magnetic field Bz (shaded region) assuming 10% uncertainties in ne and kz. Black and blue circles are the measured values of fAW for Bz = 225 Gauss and 375 Gauss.
Interestingly, the dispersion relation also predicts the observed increase in fAW with increasing Bz, a key feature that could not be explained by Doppler shift effects arising from axial flow or E × B drift. We note that other kink experiments have also reported an increase in the kink instability frequency with increasing Bz.14,42 Fluctuations with somewhat similar characteristics were observed on the Encore Tokamak at similar plasma parameters. These fluctuations were identified as drift Alfvén waves.43 However, in PHASMA, the coupling between drift waves and Alfvén waves is not expected to be significant because the expected parallel phase speed (60 km/s) of the drift wave is much larger than VA.
D. Frequency limit of the mode
If the observed wave mode has Alfvénic nature, it should exhibit significant changes for frequencies near the ion cyclotron frequency. To determine if there is a frequency limit for this mode, we compare measurements in two working gases, argon ( kHz) and hydrogen ( kHz). Figure 13 shows the dependence of fAW and normalized mode frequency on the bias current (argon as black squares and hydrogen as blue squares). In contrast to a nearly constant fAW in argon, fAW in hydrogen increases monotonically with bias current. These observations are consistent with those in other kink experiments in which no limit in mode frequency was seen when the normalized mode frequency is small, e.g., in hydrogen RSX plasmas.11,14 Whether or not these observations point to resonance at is not yet clear.
The dependence of the mode frequency fAW and the normalized mode frequency on bias current for (a) argon plasma and (b) hydrogen plasma. A frequency limit is observed in argon plasma when fAW is a significant portion of fci.
The dependence of the mode frequency fAW and the normalized mode frequency on bias current for (a) argon plasma and (b) hydrogen plasma. A frequency limit is observed in argon plasma when fAW is a significant portion of fci.
V. CONCLUSION
In summary, a magnetic flux rope is successfully created with a pulsed plasma gun operating with argon gas in the PHASMA facility. A kink instability appears when the bias current exceeds , consistent with the prediction for the threshold current value for a non-line tied axial boundary condition. Both magnetic and visible light intensity fluctuation measurements indicate that the kink instability is globally unstable. From measurements of and kz, the m = 1 fluctuations have right (left)-handed helicity when ( ), exactly as predicted for a kink instability. The predicted paramagnetic effect of a kink instability is measured experimentally, and the magnitude of the paramagnetism is consistent with predictions for a kink instability. Linear growth of the kink instability after onset is observed. However, the frequency of the observed fluctuations f is much larger than that explainable by a Doppler shift, fVZ, due to axial flow. The differential frequency of the fluctuations is attributed to another wave mode co-existing with kink instability. This wave mode is driven by the kink instability since the mode amplitude dependence on bias current tracks that of the kink instability. The new wave mode is consistent with an Alfvén wave dispersion relation. The frequency of the mode is constrained when it reaches a significant portion of fci. It is unclear if an Alfvénic resonance is occurring at . In future studies, the ion temperature and bulk flow speed will be measured with LIF44,45 to identify any possibly resonant interaction of the ions with the Alfvén wave that is coupled to the kink instability.
ACKNOWLEDGMENTS
We would like to express our deep appreciation to the late Dr. Thomas Intrator who constructed the original RSX device from which many of the components of PHASMA come. We also thank Dr. Glen Wurden of the Los Alamos National Laboratory who helped with the transfer of the RSX components to the PHASMA Laboratory. We thank the plasma physics group at the University of Wisconsin-Madison for the use of the plasma guns. This work was supported by NSF Award PHYS Nos. 1827325 and 1902111, NASA Award No. 80NSSC19M0146, and DoE Award No. DE-SC0020294.
DATA AVAILABILITY
The data that support the findings of this study are available from the corresponding authors upon reasonable requests.