We investigate the effect of magnetic helicity on the stability of buoyant magnetic cavities as found in the intergalactic medium. In these cavities, we insert helical magnetic fields and test whether or not helicity can increase their stability to shredding through the Kelvin-Helmholtz instability and, with that, their lifetime. This is compared to the case of an external vertical magnetic field which is known to reduce the growth rate of the Kelvin-Helmholtz instability. By comparing a low-helicity configuration with a high helicity one with the same magnetic energy we find that an internal helical magnetic field stabilizes the cavity. This effect increases as we increase the helicity content. Stabilizing the cavity with an external magnetic field requires instead a significantly stronger field at higher magnetic energy. We conclude that the presence of helical magnetic fields is a viable mechanism to explain the stability of intergalactic cavities on time scales longer than 100 Myr.

The full-sky Planck polarization data at 850um revealed unexpected properties of the E and B mode power spectra of dust emission in the interstellar medium (ISM). The positive cross-correlation between the total dust intensity, T, with the B modes has raised new questions about the physical mechanisms that affect dust polarization, such as the Galactic magnetic-field structure. This is key both to better understanding ISM dynamics and to accurately describing Galactic foregrounds to the polarization of the Cosmic Microwave Background (CMB). In this theoretical paper we investigate the possibility that the observed cross-correlations in the dust polarization power spectra, and specifically between T and B, can be related to a parity-odd quantity in the ISM such as the magnetic helicity. We produce synthetic dust polarization data, derived from 3D analytical toy models of density structures and helical magnetic fields, to compare with the E and B modes of observations. Focusing on the observed T-B correlation, we propose a new line of interpretation of the Planck observations based on a large-scale helical component of the Galactic magnetic field in the solar neighborhood. Our analysis shows that: I) the sign of magnetic helicity does not affect E and B modes for isotropic magnetic-field configurations; II) helical magnetic fields threading interstellar filaments cannot reproduce the Planck results; III) a weak helical left-handed magnetic field structure in the solar neighborhood may explain the T-B correlation seen in the Planck data. This work suggests a new perspective for the interpretation of the dust polarization power spectra, which strongly supports the imprint of a large-scale structure of the Galactic magnetic field in the solar neighborhood.

Building on previous results on the quadratic helicity in magnetohydrodynamics (MHD) we investigate particular minimum helicity states. Those are eigenfunctions of the curl operator and are shown to constitute solutions of the quasi-stationary incompressible ideal MHD equations. We then show that these states have indeed minimum quadratic helicity.

In this paper we seek to understand the timescale on which the photospheric motions on the Sun braid coronal magnetic field lines. This is a crucial ingredient for determining the viability of the braiding mechanism for explaining the high temperatures observed in the corona. We study the topological complexity induced in the coronal magnetic field, primarily using plasma motions extracted from magneto-convection simulations. This topological complexity is quantified using the field line winding, finite time topological entropy and passive scalar mixing. With these measures we contrast mixing efficiencies of the magneto-convection simulation, a benchmark flow known as a ``blinking vortex', and finally photospheric flows inferred from sequences of observed magnetograms using local correlation tracking. While the highly resolved magneto-convection simulations induce a strong degree of field line winding and finite time topological entropy, the values obtained from the observations from the plage region are around an order of magnitude smaller. This behavior is carried over to the finite time topological entropy. Nevertheless, the results suggest that the photospheric motions induce complex tangling of the coronal field on a timescale of hours.

We investigate the existence of magnetohydrostatic equilibria for topologically complex magnetic fields. The approach employed is to perform ideal numerical relaxation experiments. We use a newly-developed Lagrangian relaxation scheme that exactly preserves the magnetic field topology during the relaxation. Our configurations include both twisted and sheared fields, of which some fall into the category for which Parker (1972) predicted no force-free equilibrium. The first class of field considered contains no magnetic null points, and field lines connect between two perfectly conducting plates. In these cases we observe only resolved current layers of finite thickness. In further numerical experiments we confirm that magnetic null points are loci of singular currents.

We study the effect of photospheric footpoint motions on magnetic field structures containing magnetic nulls. The footpoint motions are prescribed on the photospheric boundary as a velocity field which entangles the magnetic field. We investigate the propagation of the injected energy, the conversion of energy, emergence of current layers and other consequences of the non-trivial magnetic field topology in this situation. These boundary motions lead initially to an increase in magnetic and kinetic energy. Following this, the energy input from the photosphere is partially dissipated and partially transported out of the domain through the Poynting flux. The presence of separatrix layers and magnetic null-points fundamentally alters the propagation behavior of disturbances from the photosphere into the corona. Depending on the field line topology close to the photosphere, the energy is either trapped or free to propagate into the corona.

We examine the dynamics of magnetic flux tubes containing non-trivial field line braiding (or linkage), using mathematical and computational modelling, in the context of testable predictions for the laboratory and their significance for solar coronal heating. We investigate the existence of braided force-free equilibria, and demonstrate that for a field anchored at perfectly-conducting plates, these equilibria exist and contain current sheets whose thickness scales inversely with the braid complexity - as measured for example by the topological entropy. By contrast, for a periodic domain braided exact equilibria typically do not exist, while approximate equilibria contain thin current sheets. In the presence of resistivity, reconnection is triggered at the current sheets and a turbulent relaxation ensues. We finish by discussing the properties of the turbulent relaxation and the existence of constraints that may mean that the final state is not the linear force-free field predicted by Taylor's hypothesis.

Kepler data from G-, K- and M-type stars are used to study conditions that lead to superflares with energies above $10^{34} {\rm erg}$. From the 117,661 stars included, 380 show superflares with a total of 1690 such events. We study whether parameters, like effective temperature or the rotation rate, have any effect on the superflare occurrence rate or energy. With increasing effective temperature we observe a decrease in the superflare rate, which is analogous to the previous findings of a decrease in dynamo activity with increasing effective temperature. For slowly rotating stars, we find a quadratic increase of the mean occurrence rate with the rotation rate up to a critical point, after which the rate decreases linearly. Motivated by standard dynamo theory, we study the behavior of the relative starspot coverage, approximated as the relative brightness variation. For faster rotating stars, an increased fraction of stars shows higher spot coverage, which leads to higher superflare rates. A turbulent dynamo is used to study the dependence of the Ohmic dissipation as a proxy of the flare energy on the differential rotation or shear rate. The resulting statistics of the dissipation energy as a function of dynamo number is similar to the observed flare statistics as a function of the inverse Rossby number and shows similarly strong fluctuations. This supports the idea that superflares might well be possible for solar-type G stars.

We present a new code that performs a relaxation of a magnetic field towards a force-free state (Beltrami field) using a Lagrangian numerical scheme. Beltrami fields are of interest for the dynamics of many technical and astrophysical plasmas as they are the lowest energy states that the magnetic field can reach. The numerical method strictly preserves the magnetic flux and the topology of magnetic field lines. In contrast to other implementations we use mimetic operators for the spatial derivatives in order to improve accuracy for high distortions of the grid. Compared with schemes using direct derivatives we find that the final state of the simulation approximates a force-free state with a significantly higher accuracy. We implement the scheme in a code which runs on graphical processing units (GPU), which leads to an enhanced computing speed compared to previous relaxation codes.

Magnetic field generation on scales large compared with the scale of the turbulent eddies is known to be possible via the so-called $\alpha$ effect when the turbulence is helical and if the domain is large enough for the $\alpha$ effect to dominate over turbulent diffusion. Using three-dimensional turbulence simulations, we show that the energy of the resulting mean magnetic field of the saturated state increases linearly with the product of normalized helicity and the ratio of domain scale to eddy scale, provided this product exceeds a critical value of around unity. This implies that large-scale dynamo action commences when the normalized helicity is larger than the inverse scale ratio. Our results show that the emergence of small-scale dynamo action does not have any noticeable effect on the large-scale dynamo. Recent findings by Pietarila Graham et al. (2012, Phys. Rev. E85, 066406) of a smaller minimal helicity may be an artifact due to the onset of small-scale dynamo action at large magnetic Reynolds numbers. However, the onset of large-scale dynamo action is difficult to establish when the kinetic helicity is small. Instead of random forcing, they used an ABC-flow with time-dependent phases. We show that such dynamos saturate prematurely in a way that is reminiscent of inhomogeneous dynamos with internal magnetic helicity fluxes. Furthermore, even for very low fractional helicities, such dynamos display large-scale fields that change direction, which is uncharacteristic of turbulent dynamos.

Stability and reconnection of magnetic fields play a fundamental role in natural and man-made plasma. In these applications the field's topology determines the stability of the magnetic field. Here I will describe the importance of one topology quantifier, the magnetic helicity, which impedes any free decay of the magnetic energy. Further constraints come from the fixed point index which hinders the field to relax into the Taylor state.

Magnetic field relaxation is determined by both the field's geometry and its topology. For relaxation processes, however, it turns out that its topology is a much more stringent constraint. As quantifier for the topology we use magnetic helicity and test whether it is a stronger condition than the linking of field lines. Further, we search for evidence of other topological invariants, which give rise to further restrictions in the field's relaxation. We find that magnetic helicity is the sole determinant in most cases. Nevertheless, we see evidence for restrictions not captured through magnetic helicity.

Magnetic helicity fluxes in turbulently driven alpha^2 dynamos are studied to demonstrate their ability to alleviate catastrophic quenching. A one-dimensional mean-field formalism is used to achieve magnetic Reynolds numbers of the order of 10^5. We study both diffusive magnetic helicity fluxes through the mid-plane as well as those resulting from the recently proposed alternate dynamic quenching formalism. By adding shear we make a parameter scan for the critical values of the shear and forcing parameters for which dynamo action occurs. For this $\alpha\Omega$ dynamo we find that the preferred mode is antisymmetric about the mid-plane. This is also verified in 3-D direct numerical simulations.

In turbulent dynamos the production of large-scale magnetic fields is accompanied by a separation of magnetic helicity in scale. The large- and small-scale parts increase in magnitude. The small-scale part can eventually work against the dynamo and quench it, especially at high magnetic Reynolds numbers. A one-dimensional mean-field model of a dynamo is presented where diffusive magnetic helicity fluxes within the domain are important. It turns out that this effect helps to alleviate the quenching. Here we show that internal magnetic helicity fluxes, even within one hemisphere, can be important for alleviating catastrophic quenching.

Observations have shown that the Sun's magnetic field has helical structures. The helicity content in magnetic field configurations is a crucial constraint on the dynamical evolution of the system. Since helicity is connected with the number of links we investigate configurations with interlocked magnetic flux rings and one with unlinked rings. It turns out that it is not the linking of the tubes which affects the magnetic field decay, but the content of magnetic helicity.

Two setups with interlocked magnetic flux tubes are used to study the evolution of magnetic energy and helicity on magnetohydrodynamical (MHD) systems like plasmas. In one setup the initial helicity is zero while in the other it is finite. To see if it is the actual linking or merely the helicity content that influences the dynamics of the system we also consider a setup with unlinked field lines as well as a field configuration in the shape of a trefoil knot. For helical systems the decay of magnetic energy is slowed down by the helicity which decays slowly. It turns out that it is the helicity content, rather than the actual linking, that is significant for the dynamics.

We present calculations of the relaxation of magnetic field structures that have the shape of particular knots and links. A set of helical magnetic flux configurations is considered, which we call $n$-foil knots of which the trefoil knot is the most primitive member. We also consider two nonhelical knots; namely, the Borromean rings as well as a single interlocked flux rope that also serves as the logo of the Inter-University Centre for Astronomy and Astrophysics in Pune, India. The field decay characteristics of both configurations is investigated and compared with previous calculations of helical and nonhelical triple-ring configurations. Unlike earlier nonhelical configurations, the present ones cannot trivially be reduced via flux annihilation to a single ring. For the $n$-foil knots the decay is described by power laws that range form $t^{-2/3}$ to $t^{-1/3}$, which can be as slow as the $t^{-1/3}$ behavior for helical triple-ring structures that were seen in earlier work. The two nonhelical configurations decay like $t^{-1}$, which is somewhat slower than the previously obtained $t^{-3/2}$ behavior in the decay of interlocked rings with zero magnetic helicity. We attribute the difference to the creation of local structures that contain magnetic helicity which inhibits the field decay due to the existence of a lower bound imposed by the realizability condition. We show that net magnetic helicity can be produced resistively as a result of a slight imbalance between mutually canceling helical pieces as they are being driven apart. We speculate that higher order topological invariants beyond magnetic helicity may also be responsible for slowing down the decay of the two more complicated nonhelical structures mentioned above.

Magnetic helicity fluxes are investigated in a family of gauges in which the contribution from ideal magnetohydrodynamics takes the form of a purely advective flux. Numerical simulations of magnetohydrodynamic turbulence in this advective gauge family exhibit instabilities triggered by the build-up of unphysical irrotational contributions to the magnetic vector potential. As a remedy, the vector potential is evolved in a numerically well behaved gauge, from which the advective vector potential is obtained by a gauge transformation. In the kinematic regime, the magnetic helicity density evolves similarly to a passive scalar when resistivity is small and turbulent mixing is mild, i.e. when the fluid Reynolds number is not too large. In the dynamical regime, resistive contributions to the magnetic helicity flux in the advective gauge are found to be significant owing to the development of small length scales in the irrotational part of the magnetic vector potential.

The resistive decay of chains of three interlocked magnetic flux rings is considered. Depending on the relative orientation of the magnetic field in the three rings, the late-time decay can be either fast or slow. Thus, the qualitative degree of tangledness is less important than the actual value of the linking number or, equivalently, the net magnetic helicity. Our results do not suggest that invariants of higher order than that of the magnetic helicity need to be considered to characterize the decay of the field.

We use direct numerical simulations of forced MHD turbulence with a forcing function that produces two different signs of kinetic helicity in the upper and lower parts of the domain. We show that the mean flux of magnetic helicity from the small-scale field between the two parts of the domain can be described by a Fickian diffusion law with a diffusion coefficient that is approximately independent of the magnetic Reynolds number and about one third of the estimated turbulent magnetic diffusivity. The data suggest that the turbulent diffusive magnetic helicity flux can only be expected to alleviate catastrophic quenching at Reynolds numbers of more than several thousands. We further calculate the magnetic helicity density and its flux in the domain for three different gauges. We consider the Weyl gauge, in which the electrostatic potential vanishes, the pseudo-Lorenz gauge, where the speed of light is replaced by the sound speed, and the `resistive gauge' in which the Laplacian of the magnetic vector potential acts as resistive term. We find that, in the statistically steady state, the time-averaged magnetic helicity density and the magnetic helicity flux are the same in all three gauges.

Using mean-field models with a dynamical quenching formalism we show that in finite domains magnetic helicity fluxes associated with small-scale magnetic fields are able to alleviate catastrophic quenching. We consider fluxes that result either from advection by a mean flow, the turbulent mixing down the gradient of mean small-scale magnetic helicity concentration, or the explicit removal which may be associated with the effects of coronal mass ejections in the Sun. In the absence of shear, all the small-scale magnetic helicity fluxes are found to be equally strong both for large-scale and small-scale fields. In the presence of shear there is also an additional magnetic helicity flux associated with the mean field, but this flux does not alleviate catastrophic quenching. Outside the dynamo-active region there are neither sources nor sinks of magnetic helicity, so in a steady state this flux must be constant. It is shown that unphysical behavior emerges if the small-scale magnetic helicity flux is forced to vanish within the computational domain.