A simple model for Jeans instability in the presence of a constant magnetic field arxiv:astro-ph/ v1 12 Feb 2004

Transcription

1 A siple odel for Jeans instability in the presence of a constant agnetic field arxiv:astro-ph/ v1 12 Feb 2004 A. Sandoval-Villalbazo a, L.S. García-Colín b, c and A. Arrieta-Ostos a a Departaento de Física y Mateáticas, Universidad Iberoaericana Loas de Santa Fe México D.F., México alfredo.sandoval@uia.x, anabel.arrieta@uia.x b Departaento de Física, Universidad Autónoa Metropolitana México D.F., México c El Colegio Nacional, Centro Histórico México D.F., México lgcs@xanu.ua.x March 17, 2008 Abstract Jeans instability is derived for the case of a low density self-gravitating fully ionized gas in the presence of a constant agnetic field. A new Jeans wave nuber, substantially enhanced by the agnetic field, is established under a suitable syetry assuption. Nevertheless, fro our present knowledge regarding the existence of prieval agnetic fields in the Universe, it appears that they are indeed irrelevant for our understanding of large scale structure foration. Nuerical exaples are given and analyzed. 1 Introduction Magnetic fields have a significant effect on virtually all astrophysical objects. They are observed in all scales. Close to hoe, the Earth has a bipolar agnetic field with a strength of 0.3G at the equator and 0.6G at the poles [2]. Within the interstellar ediu, agnetic fields are thought to regulate star foration via the abipolar diffusion echanis [9]. Our Galaxy has a typical interstellar agnetic field strength of 2µG in both regular ordered and rando coponents. Other spiral galaxies have been estiated to have agnetic field strengths of 5 to 10µG, with fields strengths up to 50µG found in starburst 1

2 galaxy nuclei [1]. Also agnetic fields are fundaental to the observed properties of jets and lobes in radio galaxies, and they ay be priary eleents in the generation of relativistic outflows fro accreting assive black holes [2]. Magnetic fields with typical strength of order 1µG have been easured in the intercluster ediu using a variety of techniques. Large variations in the field strength and topology are expected fro cluster to cluster, especially when coparing dynaically relaxed clusters to those that have recently undergone a erger. Magnetic fields with strengths of 10 40µG have been observed in soe locations [2]. In all cases, the agnetic fields play iportant role in the energy transport in the intercluster ediu and in gas collapse. On the other hand at the cosological level the presence or existence of agnetic fields is ore controversial. In a recent review on the subject [11] it is firly asserted that a true cosological agnetic field is one that cannot be associated with collapsing or virialized structures. Thus the particular role that they ay play in the epoch of galaxy foration is rather obscure. Although liits have been placed on the strength of cosological agnetics fields fro Faraday rotation studies, of high redshift sources, anisotropy easureents of the CMB and the light eleent abundances fro nucleosynthesis, the question reains: Is there the possibility that the Jeans ass arising fro gravitational instabilities responsible for galaxies foration be odified by the presence of a agnetic field?. In spite of the dubious background provided by our present knowledge, this question has been tackled since over fifty years ago. In fact, already Chandrasekhar & Feri [3] reached the conclusion that Jeans criteria for the onset of a hydrodynaic instability is unaffected by a agnetic field in an extended hoogeneous gas of infinite conductivity in the presence of an unifor agnetic field. However, in their calculation they assued that the ediu was agnetized and that, within the gas, there existed a fluctuating agnetic field. This proble has been retaken by several other authors in different contexts. In particular Lou [6] studied the proble of gravitational collapse in a agnetized dynaic plasa in the presence of a finite aplitude circular polarized Alfvén wave. This author does find a case in which Jeans wave nuber k J is indeed odified by the agnetic field by a ter proportional to [c c2 A ]1/2 where c 0 is the velocity of sound and c A = Bz 0 (4πρ 0 ) 1/2 is Alfvén wave speed Bz 0 being the z-coponent of the unifor agnetic field. Other attepts to show that agnetic fields do play an essential of roll in galaxy foration have been perfored, Ki, Olinto & Rosner [5], although not specifically addressing the question of a agnetic instability. In view of all these efforts we still feel that the siple question of whether or not a dilute non-agnetized plasa cloud placed in the presence of an external, unifor agnetic field in which density fluctuations are also present due to a fluctuating gravitational field, exhibits a Jeans wave nuber which is odified by the presence of the field, has not yet been fully discussed in the literature. This is the purpose of the present work. The basic and rather siple foralis is given in 2. 3 is devoted to the derivation of the dispersion relation leading to the odified for of k J and soe attepts to place the relevance of the 2

3 results within a realistic frae for existing agnetic field intensities. Soe concluding rearks are given in 4. 2 Basic Foralis We start by assuing that the dynaics of the dilute plasa is governed by Euler s equations of hydrodynaics, naely the balance equations for the fluid s ass density ρ( r, t) and its velocity u( r, t). Thus, (ρ u) ρ + (ρ u) = 0 (1) + (ρ u u) + p = f g + f M (2) In equation (2), f g is the force arising fro the gravitational field and f M = q ρ o( u B) is the Lorentz force. This iplies that the plasa is diluted enough so that the external field B will act only on the eleentary charges ±q of the plasa. Here is the ass of the point charge q. Eqs. (1-2) can be linearized by introducing density and velocity fluctuations defined by: ρ = ρ 0 + δρ (3) and, u = u 0 + δ u (4) δθ (δ u) (5) where ρ 0 is the average density. The fluid is assued to be static, so that u 0 = 0, ϕ represents the gravitational potential and the external agnetic force is that corresponding to a constant agnetic field B = (B 0 + δb) ˆk, so that the linearized equations for the density and velocity fluctuations can be written as: and (δρ) + ρ 0 δθ = 0 (6) ρ 0 (δ u) + (δp) = ρ 0 (δϕ) + q ρ 0 (δ u B ) 0 (7) where f g = (δϕ), δϕ being the fluctuating gravitational potential. Neglecting teperature fluctuations, the pressure ter in equation (7) ay be rewritten 3

4 in ters of the density fluctuations through the local equilibriu assuption naely, p = p(ρ) so that ( ) p p = ρ = c2 0 ρ 0 γ ρ T We know recall that K T, the theral copressibility satisfies the relation K T = γ/c 2 0 where γ = C p/c v and c 0 is the velocity of sound in the ediu. Eq. (7) ay now be rewritten as ρ 0 (δ u) + c2 0 γ (δρ) = ρ 0 (δϕ) + q ρ 0 (δ u B ) 0. (8) Assuing that δϕ is given by Poisson s equation, 2 (δϕ) = 4πGδρ, that B o = B oˆk (ˆk is the unit vector along the z-axis) and noticing that for this case the last ter equals B o (ĵu x îu y ), equation (8) reduces to ρ o (δθ) + c2 o γ 2 (δρ) = 4πGρ o (δρ) + q ρ o ( δ u)ˆk (9) after taking its divergence and using equation (6). Equations (6) and (9) are now two siultaneous equations for δρ and δ u which need to be solved. To do so we introduce an approxiation, naely that there exists cylindrical syetry with respect to the z axis so that z = 0. With this approxiation the last ter in equation (9) can be calculated taking the curl of equation (8). This yields ρ 0 ( δ u)ˆk = q ) B 0ρ 0 (ĵδux îδu y ˆk = q B 0ρ 0 δθ. Finally, taking the tie derivative of equation (9), using the previous equation and aking use of Eq. (6), one is finally led to the result that 3 3 (δρ) + c2 0 γ 2 + ( qb0 ( (δρ) ) ( ) (δρ) + 4πGρ 0 δt dt ) 2 (δρ) = 0 Integrating once with respect to tie and setting the integrating constant equal to zero which deprives of no generality to our arguent, we get that 4

5 2 2 (δρ) + c2 0 γ 2 (δρ) + 4πGρ 0 (δρ) ( ) 2 qb0 + (δρ) = 0 (10) Equation (10) is now a single equation for the density fluctuations δρ. Its solution readily achieved proposing that δρ is described by a plane wave δρ = Ae i( k r ωt), (11) where all sybols are self defined. Substitution of equation (11) into equation (10) yields the dispersion relation, ω 2 c2 0 γ k2 + 4πGρ 0 + ( ) 2 qb0 = 0. (12) Instabilities enhanced by the gravitational and agnetic field arise when the roots for ω in this equation are iaginary. The threshold value of k beyond which this happens is precisely Jeans wave nuber and is here given by [ kj 2 = γ ( ) ] 2 qb0 c 2 4πGρ 0 +. (13) 0 Clearly, if B 0 = 0 we recover the well known expression for k J. The question now is how relevant is the second ter in enhancing or hindering structure foration. This will be analyzed in the following section. 3 Analysis of the Dispersion Relation As indicated in the previous section, the linearized version of fluctuating nondissipative hydrodynaics, applied to a a dilute hoogeneous plasa in the presence of a unifor agnetic field, leads to a Jeans wave nuber which is, as depicted in equation (13), a superposition of the gravitational and agnetic fields. The question is if this result has any bearing on the value of Jeans ass in realistic cases. Clearly equation (13) points at three possibilities, naely, if the agnetic ter is a) negligible, b) of the sae order or c) larger than the gravitational contribution. As we recall, Jeans ass is defined as ( 2π ) 3 so that using equation (13) M J 4π 3 ρ oλ 3 J = 4π 3 ρ o M J = 32π4 3 ρ o γ C2 o 1 k J 4πGρ o + ( qbo ) 2 3/2 (14) 5

6 where ρ o = n o, n o being the particle density in the plasa. As it has been exhaustively discussed in the literature [4] [7] [8] [10] without the agnetic contribution, present values of ρ o gr/c 3, = H gr and T 10 5 K yield M J M. Taking now q/ as the proton-ass ratio which is, in absolute value, approxiately equal to 10 8 c/kg, how does the agnetic ter rates with respect to the gravitational one for different values of B 2 o?. Due to the uncertainties of the estiates for B o according to current data [11] this ay not be so easy to quantify. According to the constraints iposed on cosological agnetic fields fro the different observations entioned in the introduction, their strength sees to be confined in the range [11], B 10 6 G. Taking B = 10 6 G = T, the ter q B o 10 2 s 1 eaning that its squared value is 10 2 s 2. This would iply that M J M, contrary to observations. Thus in agreeent with Chandrasekhar and Feri [3], the agnetic ter is copletely negligible insofar as the instabilities required for the foration of large structures, say galaxies. In order that q B o s 1, B o ust be of the order of G which by the way is indeed of the order of agnitude of strengths that are envisaged to be detected at the cosological level [11]. However, this is still a thought and not a fact so that expectation will reain until these easureents are achieved. 4 Conclusions Even for the very siple case of a dilute not agnetized plasa in the presence of a hoogeneous unifor agnetic field the value of Jeans s ass as predicted by gravitational instabilities sees not to be substantially odified according to our present knowledge of the values available for cosological agnetic fields. Yet there is a possibility that uch lower values could be observed in which case the value of M J would be substantially odified. Thus, it appears that agnetic fields are probably not essential to our understanding of large-scale structures foration in the Universe. This work was perfored under grants fro CONACyT F. References [1] Beck, R., Brandenburg, A., Moss, D., Shukurov, A. & Sokoloff, D. 1996, Annu. Rev. Astron. Astrophys., 34, 155 [2] Carilli, C.L. & Taylor, G. B. 2002, Annu. Rev. Astron. Astrophys. 40, 319 [3] Chandrasekhar, S. & Feri, E. 1953, ApJ, 118, 116 6

7 [4] Jeans J. 1945, Astronoy and Cosogony, Dover Publications, Inc. p.645 [5] Ki, E., Olinto, A. & Rosner, R. 1996, ApJ, 468, 28 [6] Lou, Y.-Q. 1996, MNRAS, 279, L67 [7] Peebles, P.J.E. 1993, Principles of Physical Cosology, Princeton Univ. Press, Princeton, N.J. 2nd. Edition [8] Sandoval-Villalbazo, A. & García-Colín, L. S. 2002, Class. and Quant. Gravity, 19, 2171 [9] Spitzer, L. 1978, Physical Processes in the Interstellar Mediu. New York: Wiley [10] Weinberg, S. (1971), ApJ, 168, 175 [11] Widrow, L.M. 2002, Rev.Mod.Phys., 74,775 7