The negative effective magnetic pressure instability discovered recently in direct numerical simulations (DNSs) may play a crucial role in the formation of sunspots and active regions in the Sun and stars. This instability is caused by a negative contribution of turbulence to the effective mean Lorentz force (the sum of turbulent and non-turbulent contributions) and results in the formation of large-scale inhomogeneous magnetic structures from an initially uniform magnetic field. Earlier investigations of this instability in DNSs of stably stratified, externally forced, isothermal hydromagnetic turbulence in the regime of large plasma β are now extended into the regime of larger scale separation ratios where the number of turbulent eddies in the computational domain is about 30. Strong spontaneous formation of large-scale magnetic structures is seen even without performing any spatial averaging. These structures encompass many turbulent eddies. The characteristic time of the instability is comparable to the turbulent diffusion time, L²/η_t, where η_t is the turbulent diffusivity and L is the scale of the domain. DNSs are used to confirm that the effective magnetic pressure does indeed become negative for magnetic field strengths below the equipartition field. The dependence of the effective magnetic pressure on the field strength is characterized by fit parameters that seem to show convergence for larger values of the magnetic Reynolds number