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.

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 $\chi_{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 $\chi_{m}$, we find that the paramagnetic tori ($\chi_{m}>0$) are more dense and less magnetized than the diamagnetic ones ($\chi_{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 $\chi_{m}$, the tori become more magnetized than the Komissarov solution in the region where $\partial\chi_{m}/\partial r<0$, and less magnetized when $\partial\chi_{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.

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^{\mu}$. 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 $\chi_m$. This approximation allows us to describe diamagnetic fluids, for which $\chi_m<0$, and paramagnetic fluids where $\chi_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 $\chi_{m}\lesssim 0.005$, for all the magnetic field strength $\beta$ we considered.