We present a linear analysis of the radiative Rayleigh-Taylor (RT) instability in the presence of magnetic field for both optically thin and thick regimes. When the flow is optically thin, magnetic field not only stabilizes perturbations with short wavelengths, but also growth rate of the instability at long wavelengths is reduced compared to a nonmagnetized case. Then, we extend our analysis to the optically thick flows with a conserved total specific entropy and properties of the unstable perturbations are investigated in detail. Growth rate of the instability at short wavelengths is suppressed due to the presence of the magnetic field, however, growth rate is nearly constant at long wavelengths because of the radiation field. Since the radiative bubbles around massive protostars are subject to the RT instability, we also explore implications of our results in this context. In the nonmagnetized case, the growth time-scale of the instability for a typical bubble is found less than one thousand years which is very short compared to the typical star formation time-scale. Magnetic field with a reasonable strength significantly increases the growth time-scale to more than hundreds of thousands years. The instability, furthermore, is more efficient at large wavelengths, whereas in the non-magnetized case, growth rate at short wavelengths is more significant.

We study linear theory of the magnetized Rayleigh-Taylor instability in a system consisting of ions and neutrals. Both components are affected by a uniform vertical gravitational field. We consider ions and neutrals as two separate fluid systems where they can exchange momentum through collisions. However, ions have direct interaction with the magnetic field lines but neutrals are not affected by the field directly. The equations of our two-fluid model are linearized and by applying a set of proper boundary conditions, a general dispersion relation is derived for our two superposed fluids separated by a horizontal boundary. We found two unstable modes for a range of the wavenumbers. It seems that one of the unstable modes corresponds to the ions and the other one is for the neutrals. Both modes are reduced with increasing the collision rate of the particles and the ionization fraction. We show that if the two-fluid nature is considered, RT instability would not be suppressed and also show that the growth time of the perturbations increases. As an example, we apply our analysis to the Local clouds which seems to have arisen because of the RT instability. Assuming that the clouds are partially ionized, we find that the growth rate of these clouds increases in comparison to a fully ionized case.