We present the first demonstration of the negative effective magnetic pressure instability in direct numerical simulations of stably stratified, externally forced, isothermal hydromagnetic turbulence in the regime of large plasma beta. By the action of this instability, an initially uniform horizontal magnetic field forms flux concentrations whose scale is large compared to the turbulent scale. We further show that the magnetic energy of these large-scale structures is only weakly dependent on the magnetic Reynolds number, provided its value is large enough for the instability to be excited. Our results support earlier mean-field calculations and analytic work which identified this instability. Applications to the formation of active regions in the Sun are discussed.

The negative effective magnetic pressure instability operates on scales encompassing many turbulent eddies and is here discussed in connection with the formation of active regions near the surface layers of the Sun. This instability is related to the negative contribution of turbulence to the mean magnetic pressure that causes the formation of large-scale magnetic structures. For an isothermal layer, direct numerical simulations and mean-field simulations of this phenomenon are shown to agree in many details in that their onset occurs at the same depth. This depth increases with increasing field strength, such that the maximum growth rate of this instability is independent of the field strength, provided the magnetic structures are fully contained within the domain. A linear stability analysis is shown to support this finding. The instability also leads to a redistribution of turbulent intensity and gas pressure that could provide direct observational signatures.