The plasma of the lower solar atmosphere consists of mostly neutral particles, whereas the upper solar atmosphere is mostly ionised particles and electrons. A shock that propagates upwards in the solar atmosphere therefore undergoes a transition where the dominant fluid is either neutral or ionised. An upwards propagating shock also passes a point where the sound and Alfv\'en speed are equal. At this point the energy of the acoustic shock can separated into fast and slow components. How the energy is distributed between the two modes depends on the angle of magnetic field. The separation of neutral and ionised species in a gravitationally stratified atmosphere is investigated. The role of two-fluid effects on the structure of the shocks post-mode-conversion and the frictional heating is quantified for different levels of collisional coupling. Two-fluid numerical simulations are performed using the (P\underline{I}P) code of a wave steepening into a shock in an isothermal, partially-ionised atmosphere. The collisional coefficient is varied to investigate the regimes where the plasma and neutral species are weakly, strongly and finitely coupled. The propagation speeds of the compressional waves hosted by neutral and ionised species vary, therefore velocity drift between the two species is produced as the plasma attempts to propagate faster than the neutrals. This is most extreme for a fast-mode shock. We find that the collisional coefficient drastically changes the features present in the system, specifically the mode conversion height, type of shocks present, and the finite shock widths created by the two-fluid effects. In the finitely-coupled regime fast-mode shock widths can exceed the pressure scale height leading to a new potential observable of two-fluid effects in the lower solar atmosphere.

The magnetic network extending from the photosphere (solar radius $\simeq R_\odot$) to lower corona ($R_\odot + 10$ Mm) plays an important role in the heating mechanisms of the solar atmosphere. Here we further develop the models with realistic open magnetic flux tubes of the authors in order to model more complicated configurations. Closed magnetic loops, and combinations of closed and open magnetic flux tubes are modelled. These are embedded within a stratified atmosphere, derived from observationally motivated semi-empirical and data-driven models subject to solar gravity and capable of spanning from the photosphere up into the chromosphere and lower corona. Constructing a magnetic field comprising self-similar magnetic flux tubes, an analytic solution for the kinetic pressure and plasma density is derived. Combining flux tubes of opposite polarity it is possible to create a steady background magnetic field configuration modelling a solar atmosphere exhibiting realistic stratification. The result can be applied to SOHO/MDI and SDO/HMI and other magnetograms from the solar surface, upon which photospheric motions can be simulated to explore the mechanism of energy transport. We demonstrate this powerful and versatile method with an application to Helioseismic and Magnetic Imager data.

Slow-mode shocks are important in understanding fast magnetic reconnection, jet formation and heating in the solar atmosphere, and other astrophysical systems. The atmospheric conditions in the solar chromosphere allow both ionised and neutral particles to exist and interact. Under such conditions, fine substructures exist within slow-mode shocks due to the decoupling and recoupling of the plasma and neutral species. We study numerically the fine substructure within slow-mode shocks in a partially ionised plasma, in particular, analysing the formation of an intermediate transition within the slow-mode shock. High-resolution 1D numerical simulations are performed using the (P\underline{I}P) code using a two-fluid approach. We discover that long-lived intermediate (Alfven) shocks can form within the slow-mode shock, where there is a shock transition from above to below the Alfven speed and a reversal of the magnetic field across the shock front. The collisional coupling provides frictional heating to the neutral fluid, resulting in a Sedov-Taylor-like expansion with overshoots in the neutral velocity and neutral density. The increase in density results in a decrease of the Alfven speed and with this the plasma inflow is accelerated to above the Alfven speed within the finite width of the shock leading to the intermediate transition. This process occurs for a wide range of physical parameters and an intermediate shock is present for all investigated values of plasma-$\beta$, neutral fraction, and magnetic angle. As time advances the magnitude of the magnetic field reversal decreases since the neutral pressure cannot balance the Lorentz force. The intermediate shock is long-lived enough to be considered a physical structure, independent of the initial conditions.

Synthetic intensity maps are generated from a 3D kink-unstable flux rope simulation using several DKIST/DL-NIRSP spectral lines to make a prediction of the observational signatures of energy transport and release. The reconstructed large field-of-view intensity mosaics and single tile sit-and-stare high-cadence image sequences show detailed, fine-scale structure and exhibit signatures of wave propagation, redistribution of heat, flows and fine-scale bursts. These fine-scale bursts are present in the synthetic Doppler velocity maps and can be interpreted as evidence for small-scale magnetic reconnection at the loop boundary. The spectral lines reveal the different thermodynamic structures of the loop, with the hotter lines showing the loop interior and braiding, and the cooler lines showing the radial edges of the loop. The synthetic observations of DL-NIRSP are found to preserve the radial expansion and hence the loop radius can be measured accurately. The electron number density can be estimated using the intensity ratio of the Fe~\textsc{xiii} lines at 10747 and 10798~\AA. The estimated density from this ratio is correct to within $\pm10\%$ during the later phases of the evolution, however it is less accurate initially when line-of-sight density inhomogeneities contribute to the Fe~\textsc{xiii} intensity, resulting in an overprediction of the density by $\approx30\%$. The identified signatures are all above a conservative estimate for instrument noise and therefore will be detectable. In summary, we have used forward modelling to demonstrate that the coronal off-limb mode of DKIST/DL-NIRSP will be able to detect multiple independent signatures of a kink-unstable loop and observe small-scale transient features including loop braiding/twisting and small-scale reconnection events occurring at the radial edge of the loop.

Vortex motions are frequently observed on the solar photosphere. These motions may play a key role in the transport of energy and momentum from the lower atmosphere into the upper solar atmosphere, contributing to coronal heating. The lower solar atmosphere also consists of complex networks of flux tubes that expand and merge throughout the chromosphere and upper atmosphere. We perform numerical simulations to investigate the behaviour of vortex driven waves propagating in a pair of such flux tubes in a non-force-free equilibrium with a realistically modelled solar atmosphere. The two flux tubes are independently perturbed at their footpoints by counter-rotating vortex motions. When the flux tubes merge, the vortex motions interact both linearly and nonlinearly. The linear interactions generate many small-scale transient magnetic substructures due to the magnetic stress imposed by the vortex motions. Thus, an initially monolithic tube is separated into a complex multi-threaded tube due to the photospheric vortex motions. The wave interactions also drive a superposition that increases in amplitude until it exceeds the local Mach number and produces shocks that propagate upwards with speeds of approximately $ 50$ km s$^{-1}$. The shocks act as conduits transporting momentum and energy upwards, and heating the local plasma by more than an order of magnitude, with peak temperature approximately $60,000$ K. Therefore, we present a new mechanism for the generation of magnetic waveguides from the lower solar atmosphere to the solar corona. This wave guide appears as the result of interacting perturbations in neighbouring flux tubes. Thus, the interactions of photospheric vortex motions is a potentially significant mechanism for energy transfer from the lower to upper solar atmosphere.

We investigate the onset of 2D time-dependent magnetic reconnection that is triggered using an external velocity driver located away from, and perpendicular to, an equilibrium Harris current sheet. Previous studies have typically utilised an internal trigger to initiate reconnection, e.g. initial conditions centred on the current sheet. Numerical simulations solving the compressible, resistive magnetohydrodynamics equations were performed to investigate the reconnection onset within different atmospheric layers of the Sun, namely the corona, chromosphere and photosphere. A reconnecting state is reached for all atmospheric heights considered, with the dominant physics being highly dependent on atmospheric conditions. The coronal case achieves a sharp rise in electric field for a range of velocity drivers. For the chromosphere, we find a larger velocity amplitude is required to trigger reconnection. For the photospheric environment, the electric field is highly dependent on the inflow speed; a sharp increase in electric field is obtained only as the velocity entering the reconnection region approaches the Alfven speed. Additionally, the role of ambipolar diffusion is investigated for the chromospheric case and we find that the ambipolar diffusion alters the structure of the current density in the inflow region. The rate at which flux enters the reconnection region is controlled by the inflow velocity. This determines all aspects of the reconnection start-up process, i.e. the early onset of reconnection is dominated by the advection term in Ohms law in all atmospheric layers. A lower plasma-$\beta$ enhances reconnection and creates a large change in the electric field. A high plasma-$\beta$ hinders the reconnection, yielding a sharp rise in the electric field only when the velocity flowing into the reconnection region approaches the local Alfven speed.

The signatures of energy release and energy transport for a kink-unstable coronal flux rope are investigated via forward modelling. Synthetic intensity and Doppler maps are generated from a 3D numerical simulation. The CHIANTI database is used to compute intensities for three Hinode/EIS emission lines that cover the thermal range of the loop. The intensities and Doppler velocities at simulation resolution are spatially degraded to the Hinode/EIS pixel size (1\arcsec), convolved using a Gaussian point-spread function (3\arcsec), and exposed for a characteristic time of 50 seconds. The synthetic images generated for rasters (moving slit) and sit-and-stare (stationary slit) are analysed to find the signatures of the twisted flux and the associated instability. We find that there are several qualities of a kink-unstable coronal flux rope that can be detected observationally using Hinode/EIS, namely the growth of the loop radius, the increase in intensity towards the radial edge of the loop, and the Doppler velocity following an internal twisted magnetic field line. However, EIS cannot resolve the small, transient features present in the simulation, such as sites of small-scale reconnection (e.g. nanoflares)

The acoustic resonator is an important model for explaining the three-minute oscillations in the chromosphere above sunspot umbrae. The steep temperature gradients at the photosphere and transition region provide the cavity for the acoustic resonator, which allows waves to be both partially transmitted and partially reflected. In this paper, a new method of estimating the size and temperature profile of the chromospheric cavity above a sunspot umbra is developed. The magnetic field above umbrae is modelled numerically in 1.5D with slow magnetoacoustic wave trains travelling along magnetic fieldlines. Resonances are driven by applying the random noise of three different colours---white, pink and brown---as small velocity perturbations to the upper convection zone. Energy escapes the resonating cavity and generates wave trains moving into the corona. Line of sight (LOS) integration is also performed to determine the observable spectra through SDO/AIA. The numerical results show that the gradient of the coronal spectra is directly correlated with the chromosperic temperature configuration. As the chromospheric cavity size increases, the spectral gradient becomes shallower. When LOS integrations is performed, the resulting spectra demonstrate a broadband of excited frequencies that is correlated with the chromospheric cavity size. The broadband of excited frequencies becomes narrower as the chromospheric cavity size increases. These two results provide a potentially useful diagnostic for the chromospheric temperature profile by considering coronal velocity oscillations.