In this paper, we study the effects of polynomial $f(R)$ model on the stability of homogeneous energy density in self-gravitating spherical stellar object. For this purpose, we construct couple of evolution equations which relate the Weyl tensor with matter parameters. We explore different factors responsible for density inhomogeneities with non-dissipative dust, isotropic as well as anisotropic fluids and dissipative dust cloud. We find that shear, pressure, dissipative parameters and $f(R)$ terms affect the existence of inhomogeneous energy density.

This paper is devoted to investigate the cold plasma wave properties outside the event horizon of the Schwarzschild planar analogue. The dispersion relations are obtained from the corresponding Fourier analyzed equations for non-rotating and rotating, non-magnetized and magnetized backgrounds. These dispersion relations provide complex wave numbers. The wave numbers are shown in graphs to discuss the nature and behavior of waves and the properties of plasma lying in the vicinity of the Schwarzschild event horizon.

The behavior of isothermal plasma waves has been analyzed near the Schwarzschild horizon. We consider a non-rotating background with non-magnetized and magnetized plasmas. The general relativistic magnetohydrodynamical equations for the Schwarzschild planar analogue spacetime with an isothermal state of the plasma are formulated. The perturbed form of these equations is linearized and Fourier analyzed by introducing simple harmonic waves. The determinant of these equations in each case leads to a complex dispersion relation, which gives complex values of the wave number. This has been used to discuss the nature of the waves and their characteristics near the horizon.

We investigate the wave properties of the Kerr black hole with isothermal plasma using 3+1 ADM formalism. The corresponding Fourier analyzed perturbed GRMHD equations are used to obtain the dispersion relations. These relations lead to the real values of the components of wave vector $\textbf{k}$ which are used to evaluate the quantities like phase and group velocities etc. These have been discussed graphically in the neighborhood of the pair production region. The results obtained verify the conclusion of Mackay et al. according to which rotation of a black hole is required for negative phase velocity propagation.

The 3+1 GRMHD equations for Schwarzschild spacetime in Rindler coordinates with isothermal state of plasma are formulated. We consider the cases of non-rotating and rotating backgrounds with non-magnetized and magnetized plasmas. For these cases, the perturbed form of these equations are linearized and Fourier analyzed by introducing plane wave type solutions. The determinant of these equations in each case leads to two dispersion relations which give value of the wave number $k$. Using the wave number, we obtain information like phase and group velocities etc. which help to discuss the nature of the waves and their characteristics. These provide interesting information about the black hole magnetosphere near the horizon. There are cases of normal and anomalous dispersion. We find a case of normal dispersion of waves when the plasma admits the properties of Veselago medium. Our results agree with those of Mackay et al. according to which rotation of a black hole is required for negative phase velocity propagation.

We apply the ADM 3+1 formalism to derive the general relativistic magnetohydrodynamic equations for cold plasma in spatially flat Schwarzschild metric. Respective perturbed equations are linearized for non-magnetized and magnetized plasmas both in non-rotating and rotating backgrounds. These are then Fourier analyzed and the corresponding dispersion relations are obtained. These relations are discussed for the existence of waves with positive angular frequency in the region near the horizon. Our results support the fact that no information can be extracted from the Schwarzschild black hole. It is concluded that negative phase velocity propagates in the rotating background whether the black hole is rotating or non-rotating.