<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">JMP</journal-id><journal-title-group><journal-title>Journal of Modern Physics</journal-title></journal-title-group><issn pub-type="epub">2153-1196</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/jmp.2016.710101</article-id><article-id pub-id-type="publisher-id">JMP-67526</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Physics&amp;Mathematics</subject></subj-group></article-categories><title-group><article-title>
 
 
  Stability Analysis of Electromagnetic Ordinary and Extraordinary Modes
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>N.</surname><given-names>Noreen</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>S.</surname><given-names>Zaheer</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>H.</surname><given-names>A. Shah</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib></contrib-group><aff id="aff1"><addr-line>Forman Christian College, Chartered University, Lahore, Pakistan</addr-line></aff><aff id="aff2"><addr-line>Government College University, Lahore, Pakistan</addr-line></aff><pub-date pub-type="epub"><day>08</day><month>06</month><year>2016</year></pub-date><volume>07</volume><issue>10</issue><fpage>1120</fpage><lpage>1131</lpage><history><date date-type="received"><day>5</day>	<month>April</month>	<year>2016</year></date><date date-type="rev-recd"><day>accepted</day>	<month>16</month>	<year>June</year>	</date><date date-type="accepted"><day>21</day>	<month>June</month>	<year>2016</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p><html>
 <head></head>
 
  By using kinetic theory, we derived the general dispersion relations for ordinary mode (O-mode) and Extra-ordinary mode (X-mode) in anisotropic magnetized plasma. The effects of energy anisotropy, magnetic field to density ratio (
  <img src="Edit_c0be58c2-48b6-42e2-89ac-92b8394a5873.bmp" alt="" />) and the plasma beta 
  <img src="Edit_0d77fa39-7c76-4410-a7f0-596ea2631b94.bmp" alt="" /> on the propagation characteristics, have been analyzed. The stability analysis and the growth rates have been presented. The marginal threshold condition for oscillatory and purely growing mode has been obtained for higher harmonics and we have also calculated their growth rates in terms of plasma beta 
  <img src="Edit_3bb9fb4f-d3bf-4723-9c2d-209255cf69e7.bmp" alt="" /> and energy anisotropy 
  <img src="Edit_c936187c-0323-49a0-9130-4d3e80ed6d3d.bmp" alt="" />. The X-mode satisfies the instability condition according to difference of geometry with the O-Mode. These modes are important for spherical tokamaks, and their coupling leads to the generation of the Bernstein mode, which causes the heating effects.
 
</html></p></abstract><kwd-group><kwd>Instabilities</kwd><kwd> Growth Rate</kwd><kwd> Anisotropy</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>The ordinary mode (O-mode) is a linearly polarized electromagnetic perpendicularly propagating wave, which propagates only when wave frequency is greater than the plasma frequency. The work is related to the electromagnetic cyclotron harmonic instability for its possible role in solar and interplanetary radio emission processes where the ratio <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x10.png" xlink:type="simple"/></inline-formula> (where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x11.png" xlink:type="simple"/></inline-formula> is the electron plasma frequency and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x12.png" xlink:type="simple"/></inline-formula> is the electron cyclotron frequency) is relatively high i.e., the ratio is of the order of 10 or can be as high as 50 or even 100 near 1 a.u. It may be useful for the heating and current drive mechanism in the spherical tori like the NSTX [<xref ref-type="bibr" rid="scirp.67526-ref1">1</xref>] and MAST [<xref ref-type="bibr" rid="scirp.67526-ref2">2</xref>] where the <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x13.png" xlink:type="simple"/></inline-formula></p><p>It is found that extraordinary mode (X-mode) power is not absorbed at the cyclotron resonance but uniquely at the upper hybrid resonance, displaced to the low field side of the cyclotron resonance. O-mode power, however, is absorbed at the cyclotron resonance as well. The displacement of the upper hybrid resonance to the low field side with O-mode launch is significantly smaller than that with X-mode launch because of the lower densities produced by O-mode launch at the same microwave power level [<xref ref-type="bibr" rid="scirp.67526-ref3">3</xref>] . Hamasaki [<xref ref-type="bibr" rid="scirp.67526-ref4">4</xref>] [<xref ref-type="bibr" rid="scirp.67526-ref5">5</xref>] investigated the electromagnetic o-mode instability with perpendicularly propagating waves for a two temperature Maxwellian distribution function. Lee [<xref ref-type="bibr" rid="scirp.67526-ref6">6</xref>] studied the same mode in counterstreaming plasmas and showed that the ordinary mode became unstable as the magnetic field changed. Later Bornatici and Lee [<xref ref-type="bibr" rid="scirp.67526-ref7">7</xref>] worked on O-mode and determined that for counterstreaming plasmas an instability occurred if the streaming velocity exceeded a certain threshold value which can be below the required velocity to excite the electrostatic two-stream instability. They also concluded that whereas the perpendicular temperature stabilized the effect the parallel temperature enhanced the instability. Shivamoggi [<xref ref-type="bibr" rid="scirp.67526-ref8">8</xref>] also discussed the destabilization of the O-mode due to magnetic field and thermal effects. Ibscher et al. [<xref ref-type="bibr" rid="scirp.67526-ref9">9</xref>] investigated the nonresonant Wieble mechanism which can drive the O-mode unstable. They studied the instability on the basis of a threshold which gave the instability conditions and upper limits of the growth rate. Their problem was restricted for fundamental harmonic only. Iqbal et al. [<xref ref-type="bibr" rid="scirp.67526-ref10">10</xref>] studied the O-mode in degenerate anisotropic plasmas and proposed the excitation of a new banded type of instability which grew at some particular values of temperature anisotropy. Hadi et al. [<xref ref-type="bibr" rid="scirp.67526-ref11">11</xref>] also revised the analysis of the O-mode instability with Maxwellian parallel distribution coupled with thermal ring perpendicular distribution. They demonstrated that O-mode for thermal ring distribution may be excited for cyclotron harmonics as well as for the purely growing branch, depending on the value of the normalized ring speed. Lazar</p><p>et al. [<xref ref-type="bibr" rid="scirp.67526-ref12">12</xref>] concluded that O-mode instability was driven by an excess of parallel temperature where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x14.png" xlink:type="simple"/></inline-formula></p><p>for <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x15.png" xlink:type="simple"/></inline-formula> Vafin et al. [<xref ref-type="bibr" rid="scirp.67526-ref13">13</xref>] derived the analytical marginal instability condition for magnetized plasmas when charged particles were distributed in counter-streams with equal temperatures. They confirmed the O-mode instability at small plasma beta values, when the parallel counter-stream free energy exceeded the perpendicular bi-Maxwellian free energy. Farrell [<xref ref-type="bibr" rid="scirp.67526-ref14">14</xref>] presented a theory in which he described the direct generation of electromagnetic O-mode emission via mildly energetic electron beams in a highly dense and warm plasma.</p><p>In this manuscript, the energy anisotropic Heaviside distribution function is used for understanding the behavior of O-mode and X-mode. Such distribution function provides the detailed information about banded emission of O-mode instability. Such type of emission has been observed in space plasmas, where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x16.png" xlink:type="simple"/></inline-formula> e.g. solar wind. Satellite wave instruments commonly detect banded magnetospheric emissions between har- monics of the electron gyrofrequency in the outer magnetosphere [<xref ref-type="bibr" rid="scirp.67526-ref15">15</xref>] . This type of banded emission has been observed in the terristial magnetosphere. Frequency-banded electromagnetic waves up to 2000 Hz are observed concurrently with warm energy-banded ions in the low latitude auroral and sub-auroral zones during every large geomagnetic storm, observed by the FAST and DEMETER satellites. The appearance of the banded wave activity suggests that there may be distinct changes in the geospace system that characterize large magnetic storms [<xref ref-type="bibr" rid="scirp.67526-ref16">16</xref>] .</p><p>Coupling of the O-mode and X-mode is a necessary tool for generation of the bernstein mode which is a powerful source of heating in spherical tokamaks. Literature shows the different methods of their coupling. But their unstable regions are a major problem in the coupling. Padoba et al. [<xref ref-type="bibr" rid="scirp.67526-ref17">17</xref>] first time demonstrated the conversion from an O-mode to an X-mode by probe measurements of amplitude and phase of the wave field in the conversion region. Cairns et al. [<xref ref-type="bibr" rid="scirp.67526-ref18">18</xref>] used sheared magnetic field to calculate the linear conversion of the O-mode to the X-mode. Because electron Bernstein waves are analyzed as possible candidates for heating spherical tokamaks. Ram et al. [<xref ref-type="bibr" rid="scirp.67526-ref19">19</xref>] developed a kinetic model for studying the energy flow transfer between the X-mode, the O-mode and the EBW in the mode conversion region in the vicinity of the cold plasma upper hybrid resonance. Sodha et al. [<xref ref-type="bibr" rid="scirp.67526-ref20">20</xref>] derived the dispersion relation for modulational instabilities of a Gaussian electromagnetic beam propagating in the two modes: O-mode and X-mode, along the externally applied d.c. magnetic field, in a homogeneous magnetoplasma.</p><p>The layout of this paper is as follow. Section 2 gives information about the mathematical model of O-mode and X-mode. This section presents the stability analysis and calculates the maximum growth rate. A brief summary of results and discussions is given in Section 3. Section 4 will conclude the results.</p></sec><sec id="s2"><title>2. Mathematical Model</title><sec id="s2_1"><title>2.1. The Ordinary Mode (O-Mode)</title><p>By using kinetic model, the general dispersion relation for perpendicularly propagating O-mode with <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x17.png" xlink:type="simple"/></inline-formula> in collionless plasma is as follow [<xref ref-type="bibr" rid="scirp.67526-ref21">21</xref>]</p><disp-formula id="scirp.67526-formula1460"><label>(1)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/9-7502712x18.png"  xlink:type="simple"/></disp-formula><p>Here <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x19.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x20.png" xlink:type="simple"/></inline-formula> is distribution function.</p><p>The energy anisotropic Heaviside distribution function is [<xref ref-type="bibr" rid="scirp.67526-ref22">22</xref>] [<xref ref-type="bibr" rid="scirp.67526-ref23">23</xref>]</p><disp-formula id="scirp.67526-formula1461"><label>(2)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/9-7502712x21.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.67526-formula1462"><graphic  xlink:href="http://html.scirp.org/file/9-7502712x22.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x23.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x23.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x24.png" xlink:type="simple"/></inline-formula> are the effective temperatures in the perpendicular and parallel directions defined as follows</p><disp-formula id="scirp.67526-formula1463"><graphic  xlink:href="http://html.scirp.org/file/9-7502712x25.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.67526-formula1464"><graphic  xlink:href="http://html.scirp.org/file/9-7502712x26.png"  xlink:type="simple"/></disp-formula><p>and their corresponding integrations yields the results</p><disp-formula id="scirp.67526-formula1465"><graphic  xlink:href="http://html.scirp.org/file/9-7502712x27.png"  xlink:type="simple"/></disp-formula><p>Using Equations (1) and (2), we obtain</p><disp-formula id="scirp.67526-formula1466"><label>(3)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/9-7502712x28.png"  xlink:type="simple"/></disp-formula><p>where</p><disp-formula id="scirp.67526-formula1467"><graphic  xlink:href="http://html.scirp.org/file/9-7502712x29.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.67526-formula1468"><graphic  xlink:href="http://html.scirp.org/file/9-7502712x30.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.67526-formula1469"><graphic  xlink:href="http://html.scirp.org/file/9-7502712x31.png"  xlink:type="simple"/></disp-formula><p>For principle harmonic i.e., n = 1, we get the following linear dispersion relation</p><disp-formula id="scirp.67526-formula1470"><label>(4)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/9-7502712x32.png"  xlink:type="simple"/></disp-formula><p>where<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x33.png" xlink:type="simple"/></inline-formula>.</p><p>We note that <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x34.png" xlink:type="simple"/></inline-formula> is the condition for instability.</p><p>However, for higher harmonics, the linear dispersion relation takes the form</p><disp-formula id="scirp.67526-formula1471"><label>(5)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/9-7502712x35.png"  xlink:type="simple"/></disp-formula></sec><sec id="s2_2"><title>2.2. The Extra-Ordinary Mode (X-Mode)</title><p>The general dispersion relation of the X-mode is as</p><disp-formula id="scirp.67526-formula1472"><label>(6)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/9-7502712x36.png"  xlink:type="simple"/></disp-formula><p>By using the simple mathematical analysis, the dispersion relation of the X-mode is</p><disp-formula id="scirp.67526-formula1473"><label>(7)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/9-7502712x37.png"  xlink:type="simple"/></disp-formula><p>In terms of A and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x38.png" xlink:type="simple"/></inline-formula>, the relation can be expressed as</p><disp-formula id="scirp.67526-formula1474"><label>(8)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/9-7502712x39.png"  xlink:type="simple"/></disp-formula><p>where</p><disp-formula id="scirp.67526-formula1475"><graphic  xlink:href="http://html.scirp.org/file/9-7502712x40.png"  xlink:type="simple"/></disp-formula></sec></sec><sec id="s3"><title>3. Results and Discussion</title><p>In this section we will discuss the stability condition and calculate the growth rate for different combinations of A and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x41.png" xlink:type="simple"/></inline-formula>.</p><p>We first numerically discussed the results obtained for the O-mode from Equation (3). Lee [<xref ref-type="bibr" rid="scirp.67526-ref12">12</xref>] has calculated O-mode for three harmonics with the Maxwellian distribution function and concluded that the mode is stable for the Maxwellian distribution. Ichimaru [<xref ref-type="bibr" rid="scirp.67526-ref24">24</xref>] has discussed the O-mode for higher harmonics with nonlocal effects and confirmed the existence of Azbel-Karner resonance when the wave frequency is multiple of electron cy- clotron frequency.</p><p>The banded emission is observed in plots of A vs <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x42.png" xlink:type="simple"/></inline-formula> in the case of energy anisotropic Heaviside distribution <xref ref-type="fig" rid="fig1">Figure 1</xref>. This banded emission strongly agrees with the results of Iqbal et al. [<xref ref-type="bibr" rid="scirp.67526-ref10">10</xref>] where the anisotropic Fermi Dirac distribution function was used. The wave provides a wide range of stable and unstable regions.</p><p>In <xref ref-type="fig" rid="fig1">Figure 1</xref>, the relation of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x43.png" xlink:type="simple"/></inline-formula> and anisotropy A is plotted, it provides a marginal threshold value. The</p><p>dotted curve shows that<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x44.png" xlink:type="simple"/></inline-formula>, this curve plays the role of threshold value between stable and unstable</p><p>O-mode. Below the dotted curve the condition <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x45.png" xlink:type="simple"/></inline-formula> satisfies and mode is unstable which is presented by dashed curve. Above that dotted curve the condition is <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x45.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x46.png" xlink:type="simple"/></inline-formula> it means that there is a stable region</p><p>i.e., the solid curve. The comparison of plots defines that for small <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x47.png" xlink:type="simple"/></inline-formula> the <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x47.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x48.png" xlink:type="simple"/></inline-formula> contains large value, this is the region where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x47.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x48.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x49.png" xlink:type="simple"/></inline-formula> is large enough to provide a growth rate much larger than the oscillatory frequency <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x47.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x48.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x49.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x50.png" xlink:type="simple"/></inline-formula> so these results strongly agree with the environment i.e., solar wind. For large anisotropy, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x47.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x48.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x49.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x50.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x51.png" xlink:type="simple"/></inline-formula>, or for large<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x47.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x48.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x49.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x50.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x51.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x52.png" xlink:type="simple"/></inline-formula>, the O-mode instability is faster than the firehose instability. Larger values of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x47.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x48.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x49.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x50.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x51.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x52.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x53.png" xlink:type="simple"/></inline-formula> means low magnetic fields or more dense and hotter plasma, these conditions can come across at different altitudes in the solar wind regime.</p><p>In series of Figures 2-5, growth rates of higher harmonics have been plotted. For analytical threshold we consider complex form <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x54.png" xlink:type="simple"/></inline-formula> in plots the solid lines represent the <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x54.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x55.png" xlink:type="simple"/></inline-formula> and dashed shows the <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x54.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x55.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x56.png" xlink:type="simple"/></inline-formula> where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x54.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x55.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x56.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x57.png" xlink:type="simple"/></inline-formula> for Figures 2-5.</p><fig id="fig1"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref></label><caption><title> Marginal stability condition</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/9-7502712x58.png"/></fig><fig id="fig2"  position="float"><label><xref ref-type="fig" rid="fig2">Figure 2</xref></label><caption><title> A = 0.1,<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x60.png" xlink:type="simple"/></inline-formula></title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/9-7502712x59.png"/></fig><fig id="fig3"  position="float"><label><xref ref-type="fig" rid="fig3">Figure 3</xref></label><caption><title> A = 0.1<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x62.png" xlink:type="simple"/></inline-formula></title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/9-7502712x61.png"/></fig><fig id="fig4"  position="float"><label><xref ref-type="fig" rid="fig4">Figure 4</xref></label><caption><title> A = 0.6<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x64.png" xlink:type="simple"/></inline-formula></title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/9-7502712x63.png"/></fig><fig id="fig5"  position="float"><label><xref ref-type="fig" rid="fig5">Figure 5</xref></label><caption><title> Growth rate for different A</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/9-7502712x65.png"/></fig><p>In <xref ref-type="fig" rid="fig2">Figure 2</xref>, there is a stable form of O-mode but at A = 0.1 and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x66.png" xlink:type="simple"/></inline-formula> the harmonics start to intersect with each other.</p><p>Figures 3-5 show the real part of dispersion relation and dependence of O-mode on magnetic field. As we increase values of magnetic field, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x67.png" xlink:type="simple"/></inline-formula>it becomes unstable and the first unphysical state generates as in <xref ref-type="fig" rid="fig3">Figure 3</xref>. These results also satisfy the marginal instability condition as discussed earlier numerically. In above plots, noticeable thing is the value of A = 0.1. The parallel streaming is dominating in O-mode and playing a role to destabilize the wave. The plasma beta is greater than one so these effects satisfy the high plasma beta regimes.</p><p>On further increasing the magnetic field, unstable regions are obtained and at <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x68.png" xlink:type="simple"/></inline-formula> and A = 0.9, the wave becomes totally unstable as in <xref ref-type="fig" rid="fig4">Figure 4</xref>.</p><p>The growth rate shows that parallel streaming responsible to grow the wave. The complex part of the dispersion relation tells that the wave is growing in the gaps. <xref ref-type="fig" rid="fig5">Figure 5</xref> shows the growing parts of the first two gaps.</p><p>The O-mode instability divides in two branches for complex<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x69.png" xlink:type="simple"/></inline-formula>. First branch is oscillatory when <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x69.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x70.png" xlink:type="simple"/></inline-formula> and second branch is aperiodic or purely growing when <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x69.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x70.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x71.png" xlink:type="simple"/></inline-formula> as in <xref ref-type="fig" rid="fig6">Figure 6</xref>. Further increasing the value of parallel streaming, the aperiodic branch is obtained. For oscillatory branch the magnetic field plays a role to destabilize the wave and increase the growth of the wave.</p><p>For second branch, which is aperiodic or purely growing, the trend totally reverses . The growth rate increases with the decreasing value of A. The result proves that the anisotropy stabilizes the purely growing part. <xref ref-type="fig" rid="fig7">Figure 7</xref> shows the increasing growth rate of aperiodic mode with decreasing value of A. The noticeable thing is that this part also satisfies the condition of the firehose instability i.e., <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x72.png" xlink:type="simple"/></inline-formula>and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x72.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x73.png" xlink:type="simple"/></inline-formula>. The purely growing wave is also called non-propagating firehose instability [<xref ref-type="bibr" rid="scirp.67526-ref25">25</xref>] The study of variation of anisotropy tells us that with the increasing value of A, the growth rate is also increases that means by increasing value of A destabilizes the wave. When <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x72.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x73.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x74.png" xlink:type="simple"/></inline-formula> and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x72.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x73.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x74.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x75.png" xlink:type="simple"/></inline-formula>, then wave becomes more unstable this result proves that O-mode instability satisfies the condition of the firehose instability.</p><p>For number of harmonics, the X-mode is also unstable but for this mode perpendicular temperature is</p><fig id="fig6"  position="float"><label><xref ref-type="fig" rid="fig6">Figure 6</xref></label><caption><title> A = 0.09<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x77.png" xlink:type="simple"/></inline-formula></title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/9-7502712x76.png"/></fig><fig id="fig7"  position="float"><label><xref ref-type="fig" rid="fig7">Figure 7</xref></label><caption><title> Growth rate for aperiodic branch</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/9-7502712x78.png"/></fig><p>dominating. The wave becomes unstable for larger value of A. In <xref ref-type="fig" rid="fig8">Figure 8</xref> when A = 6.5, the wave is stable. But after that when A = 6.96, the harmonics overlap each other and wave starts to be unstable as in <xref ref-type="fig" rid="fig9">Figure 9</xref>.</p><fig id="fig8"  position="float"><label><xref ref-type="fig" rid="fig8">Figure 8</xref></label><caption><title> A = 6.5 <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x80.png" xlink:type="simple"/></inline-formula> <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x80.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x81.png" xlink:type="simple"/></inline-formula></title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/9-7502712x79.png"/></fig><fig id="fig9"  position="float"><label><xref ref-type="fig" rid="fig9">Figure 9</xref></label><caption><title> A = 6.69 <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x83.png" xlink:type="simple"/></inline-formula> <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x83.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x84.png" xlink:type="simple"/></inline-formula></title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/9-7502712x82.png"/></fig><p>On further increasing the value of A the mode becomes more unstable as in <xref ref-type="fig" rid="fig1">Figure 1</xref>0 the value of A is 8 and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x85.png" xlink:type="simple"/></inline-formula>. The solid curves show the real part of the wave and dashed curves show the growth of the said wave. The increasing value of A shows that in X-mode instability perpendicular streaming is dominating. At A = 11, it becomes totally unstable <xref ref-type="fig" rid="fig1">Figure 1</xref>1.</p><p><xref ref-type="fig" rid="fig1">Figure 1</xref>2 discusses that the growth rate increases with the increasing value of anisotropy. Anisotropy de- stabilizes the X-mode, the X-mode follows the same trend as that of the O-mode.</p><fig id="fig10"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref>0</label><caption><title> A = 8.0 <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x87.png" xlink:type="simple"/></inline-formula> <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x87.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x88.png" xlink:type="simple"/></inline-formula></title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/9-7502712x86.png"/></fig><fig id="fig11"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref>1</label><caption><title> A = 11 <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x90.png" xlink:type="simple"/></inline-formula> <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x90.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x91.png" xlink:type="simple"/></inline-formula></title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/9-7502712x89.png"/></fig><fig id="fig12"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref>2</label><caption><title> Growth rate of X-mode for different values of A</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/9-7502712x92.png"/></fig></sec><sec id="s4"><title>4. Conclusion</title><p>O-mode instability, for principle harmonic, depends upon the magnetic field even it is weaker. The instability generates due to temperature anisotropy and free energy of anisotropy converted in the magnetic induction which is the reason of growing wave. The growth rate varies directly with the value of ratio of anisotropy. Here we have calculated the marginal threshold condition in form of plasma parameters A and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x93.png" xlink:type="simple"/></inline-formula> for principle harmonic. For higher harmonics, oscillatory branch satisfies the statement and purely growing part inverts the condition. It varies inversely with the anisotropy. The oscillatory and purely growing mode both satisfies the conditions of firehose instability i.e., <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x93.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x94.png" xlink:type="simple"/></inline-formula>and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x93.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x94.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x95.png" xlink:type="simple"/></inline-formula> as <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x93.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x94.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x95.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x96.png" xlink:type="simple"/></inline-formula> in O-mode. The stability analysis of the X-mode tells that perpendicular temperature is dominating. The mode is unstable for <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x93.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x94.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x95.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x96.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x97.png" xlink:type="simple"/></inline-formula> according to geometry of the X-mode that is <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x93.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x94.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x95.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x96.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x97.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-7502712x98.png" xlink:type="simple"/></inline-formula> . Coupling of these two modes converts them into the Bernstein mode which is responsible of heating effects in tokamak. The O-X conversion is the method of achieving the Bernstein mode.</p></sec><sec id="s5"><title>Acknowledgements</title><p>Authors are thankful to the Department of Physics, FC College (A Chartered University) for financial assistance.</p></sec><sec id="s6"><title>Cite this paper</title><p>N. Noreen,S. Zaheer,H. A. Shah, (2016) Stability Analysis of Electromagnetic Ordinary and Extraordinary Modes. Journal of Modern Physics,07,1120-1131. doi: 10.4236/jmp.2016.710101</p></sec></body><back><ref-list><title>References</title><ref id="scirp.67526-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Ono, M., et al., (2004) Nuclear Fusion, 44, 452. http://dx.doi.org/10.1088/0029-5515/44/3/011</mixed-citation></ref><ref id="scirp.67526-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Akers, R.-J., et al., (2002) Physics of Plasmas, 9, 3919. http://dx.doi.org/10.1063/1.1490928</mixed-citation></ref><ref id="scirp.67526-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Whaley, D.R., et al. (1992) Nuclear Fusion, 32, 757. http://dx.doi.org/10.1088/0029-5515/32/5/I04</mixed-citation></ref><ref id="scirp.67526-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Hamasaki, S. (1968) Physics of Fluids, 11, 12.</mixed-citation></ref><ref id="scirp.67526-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Hamasaki, S. (1968) Physics of Fluids, 6, 11.</mixed-citation></ref><ref id="scirp.67526-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Lee, K.F. (1969) Physics Review, 1, 181.</mixed-citation></ref><ref id="scirp.67526-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Bornatici, M. and Lee, K.F. (1971) Physics Fluids, 13, 42.</mixed-citation></ref><ref id="scirp.67526-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Shivamoggi, B.K. (1982) Astrophysics and Space Science, 82, 481-483. http://dx.doi.org/10.1007/BF00651455</mixed-citation></ref><ref id="scirp.67526-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Ibschar, D., Lazar, M. and Schlickeiser, R. (2012) Physics of Plasmas, 19, 072116. http://dx.doi.org/10.1063/1.4736992</mixed-citation></ref><ref id="scirp.67526-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Iqbal, Z., Hussain, A., Murtaza, G. and Tsintsadze, N.L. (2014) Physics of Plasmas, 21, 032128. http://dx.doi.org/10.1063/1.4870007</mixed-citation></ref><ref id="scirp.67526-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">Hadi, F., Yoon, P.H. and Qamar, A. (2015) Physics of Plasmas, 22, 022112. http://dx.doi.org/10.1063/1.4907657</mixed-citation></ref><ref id="scirp.67526-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">Lazar, M., Schlickeiser, R., Poedts, S., Stockem, A. and Van, S. (2014) http://arxiv.org/abs/1411.1508</mixed-citation></ref><ref id="scirp.67526-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">Vafin, S., Schlickeiser, R. and Yoon, P.H. (2014) Physics of Plasmas, 21, 104504. http://dx.doi.org/10.1063/1.4897373</mixed-citation></ref><ref id="scirp.67526-ref14"><label>14</label><mixed-citation publication-type="other" xlink:type="simple">Farrell, W.M. (2001) Journal of Geophysical Research: Space Physics, 106, 15701-15709. http://dx.doi.org/10.1029/2000JA000156</mixed-citation></ref><ref id="scirp.67526-ref15"><label>15</label><mixed-citation publication-type="other" xlink:type="simple">Colpitts, C.A., Cattell, C.A., Kozyra, J.U. and Parrot, M. (2012) Journal of Geophysical Research: Space Physics, 117. http://dx.doi.org/10.1029/2011JA017329</mixed-citation></ref><ref id="scirp.67526-ref16"><label>16</label><mixed-citation publication-type="other" xlink:type="simple">LaBelle, J., Ruppert, D.R. and Treumann, R.A. (1999) Journal of Geophysical Research, 104, 293-303. http://dx.doi.org/10.1029/1998JA900050</mixed-citation></ref><ref id="scirp.67526-ref17"><label>17</label><mixed-citation publication-type="other" xlink:type="simple">Podoba, Y.Y., Laqua, H.P., Warr, G.B., Schubert, M., Otte, M., Marsen, S. and Wagner, F. (2007) Physical Review Letters, 98, 255003. http://dx.doi.org/10.1103/PhysRevLett.98.255003</mixed-citation></ref><ref id="scirp.67526-ref18"><label>18</label><mixed-citation publication-type="other" xlink:type="simple">Cairns, R.A. and Lashmore-Daviesa, C.N. (2000) Physics of Plasmas, 7, 4126-4134.</mixed-citation></ref><ref id="scirp.67526-ref19"><label>19</label><mixed-citation publication-type="other" xlink:type="simple">Ram, A.K., Bers, A. and Lashmore-Davies, C.N. (2002) Physics of Plasmas, 9, 409-418. http://dx.doi.org/10.1063/1.1429634</mixed-citation></ref><ref id="scirp.67526-ref20"><label>20</label><mixed-citation publication-type="other" xlink:type="simple">Sodha, M.S., Sharma, R.P., Maheshwari, K.P. and Kaushik, S.C. (1977) Plasma Physics, 20, 585. http://dx.doi.org/10.1088/0032-1028/20/6/009</mixed-citation></ref><ref id="scirp.67526-ref21"><label>21</label><mixed-citation publication-type="other" xlink:type="simple">Bashir, M.F. and Murtaza, G. (2012) Brazilian Journal of Physics, 42, 487-504.</mixed-citation></ref><ref id="scirp.67526-ref22"><label>22</label><mixed-citation publication-type="other" xlink:type="simple">Yoon, P.H. and Davidson, R.C. (1987) Physical Review A, 35, 2718-2721.</mixed-citation></ref><ref id="scirp.67526-ref23"><label>23</label><mixed-citation publication-type="other" xlink:type="simple">Bashir, M.F., Noreen, N., Murtaza, G. and Yoon, P.H. (2014) Plasma Physics and Controlled Fusion, 56, Article ID: 055009.</mixed-citation></ref><ref id="scirp.67526-ref24"><label>24</label><mixed-citation publication-type="other" xlink:type="simple">Ichimaru, S. (1973) Basic Principles of Plasma Physics: A Statistical Approach. Addison-Wesley Publishing Company, New York, 104.</mixed-citation></ref><ref id="scirp.67526-ref25"><label>25</label><mixed-citation publication-type="other" xlink:type="simple">Lazar, M., Poedts, S., Schlickeiser, R. and Ibscher, D. (2014) Solar Physics, 289, 369-378. http://arxiv.org/abs/1307.0768 </mixed-citation></ref></ref-list></back></article>