Advances in Mathematical Physics

Volume 2015, Article ID 916026, 13 pages

http://dx.doi.org/10.1155/2015/916026

## Exact Relativistic Magnetized Haloes around Rotating Disks

^{1}Facultad de Ciencias Básicas, Universidad Tecnológica de Bolívar, CP 131001, Cartagena, Colombia^{2}Instituto de Ciencias Nucleares, Universidad Nacional Autónoma de México, AP 70543, 04510 México, DF, Mexico^{3}Federal University of Latin-American Integration, P.O. Box 2123, 85867-670 Foz do Iguaçu, PR, Brazil^{4}Casimiro Montenegro Filho Astronomy Center, Itaipu Technological Park, 85867-900 Foz do Iguaçu, PR, Brazil

Received 16 August 2015; Revised 9 September 2015; Accepted 18 November 2015

Academic Editor: Klaus Kirsten

Copyright © 2015 Antonio C. Gutiérrez-Piñeres and Abraão J. S. Capistrano. 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

The study of the dynamics of magnetic fields in galaxies is one of important problems in formation and evolution of galaxies. In this paper, we present the exact relativistic treatment of a rotating disk surrounded by a magnetized material halo. The features of the halo and disk are described by the distributional energy-momentum tensor of a general fluid in canonical form. All the relevant quantities and the metric and electromagnetic potentials are exactly determined by an arbitrary harmonic function only. For instance, the generalized Kuzmin-disk potential is used. The particular class of solutions obtained is asymptotically flat and satisfies all the energy conditions. Moreover, the motion of a charged particle on the halo is described. As far as we know, this is the first relativistic model describing analytically the magnetized halo of a rotating disk.

#### 1. Introductory Remarks

In the observational context, many ambiguities still exist about the main constituents, geometry, and dynamics (thermodynamics) of the galactic disk-haloes. However, there are several different observations which probe the galactic and surrounding galactic magnetic field. For instance, it can be measured by nonthermal radio emission from energy equipartition that results from the interaction of magnetic energy with relativistic particles which can play a role in the formation of arms in spiral galaxies (see Krause’s “Magnetic Fields and Halos in Spiral Galaxies” [1] and “Magnetic Fields in Spiral Galaxies” [2] and references therein). For nearby galaxies, other probes are used as optical polarization, polarized emission of clouds and dust grains, maser emissions, diffuse radio polarized emission, and rotation measures of background polarization sources also. Magnetic fields can be seen as a new structural quantity and added to other constituent parameters (e.g., -alpha lines, density mass, local shocks, etc.) are important to study the formation and dynamics of galaxies. Other works propose that magnetic fields also play important role in large-structure formation [3].

It is important to stress that magnetic fields are found mainly in interstellar medium and can be found in every type of galaxies but remarkably noticed in spiral galaxies (see [4] and references therein). For instance, Milky Way has been actively studied in its three regions (central bulge, halo, and accretion disk). This motivates obtaining a proper description of magnetic fields from an oriented general principle and relativistic thin disks have been revealed as one of the best ways to study them. Relativistic disks were also extensively studied in literature in several configurations. Exact solutions that have relativistic static thin disks as their sources were first studied by Bonnor and Sackfield [5] and T. Morgan and L. Morgan [6, 7]. Subsequently, several classes of exact solutions corresponding to static [8–17] and stationary [18–22] thin disks have been obtained by different authors.

The superposition of a static or stationary thin disk with a black hole has been considered in [23–31]. Thin disks around static black holes in a magnetic field have been studied in [32]. Relativistic disks embedded in an expanding Friedmann-Lemaître-Robertson-Walker universe have been studied in [33] and perfect fluid disks with halos in [34]. Furthermore, the stability of thin disks models has been investigated using a first-order perturbation of the energy-momentum tensor in [35]. On the other hand, thin disks have been discussed as sources for Kerr-Newman fields [36, 37], magnetostatic axisymmetric fields [38, 39], and conformastatic and conformastationary metrics [40–42]. Also, models of electrovacuum static counterrotating dust disks were presented in [43], charged perfect fluid disks were studied in [44], and charged perfect fluid disks were presented as sources of static and Taub-NUT-type spacetimes in [45, 46]. Also, monopole and dipole layers in curved spacetimes were analyzed in [47] and electromagnetic sources distributed on shells in a Schwarzschild background in [48].

The thin disks with magnetic fields studied in [36–38] were obtained by means of the well-known “displace, cut, and reflect” method. Solutions that can be interpreted as thin disks of finite extension can be obtained if an appropriate coordinate system is introduced. A coordinate system that adapts naturally to a finite source and presents the required discontinuous behavior is given by the spheroidal coordinates. Some examples of finite thin disks obtained from vacuum solutions expressed in these coordinates can be found in [5, 6, 9, 11] and from electrovacuum solutions in [42].

Moreover, interestingly, magnetic fields seem to play an important role in the formation of jets (resulting from collimated bipolar outflows of relativistic particles) and the accretion disk near supermassive black holes [49]. A current revision of the status of our knowledge about the magnetic fields in our Milky Way and in nearby star-forming galaxies is summarized in [50]. Additionally, a study of the disk and halo rotation is reported in [51], whereas the possibility of magnetic fields that can be generated in the outskirts of disks is studied in [52]. Solutions for the Einstein and Einstein-Maxwell field equations which are consistently applicable to the context of astrophysics remain a topical problem. Nevertheless, the effects of magnetic fields on the physical processes in galaxies and their disk-halo interaction have been scarcely considered in the past. Similarly, the relevance of relativistic models of disks around black holes in a magnetic field is discussed in [32].

The presence of the electric field on the dark matter halo models has been considered in [53], whereas the presence of electromagnetic field in the halo-disk system has been studied in [54, 55] in which the gravitational sources are statics. In [54] we provided detailed overview of the research in the relativistic disks; accordingly we will not repeat them here.

In this paper, we do not focus on the whys and wherefores of the detailed discussion on the probes and relevance of magnetic fields on disks but we apply a standard galaxy modeling as a stationary thin disk and, correspondingly, we associate the halo with the region surrounding the disk. We present the conformastationary version of the static thin disk-halo systems studied in [54]. In addition, we study the features of the principal quantities characterizing the dynamic of the magnetized haloes corresponding to the disks presented in [56]. Therefore, we take the definition in [57] as standard, following the original terminology by Synge [58]: conformastationary are those stationary spacetimes with a conformally flat space of orbits.

Accordingly, we show that the rotating disk-haloes with isotropic pressure, stress tensor, and heat flow generalize the static disk-haloes obtained in [54]. Our results are compatible with those presented in [53] on possible features of galactic halo. Moreover, the description of the motion of charged particles on disk is deduced and is in agreement with the results of the similar analysis discussed in [59]. As far as we know, this is the first relativistic model describing analytically the relativistic magnetized halo of a rotating disk.

The paper is organized as follows. In Section 2, the distributional Einstein-Maxwell equations for haloes surrounding thin disks are obtained. In Section 3 we obtain expressions, in terms of an arbitrary harmonic function, for the most important quantities characterizing the dynamic of the disk and halo. In Section 4 we first calculate quantities for harmonic function described by the generalized Kuzmin-disk potential. Then, we analyze the obtained results and calculate the constants of motion of the disk. Moreover, the description of the motion of a charged particle on the halo is shown in Section 5. Finally, we complete the paper with a discussion of the results in Section 6.

#### 2. Exact Solutions for Relativistic Magnetized Haloes Surrounding Thin Disks

In this section we consider the conventional treatment of rotating galaxies modelled as a stationary thin disk and, correspondingly, we associate the magnetized halo with the region surrounding the disk. This is motivated from the fact that the magnetic field is remarkably noticed on spiral galaxies and can play a fundamental role in formation of arms. To do so, we formulate the distributional Einstein-Maxwell field equations assuming axial symmetry [60]. We also suppose that the derivatives of the metric and electromagnetic potential across the disk space-like hypersurface are discontinuous. To formulate the corresponding distributional form of the Einstein-Maxwell field equations, we introduce the usual cylindrical coordinates and assume that there exists an infinitesimally thin disk located at the hypersurface . Accordingly, we identify the halo surrounding the disk with the positive and negative regions around the equatorial plane , denoted here by the superscripts “,” so that the metric and the electromagnetic potential can be written, respectively, as and . Here denotes the Heaviside distribution. As a consequence, the Ricci tensor readswhere is the Dirac distribution and with and all the quantities are evaluated at . In agreement with (1) the energy-momentum tensor and the electric current density acquire the formwhere and are, respectively, the energy-momentum tensors and electric current density of halo. Moreover, and represent the part of the energy-momentum tensor and the electric current density of disk. The energy-momentum tensor in (3a) is taken to be the sum of two distributional components, the purely electromagnetic (trace-free) part and a “material” (trace) part:where is the electromagnetic energy-momentum tensorwith and is an unknown “material” energy-momentum tensor to be obtained. Accordingly, the Einstein-Maxwell equations, in geometrized units such that , are equivalent to the system of equationswhere , and semicolon “;” in the index indicates the covariant derivative with respect to the metric. The square brackets in expressions such as denote the jump of across the surface and denotes a unitary vector in the direction normal to it. To obtain a solution of the Einstein-Maxwell equations describing a system composed of a magnetized halo surrounding a rotating thin disk we will restrict ourselves to the case where the electric potential .

To solve the Einstein-Maxwell equations (6a), (6b), (6c), and (6d) we assume the conformastationary spacetime given by the line element where depends only on and and is an arbitrary real constant. We also assume that the magnetic potential is time independent. Accordingly, by computing the Einstein tensor explicitly from line element (7) and electromagnetic energy-momentum tensor (5), we obtain for the nonzero components of the energy-momentum tensor of the halo (EMTH) :Moreover, from (6c) the nonzero components of the electric current density on the halo have the formwhere all the quantities depend on and .

The discontinuity in the -direction of and defines, respectively, the surface energy-momentum tensor (SEMT) and the surface electric current density (SECD) of the disk, more preciselywhere is the “physical measure” of length in the direction normal to the surface. Accordingly, for metric (7), the nonzero components of and are given byrespectively. Note that in (11a), (11b), (11c), and (11d) and (12a) and (12b) all the quantities are evaluated on the surface of the disk .

In order to reduce the complexity of the last field equation systems we assume that the halo’s electric current density vanishes (i.e., in (9a) and (9b)), and it turns out that the magnetic potential and the metric functions and become completely determined in terms of an arbitrary harmonic function as follows (see [54, 56] for more details): with and arbitrary constants. Since the nonzero components of the EMTD and EMTH and the electric current density directly depend on the metric functions and magnetic potential, we observe that the entire solution is determined by a single harmonic function.

#### 3. Exact Relativistic Model for Magnetized Disk-Haloes

So far, by using the inverse method and the distributional formulation of the Einstein-Maxwell equations, we have obtained the separate energy-momentum tensor of the disk and halo. In addition, we have discussed a method to determine its components in terms of an arbitrary harmonic function. Now, the behavior of the energy-momentum tensors obtained must be investigated to find what conditions must be imposed on the solutions and the parameters that appear in the disk-haloes models in such a way that it can describe reasonably physical sources. We will now study the possible features of the disk by assuming that it is possible to express its energy-momentum tensor in the canonical form: where , , and all the quantities are evaluated in . Similarly, we assume that it is possible to express the energy-momentum tensor of the halo in the canonical form: where , , and all the quantities depend on and . Consequently, we can say that the disk and halo are constituted by some mass-energy distributions described by the energy-momentum tensors (14) and (15), respectively. is the four-velocity of a certain observer. Correspondingly, , , , and are then the energy density, the isotropic pressure, the heat flux, and the anisotropic tensor on the surface of the disk. Analogously, , , , and are then the energy density, the isotropic pressure, the heat flux, and the anisotropic tensor on the halo, respectively. Thus, it is straightforward to see that for the halo we havewhere the projection tensor is defined by and all the quantities depend on and , whereas for the disk we have where all the quantities are evaluated in . It is easy to note that by choosing the angular velocity to be zero in (A.7) we have then a fluid comoving in our coordinates system. Hence, we may introduce a suitable reference frame in terms of the local observers tetrads (A.3a), (A.3b), (A.3c), and (A.3d) and (A.4a), (A.4b), (A.4c), and (A.4d) in the formwith the corresponding dual tetrad Since the SECD of the disk can be also written in the canonical form , can be interpreted as the surface electric charge density and as the “current of magnetization” of the disk. A direct calculation shows that the surface electric charge density , whereas the “current of magnetization” of the disk is given by , where, as above, denotes the jump of the -derivative of the magnetic potential across the disk and, all quantities are evaluated on the disk.

By using the results obtained in the precedent section, we can write the surface energy density of the disk and the energy density of the halo can written asrespectively. Moreover, we have a barotropic equation of state on the surface of the disk, which can be given by , with , in such a way that we can use the energy conditions and the causality requirement for the speed of sound on the disk to obtain the physical range of possible values of the parameter . Analogously, the pressure of the halo is as follow: , where in such a way that the pressure of the halo depends not only on the energy density but also on the gravitational and magnetic fields through the function . The heat function of the disk is given bySimilarly, by inserting (15) into (16a), (16b), (16c), and (16d) we obtain for the heat flux of the haloThe nonzero components of the anisotropic tensor of the disk read , where Moreover, it is easy to see that the anisotropic tensor of the halo reads where, , as in the usual manner andNotice that and, consequently, the trace . We have obtained expressions for the energy, pressure, and the other quantities characterizing the dynamic of the halo. All the dynamic quantities have been expressed in terms of an arbitrary harmonic function. Finally, as we know, the electric current density of the halo is zero whereas it is easy to note that the magnetization current density on surface of the disk is

It is important to remark that is a defining constant in (23) and (24). Indeed, when the heat flux functions and vanish, a feature of the static systems. Because we used the inverse method, no “a priori” restriction is imposed on the physical properties of the material constituting the disk and halo. The nonzero components of the energy-momentum tensors of the disk and halo result from “the nature” of the chosen metric and the corresponding solutions. So, in our case, the nonzero components and are conditioned by the parameter and the metric function in such a way that when the component vanishes, whereas when vanishes. The decomposition of the energy-momentum tensor of the disk-halo system into (14) and (15) was chosen with the aim of describing the SEMT and EMTH by the more general fluid model. Hence, the heat flux appears here in a “natural” way as a function determined by the metric function and, consequently, by the “rotation.” Unfortunately, as we can see from (17c) and (16c), this function is oriented along the closed circular orbits and thus its physical interpretation is unclear. It is an issue that remains unanswered in this paper but should be addressed in the future.

#### 4. Rotating Kuzmin-Like Disk with Magnetized Haloes

As an example of application of the formalism described in the precedent sections, we now consider the magnetized haloes surrounding the rotating disks generated by a generalization of the Kuzmin-disk potential in the form [57, 61] where is the Legendre polynomial in cylindrical coordinates that was derived in the present form by a direct comparison of the Legendre polynomial expansion of the generating function with a Taylor series expansion of [62], being and arbitrary constant coefficients. The corresponding magnetic potential, obtained from (13a) and (13b), is where we have imposed in order to preserve the regularity of the axis of symmetry. Next, to introduce the corresponding discontinuity in the first-order derivatives of the metric potential and the magnetic potential required to define the disk we perform the transformation . It is worth noting that, for the two first members of the family of solutions ( and ) we haverespectively, where and whereas and ; moreover . For the two first members of the family of solutions the surface energy density of the disks can be written asrespectively. Similarly, for the two first members of the family we have the heat flux of the disksand the corresponding anisotropic tensorAs we know, other quantities are and . In the last expressions we have used the dimensionless expressions , . Finally, for the two first members of the family we have the current of magnetization aswhere and we first have assumed that the -derivative of the magnetic potential presents a finite discontinuity through the disk. In fact, as we have said above, the derivatives of and are continuous functions across the surface of the disk. We artificially introduce the discontinuity through the transformation .

It is worth noticing that the mass surface density as well as the isotropic pressure of the disk decays very rapidly (as and for and , resp.) indicating that the above solution can be associated with a disk with a finite energy-momentum distribution. In every case, the characteristic size can be adjusted through the parameters and of the solutions. Moreover, a simple calculation of the curvature invariants reveals that the solution is asymptotically flat and singularity-free.

To illustrate the results corresponding to the principal quantities describing the halo in Figure 1, we show the behavior of energy densities on the halo as a function of and . In each case, we plot (Figure 1(a)) and (Figure 1(b)) for the indicated values of the parameters. It can be seen that the energy density is everywhere positive and vanishes sufficiently fast as increases.