Negative Effective Magnetic Pressure Instability Research Articles

Magnetic bipoles in rotating turbulence with coronal envelope

Context. The formation mechanism of sunspots and starspots is not yet fully understood. It is a major open problem in astrophysics.Aims. Magnetic flux concentrations can be produced by the negative effective magnetic pressure instability (NEMPI). This instability is strongly suppressed by rotation. However, the presence of an outer coronal envelope was previously found to strengthen the flux concentrations and make them more prominent. It also allows for the formation of bipolar regions (BRs). We aim to understand the important issue of whether the presence of an outer coronal envelope also changes the excitation conditions and the rotational dependence of NEMPI.Methods. We have used direct numerical simulations and mean-field simulations. We adopted a simple two-layer model of turbulence that mimics the jump between the convective turbulent and coronal layers below and above the surface of a star, respectively. The computational domain is Cartesian and located at a certain latitude of a rotating sphere. We investigated the effects of rotation on NEMPI by changing the Coriolis number, the latitude, the strengths of the imposed magnetic field, and the box resolution.Results. Rotation has a strong impact on the process of BR formation. Even rather slow rotation is found to suppress BR formation. However, increasing the imposed magnetic field strength also makes the structures stronger and alleviates the rotational suppression somewhat. The presence of a coronal layer itself does not significantly reduce the effects of rotational suppression.

Nonlinear mean-field dynamo and prediction of solar activity

We apply a nonlinear mean-field dynamo model which includes a budget equation for the dynamics of Wolf numbers to predict solar activity. This dynamo model takes into account the algebraic and dynamic nonlinearities of the$\unicode[STIX]{x1D6FC}$effect, where the equation for the dynamic nonlinearity is derived from the conservation law for the magnetic helicity. The budget equation for the evolution of the Wolf number is based on a formation mechanism of sunspots related to the negative effective magnetic pressure instability. This instability redistributes the magnetic flux produced by the mean-field dynamo. To predict solar activity on the time scale of one month we use a method based on a combination of the numerical solution of the nonlinear mean-field dynamo equations and the artificial neural network. A comparison of the results of the prediction of the solar activity with the observed Wolf numbers demonstrates a good agreement between the forecast and observations.

Spontaneous flux concentrations from the negative effective magnetic pressure instability beneath a radiative stellar surface

Context. The formation of sunspots requires the concentration of magnetic flux near the surface. The negative effective magnetic pressure instability (NEMPI) might be a possible mechanism for accomplishing this, but it has mainly been studied in simple systems using an isothermal equation of state without a natural free surface. Aims. We study NEMPI in a stratified Cartesian mean-field model where turbulence effects are parameterized. We use an ideal equation of state and include radiation transport, which establishes selfconsistently a free surface. Methods. We use a Kramers-type opacity with adjustable exponents chosen such that the deeper layers are approximately isentropic. No convection is therefore possible in this model, allowing us to study NEMPI with radiation in isolation. We restrict ourselves to two-dimensional models. We use artificially enhanced mean-field coefficients to allow NEMPI to develop, thereby making it possible to study the reason why it is much harder to excite in the presence of radiation. Results. NEMPI yields moderately strong magnetic flux concentrations a certain distance beneath the surface where the optical depth is unity. The instability is oscillatory and in the form of upward traveling waves. This seems to be a new effect that has not been found in earlier models without radiative transport. The horizontal wavelength is about ten times smaller than what has previously been found in more idealized isothermal models. Conclusions. In our models, NEMPI saturates at field strengths too low to explain sunspots. Furthermore, the structures appear too narrow and too far beneath the surface to cause significant brightness variations at the radiative surface. We speculate that the failure to reproduce effects resembling sunspots may be related to the neglect of convection.

Magnetic concentrations in stratified turbulence: the negative effective magnetic pressure instability

In the presence of strong density stratification, hydromagnetic turbulence attains qualitatively new properties: the formation of magnetic flux concentrations. We review here the theoretical foundations of this mechanism in terms of what is now called the negative effective magnetic pressure instability. We also present direct numerical simulations of forced turbulence in strongly stratified layers and discuss the qualitative and quantitative similarities with corresponding mean-field simulations. Finally, the relevance to sunspot formation is discussed.

A new look at sunspot formation using theory and observations

AbstractSunspots are of basic interest in the study of the Sun. Their relevance ranges from them being an activity indicator of magnetic fields to being the place where coronal mass ejections and flares erupt. They are therefore also an important ingredient of space weather. Their formation, however, is still an unresolved problem in solar physics. Observations utilize just 2D surface information near the spot, but it is debatable how to infer deep structures and properties from local helioseismology. For a long time, it was believed that flux tubes rising from the bottom of the convection zone are the origin of the bipolar sunspot structure seen on the solar surface. However, this theory has been challenged, in particular recently by new surface observation, helioseismic inversions, and numerical models of convective dynamos. In this article we discuss another theoretical approach to the formation of sunspots: the negative effective magnetic pressure instability. This is a large-scale instability, in which the total (kinetic plus magnetic) turbulent pressure can be suppressed in the presence of a weak large-scale magnetic field, leading to a converging downflow, which eventually concentrates the magnetic field within it. Numerical simulations of forced stratified turbulence have been able to produce strong super-equipartition flux concentrations, similar to sunspots at the solar surface. In this framework, sunspots would only form close to the surface due to the instability constraints on stratification and rotation. Additionally, we present some ideas from local helioseismology, where we plan to use the Hankel analysis to study the pre-emergence phase of a sunspot and to constrain its deep structure and formation mechanism.

The dynamics of Wolf numbers based on nonlinear dynamos with magnetic helicity: comparisons with observations

We investigate the dynamics of solar activity using a nonlinear one-dimensional dynamo model and a phenomenological equation for the evolution of Wolf numbers. This system of equations is solved numerically. We take into account the algebraic and dynamic nonlinearities of the alpha effect. The dynamic nonlinearity is related to the evolution of a small-scale magnetic helicity, and it leads to a complicated behavior of solar activity. The evolution equation for the Wolf number is based on a mechanism of formation of magnetic spots as a result of the negative effective magnetic pressure instability (NEMPI). This phenomenon was predicted 25 years ago and has been investigated intensively in recent years through direct numerical simulations and mean-field simulations. The evolution equation for the Wolf number includes the production and decay of sunspots. Comparison between the results of numerical simulations and observational data of Wolf numbers shows a 70 % correlation over all intervals of observation (about 270 years). We determine the dependence of the maximum value of the Wolf number versus the period of the cycle and the asymmetry of the solar cycles versus the amplitude of the cycle. These dependencies are in good agreement with observations.

Bipolar region formation in stratified two-layer turbulence

This work presents an extensive study of the previously discovered formation of bipolar flux concentrations in a two-layer model. We interpret the formation process in terms of negative effective magnetic pressure instability (NEMPI), which is a possible mechanism to explain the origin of sunspots. In our simulations, we use a Cartesian domain of isothermal stratified gas that is divided into two layers. In the lower layer, turbulence is forced with transverse nonhelical random waves, whereas in the upper layer no flow is induced. A weak uniform magnetic field is imposed in the entire domain at all times. In this study we vary the stratification by changing the gravitational acceleration, magnetic Reynolds number, strength of the imposed magnetic field, and size of the domain to investigate their influence on the formation process. Bipolar magnetic structure formation takes place over a large range of parameters. The magnetic structures become more intense for higher stratification until the density contrast becomes around $100$ across the turbulent layer. For the Reynolds numbers considered, magnetic flux concentrations are generated at magnetic Prandtl number between 0.1 and 1. The magnetic field in bipolar regions increases with higher imposed field strength until the field becomes comparable to the equipartition field strength of the turbulence. A larger horizontal extent enables the flux concentrations to become stronger and more coherent. The size of the bipolar structures turns out to be independent of the domain size. In the case of bipolar region formation, we find an exponential growth of the large-scale magnetic field, which is indicative of a hydromagnetic instability. Additionally, the flux concentrations are correlated with strong large-scale downward and converging flows. These findings imply that NEMPI is responsible for magnetic flux concentrations.

Magnetic flux concentrations from turbulent stratified convection

(abridged) Context: The mechanisms that cause the formation of sunspots are still unclear. Aims: We study the self-organisation of initially uniform sub-equipartition magnetic fields by highly stratified turbulent convection. Methods: We perform simulations of magnetoconvection in Cartesian domains that are $8.5$-$24$ Mm deep and $34$-$96$ Mm wide. We impose either a vertical or a horizontal uniform magnetic field in a convection-driven turbulent flow. Results: We find that super-equipartition magnetic flux concentrations are formed near the surface with domain depths of $12.5$ and $24$ Mm. The size of the concentrations increases as the box size increases and the largest structures ($20$ Mm horizontally) are obtained in the 24 Mm deep models. The field strength in the concentrations is in the range of $3$-$5$ kG. The concentrations grow approximately linearly in time. The effective magnetic pressure measured in the simulations is positive near the surface and negative in the bulk of the convection zone. Its derivative with respect to the mean magnetic field, however, is positive in the majority of the domain, which is unfavourable for the negative effective magnetic pressure instability (NEMPI). Furthermore, we find that magnetic flux is concentrated in regions of converging flow corresponding to large-scale supergranulation convection pattern. Conclusions: The linear growth of large-scale flux concentrations implies that their dominant formation process is tangling of the large-scale field rather than an instability. One plausible mechanism explaining both the linear growth and the concentrate on of the flux in the regions of converging flow pattern is flux expulsion. Possible reasons for the absence of NEMPI are that the derivative of the effective magnetic pressure with respect to the mean magnetic field has an unfavourable sign and that there may not be sufficient scale separation.

Hydraulic effects in a radiative atmosphere with ionization

In a paper of 1978, Eugene Parker postulated the need for hydraulic downward motion to explain magnetic flux concentrations at the solar surface. A similar process has recently also been seen in simplified (e.g., isothermal) models of flux concentrations from the negative effective magnetic pressure instability. We study the effects of partial ionization near the radiative surface on the formation of such magnetic flux concentrations. We first obtain one-dimensional (1D) equilibrium solutions using either a Kramers-like opacity or the ${\rm H}^{-}$ opacity. The resulting atmospheres are then used as initial conditions in two-dimensional (2D) models where flows are driven by an imposed gradient force resembling a localized negative pressure in the form of a blob. To isolate the effects of partial ionization and radiation, we ignore turbulence and convection. In 1D models, due to partial ionization, an unstable stratification forms always near the surface. We show that the extrema in the specific entropy profiles correspond to the extrema in degree of ionization. In the 2D models without partial ionization, flux concentrations form close to the height where the blob is placed. In models with partial ionization, such flux concentrations form at the surface much above the blob. This is due to the corresponding unstable layer in specific entropy. With ${\rm H}^{-}$ opacity, flux concentrations are weaker due to the stably stratified deeper parts. We demonstrate that, together with density stratification, the imposed source of negative pressure drives the formation of flux concentrations. We find that the inclusion of partial ionization affects entropy profiles causing the strong flux concentrations to form closer to the surface. We speculate that turbulence is needed to limit the strength of flux concentrations and homogenize the specific entropy to a more nearly marginal stratification.

Magnetic flux concentrations from dynamo-generated fields

The mean-field theory of magnetized stellar convection gives rise to the two possibility of distinct instabilities: the large-scale dynamo instability, operating in the bulk of the convection zone, and a negative effective magnetic pressure instability (NEMPI) operating in the strongly stratified surface layers. The latter might be important in connection with magnetic spot formation, but the growth rate of NEMPI is suppressed with increasing rotation rates, although recent direct numerical simulations (DNS) have shown a subsequent increase in the growth rate. We examine quantitatively whether this increase in the growth rate of NEMPI can be explained by an alpha squared mean-field dynamo, and whether both NEMPI and the dynamo instability can operate at the same time. We use both DNS and mean-field simulations (MFS) to solve the underlying equations numerically either with or without an imposed horizontal field. We use the test-field method to compute relevant dynamo coefficients. DNS show that magnetic flux concentrations are still possible up to rotation rates above which the large-scale dynamo effect produces mean magnetic fields. The resulting DNS growth rates are quantitatively well reproduced with MFS. As expected, for weak or vanishing rotation, the growth rate of NEMPI increases with increasing gravity, but there is a correction term for strong gravity and large turbulent magnetic diffusivity. Magnetic flux concentrations are still possible for rotation rates above which dynamo action takes over. For the solar rotation rate, the corresponding turbulent turnover time is about 5 hours, with dynamo action commencing in the layers beneath.

Magnetic flux concentrations in a polytropic atmosphere

Strongly stratified hydromagnetic turbulence has recently been identified as a candidate for explaining the spontaneous formation of magnetic flux concentrations by the negative effective magnetic pressure instability (NEMPI). Much of this work has been done for isothermal layers for which the density scale height is constant throughout. We now study the validity of earlier conclusions about the size and growth rate of magnetic structures in the case of polytropic layers, which scale height decreases sharply towards the surface. To allow for a continuous transition from isothermal to polytropic layers, we employ a generalization of the exponential function known as the q-exponential. Now, the top of the polytropic layer shifts with the polytropic index such that the scale height at some reference height is always the same. We use both mean-field and direct numerical simulations of forced stratified turbulence to determine the resulting flux concentrations. Magnetic structures begin to form at a depth where magnetic field strength is about 3-4% the local equipartition field strength with respect to the turbulent kinetic energy. Unlike the isothermal case where stronger fields can give rise to magnetic flux concentrations at larger depths, in the polytropic one the growth rates decreases for structures deeper down. For vertical fields, magnetic structures of super-equipartition strengths are formed because such fields survive downward advection, unlike NEMPI under horizontal magnetic fields. The horizontal cross-section of such structures is approximately circular. Results based on isothermal models can be applied locally to polytropic layers. For vertical fields, magnetic flux concentrations of super-equipartition strengths form, which supports suggestions that sunspot formation might be a shallow phenomenon.

Mean-field and direct numerical simulations of magnetic flux concentrations from vertical field

Strongly stratified hydromagnetic turbulence has previously been found to produce magnetic flux concentrations if the domain is large enough compared with the size of turbulent eddies. Mean-field simulations (MFS) using parameterizations of the Reynolds and Maxwell stresses show a negative effective magnetic pressure instability and have been able to reproduce many aspects of direct numerical simulations (DNS) regarding the growth rate of this large-scale instability, shape of the resulting magnetic structures, and their height as a function of magnetic field strength. Unlike the case of an imposed horizontal field, for a vertical one, magnetic flux concentrations of equipartition strength with the turbulence can be reached. This results in magnetic spots that are reminiscent of sunspots. Here we want to find out under what conditions magnetic flux concentrations with vertical field occur and what their internal structure is. We use a combination of MFS, DNS, and implicit large-eddy simulations to characterize the resulting magnetic flux concentrations in forced isothermal turbulence with an imposed vertical magnetic field. We confirm earlier results that in the kinematic stage of the large-scale instability the horizontal wavelength of structures is about 10 times the density scale height. At later times, even larger structures are being produced in a fashion similar to inverse spectral transfer in helically driven turbulence. Using turbulence simulations, we find that magnetic flux concentrations occur for different values of the Mach number between 0.1 and 0.7. DNS and MFS show magnetic flux tubes with mean-field energies comparable to the turbulent kinetic energy. The resulting vertical magnetic flux tubes are being confined by downflows along the tubes and corresponding inflow from the sides, which keep the field concentrated.

BIPOLAR MAGNETIC STRUCTURES DRIVEN BY STRATIFIED TURBULENCE WITH A CORONAL ENVELOPE

We report the spontaneous formation of bipolar magnetic structures in direct numerical simulations of stratified forced turbulence with an outer coronal envelope. The turbulence is forced with transverse random waves only in the lower (turbulent) part of the domain. Our initial magnetic field is either uniform in the entire domain or confined to the turbulent layer. After about 1--2 turbulent diffusion times, a bipolar magnetic region of vertical field develops with two coherent circular structures that live during one turbulent diffusion time, and then decay during 0.5 turbulent diffusion times. The resulting magnetic field strengths inside the bipolar region are comparable to the equipartition value with respect to the turbulent kinetic energy. The bipolar magnetic region forms a loop-like structure in the upper coronal layer. We associate the magnetic structure formation with the negative effective magnetic pressure instability in the two-layer model.

Surface flux concentrations in a sphericalα2dynamo

In the presence of strong density stratification, turbulence can lead to a large-scale instability of a horizontal magnetic field if its strength is in a suitable range (within a few percent of the turbulent equipartition value). This instability is related to a suppression of the turbulent pressure so that the turbulence contribution to the mean magnetic pressure becomes negative. This results in the excitation of a negative effective magnetic pressure instability (NEMPI). This instability has so far only been studied for an imposed magnetic field. We want to know how NEMPI works when the mean magnetic field is generated self-consistently by an $\alpha^2$ dynamo, whether it is affected by global spherical geometry, and whether it can influence the properties of the dynamo itself. We adopt the mean-field approach which has previously been shown to provide a realistic description of NEMPI in direct numerical simulations. We assume axisymmetry and solve the mean-field equations with the Pencil-Code for an adiabatic stratification at a total density contrast in the radial direction of approximately 4 orders of magnitude. NEMPI is found to work when the dynamo-generated field is about 4% of the equipartition value, which is achieved through strong $\alpha$ quenching. This instability is excited in the top 5% of the outer radius provided the density contrast across this top layer is at least 10. NEMPI is found to occur at lower latitudes when the mean magnetic field is stronger. For weaker fields, NEMPI can make the dynamo oscillatory with poleward migration. NEMPI is a viable mechanism for producing magnetic flux concentrations in a strongly stratified spherical shell in which a magnetic field is generated by a strongly quenched $\alpha$ effect dynamo.

Competition of rotation and stratification in flux concentrations

In a strongly stratified turbulent layer, a uniform horizontal magnetic field can become unstable to spontaneously form local flux concentrations due to a negative contribution of turbulence to the large-scale (mean-field) magnetic pressure. This mechanism, called the negative effective magnetic pressure instability (NEMPI), is of interest in connection with dynamo scenarios where most of the magnetic field resides in the bulk of the convection zone, and not at the bottom. Recent work using the mean-field hydromagnetic equations has shown that NEMPI becomes suppressed at rather low rotation rates with Coriolis numbers as low as 0.1.}{Here we extend these earlier investigations by studying the effects of rotation both on the development of NEMPI and on the effective magnetic pressure. We also quantify the kinetic helicity from direct numerical simulations (DNS) and compare with earlier work.}{To calculate the rotational effect on the effective magnetic pressure we consider both DNS and analytical studies using the $\tau$ approach. To study the effects of rotation on the development of NEMPI we use both DNS and mean-field calculations of the 3D hydromagnetic equations in a Cartesian domain.}{We find that the growth rates of NEMPI from earlier mean-field calculations are well reproduced with DNS, provided the Coriolis number is below about 0.06. In that case, kinetic and magnetic helicities are found to be weak. For faster rotation, dynamo action becomes possible. However, there is an intermediate range of rotation rates where dynamo action on its own is not yet possible, but the rotational suppression of NEMPI is being alleviated.}{Production of magnetic flux concentrations through the suppression of turbulent pressure appears to be possible only in the upper-most layers of the Sun, where the convective turnover time is less than 2 hours.}

Non-uniformity effects in the negative effective magnetic pressure instability

In direct numerical simulations of strongly stratified turbulence we have previously studied the development of large scale magnetic structures starting from a uniform background field. This is caused by an instability resulting from a negative contribution of small-scale turbulence to the effective (mean-field) magnetic pressure, and was qualitatively reproduced in mean-field simulations (MFS) where this pressure reduction was modeled as a function of the mean magnetic field normalized by the equipartition field. We now investigate the effect of mean current density on the turbulent pressure reduction. In our MFS, such currents are associated with sharp gradients of the growing structures. We find that an enhanced mean current density increases the suppression of the turbulent pressure.

Rotational effects on the negative magnetic pressure instability

The surface layers of the Sun are strongly stratified. In the presence of turbulence with a weak mean magnetic field, a large-scale instability resulting in the formation of non-uniform magnetic structures, can be excited over the scale of many turbulent eddies or convection cells. This instability is caused by a negative contribution of turbulence to the effective (mean-field) magnetic pressure and has previously been discussed in connection with the formation of active regions and perhaps sunspots. We want to understand the effects of rotation on this instability in both two and three dimensions. We use mean-field magnetohydrodynamics in a parameter regime in which the properties of the negative effective magnetic pressure instability have previously been found to be in agreement with those of direct numerical simulations. We find that the instability is suppressed already for relatively slow rotation with Coriolis numbers (i.e. inverse Rossby numbers) around 0.2. The suppression is strongest at the equator. In the nonlinear regime, we find traveling wave solutions with propagation in the prograde direction at the equator with additional poleward migration away from the equator. The prograde rotation of the magnetic pattern near the equator is argued to be a possible explanation for the faster rotation speed of magnetic tracers found on the Sun. In the bulk of the domain, kinetic and current helicities are negative in the northern hemisphere and positive in the southern.

Flux concentrations in turbulent convection

AbstractWe present preliminary results from high resolution magneto-convection simulations where we find the formation of flux concentrations from an initially uniform magnetic field. The structures appear in roughly ten convective turnover times and live close to a turbulent diffusion time. The time scales are compatible with the negative effective magnetic pressure instability (NEMPI), although structure formation is not restricted to regions where the effective magnetic pressure is negative.

Role of longitudinal activity complexes for solar and stellar dynamos

AbstractIn this paper we first discuss observational evidence of longitudinal concentrations of magnetic activity in the Sun and rapidly rotating late-type stars with outer convective envelopes. Scenarios arising from the idea of rotationally influenced anisotropic convective turbulence being the key physical process generating these structures are then presented and discussed - such effects include the turbulent dynamo mechanism, negative effective magnetic pressure instability (NEMPI) and hydrodynamical vortex instability. Finally, we discuss non-axisymmetric stellar mean-field dynamo models, the results obtained with them, and compare those with the observational information gathered up so far. We also present results from a pure α2 mean-field dynamo model, which show that time-dependent behavior of the dynamo solutions can occur both in the form of an azimuthal dynamo wave and/or oscillatory behavior related to the alternating energy levels of the active longitudes.

Active Region Formation through the Negative Effective Magnetic Pressure Instability

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.