Based on a decomposition of the magnetic field into potential and nonpotential components, magnetic energy and relative helicity can both also be decomposed into two quantities: potential and free energies, and volume-threading and current-carrying helicities. In this study, we perform a coupled analysis of their behaviors in a set of parametric 3D magnetohydrodynamic (MHD) simulations of solar-like eruptions. We present the general formulations for the time-varying components of energy and helicity in resistive MHD. We calculated them numerically with a specific gauge, and compared their behaviors in the numerical simulations, which differ from one another by their imposed boundary-driving motions. Thus, we investigated the impact of different active regions surface flows on the development of the energy and helicity-related quantities. Despite general similarities in their overall behaviors, helicities and energies display different evolutions that cannot be explained in a unique framework. While the energy fluxes are similar in all simulations, the physical mechanisms that govern the evolution of the helicities are markedly distinct from one simulation to another: the evolution of volume-threading helicity can be governed by boundary fluxes or helicity transfer, depending on the simulation. The eruption takes place for the same value of the ratio of the current-carrying helicity to the total helicity in all simulations. However, our study highlights that this threshold can be reached in different ways, with different helicity-related processes dominating for different photospheric flows. This means that the details of the pre-eruptive dynamics do not influence the eruption-onset helicity-related threshold. Nevertheless, the helicity-flux dynamics may be more or less efficient in changing the time required to reach the onset of the eruption.

We compare the coronal magnetic energy and helicity of two solar active regions (ARs), prolific in major eruptive (AR~11158) and confined (AR~12192) flaring, and analyze the potential of deduced proxies to forecast upcoming flares. Based on nonlinear force-free (NLFF) coronal magnetic field models with a high degree of solenoidality, and applying three different computational methods to investigate the coronal magnetic helicity, we are able to draw conclusions with a high level of confidence. Based on real observations of two solar ARs we checked trends regarding the potential eruptivity of the active-region corona, as suggested earlier in works that were based on numerical simulations, or solar observations. Our results support that the ratio of current-carrying to total helicity, $|H_\mathrm{J}|/|H_\mathrm{V}|$, shows a strong ability to indicate the eruptive potential of a solar AR. However, $|H_\mathrm{J}|/|H_\mathrm{V}|$ seems not to be indicative for the magnitude or type of an upcoming flare (confined or eruptive). Interpreted in context with earlier observational studies, our findings furthermore support that the total relative helicity normalized to the magnetic flux at the NLFF model's lower boundary, $H_\mathrm{V}/\phi^2$, represents no indicator for the eruptivity.

Context. In September 2017 the largest X-class flare of Solar Cycle 24 occurred from the most active region (AR) of this cycle, AR 12673. The AR attracted much interest because of its unique morphological and evolution characteristics. Among the parameters examined in the AR was magnetic helicity, but either only approximately, and/or intermittently. Aims. This work is interested in studying the evolution of the relative magnetic helicity and of the two components of its decomposition, the non-potential, and the volume-threading one, in the time interval around the highest activity of AR 12673. Special emphasis is given on the study of the ratio of the non-potential to total helicity, that was recently proposed as an indicator of ARs eruptivity. Methods. For these, we first approximate the coronal magnetic field of the AR with two different optimization-based extrapolation procedures, and choose the one that produces the most reliable helicity value at each instant. Moreover, in one of these methods, we weight the optimization by the uncertainty estimates derived from the Helioseismic and Magnetic Imager (HMI) instrument, for the first time. We then follow an accurate method to compute all quantities of interest. Results. The first observational determination of the evolution of the non-potential to total helicity ratio seems to confirm the quality it has in indicating eruptivity. This ratio increases before the major flares of AR 12673, and afterwards it relaxes to smaller values. Additionally, the evolution patterns of the various helicity, and energy budgets of AR 12673 are discussed and compared with other works.

Context. Magnetic helicity is an important quantity in studies of magnetized plasmas as it provides a measure of the geometrical complexity of the magnetic field in a given volume. A more detailed description of the spatial distribution of magnetic helicity is given by the field line helicity that expresses the amount of helicity associated to individual field lines, rather than in the full analysed volume. Aims. Magnetic helicity is not a gauge-invariant quantity in general, unless it is computed with respect to a reference field, yielding the so called relative magnetic helicity. The field line helicity corresponding to the relative magnetic helicity has only been examined under specific conditions so far. This work aims to define the field line helicity corresponding to relative magnetic helicity in the most general way. In addition to its general form, we provide the expression for the relative magnetic field line helicity in a few commonly used gauges, and reproduce known results as a limit of our general formulation. Methods. By starting from the definition of relative magnetic helicity, we derive the corresponding field line helicity, and we note the assumptions it is based on. Results. We check that the developed quantity reproduces relative magnetic helicity by using three different numerical simulations. For these cases we also show the morphology of field line helicity in the volume, and on the photospheric plane. As an application to solar situations, we compare the morphology of field line helicity on the photosphere with that of the connectivity-based helicity flux density in two reconstructions of an active region's magnetic field. We discuss how the derived relative magnetic field line helicity has a wide range of applications, notably in solar physics and magnetic reconnection studies.

Relative magnetic helicity is a gauge invariant quantity suitable for the study of the magnetic helicity content of heliospheric plasmas. Relative magnetic helicity can be decomposed uniquely into two gauge invariant quantities, the magnetic helicity of the non-potential component of the field, and a complementary volume-threading helicity. Recent analysis of numerical experiments simulating the generation of solar eruptions have shown that the ratio of the non-potential helicity to the total relative helicity is a clear marker of the eruptivity of the magnetic system, and that the high value of that quantity could be a sufficient condition for the onset of the instability generating the eruptions. The present study introduces the first analytical examination of the time variations of these non-potential and volume-threading helicities. The validity of the analytical formulas derived are confirmed with analysis of three-dimensional (3D) magnetohydrodynamics (MHD) simulations of solar coronal dynamics. Both the analytical investigation, and the numerical application show that, unlike magnetic helicity, the non-potential and the volume-threading helicities are not conserved quantities, even in the ideal MHD regime. A term corresponding to the transformation between the non-potential and volume-threading helicities frequently dominates their dynamics. This finding has an important consequence for their estimation in the solar corona: unlike with relative helicity, their volume coronal evolution cannot be ascertained by the flux of these quantities through the volume's boundaries. Only techniques extrapolating the 3D coronal field will enable both the proper study of the non-potential and volume-threading helicities, and the observational analysis of helicity-based solar-eruptivity proxies.

Magnetic helicity is a quantity of great importance in solar studies because it is conserved in ideal magneto-hydrodynamics. While many methods to compute magnetic helicity in Cartesian finite volumes exist, in spherical coordinates, the natural coordinate system for solar applications, helicity is only treated approximately. We present here a method to properly compute relative magnetic helicity in spherical geometry. The volumes considered are finite, of shell or wedge shape, and the three-dimensional magnetic field is considered fully known throughout the studied domain. Testing of the method with well-known, semi-analytic, force-free magnetic-field models reveals that it has excellent accuracy. Further application to a set of nonlinear force-free reconstructions of the magnetic field of solar active regions, and comparison with an approximate method used in the past, indicates that the proposed methodology can be significantly more accurate, thus making our method a promising tool in helicity studies that employ the spherical geometry. Additionally, the range of applicability of the approximate method is determined and discussed.

We study the writhe, twist and magnetic helicity of different magnetic flux ropes, based on models of the solar coronal magnetic field structure. These include an analytical force-free Titov--D\'emoulin equilibrium solution, non force-free magnetohydrodynamic simulations, and nonlinear force-free magnetic field models. The geometrical boundary of the magnetic flux rope is determined by the quasi-separatrix layer and the bottom surface, and the axis curve of the flux rope is determined by its overall orientation. The twist is computed by the Berger--Prior formula that is suitable for arbitrary geometry and both force-free and non-force-free models. The magnetic helicity is estimated by the twist multiplied by the square of the axial magnetic flux. We compare the obtained values with those derived by a finite volume helicity estimation method. We find that the magnetic helicity obtained with the twist method agrees with the helicity carried by the purely current-carrying part of the field within uncertainties for most test cases. It is also found that the current-carrying part of the model field is relatively significant at the very location of the magnetic flux rope. This qualitatively explains the agreement between the magnetic helicity computed by the twist method and the helicity contributed purely by the current-carrying magnetic field.

Magnetic helicity is a conserved quantity of ideal magneto-hydrodynamics characterized by an inverse turbulent cascade. Accordingly, it is often invoked as one of the basic physical quantities driving the generation and structuring of magnetic fields in a variety of astrophysical and laboratory plasmas. We provide here the first systematic comparison of six existing methods for the estimation of the helicity of magnetic fields known in a finite volume. All such methods are reviewed, benchmarked, and compared with each other, and specifically tested for accuracy and sensitivity to errors. To that purpose, we consider four groups of numerical tests, ranging from solutions of the three-dimensional, force-free equilibrium, to magneto-hydrodynamical numerical simulations. Almost all methods are found to produce the same value of magnetic helicity within few percent in all tests. In the more solar-relevant and realistic of the tests employed here, the simulation of an eruptive flux rope, the spread in the computed values obtained by all but one method is only 3%, indicating the reliability and mutual consistency of such methods in appropriate parameter ranges. However, methods show differences in the sensitivity to numerical resolution and to errors in the solenoidal property of the input fields. In addition to finite volume methods, we also briefly discuss a method that estimates helicity from the field lines' twist, and one that exploits the field's value at one boundary and a coronal minimal connectivity instead of a pre-defined three-dimensional magnetic-field solution.

Coronal heating through the explosive release of magnetic energy remains an open problem in solar physics. Recent hydrodynamical models attempt an investigation by placing swarms of 'nanoflares' at random sites and times in modeled one-dimensional coronal loops. We investigate the problem in three dimensions, using extrapolated coronal magnetic fields of observed solar active regions. We apply a nonlinear force-free field extrapolation above an observed photospheric magnetogram of NOAA active region (AR) 11158. We then determine the locations, energy contents, and volumes of 'unstable' areas, namely areas prone to releasing magnetic energy due to locally accumulated electric current density. Statistical distributions of these volumes and their fractal dimension are inferred, investigating also their dependence on spatial resolution. Further adopting a simple resistivity model, we infer the properties of the fractally distributed electric fields in these volumes. Next, we monitor the evolution of 10^5 particles (electrons and ions) obeying an initial Maxwellian distribution with a temperature of 10 eV, by following their trajectories and energization when subjected to the resulting electric fields. For computational convenience, the length element of the magnetic-field extrapolation is 1 arcsec, much coarser than the particles collisional mean free path in the low corona. The presence of collisions traps the bulk of the plasma around the unstable volumes, or current sheets (UCS), with only a tail of the distribution gaining substantial energy. Assuming that the distance between UCS is similar to the collisional mean free path we find that the low active-region corona is heated to 100-200 eV, corresponding to temperatures exceeding 2 MK, within tens of seconds for electrons and thousands of seconds for ions. Fractally distributed, nanoflare-triggening fragmented UCS ...

We assess the validity of the free magnetic energy - relative magnetic helicity diagram for solar magnetic structures. We used two different methods of calculating the free magnetic energy and the relative magnetic helicity budgets: a classical, volume-calculation nonlinear force-free (NLFF) method applied to finite coronal magnetic structures and a surface-calculation NLFF derivation that relies on a single photospheric or chromospheric vector magnetogram. Both methods were applied to two different data sets, namely synthetic active-region cases obtained by three-dimensional magneto-hydrodynamic (MHD) simulations and observed active-region cases, which include both eruptive and noneruptive magnetic structures. The derived energy--helicity diagram shows a consistent monotonic scaling between relative helicity and free energy with a scaling index 0.84$\pm$0.05 for both data sets and calculation methods. It also confirms the segregation between noneruptive and eruptive active regions and the existence of thresholds in both free energy and relative helicity for active regions to enter eruptive territory. We consider the previously reported energy-helicity diagram of solar magnetic structures as adequately validated and envision a significant role of the uncovered scaling in future studies of solar magnetism.

In earlier works we introduced and tested a nonlinear force-free (NLFF) method designed to self-consistently calculate the free magnetic energy and the relative magnetic helicity budgets of the corona of observed solar magnetic structures. The method requires, in principle, only a single, photospheric or low-chromospheric, vector magnetogram of a quiet-Sun patch or an active region and performs calculations in the absence of three-dimensional magnetic and velocity-field information. In this work we strictly validate this method using three-dimensional coronal magnetic fields. Benchmarking employs both synthetic, three-dimensional magnetohydrodynamic simulations and nonlinear force-free field extrapolations of the active-region solar corona. We find that our time-efficient NLFF method provides budgets that differ from those of more demanding semi-analytical methods by a factor of ~3, at most. This difference is expected from the physical concept and the construction of the method. Temporal correlations show more discrepancies that, however, are soundly improved for more complex, massive active regions, reaching correlation coefficients of the order of, or exceeding, 0.9. In conclusion, we argue that our NLFF method can be reliably used for a routine and fast calculation of free magnetic energy and relative magnetic helicity budgets in targeted parts of the solar magnetized corona. As explained here and in previous works, this is an asset that can lead to valuable insight into the physics and the triggering of solar eruptions.

We study the expected variability patterns of blazars within the two-zone acceleration model putting special emphasis on flare shapes and spectral lags. We solve semi-analytically the kinetic equations which describe the particle evolution in the acceleration and radiation zone. We then perturb the solutions by introducing Lorentzian variations in its key parameters and examine the flaring behavior of the system. We apply the above to the X-ray observations of blazar 1ES 1218+304 which exhibited a hard lag behavior during a flaring episode and discuss possibilities of producing it within the context of our model. The steady-state radio to X-rays emission of 1ES 1218+304 can be reproduced with parameters which lie well within the ones generally accepted from blazar modeling. Additionally, we find that the best way to explain its flaring behavior is by varying the rate of particles injected in the acceleration zone.

Aims: We investigate the one-zone SSC model of TeV blazars in the presence of electron acceleration. In this picture electrons reach a maximum energy where acceleration saturates from a combination of synchrotron and inverse Compton scattering losses. Methods: We solve the spatially averaged kinetic equations which describe the simultaneous evolution of particles and photons, obtaining the multi-wavelength spectrum as a function of time. Results: We apply the model to the rapid flare of Mrk 501 of July 9, 2005 as this was observed by the MAGIC telescope and obtain the relevant parameters for the pre-flare quasi steady state and the ones during the flare. We show that a hard lag flare can be obtained with parameters which lie well within the range already accepted for this source. Especially the choice of a high value of the Doppler factor seems to be necessary.

We study the acceleration and radiation of charged particles in the shock waves of supernova remnants using a recent version of the "box model". According to this, particles are accelerated in an energy-dependent region around the shock by the first order Fermi mechanism and lose energy through radiation. The particle distribution function is obtained from a spatially averaged kinetic equation that treats the energy losses self-consistently. There exists also a second population that consists of those particles that escape behind the shock where they also radiate. The energy distribution of this population is calculated in a similar manner. The application of the model to the supernova remnant RX J1713.7-3946, which was recently confirmed as a TeV source by H.E.S.S., shows that the X-ray emission can be attributed to electron synchrotron radiation while in gamma-rays there are contributions from both electrons and protons, with protons playing the dominant role. Additionally, there are strong indications that particles diffuse in turbulence that has a Kolmogorov spectrum.