ISSN: 2320-2459

**Nafis Ahmad ^{1*}, Moiz Ahmad^{2}**

^{1}Department of Physics, College of Science, UAE University, PO Box 15551, Al-Ain, United Arab Emirates

^{2}Department of Physics, Integral University, Lucknow 226 026, India

- *Corresponding Author:
- Nafis A

Department of Physics

College of Science, UAE University

PO Box 15551, Al-Ain, United Arab Emirates

**Tel:**+971 3 767 333

**E-mail:**nafis.jmi@gmail.com

**Received Date:** 23/03/2016; ** Accepted Date:** 06/05/2016; **Published Date:** 08/05/2016

**Visit for more related articles at** Research & Reviews: Journal of Pure and Applied Physics

Three wave parametric instability of a large amplitude lower hybrid wave in two-ion species plasma is studied analytically. The pump decays into a low frequency mode or quasi mode in the ion cyclotron range of frequency and a lower hybrid wave sideband. The dominant channel of decay is the one for which the low frequency mode is an ion cyclotron wave with frequency close to the ion cyclotron frequency of either ion species. For typical D-T plasma of a tokamak, the growth rate for ω close to deuterium cyclotron frequency, increases with the wave number of the ion cyclotron wave and decreases as the ratio of deuterium to tritium density increases.

Parametric Instabilities; Two-ion species plasma; Lower hybrid wave; Tokamak; Radio frequency heating

Parametric instabilities have been an important feature of large amplitude lower hybrid waves in tokamak as well as smaller plasma devices. During the first decades of nuclear fusion research, high power lower hybrid waves in the 500 MHz-1 GHz frequency range were candidates for heating magnetized plasmas to thermonuclear temperatures. More recently, such waves were found to be mostly attractive at higher frequency (1-5 GHz) for driving non-inductive currents in tokamak plasmas, thus opening the possibility of running tokamak in steady state, and leading to the conceptual design of a steady state tokamak reactor. Radio frequency (RF) heating and current drive of tokamak in the lower hybrid range of frequency has been an active field of research for over three decades. In this frequency range, toroidal plasma current is sustained by replenishing the collisional momentum loss of current carrying suprathermal electrons while ions are heated directly through the Landau damping or via the excitation of parametric instabilities [1-5]. The prominent channels of three wave parametric decay involve an ion-cyclotron mode or quasi mode and a lower hybrid sideband. In the four wave coupling oscillating two stream instability is the most dominant process. Resonant decay into two lower hybrid waves is also a completing process [6-11].

The current drive experiments on a number of tokamak machines have been successfully performed. Fully non inductive discharges of up to 3.6 MA in JT-60U [12], 3 MA in the Joint European Torus (JET) [13] and 0.5 MA during 6 min in Tore Supra [14] have been achieved with lower hybrid waves. The projected ITER lower hybrid systems will inject approximately 50 MW of power and they are optimized for the off-axis current drive functions [15]. Liu et al. have investigated the Lower hybrid wave (LHW) heating in deuterium plasma in the HT-7 tokamak [16]. The ion and electron temperatures were increased by 0.4 keV and 0.3 keV respectively for P_{LH}∼300kW. The ion heating has been observed in the HL-1M tokamak when the plasma density exceeded 3.5 × 10^{13} cm^{-3} [17]. LHW experiments have been performed in HT-6M tokamak by Li et al., they have reported the quasi-steady state H-mode with high plasma density by the injection of lower hybrid heating (LHH) and lower hybrid current drive (LHCD) with a power threshold of 50 kW [18]. A good confinement was obtained in ASDEX by combining LHCD and neutral beam injection (NBI) [19]. A full wave current drive with a very high density (4.5 × 10^{19} m^{-3}) was obtained in the high-field Frascati Tokamak Upgrade (FTU) [20]. Fisch [21,22] has given an elegant review of lower hybrid wave heating and current drive in tokamak. Recently, Ahmad [23] has pointed out the parametric excitation of ion-ion hybrid mode by a lower hybrid wave in D-T plasma where growth is faster for higher deuterium to tritium density ratio.

Observations in space indicate that LH waves are among the most important waves in the earth’s magnetosphere. In fact, LH wave play a central role in the process of collision less energy and momentum transport in space plasmas. Recently there has been much interest [24,25] in magneto sonic and other waves in multi-ion plasmas.

Two ion species plasma in tokamak has the possibility of extra channels of parametric decay. In this paper, we study the parametric decay of a large amplitude lower hybrid wave into an ion cyclotron wave and lower hybrid sideband wave in two ion species plasma such as those considered for nuclear fusion. The channels of decay are: i) resonant decay into an ion cyclotron mode near the ion cyclotron frequency of the either ion species, ii) nonlinear cyclotron damping on either ion species. The lower hybrid pump wave imparts an oscillatory velocity to electrons. The latter beats with the density perturbation due to the low frequency mode to produce a nonlinear density perturbation, driving the sideband lower hybrid wave , where The sideband couples with the pump to produce a ponder motive force on electrons that drive the low frequency mode.

In Sec. 2, we obtain the linear response of electrons to the pump and sideband waves and the low frequency mode. In Section 3, we study the nonlinear coupling and obtain the nonlinear dispersion relation and the growth rate. The results have been discussed in Section 4.

Consider a two-ion species plasma (e.g., tokamak) of equilibrium electron density , in a static magnetic field The ion species are characterized by mass m1, m2, density and charge Z1e, Z2e such that A large amplitude lower hybrid pump wave propagates through the plasma in the x-z plane with electrostatic potential

(1)

(2)

(3)

where is the mass of electrons, and -e and m are the electron charge and mass.

The pump wave imparts oscillatory velocity to electrons which is governed by the equation of motion

where we have ignored the pressure term. Linearizing this equation and replacing

we obtain

(4)

(5)

Using the velocity perturbation in the linearized equation of continuity

we obtain the density perturbation

(6)

The pump wave decays into an electrostatic ion cyclotron mode / quasimode of potential φ and a lower hybrid sideband wave of potential

(7)

(8)

where

The linear response of electrons to the sideband at is same as given by Equations (4–5), with replaced by (9)

(10)

The linear density perturbations of electrons and ions at can be written in terms of electron and ion susceptibilities (11)

where

.

where are the ions thermal speeds, are the ion-cyclotron frequencies and are the modified Bessel functions of order n and arguments b1 and b2, and

The sideband couples with the pump to produce a low frequency ponder motive force on the electrons. has two components, perpendicular and parallel to the magnetic field. The response of electrons to is strongly suppressed by the magnetic field and is usually weak. In the parallel direction, the electrons can effectively respond to hence, the low frequency nonlinearity arises mainly through

The parallel pondermotive force on electrons at can be written as

(14)

The pondermotive potential turns out to be

(15)

The electron density fluctuations in response to and the self-consistent potential φ at can be written as

(16)

The density perturbations of the two-ion species can be written as

(17)

Using Eqs. (16) and (17) in the Poisson’s equation

we obtain,

(18)

The density perturbation at couples with the oscillatory velocity of electrons to produce a nonlinear density perturbation at the sideband Solving the equation of continuity for the nonlinear density perturbation at the sideband,

we obtain

(19)

The linear density perturbation of electrons due to the self consistent potential (20)

The ion density perturbations can be written as

(21)

where

Here we have ignored the nonlinearity arising through ions as it is suppressed by their large mass.

Using these electron and ion density perturbations in the Poisson’s equation

we obtain

(22)

Equations (18) and (22) yield the nonlinear dispersion relation

ε ε=μ, (23)

where

(24)

We solve Eq. (23) in the case of resonant decay.

In case and in the absence of the pump, satisfy the linear dispersion relations corresponding to ion cyclotron and lower hybrid waves respectively, the decay process is termed as resonant decay.

In the vicinity of when the linear dispersion relation for the low frequency mode (=0) takes the form

giving

(25)

where

The linear dispersion relation for the sideband gives

(26)

In the presence of nonlinear coupling we write

and expand ε and εs around ω_{r} and ω_{s}r as

(27)

or

(28)

where we have suppressed the subscript r on ω_{r} and ω_{s} for the sake of brevity. Then the nonlinear dispersion relation gives

(29)

Equation (29) simplifies to

(30)

where δ is the angle between

In order to have numerical appreciation of the growth rate, we have carried computations of growth rate for the following parameters: (D-T plasma), We have plotted in **Figure 1**, the variation of normalized growth rate, as a function of (where

At high power, lower hybrid waves are prone to parametric decay because of their large electrostatic component in the direction perpendicular to the equilibrium magnetic field. This was particularly the case in experiments where the wave frequency was chosen for ion heating through the lower hybrid resonance. A high fraction of the injected power was found to decay nonlinearly into lower frequency daughter waves, and this phenomenon often prevented the power to reach the plasma core and heat the plasma efficiently. At the higher frequencies required for current drive, the question of parametric decay seems to be less severe, but it remains open, especially at high density, and it is the subject of extensive theoretical and experimental investigations. Lower hybrid wave is parametrically unstable when the oscillatory velocity of electrons is of the order or greater than the sound velocity. For tokamak parameters, the powers required are of the order of ≥10 MW. The decay into an ion cyclotron wave and a lower hybrid wave possesses large growth rate. The coupling between high and low frequency modes is provided primarily by the parallel ponder motive force on electrons involving drift.

For a typical D-T plasma of a tokamak, the growth rate for ion cyclotron wave with frequency ω close to deuterium cyclotron frequency increases with the wave number of the ion cyclotron wave. However, the growth rate decreases as the ratio of deuterium to tritium density increases.

- Kindel JM et al.Parametric Instabilities and Anomalous Heating of Plasmas near the Lower Hybrid Frequency. Phys. Rev. Lett. 1972;29:995.
- PorkolabM.Fusion. Phys. Fluids.1977;20
**:**2058. - Pamela J and JET EFDA ContributorsStability and Transport in Magnetic Confinement Systems.Nucl. Fusion.2003;43
**:**1540. - Ide S et al.Sustainment and modification of reversed magnetic shear by LHCD on JT-60U. Plasma Phys. Control. Fusion. 1996;38:1645.
- Hoang GT et al.Improved confinement in high li lower hybrid driven steady state plasmas in TORE SUPRA. Nucl. Fusion. 1994;34:75.
- Ott E. Decay instability of lower hybrid waves. Phys. Fluids. 1975;18:566.
- Reiman A. Parametric excitation of Alfvén wave by a lower hybrid wave. Phys. Fluids. 1980;23:1254.
- Tripathi VK and Sivaramaprasad PV.Oscillating two-stream instability in a magnetized electron-positron-ion plasma.Plasma Physics. 1989;41:13.
- Panchenko VG et al. Fluorescent Tracers of Cell Morphology and Fluid Flow. Physica Scripta. 1985;13:594.
- Kuo SP and Lee MC.Oscillating twoÃ¢ÂÂstream instability of a ducted whistler pump.Phys. Fluids. 1984;
**27:**1434. - TripathiVK et al. Electromagnetic particle simulation of the effect of toroidicity on linear mode conversion and absorption of lower hybrid waves. Phys. Fluids. 1977;20:1525.
- Naito O and the JT-60 TeamSteady state plasma performance on JT-60U. Plasma Phys
*.*Control.Fusion. 1993;35:B215. - SoldnerFX. Plasma Physics and Controlled Nuclear fusion Research, Seville, 1994;1:423.
- Van HD et al.One minute pulse operation in the TORE SUPRA tokamak. Nucl.Fusion. 2004;44:L11.
- PorkolabMet al.Quantum Stress Tensor in Schwarzschild Space-Time. Phys. Rev. Lett. 1984;53:450.
- LiuSXet al. Studies of high confinement plasma with lower hybrid current drive in the HT-7 superconducting tokamak. Plasma Phys. Controlled Fusion. 2001;43:805.
- Liu Y et al. International Atomic Energy. 1998.
- Jiangang Li et al.The plasma current profile during current reversal in AC operation of the CT-6B tokamak.Nuclear Fusion. 2000;40.
- Soldner FXet al. Electron and ion heat transport with lower hybrid current drive and neutral beam injection heating in ASDEX. Nucl. Fusion.1993;33:333.
- Pericoli-Ridolfini V. et al. High Plasma Density Lower-Hybrid Current Drive in the FTU Tokamak. Phys. Rev. Lett.1999;82:93.
- Fisch NJ. Confining a Tokamak Plasma with rf-Driven Currents.Phys. Rev. Lett.1978;41:873.
- FischNJ.Theory of current drive in plasmas.Rev. Mod. Phys
*.*1987;59:175. - Nafis A.Parametric excitation of ion-ion hybrid mode by a lower hybrid wave in a D-T plasma. Eur. Phys. J. D 2014;68:371.
- Boldyrev S. Nonlinear magnetosonic waves in a multi-ion-species plasma. Phys. Plasmas. 1998;5:1315.
- Irie S and Ohsawa Y. Nonstochastic Particle Acceleration in Collisionless Shock Waves. Phys. Plasmas. 2003;10:1253.