The Flare-energy Distributions Generated by Kink-unstable Ensembles of Zero-net-current Coronal Loops

It has been proposed that the million degree temperature of the corona is due to the combined effect of barely-detectable energy releases, so called nanoflares, that occur throughout the solar atmosphere. Alas, the nanoflare density and brightness implied by this hypothesis means that conclusive verification is beyond present observational abilities. Nevertheless, we investigate the plausibility of the nanoflare hypothesis by constructing a magnetohydrodynamic (MHD) model that can derive the energy of a nanoflare from the nature of an ideal kink instability. The set of energy-releasing instabilities is captured by an instability threshold for linear kink modes. Each point on the threshold is associated with a unique energy release and so we can predict a distribution of nanoflare energies. When the linear instability threshold is crossed, the instability enters a nonlinear phase as it is driven by current sheet reconnection. As the ensuing flare erupts and declines, the field transitions to a lower energy state, which is modelled by relaxation theory, i.e., helicity is conserved and the ratio of current to field becomes invariant within the loop. We apply the model so that all the loops within an ensemble achieve instability followed by energy-releasing relaxation. The result is a nanoflare energy distribution. Furthermore, we produce different distributions by varying the loop aspect ratio, the nature of the path to instability taken by each loop and also the level of radial expansion that may accompany loop relaxation. The heating rate obtained is just sufficient for coronal heating. In addition, we also show that kink instability cannot be associated with a critical magnetic twist value for every point along the instability threshold

Review article: MHD wave propagation near coronal null points of magnetic fields

We present a comprehensive review of MHD wave behaviour in the neighbourhood of coronal null points: locations where the magnetic field, and hence the local Alfven speed, is zero. The behaviour of all three MHD wave modes, i.e. the Alfven wave and the fast and slow magnetoacoustic waves, has been investigated in the neighbourhood of 2D, 2.5D and (to a certain extent) 3D magnetic null points, for a variety of assumptions, configurations and geometries. In general, it is found that the fast magnetoacoustic wave behaviour is dictated by the Alfven-speed profile. In a $\beta=0$ plasma, the fast wave is focused towards the null point by a refraction effect and all the wave energy, and thus current density, accumulates close to the null point. Thus, null points will be locations for preferential heating by fast waves. Independently, the Alfven wave is found to propagate along magnetic fieldlines and is confined to the fieldlines it is generated on. As the wave approaches the null point, it spreads out due to the diverging fieldlines. Eventually, the Alfven wave accumulates along the separatrices (in 2D) or along the spine or fan-plane (in 3D). Hence, Alfven wave energy will be preferentially dissipated at these locations. It is clear that the magnetic field plays a fundamental role in the propagation and properties of MHD waves in the neighbourhood of coronal null points. This topic is a fundamental plasma process and results so far have also lead to critical insights into reconnection, mode-coupling, quasi-periodic pulsations and phase-mixing.Comment: 34 pages, 5 figures, invited review in Space Science Reviews => Note this is a 2011 paper, not a 2010 pape

A Nonlinear Force-Free Magnetic Field Approximation Suitable for Fast Forward-Fitting to Coronal Loops. I. Theory

We derive an analytical approximation of nonlinear force-free magnetic field solutions (NLFFF) that can efficiently be used for fast forward-fitting to solar magnetic data, constrained either by observed line-of-sight magnetograms and stereoscopically triangulated coronal loops, or by 3D vector-magnetograph data. The derived NLFFF solutions provide the magnetic field components $B_x({\bf x})$, $B_y({\bf x})$, $B_z({\bf x})$, the force-free parameter $\alpha({\bf x})$, the electric current density ${\bf j}({\bf x})$, and are accurate to second-order (of the nonlinear force-free $\alpha$-parameter). The explicit expressions of a force-free field can easily be applied to modeling or forward-fitting of many coronal phenomena.Comment: Solar Physics (in press), 26 pages, 11 figure

3D MHD Coronal Oscillations About a Magnetic Null Point: Application of WKB Theory

This paper is a demonstration of how the WKB approximation can be used to help solve the linearised 3D MHD equations. Using Charpit's Method and a Runge-Kutta numerical scheme, we have demonstrated this technique for a potential 3D magnetic null point, ${\bf{B}}=(x,\epsilon y -(\epsilon +1)z)$. Under our cold plasma assumption, we have considered two types of wave propagation: fast magnetoacoustic and Alfv\'en waves. We find that the fast magnetoacoustic wave experiences refraction towards the magnetic null point, and that the effect of this refraction depends upon the Alfv\'en speed profile. The wave, and thus the wave energy, accumulates at the null point. We have found that current build up is exponential and the exponent is dependent upon $\epsilon$. Thus, for the fast wave there is preferential heating at the null point. For the Alfv\'en wave, we find that the wave propagates along the fieldlines. For an Alfv\'en wave generated along the fan-plane, the wave accumulates along the spine. For an Alfv\'en wave generated across the spine, the value of $\epsilon$ determines where the wave accumulation will occur: fan-plane ($\epsilon=1$), along the $x-$axis ($0<\epsilon <1$) or along the $y-$axis ($\epsilon>1$). We have shown analytically that currents build up exponentially, leading to preferential heating in these areas. The work described here highlights the importance of understanding the magnetic topology of the coronal magnetic field for the location of wave heating.Comment: 26 pages, 12 figure

Deterministically Driven Avalanche Models of Solar Flares

We develop and discuss the properties of a new class of lattice-based avalanche models of solar flares. These models are readily amenable to a relatively unambiguous physical interpretation in terms of slow twisting of a coronal loop. They share similarities with other avalanche models, such as the classical stick--slip self-organized critical model of earthquakes, in that they are driven globally by a fully deterministic energy loading process. The model design leads to a systematic deficit of small scale avalanches. In some portions of model space, mid-size and large avalanching behavior is scale-free, being characterized by event size distributions that have the form of power-laws with index values, which, in some parameter regimes, compare favorably to those inferred from solar EUV and X-ray flare data. For models using conservative or near-conservative redistribution rules, a population of large, quasiperiodic avalanches can also appear. Although without direct counterparts in the observational global statistics of flare energy release, this latter behavior may be relevant to recurrent flaring in individual coronal loops. This class of models could provide a basis for the prediction of large solar flares.Comment: 24 pages, 11 figures, 2 tables, accepted for publication in Solar Physic

Particle acceleration in three-dimensional tearing configurations

In three-dimensional electromagnetic configurations that result from unstable resistive tearing modes particles can efficiently be accelerated to relativistic energies. To prove this resistive magnetohydrodynamic simulations are used as input configurations for successive test particle simulations. The simulations show the capability of three-dimensional non-linearly evolved tearing modes to accelerate particles perpendicular to the plane of the reconnecting magnetic field components. The simulations differ considerably from analytical approaches by involving a realistic three-dimensional electric field with a non-homogenous component parallel to the current direction. The resulting particle spectra exhibit strong pitch-angle anisotropies. Typically, about 5-8 % of an initially Maxwellian distribution is accelerated to the maximum energy levels given by the macroscopic generalized electric potential structure. Results are shown for both, non-relativistic particle acceleration that is of interest, e.g., in the context of auroral arcs and solar flares, and relativistic particle energization that is relevant, e.g., in the context of active galactic nuclei.Comment: Physics of Plasmas, in prin

Recent Advances in Understanding Particle Acceleration Processes in Solar Flares

We review basic theoretical concepts in particle acceleration, with particular emphasis on processes likely to occur in regions of magnetic reconnection. Several new developments are discussed, including detailed studies of reconnection in three-dimensional magnetic field configurations (e.g., current sheets, collapsing traps, separatrix regions) and stochastic acceleration in a turbulent environment. Fluid, test-particle, and particle-in-cell approaches are used and results compared. While these studies show considerable promise in accounting for the various observational manifestations of solar flares, they are limited by a number of factors, mostly relating to available computational power. Not the least of these issues is the need to explicitly incorporate the electrodynamic feedback of the accelerated particles themselves on the environment in which they are accelerated. A brief prognosis for future advancement is offered.Comment: This is a chapter in a monograph on the physics of solar flares, inspired by RHESSI observations. The individual articles are to appear in Space Science Reviews (2011

A Quantitative Model of Energy Release and Heating by Time-dependent, Localized Reconnection in a Flare with a Thermal Loop-top X-ray Source

We present a quantitative model of the magnetic energy stored and then released through magnetic reconnection for a flare on 26 Feb 2004. This flare, well observed by RHESSI and TRACE, shows evidence of non-thermal electrons only for a brief, early phase. Throughout the main period of energy release there is a super-hot (T>30 MK) plasma emitting thermal bremsstrahlung atop the flare loops. Our model describes the heating and compression of such a source by localized, transient magnetic reconnection. It is a three-dimensional generalization of the Petschek model whereby Alfven-speed retraction following reconnection drives supersonic inflows parallel to the field lines, which form shocks heating, compressing, and confining a loop-top plasma plug. The confining inflows provide longer life than a freely-expanding or conductively-cooling plasma of similar size and temperature. Superposition of successive transient episodes of localized reconnection across a current sheet produces an apparently persistent, localized source of high-temperature emission. The temperature of the source decreases smoothly on a time scale consistent with observations, far longer than the cooling time of a single plug. Built from a disordered collection of small plugs, the source need not have the coherent jet-like structure predicted by steady-state reconnection models. This new model predicts temperatures and emission measure consistent with the observations of 26 Feb 2004. Furthermore, the total energy released by the flare is found to be roughly consistent with that predicted by the model. Only a small fraction of the energy released appears in the super-hot source at any one time, but roughly a quarter of the flare energy is thermalized by the reconnection shocks over the course of the flare. All energy is presumed to ultimately appear in the lower-temperature T<20 MK, post-flare loops

Negatively Charged Excitons and Photoluminescence in Asymmetric Quantum Well

We study photoluminescence (PL) of charged excitons ($X^-$) in narrow asymmetric quantum wells in high magnetic fields B. The binding of all $X^-$ states strongly depends on the separation $\delta$ of electron and hole layers. The most sensitive is the ``bright'' singlet, whose binding energy decreases quickly with increasing $\delta$ even at relatively small B. As a result, the value of B at which the singlet--triplet crossing occurs in the $X^-$ spectrum also depends on $\delta$ and decreases from 35 T in a symmetric 10 nm GaAs well to 16 T for $\delta=0.5$ nm. Since the critical values of $\delta$ at which different $X^-$ states unbind are surprisingly small compared to the well width, the observation of strongly bound $X^-$ states in an experimental PL spectrum implies virtually no layer displacement in the sample. This casts doubt on the interpretation of PL spectra of heterojunctions in terms of $X^-$ recombination

Flux-rope twist in eruptive flares and CMEs : due to zipper and main-phase reconnection

Funding: UK Science and Technology Facilities CouncilThe nature of three-dimensional reconnection when a twisted flux tube erupts during an eruptive flare or coronal mass ejection is considered. The reconnection has two phases: first of all, 3D “zipper reconnection” propagates along the initial coronal arcade, parallel to the polarity inversion line (PIL); then subsequent quasi-2D “main phase reconnection” in the low corona around a flux rope during its eruption produces coronal loops and chromospheric ribbons that propagate away from the PIL in a direction normal to it. One scenario starts with a sheared arcade: the zipper reconnection creates a twisted flux rope of roughly one turn (2π radians of twist), and then main phase reconnection builds up the bulk of the erupting flux rope with a relatively uniform twist of a few turns. A second scenario starts with a pre-existing flux rope under the arcade. Here the zipper phase can create a core with many turns that depend on the ratio of the magnetic fluxes in the newly formed flare ribbons and the new flux rope. Main phase reconnection then adds a layer of roughly uniform twist to the twisted central core. Both phases and scenarios are modeled in a simple way that assumes the initial magnetic flux is fragmented along the PIL. The model uses conservation of magnetic helicity and flux, together with equipartition of magnetic helicity, to deduce the twist of the erupting flux rope in terms the geometry of the initial configuration. Interplanetary observations show some flux ropes have a fairly uniform twist, which could be produced when the zipper phase and any pre-existing flux rope possess small or moderate twist (up to one or two turns). Other interplanetary flux ropes have highly twisted cores (up to five turns), which could be produced when there is a pre-existing flux rope and an active zipper phase that creates substantial extra twist.PostprintPublisher PDFPeer reviewe