Mathematicians Link Magnetic Helicity to Free Energy in Solar Fields

Magnetic fields threading the solar atmosphere carry an invisible quantity that physicists prize almost as much as the fields themselves: helicity, a measure of how tangled, linked and twisted the field lines are. A new theoretical study, published as a Brief Report in the journal Discover Space, extends one of the most elegant results in magnetohydrodynamics to a class of magnetic fields that the classical theory could not handle. The author, David MacTaggart, derives a bound showing that the magnitude of relative magnetic helicity is constrained by the free magnetic energy of the field, the portion of energy available for release above the lowest-energy configuration.

The classical result concerns magnetic helicity defined for closed fields, those whose lines never cross the boundary of the domain they occupy. In that setting, helicity is an integral of the dot product of the magnetic field and its vector potential, and it is a conserved quantity of ideal magnetohydrodynamics. Work going back to seminal papers of the field established that the absolute value of this helicity cannot exceed a constant times the total magnetic energy, where the constant depends only on the geometry of the domain, not on the magnetic field itself. As helicity grows in magnitude, the magnetic energy must grow with it, whether by strengthening the field, by changing its topology, or by both at once.

The difficulty is that closed fields are a poor model for many real astrophysical situations. In a solar active region, the field lines pass through the photosphere, the visible solar surface, which acts as a lower boundary. Such fields are described as open, and for them the classical helicity is neither conserved nor gauge-invariant, meaning it depends on the arbitrary choice of vector potential. The remedy, developed decades ago by Berger and Field and by Finn and Antonsen, is relative helicity. It compares the actual field with a reference field that shares the same boundary conditions, and the resulting quantity inherits the conservation and gauge-invariance properties of the classical definition. The reference field is normally taken to be a potential field, which for a given set of boundary conditions is the unique field with the minimum possible magnetic energy.

With that choice of reference, the difference between the actual field and the potential field carries precisely the free magnetic energy, the energy that can in principle be released, for example through an instability. A first attempt to bound relative helicity by applying the Cauchy-Schwarz inequality directly does produce an inequality involving this free energy, but it fails the crucial test: the coefficient multiplying the energy depends on the magnetic field itself, through its vector potential. The classical bound’s power comes precisely from the fact that its constant depends only on the shape of the domain. A bound whose coefficient changes with the field tells you far less about what topology and energy can do.

The way around this problem is a gauge-invariant decomposition of relative helicity into two parts, a self helicity and a mutual helicity, which the paper presents as a form of the self-mutual decomposition long used for classical helicity. The self helicity involves the difference between the actual field and the reference field, a difference that is tangent to the boundary and therefore closed. This means the self helicity has exactly the form of classical helicity, and the classical machinery applies. Writing the associated vector potential through the Biot-Savart operator, which is self-adjoint, the problem reduces to an eigenvalue question. The self helicity is bounded by twice the largest-magnitude eigenvalue of that operator times the free magnetic energy, and that eigenvalue is fixed by the shape of the domain alone.

The eigenvalue in question is related to the lowest, in magnitude, eigenvalue of the force-free field equation for the domain. Analytical results for it are few, the paper notes, but it can be determined numerically both for simply connected and for multiply connected domains. Because the self-helicity bound reproduces the classical constant, the classical helicity-energy bound emerges as a special case of the new result, recovered when the field becomes everywhere tangent to the boundary, that is, when an open field becomes closed.

The mutual helicity, which measures the linkage between the reference field and the current-carrying part of the field, is handled with the Cauchy-Schwarz inequality. At first glance this looks no better than the naive attempt, since the coefficient again involves a vector potential. The crucial difference is that it involves only the vector potential of the reference field, not of the actual field. For a fixed domain with fixed boundary conditions, that reference potential can be chosen once and for all, independently of whatever magnetic field fills the domain. Choosing the Coulomb gauge, in which the reference potential is divergence-free and tangent to no normal direction on the boundary, minimizes its energy and makes the coefficient of the mutual bound a calculable constant of the geometry and boundary data alone. The paper notes that this gauge appears in diverse areas of physics and has already been used in relative helicity calculations for force-free extrapolations of solar active regions.

Combining the two pieces yields the central result: the magnitude of the relative helicity is bounded by a constant times the free magnetic energy plus a second constant times its square root. Both constants depend only on the domain geometry and the fixed distribution of the normal field component on the boundary. Because the bound is built from the self-mutual decomposition, it carries structural information as well. There is a special value of the free energy, equal to the square of the ratio of the two constants, at which the maximal possible helicity magnitude represents an equal balance of self and mutual contributions. Below that value, the maximum is dominated by the mutual helicity; above it, by the self helicity, which the paper associates with twisted field structure. A particular field may of course deviate from this pattern, but then its helicity magnitude simply falls short of the maximum the bound allows.

The practical significance lies in what the bound might reveal about solar activity. In an extrapolation of the magnetic field of an active region, the bound could help indicate the probable dynamics of the region. If the free energy lies below the special value and the measured relative helicity sits close to its maximum, the configuration may be unlikely to produce an eruption, since eruptions typically involve structures related to self helicity. If the free energy lies above the special value and the helicity is again near its maximum, the field may be highly twisted and unstable. The special energy value itself may serve as an indicator of when a given field becomes eruptive, since instabilities are typically associated with the twisted fields described by the self-helicity term. The paper is careful to frame these as possibilities for analysis rather than demonstrated forecasts; the study is theoretical, no datasets were generated or analyzed, and the author notes that the bound may be determined numerically with existing methods.

What the result delivers, in the end, is a bridge between two quantities that solar physicists measure and model constantly: the topological complexity of an active region’s field and the energy it has available to release. For any fixed domain and boundary conditions, the bound describes how relative helicity can be maximized for a given free energy, and it does so with coefficients that a modeler can compute from geometry alone, before knowing anything about the field inside. That separation, between what the domain dictates and what the field supplies, is what made the classical bound a cornerstone of ideal magnetohydrodynamics, and the new result carries that separation, for the first time, into the open fields that dominate the solar atmosphere.

Subject of Research: A theoretical derivation of a bound relating relative magnetic helicity to free magnetic energy in open magnetic fields

Article Title: A bound for relative magnetic helicity in terms of free magnetic energy

Article References: MacTaggart, D. (2025). A bound for relative magnetic helicity in terms of free magnetic energy. Discover Space, 129(1), Article 19. https://doi.org/10.1007/s11038-025-09578-8

Image Credits: AI Generated

DOI: 10.1007/s11038-025-09578-8

Keywords: magnetic helicity, relative helicity, free magnetic energy, magnetohydrodynamics, solar active regions, solar physics, force-free fields, Biot-Savart operator, Coulomb gauge, topology, theoretical physics, eruptions