Journal of AstrophysicsVolume 2013 (2013), Article ID 838534, 7 pageshttp://dx.doi.org/10.1155/2013/838534

Research Article

Ion-Acoustic Instabilities in a Multi-Ion Plasma

Noble P. Abraham,¹ Sijo Sebastian,¹ G. Sreekala,¹ R. Jayapal,¹ C. P. Anilkumar,² and Venugopal Chandu¹

¹School of Pure & Applied Physics, Mahatma Gandhi University, Priyadarshini Hills, Kottayam, Kerala 686 560, India²Equatorial Geophysical Research Laboratory, Indian Institute of Geomagnetism, Krishnapuram, Tirunelveli, Tamil Nadu 627 011, India

Received 1 April 2013; Accepted 23 June 2013

Academic Editors: M. S. Dimitrijevic, A. Meli, and S. Naik

Copyright © 2013 Noble P. Abraham et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

We have, in this paper, studied the stability of the ion-acoustic wave in a plasma composed of hydrogen, positively and negatively charged oxygen ions, and electrons, which approximates very well the plasma environment around a comet. Modelling each cometary component (, , and ) by a ring distribution, we find that ion-acoustic waves can be generated at frequencies comparable to the hydrogen ion plasma frequency. The dispersion relation has been solved both analytically and numerically. We find that the ratio of the ring speed () to the thermal spread (ts) modifies the dispersion characteristics of the ion-acoustic wave. The contrasting behaviour of the phase velocity of the ion-acoustic wave in the presence of ions for ts (and vice versa) can be used to detect the presence of negatively charged oxygen ions and also their thermalization.

1. Introduction

Low-frequency electrostatic or longitudinal ion density waves are one of the most fundamental of oscillations in a plasma [1, 2]. In the long-wavelength limit, the ions provide the inertia with the electrons as the source of the restoring force [1]. Ion-acoustic waves also exhibit strong nonlinear properties and are highly Landau damped unless , where and are, respectively, the ion and electron temperatures [3–5]. These waves have been observed in both space and laboratory plasmas; they have thus been extensively studied in many types of high-temperature laboratory plasmas [4, 6]. The waves have been invoked to explain wave characteristics observed in Earth’s ionosphere [7] and transport in the solar wind, corona, chromosphere [8], and comets [9].

In general a cometary environment contains new born hydrogen and heavier ions, with relative densities depending on the distance from the nucleus. Previous studies have concentrated on positively charged oxygen as the heavier ion species [10]. However, Giotto’s observations of the inner coma of comet Halley showed that a new component, namely, negatively charged cometary ions was present, in addition to the usual thermal electrons and ions, fast cometary pickup ions, and so forth, [11]. These negative ions were observed in three broad mass peaks at 7–19, 22–65, and 85–110 amu with being identified unambiguously [11].

A popular model of a cometary environment is the solar wind plasma environment permeated by dilute, drifting ring distribution of electrons and ions with finite thermal spreads [10]. Instabilities driven by an electron velocity ring distributions have been studied by many authors [12–14]. However, ion ring distributions are more important because of the greater amount of free energy available [15].

Instabilities driven by ion ring distributions have also been studied by a number of authors: close to the ion cyclotron frequency, electrostatic ion cyclotron waves propagating either perpendicularly or nearly perpendicularly to the magnetic field can be excited [16, 17]. At still higher frequencies, the magnetic effects on the ions can be neglected and lower-hybrid instabilities driven by the ring ions can occur [17–21]. And at even higher frequencies, the electrons too become unmagnetized and ion-acoustic-like instabilities result.

The frequency of the ion-acoustic wave is comparable with the ion plasma frequency and propagates parallel to the magnetic field. Generally, a combination of warm electrons and cold ions () with the electrons drifting relative to the ions is the condition required to excite the ion-acoustic instability.

The ion-acoustic wave is one of the more easily observed waves in the plasma environments of comets. For example, ion-acoustic waves in the frequency range of 1.0–1.5 kHz were observed by the ICE spacecraft sent to observe the comet Giacobini-Zinner [9, Figure 3]. Again ion sound waves, with a frequency slightly less than 1 kHz, were detected by the spacecraft Sakigake which observed comet Halley [22, Figure 6].

We have, therefore, studied the stability of the ion-acoustic wave in a five-component plasma: solar wind protons, electrons, cometary hydrogen, and positively and negatively charged oxygen ions. We find that the ratio of the ring speed () to the thermal spread () affects the dispersion characteristics of the ion-acoustic wave. The phase velocity of the ion-acoustic wave depends sensitively on this ratio in the presence of ions. This variation for (and vice versa) is proposed as a tool for detecting the presence of ions and also their thermalization.

2. The Dispersion Relation

As stated above, we intend to study the stability of the ion-acoustic wave in a five-component plasma. The five components are solar wind hydrogen and electrons and the ions of cometary origin picked up by the solar wind. These cometary ions are hydrogen and positively and negatively charged oxygen ions. The solar wind components are modelled by the Maxwellian distributions while the cometary ions are described by ring distributions.

The contributions to the dispersion relation of the plasma components described by the Maxwellian distributions are well known [23]; the same when they are described by ring distributions are given in [15]. The plasma under consideration contains components modelled by both these forms. Hence, combining these contributions, we can write down the dispersion relation for waves of frequency and wave vector , as

In (1), and denote, respectively, the plasma frequency and the thermal velocity for particles of species . Also , , , and denote, respectively, the density, charge, mass, and temperature of species . () denotes the plasma frequency for the ring ions of species while and are, respectively, the ring ion speed and thermal spread for these ions. The plasma dispersion function , arising from the integration [24], has arguments of (for hydrogen ions of solar wind origin) and (for electrons). is the drift velocity of the electrons parallel to the magnetic field.

The charge neutrality condition for the plasma under consideration yields

with , , and , where , in general, denotes the density of each species.

We derive expressions for the growth rate of the ion-acoustic wave under two limiting conditions, namely, that of a large thermal spreads () and small thermal spreads (). The contributions of the solar wind hydrogen ions and electrons are the same for both cases; the contributions of the ring ions vary depending on their thermal spreads.

We use the asymptotic expansion of the plasma dispersion function for the solar wind hydrogen ions and the power series expansion of for the solar wind electrons. Substituting the appropriate expansions in (1), we finally arrive at the dispersion relation for ion-acoustic waves in an electron-ion plasma with three species of ring ions, as

In (3), and is the Debye length.

2.1. Case (i): Large Thermal Spreads

In this subsection we derive expressions for the real frequency and growth/damping rate for ion-acoustic waves when the thermal spreads of the ring ions are large. As a simplifying assumption we also let . With the assumption of large thermal spreads , (3) can be simplified to yield an expression for the real frequency as

and the growth/damping rate as

where .

In the above, is the speed of sound defined by , while

2.2. Case (ii): Small Thermal Spreads

We next consider expressions for the real frequency and the growth/damping rate when the thermal spread is small; that is, . These can be extracted from (3) with the additional assumption ; the final expressions are

Thus (4) and (7) are expressions for the frequencies of the ion-acoustic wave under the two limiting cases; (5) and (8) are the corresponding expressions for the growth/damping rate.

3. Results and Discussion

As a check on our results, we note from (4) and (7) the expressions for real frequency that they reduce to the corresponding ones in an electron-ion plasma when the ring ions are absent [25]. Similarly, the expressions for the growth/damping rate (5) and (8) also reduce to that in an electron-ion plasma when the ring ions are absent [25]. However what is interesting from (4) and (7) is that the real frequency is either greater (for ) or smaller (for ) than in a single ion plasma (in the absence of ring ions). The electron drift velocity, parallel to the magnetic field, is the source of energy for the instability. Thus, for a given wave vector and ion densities, a greater drift velocity of the electrons is required to drive the wave unstable in the case where the thermal spread of the ring ions is large.

We use parameters similar to that in [10]. The solar wind hydrogen density , while the hydrogen and electron temperatures are , . The densities of the ring ions are , , and .

Ion-acoustic waves were observed in the frequency range from 1.0 to 1.5 kHz by the ICE spacecraft at the comet Giocobini-Zinner [9]. Thus to calculate the wave vector , a frequency of  kHz was used in the cold plasma dispersion relation and the wave vector was calculated. was then varied around this value in subsequent calculations.

Figure 1 is a plot of the frequency versus wave vector , with the frequency being obtained by solving (4) and (7). Curve (a) has a drift speed and a thermal spread of () while for curve (b) , (). As is evident from the figure the waves are at a higher frequency for . Since the calculated frequencies are in better agreement with the observed range of frequencies for , it is safe to assume that thermalization of the ring ions has not yet occurred. The pickup of new born ions by the solar wind has been described by a three-step process: formation of a ring beam distribution, pitch angle scattering of this initial distribution into a shell, and a slower velocity diffusion that spreads out the shell [26].

Figure 1: Plot of frequency versus wave vector . Curve (a) is for ; drift speed and thermal spread ; curve (b) is for ; and .

We first consider the growth/damping rates when the thermal spread is large.

Figure 2 is a plot of the growth rate (5) versus wave vector for , , , , and an electron drift speed . Curve (a) is for a single ion plasma without ring distributions. Curve (b) has only one type of ring distribution, namely, the ring. We find that the growth rate is enhanced. Curve (c) is for the ring replacing the ring—there is a slight decrease in the growth rate. When this single ring is the ring, the growth rate is the lowest (curve (d)). The growth rate is also shown when all species of rings are present (curve (e)). The ion-acoustic wave is thus driven unstable by the electrons drifting parallel to the magnetic field (addition of negative ions reduces the density of electrons required to maintain charge neutrality); the ring ions add to the instability, with the lighter ion hydrogen ring being the most effective in enhancing the instability of the wave.

Figure 2: Plot of the growth rate (5) versus wave vector for , , , , and an electron drift speed .  Curve (a) is for a single ion plasma without any ring ions. Curve (b) has a ring, curve (c) has an ring, curve (d) has an ring, and for curve (e) all species of ring ions are present.

Figures 3 and 4 depict the variation of the growth rate (5) versus () (Figure 3) or versus () (Figure 4), with the temperatures and electron drift speed being the same as in Figure 2. In both figures curve (a) represents the growth rate in a single ion plasma (ring ions are not present); curve (b) when the ring alone is present, curve (c) when the ring alone is present, and curve (d) when an ring alone is present. Curve (e) is for the situation when all species of ring ions are present. While the general nature of the curves is similar to that in Figure 2, the heavier ion rings extend the range of and over which the wave is unstable as compared to the lighter ion ring alone being present.

Figure 3: Plot of the growth rate (5) versus (,  ). Other relevant parameters are the same as in Figure 2. Curve (a) is for a single ion plasma without any ring ions. Curve (b) is for a ring, curve (c) has an ring, and curve (d) is for an ring. Curve (e) depicts the case when all species of ring ions are present.

Figure 4: Plot of the growth rate (5) versus (, ). Other relevant parameters are the same as in Figure 2. Curve (a) is for a single ion plasma without any ring ions. Curve (b) has a ring, curve (c) has an ring, curve (d) has an ring, and for curve (e) all species of ring ions are present.

We next consider the case of low thermal spreads ().

Figure 5 is a plot of the growth rate (8) versus wave vector for , , , , and an electron drift speed . Similar to Figure 2, curve (a) is for a single ion plasma without any ring ions. Curve (b) is for the ring and has the lowest growth rate. Curve (c) is for the ring replacing the ring—there is an enhancement in the growth rate, than the single ion plasma. When this single ring is the ring, the growth rate is slightly lowered (curve (d)). Curve (e) shows the growth rate when all the species of rings are present. The ion-acoustic wave is thus again driven unstable by the electrons drifting parallel to the magnetic field (addition of negative ions reduces the density of electrons required to maintain charge neutrality); the ring ions generally reduce the growth with the lighter ion hydrogen ring being the most effective in damping the wave.

Figure 5: Plot of the growth rate (8) versus wave vector for , , , , and an electron drift speed . Curve (a) is for a single ion plasma without any ring ions. Curve (b) has a ring, curve (c) has an ring, and curve (d) has an ring. Curve (e) has all species of ring ions.

Figures 6 and 7 depict the variation of the growth rate (8) versus () (Figure 6) or versus () (Figure 7), with the temperatures and drift speed being the same as in Figure 5. In both figures, labelling of the curves is similar to that in Figures 3 and 4. While the general nature of the curves is similar to that in Figure 5, the heavier ion rings have much larger ranges of and over which the wave is unstable.

Figure 6: Plot of the growth rate (8) versus (, ). Other relevant parameters are the same as in Figure 5. Curve (a) is for a single ion plasma without any ring ions. Curve (b) has a ring, curve (c) has an ring, curve (d) has an ring, and curve (e) has all species of ring ions.

Figure 7: Plot of the growth rate (8) versus (, ). Other relevant parameters are the same as in Figure 5. Curve (a) is for a single ion plasma without any ring distributions. Curve (b) has a ring, curve (c) has an ring, curve (d) has an ring, and curve (e) has all species of ring ions.

In a recent study Misra et al. [27] investigated the propagation characteristics of the dust ion-acoustic (DIA) waves in a collisional, negative ion plasma with immobile, charged dust grains. In the linear regime, they found two modes—the “fast” and “slow” modes. They also found that the Landau damping effect on the fast mode was negligible while the slow mode was stable.

Our dispersion diagram (Figure 1) also reveals the existence of two modes—a “slow” mode (curve (a) of Figure 1 with ) and a “fast” mode (curve (b), Figure 1 with ). Also as mentioned in Section 1, ion-acoustic waves were observed in the frequency range of 1.0–1.5 kHz at comet Giacobini-Zinner [9] and at a frequency slightly lesser than 1.0 kHz at comet Halley [22]. The waves are also expected to saturate at frequencies below the local ion plasma frequency [28]. A comparison of the growth rate for the fast (Figures 2, 3, and 4) and the slow (Figures 5, 6, and 7) modes reveals that the growth rate of the fast mode is much larger when compared to that of the slow mode. Thus, while the slow mode in [27] was freely propagating, the slow mode in the plasma under consideration is weakly unstable. It may be noted here that this study uses the full kinetic treatment and hence can also be considered as a generalisation of the fluid treatment in [27].

Finally Figure 8 depicts the variation of the real frequency versus the wave vector . The panel at the top is for the fast wave (, , ); curve (a) is for the case when there are no heavier ions (, ) while curve (b) depicts the situation where negatively charged oxygen is also added (, , and ). We find that the addition of increases the frequency and hence the phase velocity of the ion-acoustic wave. The lower panel is for the slow wave (, , )—curve (c) is for the case where only the lighter ions are present (, ) and curve (d) depicts the situation when is also added (, , and ). The decrease in the frequency and thereby the phase velocity is distinctly evident from the figure.

Figure 8: Plot of frequency versus wave vector . The panel at the top is for the fast wave (; , ) of Figure 1: curve (a) is for the case when there are no heavier ions (, ) while curve (b) depicts the situation when negatively charged oxygen is also added (, , and ). The lower panel is for the slow wave (; , ) of Figure 1: curve (c) is for the case where only the lighter ions are present (,  ) and curve (d) depicts the situation when is also added (, , and ). Other relevant parameters are the same as in Figure 1.

Variation of the phase velocity of ion-acoustic waves has been proposed as a diagnostic tool for the detection of charged dust grains in laboratory plasmas [28, 29]. We thus propose to extend this idea—the contrasting behaviour of the phase velocity of the ion-acoustic wave in the presence of ions for (and vice versa) can be used not only to detect the presence of negatively charged oxygen ions but also their thermalization as well.

4. Conclusions

We have studied the stability of the ion-acoustic waves in a five-component plasma, which approximates very well the plasma environment around a comet. The five components are solar wind hydrogen and electrons and the ions of cometary origin. The cometary ions are hydrogen and positively and negatively charged oxygen ions. The solar wind components are modelled by the Maxwellian distributions while the cometary ions are described by ring distributions. We find that the ion-acoustic waves can be generated at frequencies comparable to the observed frequency. The ratio of the ring speed to its thermal spread modifies the dispersion characteristics of the ion-acoustic wave. When the thermal spread is large compared to the ring speed (), the ring ions enhance the instability of the ion-acoustic wave. The growth rate is lowered when the thermal spread is smaller in comparison with the ring speed (). We also find that heavier ring ions extend the range of and over which the wave is unstable when compared to the lighter ring ions. The contrasting behaviour of the phase velocity of the ion-acoustic wave in the presence of ions for (and vice versa) can be used to detect the presence of negatively charged oxygen ions and also their thermalization.

Acknowledgments

The authors thank the referees for their valuable comments. Financial assistance from the University Grants Commission (SAP) and Department of Science and Technology (FIST and PURSE Programs) is gratefully acknowledged.

References

J. Castro, P. McQuillen, and T. C. Killian, “Ion acoustic waves in ultracold neutral plasmas,” Physical Review Letters, vol. 105, no. 6, Article ID 065004, 2010.

L. Tonks and I. Langmuir, “Oscillations in ionized gases,” Physical Review, vol. 33, no. 2, pp. 195–210, 1929.

T. H. Stix, Waves in Plasmas, American Institute of Physics, New York, NY, USA, 2nd edition, 1992.

Y. Nakamura, H. Bailung, and P. K. Shukla, “Observation of ion-acoustic shocks in a dusty plasma,” Physical Review Letters, vol. 83, no. 8, pp. 1602–1605, 1999.

Z. Liu, L. Liu, and J. Du, “A nonextensive approach for the instability of current-driven ion-acoustic waves in space plasmas,” Physics of Plasmas, vol. 16, no. 7, Article ID 072111, 5 pages, 2009.

M. Yamada and M. Raether, “Saturation of the ion-acoustic instability in a weakly ionized plasma,” Physical Review Letters, vol. 32, no. 3, pp. 99–102, 1974.

M. E. Koepke, “Contributions of Q-machine experiments to understanding auroral particle acceleration processes,” Physics of Plasmas, vol. 9, no. 5, pp. 2420–2427, 2002.

S. R. Cranmer, A. A. Van Ballegooijen, and R. J. Edgar, “Self-consistent coronal heating and solar wind acceleration from anisotropic magnetohydrodynamic turbulence,” Astrophysical Journal, Supplement Series, vol. 171, no. 2, pp. 520–551, 2007.

F. L. Scarf, F. V. Coroniti, C. F. Kennel, D. A. Gurnett, W.-H. Ip, and E. J. Smith, “Plasma wave observations at comet Giacobini-Zinner,” Science, vol. 232, no. 4748, pp. 377–381, 1986.

A. L. Brinca and B. T. Tsurutani, “Unusual characteristics of electromagnetic waves excited by cometary newborn ions with large perpendicular energies,” Astronomy & Astrophysics, vol. 187, no. 1-2, pp. 311–319, 1987.

P. Chaizy, H. Rème, J. A. Sauvaud et al., “Negative ions in the coma of comet Halley,” Nature, vol. 349, no. 6308, pp. 393–396, 1991.

T. J. Tataronis and F. Crawford, “Cyclotron harmonic wave propagation and instabilities,” Journal of Plasma Physics, vol. 4, no. 2, pp. 231–264, 1970.

M. Ashour-Abdalla and C. F. Kennel, “Nonconvective and convective electron cyclotron harmonic instabilities,” Journal of Geophysical Research, vol. 83, p. 1531, 1978.

P. Sprangle, J. L. Vomvoridis, and W. M. Manheimer, “A classical electron cyclotron quasioptical maser,” Applied Physics Letters, vol. 38, no. 5, pp. 310–313, 1981.

K. Akimoto, K. Papadopoulos, and D. Winske, “Ion-acoustic instabilities driven by an ion velocity ring,” Journal of Plasma Physics, vol. 34, no. 3, pp. 467–479, 1985.

J. A. Byers and M. Grewal, “Perpendicularly propagating plasma cyclotron instabilities simulated with a one-dimensional computer model,” Physics of Fluids, vol. 13, no. 7, pp. 1819–1830, 1970.

J. K. Lee and C. K. Birdsall, “Velocity space ring-plasma instability, magnetized, part I: theory,” Physics of Fluids, vol. 22, no. 7, pp. 1306–1314, 1979.

S. Seiler, M. Yamada, and H. Ikezi, “Lower hybrid instability driven by a spiraling ion beam,” Physical Review Letters, vol. 37, no. 11, pp. 700–703, 1976.

H. E. Mynick, M. J. Gerver, and C. K. Birdsall, “Stability regions and growth rates for a two-ion component plasma, unmagnetized,” Physics of Fluids, vol. 20, no. 4, pp. 606–612, 1977.

C. Cattel and M. Hudson, “Flute mode waves near ω LH excited by ion rings in velocity space,” Geophysical Research Letters, vol. 9, no. 10, pp. 1167–1170, 1982.

K. Akimoto, K. Papadopoulos, and D. Winske, “Lower-hybrid instabilities driven by an ion velocity ring,” Journal of Plasma Physics, vol. 34, no. 3, pp. 445–465, 1985.

F. Scarf, “Plasma wave observations at comets Giacobini-Zinner and Halley,” in Plasma Waves and Instabilities at Comets and in Magnetospheres, B. T. Tsurutani and H. Oya, Eds., American Geophysical Union, Washington, DC, USA, 1989.

R. C. Davidson, “Kinetic waves and instabilities in a uniform plasma,” in Basic Plasma Physics, A. Galeev and R. Sudan, Eds., North-Holland, New York, NY, USA, 1989.

B. D. Fried and S. D. Conte, The Plasma Dispersion Function, Academic Press, New York, NY, USA, 1961.

D. A. Gurnnet and A. Bhattacharjee, Introduction to Plasma Physics: with Space and Laboratory Applications, Cambridge University Press, Cambridge, UK, 1st edition, 2005.

J. D. Gaffey Jr, D. Winske, and C. S. Wu, “Time scales for formation and spreading of velocity shells of pickup ions in the solar wind,” Journal of Geophysical Research, vol. 93, no. A6, pp. 5470–5486, 1988.

A. P. Misra, N. C. Adhikary, and P. K. Shukla, “Ion-acoustic solitary waves and shocks in a collisional dusty negative-ion plasma,” Physical Review E, vol. 86, no. 5, Article ID 056406, 10 pages, 2012.

S. H. Kim and R. L. Merlino, “Charging of dust grains in a plasma with negative ions,” Physics of Plasmas, vol. 13, no. 5, Article ID 052118, 7 pages, 2006.

M. Rosenberg and R. L. Merlino, “Ion-acoustic instability in a dusty negative ion plasma,” Planetary and Space Science, vol. 55, no. 10, pp. 1464–1469, 2007.