Issue 
A&A
Volume 619, November 2018



Article Number  A57  
Number of page(s)  7  
Section  Astrophysical processes  
DOI  https://doi.org/10.1051/00046361/201833736  
Published online  07 November 2018 
Analytic solution of a magnetized tori with magnetic polarization around Kerr black holes
Grupo de Investigación en Relatividad y Gravitación, Escuela de Física, Universidad Industrial de Santander, A. A. 678, Bucaramanga, 680002 Colombia
email: fadulora@uis.edu.co
Received:
28
June
2018
Accepted:
3
August
2018
We present the first family of magnetically polarized equilibrium tori around a Kerr black hole. The models were obtained in the test fluid approximation by assuming that the tori is a linear media, making it is possible to characterize the magnetic polarization of the fluid through the magnetic susceptibility χ_{m}. The magnetohydrodynamic (MHD) structure of the models was solved by following the Komissarov approach, but with the aim of including the magnetic polarization of the fluid, the integrability condition for the magnetic counterpart was modified. We build two kinds of magnetized tori depending on whether the magnetic susceptibility is constant in space or not. In the models with constant χ_{m}, we find that the paramagnetic tori ( χ_{m} > 0) are more dense and less magnetized than the diamagnetic ones ( χ_{m} < 0) in the region between the inner edge, r_{in}, and the center of the disk, r_{c}; however, we find the opposite behavior for r > r_{c}. Now, in the models with nonconstant χ_{m}, the tori become more magnetized than the Komissarov solution in the region where ∂χ_{m}/∂r < 0, and less magnetized when ∂χ_{m}/∂r > 0. Nevertheless, it is worth mentioning that in all solutions presented in this paper the magnetic pressure is greater than the hydrodynamic pressure. These new equilibrium tori can be useful for studying the accretion of a magnetic media onto a rotating black hole.
Key words: accretion / accretion disks / magnetohydrodynamics (MHD) / black hole physics
© ESO 2018
1. Introduction
The accretion of a fluid onto a rotating black hole is believed to be the most powerful Xray source in the universe, causing it to be an active research topic in astrophysics (Frank et al. 2002). Now, due to the strong gravitational field of the central object, it is necessary to consider a general relativistic approach to properly describe the dynamics of the fluid in the vicinity of the black hole. Additionally, one of the basic features in many theoretical models and numerical simulations of accretion disks is the presence of a strong magnetic field. This field has turned out to be very important because it interacts with the differential rotation of the disk and generates turbulence via the magnetorotational instability (Balbus & Hawley 1991). This turbulence provides the necessary viscous stress to transfer angular momentum, dissipate energy, and therefore generate the accretion process (Balbus & Hawley 1998). Moreover, it is commonly assumed that strong poloidal magnetic fields are necessary to collimate and accelerate relativistic jets (Blandford & Znajek 1977). However, the origin and strength of those fields in the accretion disks are still a subject of study. Recent research has shown that a magnetically supported disk may be obtained from a thermally unstable hot disk (Fragile & Meier 2009). Another possibility is that the magnetic field increases gradually while it is dragged from the outer region to the inner region, near the black hole (McKinney et al. 2012).
It is well known that all substances contain spinning electrons that orbit around the nucleus, meaning we can model the microscopic structure of matter as an assembly of small dipoles, characterized macroscopically by the magnetization vector (Lorrain & Corson 1970). This vector enters in the Maxwell equations as a source of magnetic field, and could be important to understand the problems concerning the origin of the magnetic fields in astrophysical scenarios. Magnetization has already been considered for studying the equilibrium structure of neutron stars. For instance, in Blandford & Hernquist (1982) the authors compute the magnetic susceptibility of the star crust and find that magnetization does not considerably change the surface properties, but may be connected to observable effects. Indeed, Suh & Mathews (2010) consider the possibility that the soft gammaray repeaters and the anomalous Xray pulsars might be observational evidence for a diamagnetic phase transition that results in a domain formation. Furthermore, as mentioned by Wang et al. (2016), neutron stars are also important for testing the Haasvan Alphen effect in which the magnetic susceptibility oscillates when the applied magnetic field is increased (De Haas & Van Alphen 1930).
In the theory of galactic systems, NavarroNoguera et al. (2018) presented a static and axially symmetric selfgravitating thin disk in which the magnetic field is generated by a magnetization vector that is normal to the disk plane, presents a maximum at the disk center, and goes to zero at infinity. In accretion disk theory the magnetic polarization of the fluid has not yet been considered. Nevertheless, the first advances in this direction were made by Pimentel et al. (2018), where we presented the theoretical and numerical background to describe the evolution of a magnetically polarized fluid in a gravitational field. Among some results of this work, we can mention that the propagation speed of the fastest waves in the onedimensional (1D) Riemann problems is greater in diamagnetic materials than in paramagnetic ones, and that the magnetic field and the relativistic character of the flows may increase considerably with the magnetic susceptibility.
Analytical models of accretion disks are of great interest in numerical simulations because they are used as initial data to study the nonlinear evolution of the fluid (Font 2008; Abramowicz & Fragile 2013). Among the analytical models, the Polish doughnuts (Abramowicz et al. 1978) are the most used in numerical simulations due to their simplicity and the low computational expense they require (Abramowicz & Fragile 2013). These equilibrium models consist of a nonselfgravitating barotropic tori orbiting around a Kerr black hole. Later, Komissarov (2006) followed the method of Abramowicz to compute a magnetized tori with a purely toroidal magnetic field. These models have recently been used by Fragile & Sdowski (2017) to show that a strong toroidal magnetic field cannot be maintained during the disk evolution. Moreover, Bugli et al. (2018) used the Komissarov models to simulate the interplay of the Papaloizou–Pringle instability and the magneto rotational instability. Additionally, it is worth mentioning that GimenoSoler & Font (2017), following the Komissarov procedure, found new magnetized equilibrium tori with a nonconstant angular momentum distribution in the disk.
The analytic solutions for magnetized tori presented in the literature do not take into account the contribution of the magnetic dipoles in the torus. We therefore generalize in this work the KomisSarov models by including the magnetic polarization of the matter. For this purpose, we consider the energymomentum tensor for a magnetically polarized fluid that is given in Maugin (1978), Huang et al. (2010), Chatterjee et al. (2015), and follow the Komissarov approach (Komissarov 2006) to solve the hydrodynamic structure of the tori. These new models may be useful as initial data to study the evolution of a magnetic media in the gravitational field of a Kerr black hole. The organization of this paper is the following: in Sect. 2 we present the Euler equations that describe the equilibrium structure of a magnetically polarized tori endowed with a toroidal magnetic field around a rotating black hole. To write these equations we concentrate on the special case in which the magnetization vector is parallel to the magnetic field, so that we can relate both vectors through the magnetic susceptibility, χ_{m}. In Sect. 3 we present the conditions to write the Euler equations as an exact differential; then, as in Komissarov (2006), we assume a barotropic tori in which the angular velocity is a function of the specific angular momentum only. However, the condition related to the magnetic counterpart is modified in order to include the magnetic polarization of the fluid. In Sect. 4 we present two kinds of magnetized tori: when the magnetic susceptibility is constant in space (Sect. 4.1), and when it changes with coordinates (Sect. 4.2). Finally, the main results are presented in Sect. 5. In this paper we use the signature (−, +, +, +) and geometrized units, for which G = c = 1.
2. Magnetohydrodynamic equations with magnetic polarization
The dynamics of an ideal fluid with magnetic polarization in a magnetic field is described by the conservation laws
and by the relevant Maxwell equations,
where ρ is the rest mass density and b^{μ } is the magnetic field, both measured in a reference frame that moves with the same fourvelocity as the fluid, u^{μ}. The magnetic polarization can then be characterized macroscopically through the magnetization vector m^{μ}, which is defined as the magnetic dipole moment per unit volume. From now on, we concentrate on the physically important case in which m^{μ} and b^{μ} are related by means of the linear constitutive equation, m^{μ} = χb^{μ}, where χ = χ_{m}/(1 + χ_{m}), and χ_{m} is the magnetic susceptibility. When χ_{m} < 0 the fluid is diamagnetic and when χ_{m} > 0 the fluid is paramagnetic. In the first case the polarization results from induced orbital dipole moments in a magnetic field (Griffiths 2005), and in the second case the polarization is generated by magnetic torques in substances whose atoms have a nonzero spin dipole moment. The total energymomentum tensor, T^{αβ}, for a magnetically polarized fluid was computed in Maugin (1978) and more recently in Chatterjee et al. (2015) by following a different approach. The resulting tensor takes the form
in the linear media approximation. Here, w is the enthalpy density, p is the thermodynamic pressure, b^{2} = b^{μ}b_{μ}, and g^{μν} is the metric tensor.
Following previous works (Komissarov 2006; Wielgus et al. 2015; GimenoSoler & Font 2017), we now assume the test fluid approximation, and the gravitational field as given by the Kerr metric in Boyer–Lindquist coordinates (t, ϕ, r, θ). We also consider that the fluid is axisymmetric and stationary, so the physical variables do not depend either on the azimuthal angle ϕ or on the time t. Finally, we restrict the movement of the fluid in such a way that u^{r} = u^{θ} = 0, and we restrict the magnetic field topology to a purely toroidal one, so b^{r} = b^{θ} = 0. With these assumptions, the baryon number conservation Eq. (2) and the relevant Maxwell Eq. (3) are identically satisfied, and the equilibrium structure of the tori is obtained from the Euler equations , where is the projection tensor. This contraction leads to the following expression.
where i = r, θ.
Nevertheless, it is useful to write this last equation in terms of the angular velocity
and the specific angular momentum
in such a way that the Euler equations take the form
where p_{m} = b^{2}/2 and . It is important to mention that when χ = 0, Eq. (8) reduces to Eq. (14) in Komissarov (2006), where the tori do not have magnetic polarization.
3. Integrability conditions
In order to write the first three terms in Eq. (8) as an exact differential, we follow the procedure used by Komissarov (2006) where it is supposed that the fluid obeys a barotropic equation of state w = w(p), and that the surfaces of the Ω and l constants coincide in such a way that Ω = Ω(l) (Abramowicz et al. 1978). With these two assumptions, Eq. (8) takes the form
where
We note that in this equation, we take χ to be an arbitrary function of the coordinates. Nevertheless, with the aim of writing ℐ as a exact differential, we assume that χ = χ(ℒ). In this way, dχ = (∂χ/∂ℒ)dℒ, and the last equation reduces to
where χ′=∂χ/∂ℒ. Now, if we can find a function z = z(p_{m}, ℒ) that satisfies the derivatives
such that ℐ = (∂z/∂p_{m})dp_{m} + (∂z/∂ℒ)dℒ, then Eq. (11) could also be written as an exact differential ℐ = dz.
By integrating Eq. (12) we obtain the general solution
if the condition
is satisfied, C_{0} being a constant of integration, and ℱ(p_{m}) a function of the magnetic pressure only. Thus, the problem reduces to find a relation between w, p_{m}, and ℒ that satisfies Eq. (14). Fortunately, it is possible to transform this condition into a partial differential equation for the enthalpy by differentiating it, first with respect to p_{m}, and then with respect to ℒ,
We solve this equation with the method of the characteristics and find that any function of the form
is a solution to the equation. In this last expression,
and f is an arbitrary function of its argument, so we can obtain different solutions depending on the particular choice of f.
In order to find magnetically polarized tori that reduce to the Komissarov solution when χ = 0, we now assume f as the function
where and η are arbitrary constants. We can then obtain the magnetic pressure from Eq. (16) in such a way that
with . Additionally, we need to choose a particular relation between χ and ℒ to compute ϕ from Eq. (17) and completely determine p_{m}. Subsequently, with the aim of analyzing disks with different magnetic polarization states, we adopt the form
where χ_{0}, χ_{1}, and α are constants. Therefore, depending on the choice of these constants we obtain a different magnetic susceptibility. When χ_{1} = 0, this susceptibility is constant, when χ_{1} ≠ 0 it becomes a function of the spatial coordinates, and when χ_{1} = χ_{0} = 0 it is zero and then we reduce to the Komissarov solution.
Once the function χ has been defined, the magnetic pressure in Eq. (19) takes the simple form
with,
Finally, with the enthalpy computed from Eq. (21), we solve the integrals in Eq. (13) to obtain
as a particular solution of Eq. (12) in which χ is given by Eq. (20).
With the function z(p_{m}, ℒ) that allows us to write ℐ = dz, the Euler equations in (9) can be solved in the form
where the integrals are computed, as in Komissarov (2006), by assuming a torus with constant angular momentum l = l_{0}, and a polytropic equation of state for the gas pressure,
where K and κ are constants. Therefore, Eq. (24) becomes
On the other hand, it is usual to introduce the relativistic effective potential (Abramowicz et al. 1978),
which reduces to W = lnu_{t} in the case of a constant angular momentum torus. With W, the Euler equations in the form (26) can be expressed as
with W_{in} = lnu_{tin} being the potential defined at the inner edge of the disk, where in turn p_{in} = 0 and p_{min} = 0. It is worth mentioning that when the magnetic polarization of the fluid is zero (χ = 0), Eqs. (21) and (28) in this paper reduce to Eqs. (29) and (30) in Komissarov (2006). Therefore, the new magnetized tori with magnetic polarization reduce to the Komissarov solution in the case where χ = 0.
4. Results
The new analytic solutions obtained in the previous sections have the following free parameters: the potential at the inner edge of the disk W_{in}, the angular momentum of the fluid l_{0}, the exponents of the pressures κ and η, the magnetic polarization parameters χ_{0}, χ_{1} and α, and finally, the enthalpy and the magnetization parameter at the disk center, w_{c} and β_{c} = p_{c}/p_{mc}, respectively. The disk center, r_{c}, is defined as one of the points where the gradient pressure vanishes, so the angular momentum l_{0} equals the Keplerian angular momentum,
where a is the angular momentum of the black hole, and the upper signs correspond to prograde motion. We highlight the fact that in these formulae we have chosen M = 1, M being the mass of the black hole. By doing l_{k} = l_{0} we therefore compute r_{c}.
The procedure to compute the physical variables consists in finding the gas pressure and the magnetic pressure at the disk center through the equations
from which we can compute the constants K and K_{m}. Next, the enthalpy as a coordinate function can be found from Eq. (28). With the enthalpy, the magnetic pressure is determined by Eq. (21), and the gas pressure by Eq. (25). The magnetic field components can be obtained from the constraint u^{μ}b_{μ} = 0, and the specific angular momentum l = −u_{ϕ}/u_{t}, which lead to the following expressions.
with . Finally, the fourvelocity is computed from,
where .
Since we are interested in analyzing the contribution of the magnetic polarization to the physical variables of the fluid, we now compare our new disk solutions with the model A proposed by Komissarov (2006), in which the disk has a constant angular momentum of l_{0} = 2.8 and an inneredge effective potential of W_{in} = −0.030. These parameters, along with the black hole’s spin parameter a = 0.9, correspond to a disk with finite outer radius, whose center is located at r_{c} = 4.62. Finally, in the Komissarov model κ = η = 4/3, β_{c} = 0.1, and w_{c} = 1. The three remaining parameters χ_{0}, χ_{1}, and α determine the magnetic polarization state of the disk, and are analyzed in the following two sections. The first of these sections is dedicated to magnetized tori with constant magnetic susceptibility, while in the second section we discuss the effect of a nonconstant χ_{m} on the equilibrium state of the disks.
4.1. Magnetized disks with constant magnetic susceptibility
In many cases, it is reasonable to consider that the spatial variations of the magnetic susceptibility are very small, so we can determine the magnetic polarization of a material through a single constant. In our models, a disk with such an approximation can be obtained by taking χ_{1} = 0, which in turn implies that χ = χ_{0}. We construct two diamagnetic disks with magnetic susceptibilities of χ_{m} = −0.2, −0.4, which correspond to χ_{0} ≈ −0.25, −0.67, and two paramagnetic disks with χ_{m} = 0.2, 0.4, which imply that χ_{0} ≈ 0.17, 0.29. These models will be compared with the Komissarov solution, which corresponds to the case χ_{0} = 0 (χ_{m} = 0).
In Fig. 1 we present the spatial behavior of the rest mass density for different values of magnetic susceptibility. The topleft panel is the radial profile of ρ in the equatorial plane θ = π/2, while the topcenter and topright panels show the angular behavior of the same physical variable at r = 3.0 (close to the inner edge of the diks) and r = 10, respectively. In these three panels, the black curve corresponds to the Komissarov solution. The bottom row shows the spatial distribution of the rest mass density in the meridional plane of the disk when the fluid is diamagnetic with χ_{m} = −0.4 and paramagnetic with χ_{m} = 0.4. We also plot the Komissarov disk (χ_{m} = 0.0) for comparison purposes. From this figure, we note that the magnetic susceptibility considerably changes the way in which the matter is distributed within the disk. More specifically, the paramagnetic tori are more dense than the tori without magnetic polarization in the region between the inner edge, r_{in}, and the center of the disk, r_{c}. On the contrary, the diamagnetic disks are less dense than the Komissarov tori in this region. Nevertheless, in r > r_{c} the rest mass density is greater in the diamagnetic disks than in those with paramagnetic properties. Additionally, we can see from the contour plots that the paramagnetic tori are more compact than the Komissarov solution, and even more compact than the diamagnetic ones.
Fig. 1.
Spatial distribution of the rest mass density ρ for different values of magnetic susceptibility χ_{m}. The topleft panel is the radial dependence of ρ on the equatorial plane defined by θ = π/2. The topcenter and topright panels correspond to the angular dependence of ρ in the radial positions r = 3.0 and r = 10.0, respectively. In the top row the black curves correspond to the Komissarov solution, obtained by taking χ_{m} = 0. The bottom row shows the density distribution of ρ in the zR plane, with R = r sin θ and z = r cos θ, for χ_{m} = −0.4, 0.0, 0.4. 
On the other hand, in Fig. 2 we show the magnetization parameter β for different values of magnetic susceptibility. The left panel describes the behavior of β on the equatorial plane, while the center and right panels describe the magnetization parameter at r = 3.0 and r = 10.0, respectively. Additionally, the dotted vertical lines are the edges of the disk, whose positions can be computed from W = W_{in}. This figure shows that when the diamagnetic character of the disks increases (a more negative value of χ_{m}), they become more magnetized in the region r < r_{c} and less magnetized in r > r_{c}, as compared to the Komissarov solution. The paramagnetic tori exhibit the opposite behavior: they are less magnetized in r < r_{c} and more magnetized in r > r_{c} than the torus without magnetic polarization.
Fig. 2.
Magnetization parameter β as a function of the spatial coordinates for different values of magnetic susceptibility χ_{m}. The left panel describes the radial dependence of β on the equatorial plane (θ = π/2). The middle and right panels show the angular dependence of β at the position r = 3.0 and r = 10.0, respectively. The vertical dotted lines represent the boundaries of the disk. 
4.2. Magnetized disks with nonconstant magnetic susceptibility
As we mentioned at the end of Sect. 2, spatial variations of χ_{m} may affect the stationary state of the disk, so it is also interesting to consider a fluid with a nonconstant magnetic susceptibility. This kind of system has its motivation in some physical phenomena; for instance, the temperature dependence of χ_{m} in paramagnetic materials (Gabold et al. 2018), the rapid decrease of the magnetic susceptibility in a ferromagneticparamagnetic phase transition (Arajs & Colvin 1964; Chevalier et al. 2002), and the interaction between the spin of a particle and the gravitational field. This last effect is described by the MathissonPapapetrouDixon equations (Papapetrou 1951; Dixon 1970), and has been studied in the kerr spacetime, for instance in Saijo et al. (1998) and Semerák (1999).
With the aim of analyzing the stationary state of a fluid with a nonconstant magnetic susceptibility, we construct two diamagnetic disk models M1 and M2, and two paramagnetic ones M3 and M4. In Table 1 we present the parameters χ_{0}, χ_{1}, and α that define the magnetic susceptibility (20) in each model. The radial and angular behavior of χ_{m} in the disk is showed in the left column of Fig. 3. We note that in the four models that we propose, the magnetic susceptibility changes considerably quickly with the radial coordinate, especially in the region between the inner edge and the center of the disk, while χ_{m} remains approximately constant in r ≳ 10. In the diamagnetic models M1 and M2, the rest mass density is reduced in the inner region of the disk (r ≲ 5) and it is higher in the outer region, in comparison with the Komissarov solution, which is labeled with letter A in Fig. 3. On the contrary, the mass density in the paramagnetic models M3 and M4 presents the opposite behavior to the diamagnetic models, with those nearer to the black hole being more dense. Nevertheless, by comparison to the profiles of ρ in Fig. 1, we note that a nonconstant magnetic susceptibility does not affect the qualitative behavior of the rest mass density.
Parameters defining the magnetic susceptibility of the disk in each model.
Fig. 3.
Magnetic susceptibility χ_{m} (left column), rest mass density ρ (middle column), and magnetization parameter β (right column) as a function of the spatial coordinates (r, θ) for the models M1, M2, M3, M4, and for the Komissarov solution, labeled by the letter A. Each one of these models is characterized by a different magnetic susceptibility function that is shown in the first column. The top row shows the radial behavior of χ_{m}, ρ, and β on the equatorial plane θ = π/2, and the bottom row shows the angular dependence of the same physical variables at the center of the disk r = r_{c}, except for the mass density, for which the behavior is showed at r = 10.0. The vertical dotted lines correspond to the boundaries of the disk. 
Now, from the right column of Fig. 3 it is clear that the gradients of χ_{m} determine the magnetization state of the disk, especially in the region between the inner edge and the center of the disk, where the spatial changes of χ_{m} are large. In this region, when ∂χ_{m}/∂r < 0, as in models M1 and M3, the disks become more magnetized than in the Komissarov solution, while in the models M2 and M4, for which ∂χ_{m}/∂r > 0, the disks considerably reduce their magnetization. Therefore, the stationary state of our tori implies that when the magnetic susceptibility decreases with the radial coordinate, the disk becomes more magnetized and vice versa.
5. Conclusions
In this paper we obtain for the first time a stationary and axisymmetric family of solutions for magnetized tori with magnetic polarization around a Kerr black hole. The magnetized tori was built by assuming a barotropic equation of state, a constant angular momentum, and a purely toroidal magnetic field. These features were previously considered by Komissarov (2006) to construct tori with dynamically important magnetic fields. In addition, the polarization of the tori was introduced by assuming the usual linear constitutive relation in which the magnetization vector and the applied magnetic field are linearly related by the magnetic susceptibility χ_{m}. We show that our models reduce to the Komissarov solution when we set the magnetic susceptibility to zero, making it possible to determine the differences in the physical variables of a tori without magnetic polarization and a paramagnetic or diamagnetic tori.
With the aim of obtaining an analytical solution of the relativistic Euler equations, we assumed that χ = χ(ℒ(r, θ)), and in particular, that χ takes the form Eq. (20). In this way, the polarization state of the tori is completely defined by choosing the constants α, χ_{0}, and χ_{1}. In this paper we present two kinds of magnetized tori, one with constant magnetic susceptibility and one where this susceptibility varies. In the models of Sect. 4.1, where χ_{m} is constant, we find that a paramagnetic torus ( χ_{m} < 0) is more compact than a diamagnetic one ( χ_{m} < 0), because the matter is more concentrated in the region between the inner edge and the center of the disk. However, in r < r_{c} the diamagnetic tori are more magnetized than those obtained by Komissarov, and therefore, more magnetized than the paramagnetic ones. The opposite behavior for the magnetization parameter is obtained in the region r > r_{c}, where the paramagnetic tori are more magnetized than the case with χ_{m} = 0.
Finally, in the models of Sect. 4.2, where the magnetic susceptibility is nonconstant in the torus, we note that the rest mass density does not change its qualitative behavior, as compared with the models of constant χ_{m}. Therefore, the way in which the matter is distributed in the tori depends on whether the fluid is diamagnetic or paramagnetic, the paramagnetic tori being more compact than the diamagnetic ones. Nevertheless, we find that the magnetization state of the tori also depends on the spatial changes of the magnetic susceptibility. In particular, when ∂χ_{m}/∂r < 0 (models M1 and M3), the tori become more magnetized in r < r_{c} than the Komissarov solution, and when ∂χ_{m}/∂r > 0 (models M2 and M4), the tori lose considerable amounts of their magnetization in the same region of the disk. All these effects are more appreciable when the changes in the magnetic susceptibility are large, as in the inner region of the tori in our models. The next step in this direction is to add this solution as initial data, in the CAFE code (LoraClavijo et al. 2015), in order to carry out numerical simulations and see the dynamics of the accretion disk once it is being accreted onto the black hole. It is worth mentioning that in the last version of CAFE the magnetic polarized matter terms were implemented following the characteristic approach presented in our recent work (Pimentel et al. 2018).
Acknowledgments
O. M. P. would like to thank the financial support from COLCIENCIAS under the program Becas Doctorados Nacionales 647 and Universidad Industrial de Santander. F. D. L.C and G. A. G. were supported in part by VIEUIS, under Grant No. 2314 and by COLCIENCIAS, Colombia, under Grant No. 8863.
References
 Abramowicz, M. A., & Fragile, P. C. 2013, Liv. Rev. Rel., 16, 1 [Google Scholar]
 Abramowicz, M., Jaroszynski, M., & Sikora, M. 1978, A&A, 63, 221 [NASA ADS] [Google Scholar]
 Arajs, S., & Colvin, R. V. 1964, J. Appl. Phys., 35, 2424 [NASA ADS] [CrossRef] [Google Scholar]
 Balbus, S. A., & Hawley, J. F. 1991, ApJ, 376, 214 [NASA ADS] [CrossRef] [Google Scholar]
 Balbus, S. A., & Hawley, J. F. 1998, Rev. Mod. Phys., 70, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Blandford, R. D., & Hernquist, L. 1982, J. Phys. C Sol. State Phys., 15, 6233 [NASA ADS] [CrossRef] [Google Scholar]
 Blandford, R. D., & Znajek, R. L. 1977, MNRAS, 179, 433 [NASA ADS] [CrossRef] [Google Scholar]
 Bugli, M., Guilet, J., Müller, E., et al. 2018, MNRAS, 475, 108 [NASA ADS] [CrossRef] [Google Scholar]
 Chatterjee, D., Elghozi, T., Novak, J., & Oertel, M. 2015, MNRAS, 447, 3785 [NASA ADS] [CrossRef] [Google Scholar]
 Chevalier, B., Kahn, M. L., Bobet, J.L., Pasturel, M., & Etourneau, J. 2002, J. Phys. Condens. Matter, 14, L365 [CrossRef] [Google Scholar]
 De Haas, W., & Van Alphen, P. 1930, Proc. Netherlands Roy. Acad. Sci, 33, 170 [Google Scholar]
 Dixon, W. G. 1970, Proc. R. Soc. Lond. A, 314, 499 [NASA ADS] [CrossRef] [Google Scholar]
 Font, J. A. 2008, Liv. Rev. Rel., 11, 7 [Google Scholar]
 Fragile, P. C., & Meier, D. L. 2009, ApJ, 693, 771 [NASA ADS] [CrossRef] [Google Scholar]
 Fragile, P. C., & Sdowski, A. 2017, MNRAS, 467, 1838 [NASA ADS] [CrossRef] [Google Scholar]
 Frank, J., King, A., & Raine, D. J. 2002, Accretion Power in Astrophysics: Third Edition, 398. [CrossRef] [Google Scholar]
 Gabold, H., Luan, Z., Paul, N., et al. 2018, Sci. Rep., 8, 4835 [NASA ADS] [CrossRef] [Google Scholar]
 GimenoSoler, S., & Font, J. A. 2017, A&A, 607, A68 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Griffiths, D. J., 2005, Introduction to Electrodynamics. [Google Scholar]
 Huang, X.G., Huang, M., Rischke, D. H., & Sedrakian, A. 2010, Phys. Rev. D, 81, 045015 [NASA ADS] [CrossRef] [Google Scholar]
 Komissarov, S. S. 2006, MNRAS, 368, 993 [NASA ADS] [CrossRef] [Google Scholar]
 LoraClavijo, F. D., CruzOsorio, A., & Guzmán, F. S. 2015, ApJS, 218, 24 [NASA ADS] [CrossRef] [Google Scholar]
 Lorrain, P., & Corson, D. 1970, Electromagnetic Fields and Waves (Freeman). [Google Scholar]
 Maugin, G. A. 1978, J. Math. Phys., 19, 1198 [NASA ADS] [CrossRef] [Google Scholar]
 McKinney, J. C., Tchekhovskoy, A., & Blandford, R. D. 2012, MNRAS, 423, 3083 [NASA ADS] [CrossRef] [Google Scholar]
 NavarroNoguera, A., LoraClavijo, F., & González, G. A. 2018, Gen. Relativ. Gravitation, 50, 76 [NASA ADS] [CrossRef] [Google Scholar]
 Papapetrou, A. 1951, Proc. R. Soc. London, Ser. A, 209, 248 [NASA ADS] [CrossRef] [Google Scholar]
 Pimentel, O. M., LoraClavijo, F. D., & González, G. A. 2018, ApJ, 861, 115 [NASA ADS] [CrossRef] [Google Scholar]
 Saijo, M., Maeda, K.I., Shibata, M., & Mino, Y. 1998, Phys. Rev. D, 58, 064005 [NASA ADS] [CrossRef] [Google Scholar]
 Semerák, O. 1999, MNRAS, 308, 863 [NASA ADS] [CrossRef] [Google Scholar]
 Suh, I.S., & Mathews, G. J. 2010, ApJ, 717, 843 [NASA ADS] [CrossRef] [Google Scholar]
 Wang, Z., Lü, G., Zhu, C., & Wu, B. 2016, PASP, 128, 104201 [NASA ADS] [CrossRef] [Google Scholar]
 Wielgus, M., Fragile, P. C., Wang, Z., & Wilson, J. 2015, MNRAS, 447, 3593 [NASA ADS] [CrossRef] [Google Scholar]
All Tables
All Figures
Fig. 1.
Spatial distribution of the rest mass density ρ for different values of magnetic susceptibility χ_{m}. The topleft panel is the radial dependence of ρ on the equatorial plane defined by θ = π/2. The topcenter and topright panels correspond to the angular dependence of ρ in the radial positions r = 3.0 and r = 10.0, respectively. In the top row the black curves correspond to the Komissarov solution, obtained by taking χ_{m} = 0. The bottom row shows the density distribution of ρ in the zR plane, with R = r sin θ and z = r cos θ, for χ_{m} = −0.4, 0.0, 0.4. 

In the text 
Fig. 2.
Magnetization parameter β as a function of the spatial coordinates for different values of magnetic susceptibility χ_{m}. The left panel describes the radial dependence of β on the equatorial plane (θ = π/2). The middle and right panels show the angular dependence of β at the position r = 3.0 and r = 10.0, respectively. The vertical dotted lines represent the boundaries of the disk. 

In the text 
Fig. 3.
Magnetic susceptibility χ_{m} (left column), rest mass density ρ (middle column), and magnetization parameter β (right column) as a function of the spatial coordinates (r, θ) for the models M1, M2, M3, M4, and for the Komissarov solution, labeled by the letter A. Each one of these models is characterized by a different magnetic susceptibility function that is shown in the first column. The top row shows the radial behavior of χ_{m}, ρ, and β on the equatorial plane θ = π/2, and the bottom row shows the angular dependence of the same physical variables at the center of the disk r = r_{c}, except for the mass density, for which the behavior is showed at r = 10.0. The vertical dotted lines correspond to the boundaries of the disk. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext 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 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.