Volume 504, Number 2, September III 2009
|Page(s)||321 - 329|
|Published online||09 July 2009|
Spectrum of plasma modes and relevant transport processes in astrophysical disks
Massachusetts Institute of Technology, Cambridge, Massachusetts, USA
Received 27 October 2008 / Accepted 29 June 2009
A simple plasma disk structure imbedded in a magnetic field and in the (prevalent) gravity of a central object is shown to be subject to the excitation of significant axisymmetric and tridimensional modes. The key factors involved in the relevant instability are the plasma pressure vertical gradient and the rotation frequency (around the central object) radial gradient. A modestly peaked vertical profile of the plasma temperature is shown to drive a ``thermo-rotational instability'' with considerable growth rates. Unstable modes are found as well for ``flat'' temperature profiles. The tri-dimensional tightly wound spirals that are found have properties that, unlike the familiar galactic spirals, depend on the vertical profiles of their amplitudes. Both radially standing and convective (quasi-modes) spirals are identified. Within the considered spectrum, unstable modes are shown to produce opposing fluxes of particles and thermal energies in the vertical direction. Thus disks with relatively flat vertical temperature profiles can expel particles (winds) from the equatorial plane while transporting thermal energy inward. An effective ``diffusion'' coefficient, for energy and angular momentum, is derived from the structure of radially convective spiral modes and shown to be consistent with significant radial transport rates.
Key words: accretion, accretion disks - black hole physics - magnetohydrodynamics (MHD) - instabilities - magnetic fields - gravitation
Identifying the plasma collective modes that can be excited in plasma disk structures (Pringle 1981; Blandford 1976; Lovelace 1976; Coppi & Rousseau 2006) surrounding compact objects such as black holes can be important to explain experimental observations associated with objects of this kind. The geometry of these disk structures (Coppi & Coppi 2001) and their physical parameters, that include in particular the radial gradient of the rotation frequency, the vertical gradients of the particle density and temperature and the effects of the magnetic field in which they are imbedded (Coppi 2008a), determine the characteristics and the spectrum of the modes that can be excited. Novel transport processes of particles and thermal energy in the vertical direction as well as radial transport of angular momentum can be produced by the identified modes. In fact, these processes and the tri-dimensional structure of the spiral modes (Coppi 2008b) that are found can be correlated with important observations connected with black holes (Coppi & Rebusco 2008a). We point out that the properties of the tightly wound spirals that are introduced depend on the vertical profiles of their amplitudes, a feature that is not shared with the better known theory of galactic spirals (Bertin 2000). We note also that the theory presented here for both axisymmetric and tri-dimensional modes include, among others, the physical ingredients of basic modes such as the so called MRI instability (Velikhov 1959; Chandrasekhar 1960; Balbus & Hawley 1991), that is appropriate for different geometries such as cylinders, or the modes that can be excited in well confined plasmas (Coppi & Spight 1978) and can explain the experimentally observed inward transport of particles in connection with the outward transport of plasma thermal energy.
In Sect. 2 the characteristic parameters of the simplest plasma disk structure, imbedded in a vertical magnetic field and surrounding a compact object, are defined in view of the analyses presented in the next sections. In Sect. 3 the basic equations are given for the spectrum of axisymmetric and tri-dimensional normal modes that can be excited in a thin plasma disk. In Sect. 4 the realistic conditions on the density and temperature vertical profiles are examined under which nearly isobaric, weakly compressible modes are found. In Sect. 5 the theory of axisymmetric modes is given and the conditions for their marginal stability are derived. The radial gradient of the rotation frequency and the vertical gradient of the plasma pressure are recognized as the key factors for the excitation of these modes. In particular, when the ratio of the relative gradient of the plasma temperature exceeds the relative gradient of the density by 2/3 a new form of instability (the ``thermo-rotational instability'') emerges (Coppi 2008a). This is in addition to that of the ballooning modes analyzed by Coppi & Keyes (2003) that are found when the ratio of the gradients mentioned earlier is 2/3 (polytropic profile). In Sect. 6 the transport produced in the vertical direction by the unstable modes described in Sect. 5 is evaluated by the relevant quasi-linear theory. Then the suggestion is made that the outward particle transport (away from the equatorial plane) occurring in the presence of relatively flat temperature profiles be the seed for the observed winds emanating from disks structures surrounding black holes (Elvis 2000). In this case an inward (toward the equatorial plane) flux of thermal energy is produced. The mode growth rates of the unstable modes on which the relevant transport coefficients depend are evaluated by completing the analysis given in Sect. 5. In Sect. 7 the characteristics of the weakly damped oscillatory modes, can be found when the temperature vertical profile is relatively flat, are examined. In Sect. 8 the theory of tri-dimensional spiral modes that are localized both in the vertical and the radial direction is given. These modes co-rotate with the plasma at the radius around which they are localized. The excitation of modes of this kind in the proximity of a black hole is, in fact, suitable to explain the observed high frequency quasi periodic oscillations (QPO's) as proposed by Coppi & Rebusco (2008a). In Sect. 9 a new class of spiral quasi-modes that are convective in the radial direction and localized in the vertical direction are introduced. These modes are weakly damped and oscillatory in time and can be associated with the coupling to modes that are linearly unstable. An important property of the convective modes which are described is that they can transfer their angular momentum and energy at relatively fast rates in the radial direction. Thus, they may be considered as possible candidates to provide the plasma ``viscosity'' (Pringle 1981) associated with the mass accretion rates inferred from the luminosities of relevant objects. In Sects. 10 concluding remarks concerning the results presented in the previous sections are given.
where and p is the total plasma pressure. The radial equilibrium equation, to lowest order in the ratio H02 /R02, reduces to
where is the rotation frequency, is the Keplerian frequency, M* is the mass of the central object and is the toroidal velocity. The relevant vertical equilibrium equation is
where is the mass density. Then Eq. (2) gives for , where . We shall consider a variety of temperature profiles including the flat profile corresponding to over the height of the disk describes the entire density profile. The scale distance for the pressure gradient is defined by for .
The electric field that is present in the equilibrium configuration we have considered can be identified by referring to the electron momentum conservation equation
where is the electron pressure, and . Realistically, .
Finally, we note that there can be different current carrying stationary plasma axisymmetric disk structures from that analyzed in this section. In particular, we may mention the configurations involving a sequence of plasma rings identified by Coppi & Rousseau (2006). We observe also that the modes analyzed in the next sections may be considered as ``seeds'' for new stationary configurations.
We consider normal mode perturbations from the indicated initial state represented by
in an interval around R0, such that , where is the mode growth rate, the frequency, the toroidal mode number and the perturbed toroidal velocity. The basic linearized equations that describe the departure from the initial state include
and the total momentum conservation equation
Here, we have taken into account that the initial magnetic field Bz is considered to be varying over scale distances of the order of R0 and have used standard symbols.
Furthermore, it is reasonable to assume that the collisional mean free path is short relative to the distance
over which the mode is localized vertically and to the mode radial wavelengths. Thus, the thermal conductivity can be neglected and the adiabatic equation of state can be adopted, that is,
where is the adiabatic index. Then, given Eq. (6) we obtain
and we choose to consider modes that co-rotate with the plasma at R = R0.
Therefore, we take
where is not linear in , and Eq. (6) reduces to
for where the subscript p indicates the relevant poloidal component. Then, the perturbed pressure and density are given by
and we consider .
The perturbed magnetic field components are derived from the combination of Eqs. (7) and (8) that gives
As we shall see, the considered limit
leads to find that
and in this case the modes of interest can be characterized as ``quasi-isobaric'' and ``weakly compressible'' in the sense that
To prove this point we take the z-component of Eq. (9)
and consider the gravitational term as important in Eq. (22), that is
where . Then and since the ordering (23) implies that
Here we limit our considerations to modes for which , as we shall show that and we restrict our analysis to for the sake of simplicity. Then we have
and that Eq. (21) is valid for . Moreover
and is obtained as a function of from Eq. (21) that becomes
We see that the sign of C0 is important in relating that of to . Finally, we can verify that, for , the toroidal pressure gradient can be neglected in the -component of Eq. (9). Then this reduces to
In fact, we limit our analysis to the case where
Finally, it is convenient to refer to the equation, derived from Eq. (9), that is
By using Eqs. (18) and (26) this can be rewritten as
Thus we have a complete set of equation given by Eq. (30) combined with
that relate and to .
2008a) with , that are purely growing. Thus, and we look for normal modes of the form
that are a special case of those represented by Eq. (6) and where , , , and is the Alfvén velocity. In addition is an even or odd function of z that is localized over a distance represented, for instance, by . The unstable modes that are found have
corresponding to the fact that the radial gradient of the rotation frequency is a key driving factor for the relevant instability. Clearly, is the representative bending factor of the magnetic field lines.
Then Eq. (31) reduces simply to
Here we consider and therefore with , that is . Then since and Eq. (29) becomes simplify
where and Eq. (34) can be reduced to
where . Since we look for localized solution soltuion in z it is convenient to take the Fournier transform of Eqs. (29) and (35) that gives
and the problem is reduced to solving a second order differential equation.
We note that at marginal stability
and Eq. (39) reduces to
Thus we obtain the known ballooning solution
for C0=0. Clearly unstable modes of this kind can be found for C0<0, that is for
and correspond to the ``thermo rotational instability'' presented in Coppi (2008a).
Since flat temperature profiles corresponding to C0<0 are an important class. to be considered, higher modes, for which
Ek0>3, have to be analyzed. For these modes the influence of the radial gradient of the rotation frequency is stronger. Then it is convenient to use the variable
and to rewrite Eq. (40) as
It is also useful to derive a quadratic form obtained by multiplying Eq. (40) by and integrating it from to . This gives
The next eigensolution is odd with
In this case
Proceeding to the next eigenfunction we find
For C0<-14/5 a still higher eigensolution has to be considered.
The considered modes can produce particle density transport, in the vertical direction, that is of contrary sign to that of the temperature transport and modify the density and temperature profiles in such a way as to lead
toward 2/3, corresponding to a polytropic. Thus, if
a particle inflow toward the equatorial plane is induced. When
including the case where
or where the surface of the disk can be hotter than the interior, the particle transport is away, from the equatorial plane. These arguments are based on the quasilinear analysis that gives the vertical particle flux produced by unstable modes as
where indicates an average over a radial distance such that . The corresponding temperature flux is .
The outflows produced by these modes when can be considered as candidates to explain the origin of the particle fluxes (winds) that have been observed to emanate from disk structures such as those at the core of AGN's (Elvis 2000).
The transport process described by Eq. (51) is similar to that proposed for the theoretical explanation of the experimentally observed particle inflow in magnetically confined toroidal plasmas that is associated (Coppi Spight 1978) with the outflow of electron thermal energy related to the ratio of the gradients of the radial electron temperature and the particle density.
To illustrate this case we refer, for simplicity, to a cylindrical or plane plasma confinement configuration where the magnetic field is in the z direction. Electrostatic modes can be excited with longitudinal phase velocities
being the electron thermal velocity. Considering the limit where the electron longitudinal thermal conductivity is relatively large but finite, it is possible to show that the radial transport velocity for the electrons is equal to and of opposite sign of that of the electron thermal energy. In particular, when
that is a condition most commonly verified in well confined plasma experiments, the thermal energy flux is outward and the particle flux is inward. We note that the experimental evidence of the particle inflow process in well confined plasmas mentioned above has expanded over the years and by now it is well documented.
It has been suggested (by the referee) that the condition be compared to the instability condition of the classical Schwarzchild instability for stellar interiors where g is the gravity acceleration and is the heat capacity at constant pressure. In fact in both cases a transport of temperature associated with an opposing flow of matter is involved.
We refer to Eq. (39) and note that the quadratic form associated with it is
In the realistic limit where this form can be used as a variational principle to evaluate . Considering this limit the ``ordering of maximum information'' for which all the remaining terms in Eq. (52) are retained is , , .
The lowest (even) eigensolution, Eq. (41), of this gives
It is clear that the instability requires relatively peaked temperature profiles, that is corresponding to the thermo-rotational instability. In particular,
In the case where
showing that the dependence of and on is weak, in this case.
The next (odd) eigenfunction, Eq. (47), corresponds to the ballooning mode investigated already by Coppi & Keyes (2003). In this case the dispersion relation is
Ek0 = 5 + C00. Clearly the mode can be unstable even if C00 < 0 that is for
We observe that, relative to the case investigated by Coppi & Keyes (2003) the relevant growth rate is increased by the term .
The next (even) ``harmonic'' has the form
which allows the instability to exist up to . We note that when , as . When a higher eigensolution has to be considered.
A class of modes that are oscillatory and weakly damped can be found when
is well below 2/3, that is when C0 is negative and finite. Clearly, these modes can be significant only if their amplitudes can be sustained by a non-linear coupling to modes that are self excited. Unlike the unstable modes discussed in the previous sections those considered here have
kR2>k02 and in particular
This means that the restoring effect of magnetic field line bending prevails over that of the radial gradient of the rotation frequency. Then instead of Eq. (39) we consider
Because of the singularity appearing in Eq. (61), will acquire a relatively small imaginary part, that is . Therefore, it is convenient to refer to the quadratic form that can be obtained from Eq. (61), observing that is complex in this case, in order to ascertain the sign of . This is
It is also convenient to introduce dimensionless variables defining
Thus , A2>1 and Eq. (61) becomes
We note that, for Ek>1, Eq. (61) gives
The fact that is related to A2 being larger than unity. In particular, . The case where and remains to be analyzed.
Referring to Eq. (13), we consider modes that are represented by
where , , varies over a radial scale distance that is considerably larger than and we take to be relatively low. These modes are of the spiral type and we discuss at first the case where and is given approximately by the dispersion relations for axisymmetric modes analyzed in Sect. 4, for
where is the radial width over which the mode is localized. In fact, this is a rather easy condition to satisfy as , and for where , is the sound velocity and .
We observe that the low-
spiral modes considered here are directly connected to the axisymmetric modes discussed in Sects. 5 and 7. Referring to the radial displacement
we note that its expression, consistent with that of
given by Eq. (65), can be expanded as follows
where , and
The expansion (67) implies that we take
i.e. , and we consider the case where and both are separable function of z and R - R0.
Then we refer to Eq. (31) and observe that when condition (66) is satisfied it reduces to
where , , and . We notice that all the terms proportional to are grouped on the r.h.s. of Eq. (69) and that, for , this reduces to
We observe also that the slow dependence of on can be expressed through if
It is important to point out that, in order that F0 be a localized function of the sign of kR and have to be related by the requirement that
and, since , , this corresponds to trailing spirals.
Then we define
and condition (66) becomes, approximately,
with . The condition can be easily satisfied as well for , . In particular, taking into account that the vertical profiles of the modes represented by Eq. (67) are given by the theory of the corresponding axisymmetric modes, the spiral modes we consider become the form
where represents the eigensolutions of Eq. (31).
We note that the ``vertical ballooning'' structure of the considered modes which depends on the density and temperature profiles of the equilibrium configuration, involves an analytical description that may be considered intermediate between a global and a shearing-box type of analysis (Umurhan 2008). It is evident from the derivation of the growth rates given in Sect. 7 that the radial gradient of the rotation frequency is the key factor for the excitation of the considered modes. In addition, a moderately peaked temperature profile corresponding to enhances the instability found when and add new unstable modes all of which require the radial gradient of the rotation frequency.
Finally, we may argue that, given the radial localization acquired by spiral modes, these should achieve higher amplitudes than axisymmetric modes which can be standing over considerably larger scale distances. Clearly, modes with higher values of have narrower radial profiles. Unlike axisymmetric modes, spiral modes have angular momentum (Bertin 2000). We note also that the radial localization distance for spiral modes is related to the vertical pressure gradient and to the rotation frequency through the growth rate . This constitutes an intrinsic difference from the spiral modes found for galactic disks whose features are not considered to depend on the characteristics of their vertical profiles.
A first report of the theory of these spiral modes and of those analyzed in the next section was given by Coppi (2008b). The excitation of modes of this kind in the proximity of a black hole has been proposed as an explanation of high frequency quasi-periodic-oscillations (QPO's) by Coppi & Rebusco (2008a).
Another kind of spirals is of special interest in view of the role they can have in the outward transport of angular momentum as needed for mass accretion by the central object. These spiral are oscillatory, in the frame where they corotate with the plasma at a given radius R = R0, and are of the convective type, radially, over scale distances that can be considerably larger than the ``microscopic'' scale distances such as that have been introduced. In fact, they may characterize the ``equilibrium'' state of a large class of accretion disks providing a means for angular momentum flow away from the central object.
Clearly, the modes here are derived from the axisymmetric oscillatory modes discussed in Sect. 8. The vertical structure of these modes can be given by the lowest order solution of Eq. (61), that is
when corresponding to relatively flat temperature profiles. The complete form of this mode can be represented as
over scale distances such that
where , and can be derived from Eq. (61) after replacing by .
Clearly, this mode is oscillatory (in R) over shorter scale distances than that required by condition (79) as
Convective modes localized over shorter scale distances can be constructed out of this class of modes providing a means to transport energy and angular momentum associated with the mode away from R0, provided . Clearly, this indicates a trailing spiral configuration as in the case of localized modes.
In particular, the type of mode packet centered around
that can be obtained by a superposition of the modes represented by Eq. (78) is exemplified by the integral
where . Then, for
Therefore, we may define an effective transport coefficient representing the ``convectivity'' of the considered modes that is defined as
It is evident that can have relatively large values and exceed the Shakura-Sunyaev coefficient (1973) for the outward transport of angular momentum, , when .
Then a reasonable scenario to be envisioned is one where the eigensolutions of Eq. (39) can be unstable and the corresponding (normal mode) spirals are radially localized around R - R0. The amplification of convective modes can be tied to the coupling to linearly unstable modes providing a natural saturation process for the latter kind of mode. Moreover, we may argue that the minimum radius R = R0 from which a significant rate of outward angular momentum transport can be produced by the excited spirals could represent the radius around which the accreting plasma accumulates. Finally, we note that when the vertical structure of oscillatory modes has to be properly analyzed.
- The identification of tridimensional spiral modes for which the vertical structure (unlike the case of galactic spirals) plays a key role in their properties.
- The vertical transport of both particles and thermal energy, in opposing vertical directions, produced both by axisymmetric and non-axisymmetric (e.g. spiral) modes which can lead to an outflow of plasma from the disk structure.
- The coexistence of both standing and convective spiral modes, providing a means to transport angular momentum radially with significant rates. These may be consistent with those needed to justify the rates of mass accretion inferred from the observed luminosities of significant astrophysical objects.
Clearly, the analysis presented here needs to be extended by appropriate computational methods, beyond the limits of the quasi-linear approximation used in Sect. 6 in order to investigate the transport processes associated with the considered modes. For this however, the shortcomings of currently adopted computational approaches such as those pointed out by Regev & Umurhan (2008) need to be resolved.
It is a pleasure to thank C. Crabtree and P. Rebusco for their comments and contributions to different subjects and issues dealt with in this paper and the (anonymous) referee for his suggestions . This work was sponsored in part by the US Department of Energy.
- Balbus, S. A., & Hawley, J. F. 1991, ApJ, 376, 214 [NASA ADS] [CrossRef] (In the text)
- Bertin, G. 2000, Dynamics of Galaxies (Publ. Cambridge Un. Press) (In the text)
- Blandford, R. D. 1976, MNRAS, 176, 465 [NASA ADS] (In the text)
- Chandrasekhar, S. 1960, Proc. Natl. Acad. Sci., 46, 253 [NASA ADS] [CrossRef] (In the text)
- Coppi, B. 2008a, Europhys. Lett., 82, 19001 [CrossRef] (In the text)
- Coppi, B. 2008b, Paper P1 (E.P.S. Conference on Plasma Physics, Crete, Greece), 177 (In the text)
- Coppi, B., & Spight, C. 1978, Phys. Rev. Lett., 41, 551 [NASA ADS] [CrossRef] (In the text)
- Coppi, B., & Coppi, P. S. 2001, Phys. Rev. Lett., 87, 051101 [NASA ADS] [CrossRef] (In the text)
- Coppi, B., & Keyes, E. A. 2003, ApJ, 595, 1000 [NASA ADS] [CrossRef] (In the text)
- Coppi, B., & Rousseau, F. 2006, ApJ, 641, 458 [NASA ADS] [CrossRef] (In the text)
- Coppi, B., & Rebusco, P. 2008, P., Paper P5 (E.P.S. Conference on Plasma Physics, Crete, Greece), 154
- Elvis, M. 2000, ApJ, 545, 689 [CrossRef] (In the text)
- Lovelace, R. V. E. 1976, Nature, 262, 649 [NASA ADS] [CrossRef] (In the text)
- Pringle, J. E. 1981, ARA&A, 19, 137 [NASA ADS] [CrossRef] (In the text)
- Regev, O., & Umurhan, O. H. 2008, A&A, 481, 1 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Shakura, N. I., & Sunyaev, R. A. 1973, A&A, 24, 337 [NASA ADS] (In the text)
- Velikhov, E. P. 1959, Sov. Phys. JETP, 36, 995 (In the text)
- Umurhan, O. H. 2008, A&A 489, 953 (In the text)
Copyright ESO 2009
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.