We have recently passed the 75th anniversary of one of the most important results in solar and space physics: Hannes Alfv\'en's discovery of Alfv\'en waves and the Alfv\'en speed. To celebrate the anniversary, this article recounts some major episodes in the history of MHD waves. Following an initially cool reception, Alfv\'en's ideas were propelled into the spotlight by Fermi's work on cosmic rays, the new mystery of coronal heating and, as scientific perception of interplanetary space shifted dramatically and the space race started, detection of Alfv\'en waves in the solar wind. From then on, interest in MHD waves boomed, laying the foundations for modern remote observations of MHD waves in the Sun, coronal seismology and some of today's leading theories of coronal heating and solar wind acceleration. In 1970, Alfv\'en received the Nobel Prize for his work in MHD, including these discoveries. The article concludes with some reflection about what the history implies about the way we do science, especially the advantages and pitfalls of idealised mathematical models.

Alfv\'enic waves have been proposed as an important energy transport mechanism in coronal loops, capable of delivering energy to both the corona and chromosphere and giving rise to many observed features, of flaring and quiescent regions. In previous work, we established that resistive dissipation of waves (ambipolar diffusion) can drive strong chromospheric heating and evaporation, capable of producing flaring signatures. However, that model was based on a simplified assumption that the waves propagate instantly to the chromosphere, an assumption which the current work removes. Via a ray tracing method, we have implemented traveling waves in a field-aligned hydrodynamic simulation that dissipate locally as they propagate along the field line. We compare this method to and validate against the magnetohydrodynamics code Lare3D. We then examine the importance of travel times to the dynamics of the loop evolution, finding that (1) the ionization level of the plasma plays a critical role in determining the location and rate at which waves dissipate; (2) long duration waves effectively bore a hole into the chromosphere, allowing subsequent waves to penetrate deeper than previously expected, unlike an electron beam whose energy deposition rises in height as evaporation reduces the mean-free paths of the electrons; (3) the dissipation of these waves drives a pressure front that propagates to deeper depths, unlike energy deposition by an electron beam.

Solar flares are driven by the release of magnetic energy from reconnection events in the solar corona, whereafter energy is transported to the chromosphere, heating the plasma and causing the characteristic radiative losses. In the collisional thick-target model, electrons accelerated to energies exceeding 10 keV traverse the corona and impact the chromosphere, where they deposit their energy through collisions with the much denser plasma in the lower atmosphere. While there are undoubtedly high energy non-thermal electrons accelerated in flares, it is unclear whether these electron beams are the sole mechanism of energy transport, or whether they only dominate in certain phases of the flare's evolution. Alfv\'enic waves are generated during the post-reconnection relaxation of magnetic field lines, so it is important to examine their role in energy transport.

Sunquakes are the surface signatures of acoustic waves in the Sun's interior that are produced by some but not all flares and coronal mass ejections (CMEs). This paper explores a mechanism for sunquake generation by the magnetic field changes that occur during flares and CMEs, using MHD simulations with a semiempirical FAL-C atmosphere to demonstrate the generation of acoustic waves in the interior in response to changing magnetic tilt in the corona. We find that Alfven-sound resonance combined with the ponderomotive force produces acoustic waves in the interior with sufficient energy to match sunquake observations when the magnetic field angle changes by the order of 10 degrees in a region where the coronal field strength is a few hundred gauss or more. The most energetic sunquakes are produced when the coronal field is strong, while the variation of magnetic field strength with height and the time scale of the tilt change are of secondary importance.

We use radiation hydrodynamic simulations to examine two models of solar flare chromospheric heating: Alfv\'en wave dissipation and electron beam collisional losses. Both mechanisms are capable of strong chromospheric heating, and we show that the distinctive atmospheric evolution in the mid-to-upper chromosphere results in Mg II k-line emission that should be observably different between wave-heated and beam-heated simulations. We also present Ca II 8542A profiles which are formed slightly deeper in the chromosphere. The Mg II k-line profiles from our wave-heated simulation are quite different from those from a beam-heated model and are more consistent with IRIS observations. The predicted differences between the Ca II 8542A in the two models are small. We conclude that careful observational and theoretical study of lines formed in the mid-to-upper chromosphere holds genuine promise for distinguishing between competing models for chromospheric heating in flares.

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.

We have developed a numerical model of flare heating due to the dissipation of Alfv\'enic waves propagating from the corona to the chromosphere. With this model, we present an investigation of the key parameters of these waves on the energy transport, heating, and subsequent dynamics. For sufficiently high frequencies and perpendicular wave numbers, the waves dissipate significantly in the upper chromosphere, strongly heating it to flare temperatures. This heating can then drive strong chromospheric evaporation, bringing hot and dense plasma to the corona. We therefore find three important conclusions: (1) Alfv\'enic waves, propagating from the corona to the chromosphere, are capable of heating the upper chromosphere and the corona, (2) the atmospheric response to heating due to the dissipation of Alfv\'enic waves can be strikingly similar to heating by an electron beam, and (3) this heating can produce explosive evaporation.

Context: Transverse loop oscillations and loop contractions are commonly associated with solar flares, but the two types of motion have traditionally been regarded as separate phenomena. Aims: We present an observation of coronal loops contracting and oscillating following onset of a flare. We aim to explain why both behaviours are seen together and why only some of the loops oscillate. Methods: A time sequence of SDO/AIA 171 \r{A} images is analysed to identify positions of coronal loops following the onset of M6.4 flare SOL2012-03-09T03:53. We focus on five loops in particular, all of which contract during the flare, with three of them oscillating as well. A simple model is then developed for contraction and oscillation of a coronal loop. Results: We propose that coronal loop contractions and oscillations can occur in a single response to removal of magnetic energy from the corona. Our model reproduces the various types of loop motion observed and explains why the highest loops oscillate during their contraction while no oscillation is detected for the shortest contracting loops. The proposed framework suggests that loop motions can be used as a diagnostic for the removal of coronal magnetic energy by flares, while rapid decrease of coronal magnetic energy is a newly-identified excitation mechanism for transverse loop oscillations.

Predicting the final state of turbulent plasma relaxation is an important challenge, both in astrophysical plasmas such as the Sun's corona and in controlled thermonuclear fusion. Recent numerical simulations of plasma relaxation with braided magnetic fields identified the possibility of a novel constraint, arising from the topological degree of the magnetic field-line mapping. This constraint implies that the final relaxed state is drastically different for an initial configuration with topological degree 1 (which allows a Taylor relaxation) and one with degree 2 (which does not reach a Taylor state). Here we test this transition in numerical resistive-magnetohydrodynamic simulations, by embedding a braided magnetic field in a linear force-free background. Varying the background force-free field parameter generates a sequence of initial conditions with a transition between topological degree 1 and 2. For degree 1, the relaxation produces a single twisted flux tube, while for degree 2 we obtain two flux tubes. For predicting the exact point of transition, it is not the topological degree of the whole domain that is relevant, but only that of the turbulent region.

We investigate the evolution of field line helicity for magnetic fields that connect two boundaries without null points, with emphasis on localized finite-B magnetic reconnection. Total (relative) magnetic helicity is already recognized as an important topological constraint on magnetohydrodynamic processes. Field line helicity offers further advantages because it preserves all topological information and can distinguish between different magnetic fields with the same total helicity. Magnetic reconnection changes field connectivity and field line helicity reflects these changes; the goal of this paper is to characterize that evolution. We start by deriving the evolution equation for field line helicity and examining its terms, also obtaining a simplified form for cases where dynamics are localized within the domain. The main result, which we support using kinematic examples, is that during localized reconnection in a complex magnetic field, the evolution of field line helicity is dominated by a work-like term that is evaluated at the field line endpoints, namely the scalar product of the generalized field line velocity and the vector potential. Furthermore, the flux integral of this term over certain areas is very small compared to the integral of the unsigned quantity, which indicates that changes of field line helicity happen in a well-organized pairwise manner. It follows that reconnection is very efficient at redistributing helicity in complex magnetic fields despite having little effect on the total helicity.

Recent results showed that the magnetic field of M-dwarf (dM) stars, currently the main targets in searches for terrestrial planets, is very different from the solar one, both in topology as well as in intensity. In particular, the magnetised environment surrounding a planet orbiting in the habitable zone (HZ) of dM stars can differ substantially to the one encountered around the Earth. These extreme magnetic fields can compress planetary magnetospheres to such an extent that a significant fraction of the planet's atmosphere may be exposed to erosion by the stellar wind. Using observed surface magnetic maps for a sample of 15 dM stars, we investigate the minimum degree of planetary magnetospheric compression caused by the intense stellar magnetic fields. We show that hypothetical Earth-like planets with similar terrestrial magnetisation (~1G) orbiting at the inner (outer) edge of the HZ of these stars would present magnetospheres that extend at most up to 6.1 (11.7) planetary radii. To be able to sustain an Earth-sized magnetosphere, the terrestrial planet would either need to orbit significantly farther out than the traditional limits of the HZ; or else, if it were orbiting within the life-bearing region, it would require a minimum magnetic field ranging from a few G to up to a few thousand G.

We study the relationship between implosive motions in a solar flare, and the energy redistribution in the form of oscillatory structures and particle acceleration. The flare SOL2012-03-09T03:53 (M6.4) shows clear evidence for an irreversible (stepwise) coronal implosion. Extreme-ultraviolet (EUV) images show at least four groups of coronal loops at different heights overlying the flaring core undergoing fast contraction during the impulsive phase of the flare. These contractions start around a minute after the flare onset, and the rate of contraction is closely associated with the intensity of the hard X-ray (HXR) and microwave emissions. They also seem to have a close relationship with the dimming associated with the formation of the Coronal Mass Ejection (CME) and a global EUV wave. Several studies now have detected contracting motions in the corona during solar flares that can be interpreted as the implosion necessary to release energy. Our results confirm this, and tighten the association with the flare impulsive phase. We add to the phenomenology by noting the presence of oscillatory variations revealed by GOES soft X-rays (SXR) and spatially-integrated EUV emission at 94 and 335 {\AA}. We identify pulsations of $\approx 60$ seconds in SXR and EUV data, which we interpret as persistent, semi-regular compressions of the flaring core region which modulate the plasma temperature and emission measure. The loop oscillations, observed over a large region, also allow us to provide rough estimates of the energy temporarily stored in the eigenmodes of the active-region structure as it approaches its new equilibrium.

We investigate the effect of the magnetic fields of M dwarf (dM) stars on potentially habitable Earth-like planets. These fields can reduce the size of planetary magnetospheres to such an extent that a significant fraction of the planet's atmosphere may be exposed to erosion by the stellar wind. We used a sample of 15 active dM stars, for which surface magnetic-field maps were reconstructed, to determine the magnetic pressure at the planet orbit and hence the largest size of its magnetosphere, which would only be decreased by considering the stellar wind. Our method provides a fast means to assess which planets are most affected by the stellar magnetic field. We show that hypothetical Earth-like planets with similar terrestrial magnetisation (1G) orbiting at the inner (outer) edge of the habitable zone of these stars would present magnetospheres that extend at most up to 6 (11.7) planetary radii. To be able to sustain an Earth-sized magnetosphere, with the exception of only a few cases, the terrestrial planet would either (1) need to orbit significantly farther out than the traditional limits of the habitable zone; or else, (2) if it were orbiting within the habitable zone, it would require at least a magnetic field ranging from a few G to up to a few thousand G. By assuming a magnetospheric size that is more appropriate for the young-Earth (3.4Gyr ago), the required planetary magnetic fields are one order of magnitude weaker. However, in this case, the polar-cap area of the planet, which is unprotected from transport of particles to/from interplanetary space, is twice as large. As the star becomes older and, therefore, its rotation rate and magnetic field reduce, the interplanetary magnetic pressure decreases and the magnetosphere of planets probably expands. Using an empirically derived rotation-activity/magnetism relation,... (continues)

How do magnetohydrodynamic waves travel from the fully ionized corona, into and through the underlying partially ionized chromosphere, and what are the consequences for solar flares? To address these questions, we have developed a 2-fluid model (of plasma and neutrals) and used it to perform 1D simulations of Alfv\'en waves in a solar atmosphere with realistic density and temperature structure. Studies of a range of solar features (faculae, plage, penumbra and umbra) show that energy transmission from corona to chromosphere can exceed 20% of incident energy for wave periods of one second or less. Damping of waves in the chromosphere depends strongly on wave frequency: waves with periods 10 seconds or longer pass through the chromosphere with relatively little damping, however, for periods of 1 second or less, a substantial fraction (37%-100%) of wave energy entering the chromosphere is damped by ion-neutral friction in the mid and upper chromosphere, with electron resistivity playing some role in the lower chromosphere and in umbras. We therefore conclude that Alfv\'enic waves with periods of a few seconds or less are capable of heating the chromosphere during solar flares, and speculate that they could also contribute to electron acceleration or exciting sunquakes.