We investigate the transient and stationary buoyant motion of the Rayleigh-Bénard instability when the fluid layer is subjected to a vertical, steady magnetic field. For Rayleigh number, Ra, in the range 10-10, and Hartmann number, Ha, between 0 and 100, we performed three-dimensional direct numerical simulations. To predict the growth rate and the wavelength of the initial regime observed with the numerical simulations, we developed the linear stability analysis beyond marginal stability for this problem. We analyzed the pattern of the flow from linear to nonlinear regime. We observe the evolution of steady state patterns depending on [Formula: see text] and Ha. In addition, in the nonlinear regime, the averaged kinetic energy is found to depend on Ra and to be independent of Ha in the studied range.