Using a full kinetic, implicit particle-in-cell code, iPiC3D, we studied the properties of plasma kinetic turbulence, such as would be found at the interface between the solar wind and the Earth magnetosphere at low latitude during northwards periods. In this case, in the presence of a magnetic field B oriented mostly perpendicular to the velocity shear, turbulence is fed by the disruption of a Kelvin-Helmholtz vortex chain via secondary instabilities, vortex pairing and non-linear interactions. We found that the magnetic energy spectral cascade between ion and electron inertial scales, $d_i$ and $d_e$, is in agreement with satellite observations and other previous numerical simulations; however, in our case the spectrum ends with a peak beyond $d_e$ due to the occurrence of the lower hybrid drift instability. The electric energy spectrum is influenced by effects of secondary instabilities: anomalous resistivity, fed by the development of the lower hybrid drift instability, steepens the spectral decay and, depending on the alignment or anti-alignment of B and the shear vorticity, peaks due to ion-Bernstein waves may dominate the spectrum around $d_i$. These waves are generated by counter-streaming flow structures, through flux freezing also responsible for reconnection of the in-plane component of the magnetic field, which then generates electron pressure anisotropy and flattening of the field-aligned component of the electron distribution function.

Using direct numerical simulations (DNS) and mean-field simulations (MFS), the effects of non-uniformity of the magnetic field on the suppression of the turbulent pressure is investigated. This suppression of turbulent pressure can lead to an instability which, in turn, makes the mean magnetic field even more non-uniform. This large-scale instability is caused by a resulting negative contribution of small-scale turbulence to the effective (mean-field) magnetic pressure. We show that enhanced mean current density increases the suppression of the turbulent pressure. The instability leads to magnetic flux concentrations in which the normalized effective mean-field pressure is reduced to a certain value at all depths within a structure.

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.

We investigate the mean electromotive force in the kinematic framework, that is, ignoring the back-reaction of the magnetic field on the fluid velocity, under the assumption of axisymmetric turbulence determined by the presence of either rotation, density stratification, or both. We use an analogous approach for the mean passive scalar flux. As an alternative to convection, we consider forced turbulence in an isothermal layer. When using standard ansatzes, the mean magnetic transport is then determined by nine, and the mean passive scalar transport by four coefficients. We give results for all these transport coefficients. We use the test-field method and the test-scalar method, where transport coefficients are determined by solving sets of equations with properly chosen mean magnetic fields or mean scalars. These methods are adapted to mean fields which may depend on all three space coordinates. We find the anisotropy of turbulent diffusion to be moderate in spite of rapid rotation or strong density stratification. Contributions to the mean electromotive force determined by the symmetric part of the gradient tensor of the mean magnetic field, which were ignored in several earlier investigations, turn out to be important. In stratified rotating turbulence, the $\alpha$ effect is strongly anisotropic, suppressed along the rotation axis on large length scales, but strongly enhanced at intermediate length scales. Also the $\OO\times\meanJJ$ effect is enhanced at intermediate length scales. The turbulent passive scalar diffusivity is typically almost twice as large as the turbulent magnetic diffusivity. Both magnetic and passive scalar diffusion are slightly enhanced along the rotation axis, but decreased if there is gravity.

To understand the basic mechanism of the formation of magnetic flux concentrations, we determine by direct numerical simulations the turbulence contributions to the mean magnetic pressure in a strongly stratified isothermal layer with large plasma beta, where a weak uniform horizontal mean magnetic field is applied. The negative contribution of turbulence to the effective mean magnetic pressure is determined for strongly stratified forced turbulence over a range of values of magnetic Reynolds and Prandtl numbers. Small-scale dynamo action is shown to reduce the negative effect of turbulence on the effective mean magnetic pressure. However, the turbulence coefficients describing the negative effective magnetic pressure phenomenon are found to be converged for magnetic Reynolds numbers between 60 and 600, which is the largest value considered here. In all these models the turbulent intensity is arranged to be nearly independent of height, so the kinetic energy density decreases with height due to the decrease in density. In a second series of numerical experiments, the turbulent intensity increases with height such that the turbulent kinetic energy density is nearly independent of height. Turbulent magnetic diffusivity and turbulent pumping velocity are determined with the test-field method for both cases. The vertical profile of the turbulent magnetic diffusivity is found to agree with what is expected based on simple mixing length expressions. Turbulent pumping is shown to be down the gradient of turbulent magnetic diffusivity, but it is twice as large as expected. Corresponding numerical mean-field models are used to show that a large-scale instability can occur in both cases, provided the degree of scale separation is large enough and hence the turbulent magnetic diffusivity small enough.

As was demonstrated in earlier studies, turbulence can result in a negative contribution to the effective mean magnetic pressure, which, in turn, can cause a large-scale instability. In this study, hydromagnetic mean-field modelling is performed for an isothermally stratified layer in the presence of a horizontal magnetic field. The negative effective magnetic pressure instability (NEMPI) is comprehensively investigated. It is shown that, if the effect of turbulence on the mean magnetic tension force vanishes, which is consistent with results from direct numerical simulations of forced turbulence, the fastest growing eigenmodes of NEMPI are two-dimensional. The growth rate is found to depend on a parameter beta_* characterizing the turbulent contribution of the effective mean magnetic pressure for moderately strong mean magnetic fields. A fit formula is proposed that gives the growth rate as a function of turbulent kinematic viscosity, turbulent magnetic diffusivity, the density scale height, and the parameter beta_*. The strength of the imposed magnetic field does not explicitly enter provided the location of the vertical boundaries are chosen such that the maximum of the eigenmode of NEMPI fits into the domain. The formation of sunspots and solar active regions is discussed as possible applications of NEMPI.

The negative effective magnetic pressure instability discovered recently in direct numerical simulations (DNS) 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 formation of large-scale inhomogeneous magnetic structures from initial uniform magnetic field. Earlier investigations of this instability in DNS of stably stratified, externally forced, isothermal hydromagnetic turbulence in the regime of large plasma beta 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^2/eta_t, where eta_t is the turbulent diffusivity and L is the scale of the domain. DNS 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.

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.

A reduction of total mean turbulent pressure due to the presence of magnetic fields was previously shown to be a measurable effect in direct numerical simulations. However, in the studied parameter regime the formation of large-scale structures, as anticipated from earlier mean-field simulations, was not found. An analysis of the relevant mean-field parameter dependency and the parameter domain of interest is conducted in order to clarify this apparent discrepancy.

While the rising flux tube paradigm is an elegant theory, its basic assumptions, thin flux tubes at the bottom of the convection zone with field strengths two orders of magnitude above equipartition, remain numerically unverified at best. As such, in recent years the idea of a formation of sunspots near the top of the convection zone has generated some interest. The presence of turbulence can strongly enhance diffusive transport mechanisms, leading to an effective transport coefficient formalism in the mean-field formulation. The question is what happens to these coefficients when the turbulence becomes anisotropic due to a strong large-scale mean magnetic field. It has been noted in the past that this anisotropy can also lead to highly non-diffusive behaviour. In the present work we investigate the formation of large-scale magnetic structures as a result of a negative contribution of turbulence to the large-scale effective magnetic pressure in the presence of stratification. In direct numerical simulations of forced turbulence in a stratified box, we verify the existence of this effect. This phenomenon can cause formation of large-scale magnetic structures even from initially uniform large-scale magnetic field.