Gravity Induced Resonant Emission

Gravity Induced Resonant Emission

Antony Soosaleon Sunitha Amborse
1. The Astronomical Observatory, Department of Physics University of Kerala, India - 695033
2. SPAP, M.G.University, Kottayam, Kerala, India-686560,
3. Center for Fundamental Research and Computational Sciences,
Thiruvanathapuram, Kerala, India - 695043.

The gravitational drift of ions relative to the electrons induces two type of waves in magnetized plasma; ion acoustic (ia) waves and lower hybrid (lh) waves. The ia waves induced by the gravity are damped by electromagnetic (em) waves, which leads to the formation of lh waves. For higher wave vector, these lh waves results in to the resonant absorption & re-emission of em waves, called as gravity induced resonant emission (gire). A general formula has been derived for gire frequency, which is given by ; c, velocity of light; g, gravity; e, charge; B, magnetic field; M & m and & are masses and number densities of ion and electron respectively.
Key Words: Gravity, Resonant Emission, Magnetic fields, MHD theory, Ion Accoustic waves, Lower Hybrid Waves, Landau Damping, Instabilities, Wave Particle Interaction, Solar Coronal Heating.

1 Introduction

Gravity is inherent, any fundamental theory including the effect of gravity would be more realistic. The effect of gravity in magnetized plasma had been discussed by Chen[6] and Treumann & Baumjohann[19] as Rayleigh Taylor Instability, which did not yield any result. The problem with their study was that they considered the plasma as cold; for any instability, a driving force alone is not sufficient but it needs energy to support. Even otherwise too, plasma is an ionized state of matter exists only at high temperature and does not sustain with out magnetic field and therefore these two quantities are inevitable, by these reasons the plasma considered for this study is warm and magnetized. The effect of these two parameters is so volatile, hence they are set as finite for simplicity. The ionized particles in magnetic field undergo cyclotron motion, which is controlled by several factors like mass () & mass number , charge state (Ze) and the effective field strength (B) (cyclotron frequency, ). The effective magnetic field on each ions is unique; even if we assume all the ions in a plasma experience the same magnetic field, the masses are different, the charge state often changes and the direction of the cyclotron motion is different for electrons and ions. With all these variations including the temperature, a magnetized warm plasma may possess a variety of localized fields and such fields can cause localized currents; sometimes certain field resonates with the other localized fields leads to several observable phenomenon. In an ultra vacuum condition (Solar Corona; [2]), the currents by localized fields are not possible, rather the Coulomb’s interaction will be dominant, by the fact that the fields will be effective to very long distances. Here is a situation arises, that the gravitational drift of heavier particles relative to the electrons induces an electric field; this gravitationally induced electric field (gief) couples with the magnetic field, causes two physical phenomenon. 1. It generates ion acoustic (ia) waves and these waves are damped by the perpendicular propagation electromagnetic (em) waves leads to the formation of lower hybrid waves 2. The gravity induced lower hybrid (lh) waves leads to the resonant absorption and re-emission of em waves and in both the cases, the Coulomb’s force is playing a lead role. The gief may not be effective to high plasmas (, is the ratio of particle pressure to the magnetic field pressure; e.g., plasmas in Fusion devises, stellar interior and photosphere) due to the high collisions.
This study is carried out by deriving a dispersion relation (relation between the angular frequency and the wave vector ) for the em wave propagation in a warm magnetized gravitational plasma. We are considering only the perpendicular propagation ( to B), by the reason that the Lorentz force is a normal force and hence the coupling is possible only in the perpendicular direction.

2 Theoretical Description

2.1 Gravitational Plasma

The purpose of this study is to find the effect of gravitational force over the magnetized plasma, therefore we need only ionized particles, for that we consider the plasma as fully ionized (coronal plasma). As briefed in the introduction, the ions are heavier, drift more relative to the electrons. The drift velocity of ions due to the gravity is deduced from the general definition, . For gravitational drift; is replaced by , which gives where , mass; , charge and g, acceleration due to gravity respectively.

2.2 The Basic Equations Used

The dispersion relation is derived by using magneto hydro dynamical theory and the following equations are used. The momentum equation is given by


Here , number density; , velocity; , temperature; , the ratio of specific heats; E, electric field; K is the Boltzmann constant respectively. Since gief is effective only to the low plasmas, we have neglected the collision & viscosity terms. The second equation is the continuity equation which is given by


2.3 Procedure of the theory

  1. The dispersion relation is derived by solving the momentum equation by linearization method.[6]

  2. First we deduce a resultant perturbed momentum equation

  3. Then we introduce the em waves as perturbation in the momentum equation and derive perturbed velocities

  4. Substituting these velocities in the perturbed continuity equation, the dispersion relation is derived

Figure 1: Gravity induced electric field in unperturbed and perturbed magnetic field

2.4 Derivation of Dispersion Relation

Consider the gravitational force is acting in the x-axis, the wave is propagating in the y-axis and the magnetic field is applied in the z-axis as shown in the fig.1. The effect of gravity over the ions relative to the electrons causes two things; 1. a density gradient in the direction of x-axis and 2. it induces electric field due to in the -x axis. This induced field causes a drift in the absence of perturbation and this drift couples with the diamagnetic drifts leads to formation of ia oscillations. But this situation changes when there is a perturbation in the form of em waves, the direction of changes at every point and hence the direction of the drift also changes (see fig.2), which leads to the damping of ia waves and the formation of lh oscillation. For simplicity, we are considering only a two component plasma (electrons & ions), since the effect of gravity is over the density and the velocity (by the induced electric field), the theory is formulated by introducing the perturbation in the density and the velocity of the ionized particles. First, we start with the ion component of the plasma; the perturbed momentum equation for ions is obtained, by replacing with and with in the Eq.(1), we get,


The steady state momentum equation for ions is written as,


Multiplying Eq.(4) with and subtracting it from Eq.(3) and neglecting the second order terms, we get the resultant perturbed equation as


Here the resultant equation (Eq.(5)) contains two density gradient terms and . The is the perturbation in the density caused due to the propagation of em waves; while assuming the perturbation as , is replaced with . Regarding the term ; is the number density of ions, the plasma we consider is fully ionized and the gravitational drift is effective to all the ions, therefore we replace with , is the number density of ions drifted per unit length. Now the Eq.(5) becomes,


where and . Equating the x and y components of Eq.(6), we obtain the components of perturbed velocities for ions as


Similarly we can deduce the components of perturbed velocities for electrons as


Where , , is the thermal velocities of electron and ion respectively.

Figure 2: Model of Gravity Induced Ion Acoustic Oscillation

Introducing the velocity components in perturbed continuity equation, we can derive a general dispersion relation valid for all range of frequencies. But the corresponding dispersion relation is very cumulative with numerous terms and analyzing it for the entire range of em waves is very plausible. Therefore the dispersion relation is separated into two; em waves of frequencies and low frequencies . For the high frequencies, the gief is not effective, since all the terms containing the are vanishing out. This is an expected result, because is caused due to the drift of ions and therefore the possible interaction should be an ion cyclotron related one. For such a case, the maximum possible frequency will be the frequency of a coupled oscillation between the ions and electrons, i.e., lh frequency which is less than .
For deducing the dispersion relation in low frequency range; by setting the condition in the Eq.(7),(8),(9) & (10), we get the velocity components in the following form,


The perturbed equation of continuity for ions is given as


Substituting the values of and in the Eq.(15), we get


Similarly, from the perturbed continuity equation of electron, we get


Substituting the value of , from the Eq.(16) in to the Eq.(17), we can derive the dispersion relation as follows,


where ; ; , is the diamagnetic drift velocities; , is the ion acoustic velocity.

Figure 3: Two modes of propagation of EM waves, for coronal plasma.( , , )

3 Results & Discussion

The above dispersion relation is a quadratic equation in , but it is non linear in nature. There are two non linear terms in this relation, is in the coefficient of and is in the coefficient of . Generally a nonlinear term signifies the energy exchange process, therefore these two terms represent two physical processes related the energy exchange. We find both these terms contain the wave vector () and gravitational drift velocity (), which means that both these processes are happening due to the interaction between the gravity and the em waves.

3.1 Gravity induced ion acoustic waves

There are several methods to analyze this dispersion relation, but we shall start in a simple way by analyzing the term of ion acoustic waves; this term is with the coefficient of the and also independent of wave vector. Therefore we can set the condition (no perturbation) in the dispersion relation and it reduces to


For a high temperature plasma (solar corona), the term is negligible, which leads to the solution . This shows that the electrons undergoes acoustic oscillations due to the inertial (gravitational) drift of ions, even when there is no perturbation. This can be understood from fig.1 & 2, in the absence of the perturbation, the direction of the electric field is in the -x axis and therefore the drift will be in the y-axis. Since the gravitational drift is in the -y axis, for ions these drifts are equal and opposite. But for the electrons the gravitational drift is negligible and hence they experiences the drift in the y axis (see fig.2). The drift of electrons due to , causes a positive potential and by the inertia of the ions, these electrons overshoot back due to Coulomb’s force, leads to the ion acoustic oscillations. As mentioned in the introduction, in the ultra vacuum condition the Coulomb’s interaction is strong, therefore these ion acoustic oscillations are observable in the solar corona.

Figure 4: Plot of growth rate for & ( , , )

The damping of ion acoustic waves is considered as one of the proposed mechanism for the solar coronal heating problem[1, 3, 15], therefore we should find out whether such damping occurs in our case. If we look in to the dispersion relation Eq.(18), we find the ion acoustic term is in the coefficient of with the nonlinear term . This shows that certain energy sharing occurs between the em waves and the ion acoustic waves. To confirm this, we have find out the unstable mode of em waves from the dispersion relation, since the dispersion relation is a quadratic equation, the solution is given as , where a, b and c are the coefficients of respectively. To find the unstable mode of em waves, the two roots are plotted (computationally using Mat-Lab programing) against the wave vector for the coronal plasma. Since one of the root is imaginary, the plot is done for the modulus of in a logarithmic scale to cover a wide range of frequencies, is shown in the fig.3. The real root (dashed line in the graph) represents the stable mode (undistorted part) of propagating em waves and the imaginary root is the unstable mode (wave instability), where the effect of the nonlinear terms are reflected (solid line in the fig.3).

3.2 Instability & Wave particle interactions

A wave instability means, the changes occurs in the physical properties of the wave due to the interaction with the local fields (wave-wave interaction) or the particles of the medium (wave-particle). During the wave-wave interaction; the propagating wave interacts with local oscillations such as Langmuir oscillations, acoustic oscillations and etc., becomes unstable. During this process, either of the wave grow at the cost of the other or if both the waves are damped, that will lead to the heating of the plasma. In the wave particle interaction; the incoming wave interacts with the particles of the medium in such a way that, either the particles gain energy from the wave (Landau Damping process) or it looses to the wave (amplitude of the wave grow due to inverse Landau damping process). When the particles gain energy from the propagating waves, the waves are get damped (the group velocity of the wave, , decreases) or the vice verse.
A simple way to study the instability is, finding the growth rate (imaginary part of ) of the instability and analyze it for the various parameters involved in the driving mechanisms. We have already found from the computational plot that one of the root is imaginary and the condition for a root to be imaginary is, , therefore we can write the growth rate of the instability as,

3.3 Damping of Ion acoustic waves

Since we have started the discussion with ion acoustic waves, we shall see the effect ion acoustic waves over the unstable mode. The ion acoustic term is contained in the coefficient of , we analyze the growth rate for very low wave vector range. For smaller wave vector , the coefficient and which gives , i.e., the growth rate

Figure 5: Gravity Induced Resonant Emission (Plot of the growth rate for coronal plasma) ( , , )

For coronal plasma, the first term is negligible compared with ion acoustic term and hence the reduces to

. Now the growth rate contains two terms, the ion acoustic term and the wave vector terms. Here the ion acoustic term is negative, which reduces the growth of the instability. When there is no perturbation(), becomes real which leads to the earlier results , that the gravitational drift in the magnetized warm plasma generates ion acoustic waves. Here we find in the absence of perturbation, the total thermal energy is converted in to the ion acoustic waves by the gravitational drift ions and the plasma becomes stable. But when ever there is a perturbation, the couples with the of ions which results in the lh oscillation, which leads to the damping of the ion acoustic waves. This could be more clear from the fig.2, in the perturbed plasma, the direction of changes continuously at every point and the plasma becomes unstable. The growth of this unstable mode depends on the available energy; when increases, the increases at the cost of ion acoustic waves. This result is established in the fig.4, is a logarithmic plot between & , for & . In the fig.4, we find increases when (dashed line, increases); this means that the unstable part of the em wave grows at the cost of ion acoustic waves (wave-wave interaction). When , is constant ( is constant), the em wave does not grow. Thus the propagation of em waves in a magnetized warm gravitational plasma, causes the damping of ion acoustic waves. The damping of ion acoustic waves have been observed in the solar corona [21, 2], but this damping will not result in the heating of plasma; because, here the energy is transferred to another wave not to the particles.

3.4 Gravity induced lower hybrid oscillations

Now, to find out the physics of the next nonlinear term, we are analyzing the unstable mode for the entire range of wave vector; the fig.5 shows the plot of with , for a wide range of frequencies.

Figure 6: Plot of growth rate for & , (, , )

We find in this plot, the increases with up to the value of , after that it decreases and then increases. The increase of in the low wave vector range; that we have already discussed, it is due to the damping of ion acoustic waves. But beyond the value , the decreases with , i.e., decreases, signifies that em wave is losing the energy and at certain point, , where the propagating em waves are completely absorbed by the particles. But, here after again increases, which reflects the resonant absorption and the re-emission of the em waves by the ionized particles. This is a very interesting result from the computational plot of the growth rate, now let shall analyze theoretically from the expression of . For high wave vector, , the coefficient, , then

. Here we have neglected the cold plasma and the ion acoustic terms; at the resonant condition, , by setting , we can derive the formula for the resonant frequency,

, where the is the lower hybrid frequency. This shows that the perpendicular propagation of the em waves, induces a lh oscillation in the gravitational plasma, through which the ionized particles undergo the resonant absorption and re-emission. This resonant absorption and re-emission is purely a gravity induced process (see fig.6), we shall name it as, the Gravity Induced Resonance Emission (gire). The lh resonance heating has been considered as a possible mechanism for the study of coronal emissions [4, 7, 8, 9, 10, 11, 12, 13, 14, 16, 17, 18, 20], therefore it is necessary to understand the physics of gire.

3.5 Gravity Induced Resonance Emission

To understand gire, we should know the basic mechanism of lh oscillation[5]. The lh oscillation is a coupled oscillation between the ions and electrons undergoing cyclotron motion, attached by Coulomb’s force. In a magnetic field, all the ionized particles under go cyclotron motion, since the ions are inertial, the Lorentz force is less effective to the ions compared with electrons in a weak magnetic field. In such conditions, the ions are considered as less magnetic compared with electrons, which under go several rotation in the same period of ions. Therefore in a weak magnetic field, if there is a local electric field, the ions experiences more electrostatic force rather than the magnetic force. Thus the ions execute a back and forth motion, while the electrons are effective in the cyclotron motion (fig.7). Since the cyclotron motion of the ions and electrons are in opposite directions, in the every rotation they becomes closer and closer, get coupled by the Coulomb’s Force and leads to a coupled oscillation under the magnetic field. Since this is a coupled oscillation between the ion and electron cyclotron motion and hence the frequency of this oscillation lies between these two (). Now we know the basic requirement for the lh oscillation is the presence of a local electric field, so that ions would be less effective to the magnetic field, undergo a back forth motion while the electrons are in the cyclotron motion. This is happening in our plasma, the gief is the local electric field, so that the formation of lh coupling is very effective in the coronal plasma. During this lh oscillation, in each cycle of revolution, the particles are getting accelerated by the Coulomb’s Force (fig.7a). The high mobility of electrons creates a deficit of field, for both the particles and hence they acquire the field energy from the incoming em waves. This is a pure Landau damping process of em waves induced by gravity, in which both the electrons and ions accelerated in the cyclotron motion to higher and higher energy level. The landau damping occurs in both the cycle of incoming wave, for the ions it happens at the upper half of the wave and for electrons in the lower half (see fig.7c). The electrons are lighter particles, undergo several rotations in the time period of ions and absorb more energy. This increases in the energy of electrons which results in the increase of it’s Larmour radius, which gives more acceleration time for ion, result in the absorption of sufficient energy from the incoming radiation. This is just as a push and pull, the electrons are pushed by the em waves which pulls the ions by Coulomb’s force to the higher energy. When the ions are accelerated to higher energy level, they get ionized more to comply with the electrons to continue in the lh oscillation (more work is needed to calculate total energy and acceleration time). Gravity is an inherent force and hence this process is continuous one, which is happening in the solar corona may be reason for the spectral emissions. Analyzing the spectral emissions on the basis of gire is beyond the scope of this paper.

Figure 7: Model orbits in lh oscillation (a) for (b) for (c) Photon absorption

3.6 Theoretical Justification

  1. For the study of motions perpendicular to B, the fluid theory is the best approximation.

  2. We have used the MHD theory in its best representation for the study of warm plasma, i.e., by multiplying the number density , with the single particle equation.

  3. The Linearization method is a very simple and straight forward procedure for solving the differential equation.

  4. The solar corona is best known ultra vacuum natural laboratory, is suitable for the study of particles interactions.

  5. The plasma chosen for this study is a holistic one, including all the parameters.

  6. And the formula we have derived is also a holistic one, convergent of all the fundamental quantities.

4 Conclusion

In ultra vacuum condition, the fields are very effective and Coulomb’s interactions are very strong and becomes a long range force. In such situation, the particle interactions are very dominant, the electric field induced by gravitational drift of ions is sufficient to cause the perturbations in the plasma. We have analyzed the effect of gravity over ionized particles, by deriving a general dispersion relation using MHD theory. We understand that the gravitational force couples with Lorentz force leads to the ion acoustic and lower hybrid oscillations. In the absence of em perturbations, the gravity causes ion acoustic oscillations in the magnetized plasma and these waves are getting damped by the em waves in the perpendicular direction of the magnetic field, which also results to the formation of lower hybrid oscillations. For higher wave vector, these gravity induced lower hybrid oscillation leads to the resonant absorption and re-emission of em waves. We have analyzed these results for coronal plasma, but it is suitable to all low pressure plasmas.
The first author acknowledges the financial support by UGC under the UGC Research Award Scheme (2016-2018).


  • [1] Angelo, N. D.,1969, SP, 7, 321
  • [2] Aschwanden, M. J., 2005, Physics of Solar Corona, Springer, Praxis Publishing, UK
  • [3] Benz, A. O., Wentzel, D.  G., 1981, A&A, 94, 100-108
  • [4] Benz, A. O., Smith, D. F., 1987, SP 107, 299.
  • [5] Boyd, T.J. M., Sanderson, J. J., 2005, The Physics of Plasmas, Cambridge University Press
  • [6] Chen, F. F., 1981, Intro. to Plasma Physics and Controlled Fusion, Plenum Press, New York
  • [7] Goodrich, C. C., 1985, Proc-1985-Tsurutani, 153.
  • [8] Krall, N. A., Liewer, P. C., 1971, Phys. Rev, A4, 2094
  • [9] Lampe, M., Papadopoulos, K., 1977, ApJ, 212, 886.
  • [10] McBride, J. B., Ott, E., Boris, J. P., Oren, J. H., 1972, Phys. Fluids, 15 2367
  • [11] McClements, K. G., Su, J. J., Bingham, R., Dawson, J. M., Spicer, D. S., 1990, SP, 130, 229.
  • [12] McClements, K. G., Bingham, R., Su, J. J., Dawson, J. M., Spicer, D. S., 1993, ApJ, 409, 465.
  • [13] Meister, C. V., 1995, SP, 160, 65.
  • [14] Papadopoulos, K., 1979, Proc-1979-Akasofu, 289.
  • [15] Revathy, P., 1977, SP, 53, 445
  • [16] Spicer, D. S., Benz, A. O., Huba, J. D., 1981, A&A 105, 221
  • [17] Tanaka, K., Papadopoulos, K., 1983, Phys. of Fluids, 26, 1697.
  • [18] Thejappa, G., 1987, SP, 111, 45
  • [19] Treumann, R. A., Baumjohann, W., 2001, Advanced Space Plasma Physics, Imperial College Press, London
  • [20] Vlahos, L., Gergely, T. E., Papadopoulos, K., 1982, ApJ, 258, 812.
  • [21] Wong, A. Y., Motley, R. W. & D’Angelo, N., 1964, Phys. Rev., 133, A4 36
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description