-
Magnetized tori with magnetic polarization around Kerr black holes: Variable angular momentum discs
Authors:
Sergio Gimeno-Soler,
Oscar. M. Pimentel,
Fabio D. Lora-Clavijo,
Alejandro Cruz-Osorio,
José A. Font
Abstract:
Analytical models of magnetized, geometrically thick disks are relevant to understand the physical conditions of plasma around compact objects and to explore its emitting properties. This has become increasingly important in recent years in the light of the Event Horizon Telescope observations of Sgr A$^*$ and M87. Models of thick disks around black holes usually consider constant angular momentum…
▽ More
Analytical models of magnetized, geometrically thick disks are relevant to understand the physical conditions of plasma around compact objects and to explore its emitting properties. This has become increasingly important in recent years in the light of the Event Horizon Telescope observations of Sgr A$^*$ and M87. Models of thick disks around black holes usually consider constant angular momentum distributions and do not take into account the magnetic response of the fluid to applied magnetic fields. We present a generalization of our previous work on stationary models of magnetized accretion disks with magnetic polarization (Pimentel et al. 2018). This extension is achieved by accounting for non-constant specific angular momentum profiles, done through a two-parameter ansatz for those distributions. We build a large number of new equilibrium solutions of thick disks with magnetic polarization around Kerr black holes, selecting suitable parameter values within the intrinsically substantial parameter space of the models. We study the morphology and the physical properties of those solutions, finding qualitative changes with respect to the constant angular momentum tori solutions (Pimentel et al. 2018). However, the dependences found on the angular momentum distribution or on the black hole spin do not seem to be strong. Some of the new solutions, however, exhibit a local maximum of the magnetization function, absent in standard magnetized tori. Due to the enhanced development of the magneto-rotational instability as a result of magnetic susceptibility, those models might be particularly well-suited to investigate jet formation through general-relativistic MHD simulations. The new equilibrium solutions reported here can be used as initial data in numerical codes to assess the impact of magnetic susceptibility in the dynamics and observational properties of black hole-thick disk systems.
△ Less
Submitted 8 July, 2024; v1 submitted 28 March, 2023;
originally announced March 2023.
-
OSIRIS: A New Code for Ray Tracing Around Compact Objects
Authors:
Velásquez-Cadavid J. M.,
Arrieta-Villamizar J. A.,
F. D. Lora-Clavijo,
O. M. Pimentel,
Osorio-Vargas J. E
Abstract:
The radiation observed in quasars and active galactic nuclei is mainly produced by a relativistic plasma orbiting close to the black hole event horizon, where strong gravitational effects are relevant. The observational data of such systems can be compared with theoretical models to infer the black hole and plasma properties. In the comparison process, ray tracing algorithms are essential to compu…
▽ More
The radiation observed in quasars and active galactic nuclei is mainly produced by a relativistic plasma orbiting close to the black hole event horizon, where strong gravitational effects are relevant. The observational data of such systems can be compared with theoretical models to infer the black hole and plasma properties. In the comparison process, ray tracing algorithms are essential to computing the trajectories followed by the photons from the source to our telescopes. In this paper, we present OSIRIS: a new stable FORTRAN code capable of efficiently computing null geodesics around compact objects, including general relativistic effects such as gravitational lensing, redshift, and relativistic boosting. The algorithm is based on the Hamiltonian formulation and uses different integration schemes to evolve null geodesics while tracking the error in the Hamiltonian constrain to ensure physical results. We found from an error analysis that the integration schemes are all stable, and the best one maintains an error below $10^{-11}$. Particularly, to test the robustness and ability of the code to evolve geodesics in curved spacetime, we compute the shadow and Einstein rings of a Kerr black hole with different rotation parameters and obtain the image of a thin Keplerian accretion disk around a Schwarzschild black hole. Although OSIRIS is parallelized neither with MPI nor with CUDA, the computation times are of the same order as those reported by other codes with these types of parallel computing platforms.
△ Less
Submitted 31 January, 2022;
originally announced February 2022.
-
Magneto-rotational instability in magnetically polarized discs
Authors:
Oscar M. Pimentel,
P. Chris Fragile,
F. D. Lora-Clavijo,
Bridget Ierace,
Deepika Bollimpalli
Abstract:
The magneto-rotational instability (MRI) is the most likely mechanism for transportation of angular momentum and dissipation of energy within hot, ionized accretion discs. This instability is produced through the interactions of a differentially rotating plasma with an embedded magnetic field. Like all substances in nature, the plasma in an accretion disc has the potential to become magnetically p…
▽ More
The magneto-rotational instability (MRI) is the most likely mechanism for transportation of angular momentum and dissipation of energy within hot, ionized accretion discs. This instability is produced through the interactions of a differentially rotating plasma with an embedded magnetic field. Like all substances in nature, the plasma in an accretion disc has the potential to become magnetically polarized when it interacts with the magnetic field. In this paper, we study the effect of this magnetic susceptibility, parameterized by $χ_m$, on the MRI, specifically within the context of black hole accretion. We find from a linear analysis within the Newtonian limit that the minimum wavelength of the first unstable mode and the wavelength of the fastest growing mode are shorter in paramagnetic ($χ_m>0$) than in diamagnetic ($χ_m<0$) discs, all other parameters being equal. Furthermore, the magnetization parameter (ratio of gas to magnetic pressure) in the saturated state should be smaller when the magnetic susceptibility is positive than when it is negative. We confirm this latter prediction through a set of numerical simulations of magnetically polarized black hole accretion discs. We additionally find that the vertically integrated stress and mass accretion rate are somewhat larger when the disc is paramagnetic than when it is diamagnetic. If astrophysical discs are able to become magnetically polarized to any significant degree, then our results would be relevant to properly interpreting observations.
△ Less
Submitted 24 May, 2021;
originally announced May 2021.
-
On the linear and non linear evolution of the RMHD Kelvin Helmholtz instability in a magnetically polarized fluid
Authors:
Oscar M. Pimentel,
F. D. Lora-Clavijo
Abstract:
The origin of the magnetic field in some systems like AGNs or GRBs is still an open question in astrophysics. A possible mechanism to explain the magnetic field generation is the Kelvin-Helmholtz instability, since it is able to transform the kinetic energy, in a shear flow, into magnetic energy. Through the present work, we investigate the linear and non linear effects produced by the magnetic su…
▽ More
The origin of the magnetic field in some systems like AGNs or GRBs is still an open question in astrophysics. A possible mechanism to explain the magnetic field generation is the Kelvin-Helmholtz instability, since it is able to transform the kinetic energy, in a shear flow, into magnetic energy. Through the present work, we investigate the linear and non linear effects produced by the magnetic susceptibility in the development of the Kelvin-Helmholtz instability in a relativistic plasma. The magnetic field in the system is parallel to the flows and the susceptibility is assumed to be homogeneous, constant in time, and equal in both fluids. In particular, we analyze the instability in three different cases, when the fluids are diamagnetic, paramagnetic, and when the susceptibility is zero. We compute the dispersion relation in the linear regime and found that the interface between diamagnetic fluids is more stable than the original case without magnetic susceptibility. On the contrary, the paramagnetism makes the interface more unstable. This behavior is increasingly important the higher the magnetization parameter. Additionally, the interval of relativistic Mach numbers for which the interface is unstable grows with the magnetic susceptibility, i.e., it is larger for paramagnetic fluids. We check these analytical results with numerical simulations, and explore the effect of the magnetic polarization in the non linear regime. We found that the magnetic field is more amplified in paramagnetic fluids than in diamagnetic ones. Surprisingly, the effect of the susceptibility in the amplification is stronger when the magnetization parameter is smaller. These results make the KH instability a more efficient and effective amplification mechanism of seed magnetic fields when considering the susceptibility of matter.
△ Less
Submitted 25 October, 2019;
originally announced October 2019.
-
Analytic solution of a magnetized tori with magnetic polarization around Kerr black holes
Authors:
Oscar M. Pimentel,
F. D. Lora-Clavijo,
Guillermo A. Gonzalez
Abstract:
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 t…
▽ More
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 non-constant $χ_{m}$, the tori become more magnetized than the Komissarov solution in the region where $\partialχ_{m}/\partial r<0$, and less magnetized when $\partialχ_{m}/\partial 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.
△ Less
Submitted 22 August, 2018;
originally announced August 2018.
-
Numerical general relativistic MHD with magnetically polarized matter
Authors:
Oscar M. Pimentel,
F. D. Lora-Clavijo,
Guillermo A. González
Abstract:
The magnetically polarized matter in astrophysical systems may be relevant in some magnetically dominated regions. For instance, the funnel that is generated in some highly magnetized disks configurations whereby relativistic jets are thought to spread, or in pulsars where the fluids are subject to very intense magnetic fields. With the aim of dealing with magnetic media in the astrophysical conte…
▽ More
The magnetically polarized matter in astrophysical systems may be relevant in some magnetically dominated regions. For instance, the funnel that is generated in some highly magnetized disks configurations whereby relativistic jets are thought to spread, or in pulsars where the fluids are subject to very intense magnetic fields. With the aim of dealing with magnetic media in the astrophysical context, we present for the first time the conservative form of the GRMHD equations with a non-zero magnetic polarization vector $m^μ$. Then, we follow the Anile method to compute the eigenvalue structure in the case where the magnetic polarization is parallel to the magnetic field, and it is parametrized by the magnetic susceptibility $χ_m$. This approximation allows us to describe diamagnetic fluids, for which $χ_m<0$, and paramagnetic fluids where $χ_m>0$. The theoretical results were implemented in the CAFE code to study the role of the magnetic polarization in some 1D Riemann problems. We found that independently of the initial condition, the first waves that appear in the numerical solutions are faster in diamagnetic materials than in paramagnetic ones. Moreover, the constant states between the waves change notably for different magnetic susceptibilities. All these effects are more appreciable if the magnetic pressure is much bigger than the fluid pressure. Additionally, with the aim of analysing a magnetic media in a strong gravitational field, we carry out for the first time the magnetized Michel accretion of a magnetically polarized fluid. With this test, we found that the numerical solution is effectively maintained over time ($t>4000$), and that the global convergence of the code is $\gtrsim$ 2 for $χ_{m}\lesssim 0.005$, for all the magnetic field strength $β$ we considered.
△ Less
Submitted 6 June, 2018;
originally announced June 2018.