DEAP-3600 is a single-phase liquid argon detector aiming to directly detect Weakly Interacting Massive Particles (WIMPs), located at SNOLAB (Sudbury, Canada). After analyzing data taken during the first year of operation, a null result was used to place an upper bound on the WIMP-nucleon spin-independent, isoscalar cross section. This study reinterprets this result within a Non-Relativistic Effective Field Theory framework, and further examines how various possible substructures in the local dark matter halo may affect these constraints. Such substructures are hinted at by kinematic structures in the local stellar distribution observed by the Gaia satellite and other recent astronomical surveys. These include the Gaia Sausage (or Enceladus), as well as a number of distinct streams identified in recent studies. Limits are presented for the coupling strength of the effective contact interaction operators $\mathcal{O}_1$, $\mathcal{O}_3$, $\mathcal{O}_5$, $\mathcal{O}_8$, and $\mathcal{O}_{11}$, considering isoscalar, isovector, and xenonphobic scenarios, as well as the specific operators corresponding to millicharge, magnetic dipole, electric dipole, and anapole interactions. The effects of halo substructures on each of these operators are explored as well, showing that the $\mathcal{O}_5$ and $\mathcal{O}_8$ operators are particularly sensitive to the velocity distribution, even at dark matter masses above 100 GeV/$c^2$.

We present a detailed analysis of a confined filament eruption and jet associated with a C1.5 class solar flare. Multi-wavelength observations from GONG and SDO reveal the filament forming over several days following the emergence and then partial cancellation of a minority polarity spot within a decaying bipolar active region. The emergence is also associated with the formation of a 3D null point separatrix that surrounds the minority polarity. The filament eruption occurs concurrent with brightenings adjacent to and below the filament, suggestive of breakout and flare reconnection, respectively. The erupting filament material becomes partially transferred into a strong outflow jet (~ 60 km/s) along coronal loops, becoming guided back towards the surface. Utilising high resolution H$\alpha$ observations from SST/CRISP, we construct velocity maps of the outflows demonstrating their highly structured but broadly helical nature. We contrast the observations with a 3D MHD simulation of a breakout jet in a closed-field background and find close qualitative agreement. We conclude that the suggested model provides an intuitive mechanism for transferring twist/helicity in confined filament eruptions, thus validating the applicability of the breakout model not only to jets and coronal mass ejections but also to confined eruptions and flares.

We describe the design and assembly of the LUX-ZEPLIN experiment, a direct detection search for cosmic WIMP dark matter particles. The centerpiece of the experiment is a large liquid xenon time projection chamber sensitive to low energy nuclear recoils. Rejection of backgrounds is enhanced by a Xe skin veto detector and by a liquid scintillator Outer Detector loaded with gadolinium for efficient neutron capture and tagging. LZ is located in the Davis Cavern at the 4850' level of the Sanford Underground Research Facility in Lead, South Dakota, USA. We describe the major subsystems of the experiment and its key design features and requirements.

Quasi-periodic pulsations (QPPs) appear to be a common feature observed in the light curves of both solar and stellar flares. However, their quasi-periodic nature, along with the fact that they can be small in amplitude and short-lived, makes QPPs difficult to unequivocally detect. In this paper, we test the strengths and limitations of state-of-the-art methods for detecting QPPs using a series of hare-and-hounds exercises. The hare simulated a set of flares, both with and without QPPs of a variety of forms, while the hounds attempted to detect QPPs in blind tests. We use the results of these exercises to create a blueprint for anyone who wishes to detect QPPs in real solar and stellar data. We present eight clear recommendations to be kept in mind for future QPP detections, with the plethora of solar and stellar flare data from new and future satellites. These recommendations address the key pitfalls in QPP detection, including detrending, trimming data, accounting for colored noise, detecting stationary-period QPPs, detecting QPPs with nonstationary periods, and ensuring that detections are robust and false detections are minimized. We find that QPPs can be detected reliably and robustly by a variety of methods, which are clearly identified and described, if the appropriate care and due diligence are taken.

DEAP-3600 is a single-phase liquid argon (LAr) direct-detection dark matter experiment, operating 2 km underground at SNOLAB (Sudbury, Canada). The detector consists of 3279 kg of LAr contained in a spherical acrylic vessel. This paper reports on the analysis of a 758 tonne\cdot day exposure taken over a period of 231 live-days during the first year of operation. No candidate signal events are observed in the WIMP-search region of interest, which results in the leading limit on the WIMP-nucleon spin-independent cross section on a LAr target of $3.9\times10^{-45}$ cm$^{2}$ ($1.5\times10^{-44}$ cm$^{2}$) for a 100 GeV/c$^{2}$ (1 TeV/c$^{2}$) WIMP mass at 90\% C. L. In addition to a detailed background model, this analysis demonstrates the best pulse-shape discrimination in LAr at threshold, employs a Bayesian photoelectron-counting technique to improve the energy resolution and discrimination efficiency, and utilizes two position reconstruction algorithms based on PMT charge and photon arrival times.

The DEAP-3600 experiment is searching for WIMP dark matter with a 3.3 tonne single phase liquid argon (LAr) target, located 2.1 km underground at SNOLAB. The experimental signature of dark matter interactions is keV-scale $^{40}$Ar nuclear recoils (NR) producing 128 nm LAr scintillation photons observed by PMTs. The largest backgrounds in DEAP-3600 are electronic recoils (ER) induced by $\beta$ and $\gamma$-rays originating from internal and external radioactivity in the detector material. A background model of the ER interactions in DEAP-3600 was developed and is described in this work. The model is based on several components which are expected from radioisotopes in the LAr, from ex-situ material assay measurements, and from dedicated independent in-situ analyses. This prior information is used in a Bayesian fit of the ER components to a 247.2 d dataset to model the radioactivity in the surrounding detector materials. While excellent discrimination between ERs and NRs is reached with pulse shape discrimination, utilizing the large difference between fast and slow components of LAr scintillation light, detailed knowledge of the ER background and activity of detector components, sets valuable constraints on other key types of backgrounds in the detector: neutrons and alphas. In addition, the activity of $^{42}$Ar in LAr in DEAP-3600 is determined by measuring the daughter decay of $^{42}$K. This cosmogenically activated trace isotope is a relevant background at higher energies for other rare event searches using atmospheric argon e.g. DarkSide-20k, GERDA or LEGEND. The specific activity of $^{42}$Ar in the atmosphere is found to be $40.4 \pm 5.9$ $\mu$Bq/kg of argon.

Alfv\'enic waves have gained renewed interest since the existence of ubiquitous propagating kink waves were discovered in the corona. {It has long been suggested that Alfv\'enic} waves play an important role in coronal heating and the acceleration of the solar wind. To this effect, it is imperative to understand the mechanisms that enable their energy to be transferred to the plasma. Mode conversion via resonant absorption is believed to be one of the main mechanisms for kink wave damping, and is considered to play a key role in the process of energy transfer. This study examines the damping of propagating kink waves in quiescent coronal loops using the Coronal Multi-channel Polarimeter (CoMP). A coherence-based method is used to track the Doppler velocity signal of the waves, enabling us to investigate the spatial evolution of velocity perturbations. The power ratio of outward to inward propagating waves is used to estimate the associated damping lengths and quality factors. To enable accurate estimates of these quantities, {we provide the first derivation of a likelihood function suitable for fitting models to the ratio of two power spectra obtained from discrete Fourier transforms. Maximum likelihood estimation is used to fit an exponential damping model to the observed variation in power ratio as a function of frequency.} We confirm earlier indications that propagating kink waves are undergoing frequency dependent damping. Additionally, we find that the rate of damping decreases, or equivalently the damping length increases, for longer coronal loops that reach higher in the corona.

Many cool stars possess complex magnetic fields [1] that are considered to undertake a central role in the structuring and energising of their atmospheres [2]. Alfv\'enic waves are thought to make a critical contribution to energy transfer along these magnetic fields, with the potential to heat plasma and accelerate stellar winds [3] [4] [5]. Despite Alfv\'enic waves having been identified in the Sun's atmosphere, the nature of the basal wave energy flux is poorly understood. It is generally assumed that the associated Poynting flux is generated solely in the photosphere and propagates into the corona, typically through the continuous buffeting of magnetic fields by turbulent convective cells [4] [6] [7]. Here we provide evidence that the Sun's internal acoustic modes also contribute to the basal flux of Alfv\'enic waves, delivering a spatially ubiquitous input to the coronal energy balance that is sustained over the solar cycle. Alfv\'enic waves are thus a fundamental feature of the Sun's corona. Acknowledging that internal acoustic modes have a key role in injecting additional Poynting flux into the upper atmospheres of Sun-like stars has potentially significant consequences for the modelling of stellar coronae and winds.

(Abridged for ArXiv) An injection of energy towards a magnetic null point can drive reversals of current sheet polarity leading to time-dependent Oscillatory Reconnection, which is a possible explanation of how periodic phenomena can be generated when reconnection occurs in the solar atmosphere. However, the details of what controls the period of these oscillations is poorly understood, despite being crucial in assessing whether OR can account for observed periodic behaviour. This paper aims to highlight that different types of reconnection reversal are supported about null points, and that these are distinct from the oscillation on the closed-boundary, linear systems considered in the 1990s. In particular, we explore the features of a nonlinear oscillation local to the null point, and examine the effect of resistivity and perturbation energy on the period, contrasting it to the linear case. It is found that in the linear systems, the inverse Lundquist number dictates the period, provided the perturbation energy is small relative to the inverse Lundquist number defined on the boundary, regardless of the broadband structure of the initial perturbation. However, when the perturbation energy exceeds the threshold required for 'nonlinear' null collapse to occur, a complex oscillation of the magnetic field is produced which is, at best, only weakly-dependent on the resistivity. The resultant periodicity is strongly influenced by the amount of free energy, with more energetic perturbations producing higher-frequency oscillations. Crucially, with regards to typical solar-based and astrophysical-based input energies, we demonstrate that the majority far exceed the threshold for nonlinearity to develop. This substantially alters the properties and periodicity of both null collapse and subsequent OR. Therefore, nonlinear regimes of OR should be considered in solar and astrophysical contexts.

This paper reports the first results of a direct dark matter search with the DEAP-3600 single-phase liquid argon (LAr) detector. The experiment was performed 2 km underground at SNOLAB (Sudbury, Canada) utilizing a large target mass, with the LAr target contained in a spherical acrylic vessel of 3600 kg capacity. The LAr is viewed by an array of PMTs, which would register scintillation light produced by rare nuclear recoil signals induced by dark matter particle scattering. An analysis of 4.44 live days (fiducial exposure of 9.87 tonne-days) of data taken with the nearly full detector during the initial filling phase demonstrates the detector performance and the best electronic recoil rejection using pulse-shape discrimination in argon, with leakage $<1.2\times 10^{-7}$ (90% C.L.) between 16 and 33 keV$_{ee}$. No candidate signal events are observed, which results in the leading limit on WIMP-nucleon spin-independent cross section on argon, $<1.2\times 10^{-44}$ cm$^2$ for a 100 GeV/c$^2$ WIMP mass (90% C.L.).

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.

Recent instrumentation has demonstrated that the solar atmosphere supports omnipresent transverse waves, which could play a key role in energizing the solar corona. Large-scale studies are required in order to build up an understanding of the general properties of these transverse waves. To help facilitate this, we present an automated algorithm for identifying and tracking features in solar images and extracting the wave properties of any observed transverse oscillations. We test and calibrate our algorithm using a set of synthetic data which includes noise and rotational effects. The results indicate an accuracy of 1-2% for displacement amplitudes and 4-10% for wave periods and velocity amplitudes. We also apply the algorithm to data from the Atmospheric Imaging Assembly (AIA) on board the $\textit{Solar Dynamics Observatory}$ (SDO) and find good agreement with previous studies. Of note, we find that 35-41% of the observed plumes exhibit multiple wave signatures, which indicates either the superposition of waves or multiple independent wave packets observed at different times within a single structure. The automated methods described in this paper represent a significant improvement on the speed and quality of direct measurements of transverse waves within the solar atmosphere. This algorithm unlocks a wide range of statistical studies that were previously impractical.

Implosive formation of current sheets is a fundamental plasma process. Previous studies focused on the early time evolution, while here our primary aim is to explore the longer-term evolution, which may be critical for determining the efficiency of energy release. To address this problem we investigate two closely-related problems, namely: (i) 1D, pinched anti-parallel magnetic fields and (ii) 2D, null point containing fields which are locally imbalanced ('null-collapse' or 'X-point collapse'). Within the framework of resistive MHD, we simulate the full nonlinear evolution through three distinct phases: the initial implosion, its eventual halting mechanism, and subsequent evolution post-halting. In a parameter study, we find the scaling with resistivity of current sheet properties at the halting time is in good agreement - in both geometries - with that inferred from a known 1D similarity solution. We find that the halting of the implosions occurs rapidly after reaching the diffusion scale by sudden Ohmic heating of the dense plasma within the current sheet, which provides a pressure gradient sufficient to oppose further collapse and decelerate the converging flow. This back-pressure grows to exceed that required for force balance and so the post-implosion evolution is characterised by the consequences of the current sheet `bouncing' outwards. These are: (i) the launching of propagating fast MHD waves (shocks) outwards and (ii) the width-wise expansion of the current sheet itself. The expansion is only observed to stall in the 2D case, where the pressurisation is relieved by outflow in the reconnection jets. In the 2D case, we quantify the maximum amount of current sheet expansion as it scales with resistivity, and analyse the structure of the reconnection region which forms post-expansion, replete with Petschek-type slow shocks and fast termination shocks.

The Dark matter Experiment using Argon Pulse-shape discrimination (DEAP) has been designed for a direct detection search for particle dark matter using a single-phase liquid argon target. The projected cross section sensitivity for DEAP-3600 to the spin-independent scattering of Weakly Interacting Massive Particles (WIMPs) on nucleons is $10^{-46}~\rm{cm}^{2}$ for a 100 GeV/$c^2$ WIMP mass with a fiducial exposure of 3 tonne-years. This paper describes the physical properties and construction of the DEAP-3600 detector.

Null collapse is an implosive process whereby MHD waves focus their energy in the vicinity of a null point, forming a current sheet and initiating magnetic reconnection. We consider, for the first time, the case of collapsing 3D magnetic null points in nonlinear, resistive MHD using numerical simulation, exploring key physical aspects of the system as well as performing a detailed parameter study. We find that within a particular plane containing the 3D null, the plasma and current density enhancements resulting from the collapse are quantitatively and qualitatively as per the 2D case in both the linear and nonlinear collapse regimes. However, the scaling with resistivity of the 3D reconnection rate - which is a global quantity - is found to be less favourable when the magnetic null point is more rotationally symmetric, due to the action of increased magnetic back-pressure. Furthermore, we find that with increasing ambient plasma pressure the collapse can be throttled, as is the case for 2D nulls. We discuss this pressure-limiting in the context of fast reconnection in the solar atmosphere and suggest mechanisms by which it may be overcome. We also discuss the implications of the results in the context of null collapse as a trigger mechanism of Oscillatory Reconnection, a time-dependent reconnection mechanism, and also within the wider subject of wave-null point interactions. We conclude that, in general, increasingly rotationally-asymmetric nulls will be more favourable in terms of magnetic energy release via null collapse than their more symmetric counterparts.

Solar flare emission is detected in all EM bands and variations in flux density of solar energetic particles. Often the EM radiation generated in solar and stellar flares shows a pronounced oscillatory pattern, with characteristic periods ranging from a fraction of a second to several minutes. These oscillations are referred to as quasi-periodic pulsations (QPPs), to emphasise that they often contain apparent amplitude and period modulation. We review the current understanding of quasi-periodic pulsations in solar and stellar flares. In particular, we focus on the possible physical mechanisms, with an emphasis on the underlying physics that generates the resultant range of periodicities. These physical mechanisms include MHD oscillations, self-oscillatory mechanisms, oscillatory reconnection/reconnection reversal, wave-driven reconnection, two loop coalescence, MHD flow over-stability, the equivalent LCR-contour mechanism, and thermal-dynamical cycles. We also provide a histogram of all QPP events published in the literature at this time. The occurrence of QPPs puts additional constraints on the interpretation and understanding of the fundamental processes operating in flares, e.g. magnetic energy liberation and particle acceleration. Therefore, a full understanding of QPPs is essential in order to work towards an integrated model of solar and stellar flares.

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.

Here we detail the dynamic evolution of localised reconnection regions about three-dimensional (3D) magnetic null points by using numerical simulation. We demonstrate for the first time that reconnection triggered by the localised collapse of a 3D null point due to an external MHD wave involves a self-generated oscillation, whereby the current sheet and outflow jets undergo a reconnection reversal process during which back-pressure formation at the jet heads acts to prise open the collapsed field before overshooting the equilibrium into an opposite-polarity configuration. The discovery that reconnection at fully 3D nulls can proceed naturally in a time-dependent and periodic fashion is suggestive that oscillatory reconnection mechanisms may play a role in explaining periodicity in astrophysical phenomena associated with magnetic reconnection, such as the observed quasi-periodicity of solar and stellar flare emission. Furthermore, we find a consequence of oscillatory reconnection is the generation of a plethora of freely-propagating MHD waves which escape the vicinity of the reconnection region

Determining the mechanisms responsible for the heating of the coronal plasma and maintaining and accelerating the solar wind are long standing goals in solar physics. There is a clear need to constrain the energy, mass and momentum flux through the solar corona and advance our knowledge of the physical process contributing to these fluxes. Furthermore, the accurate forecasting of Space Weather conditions at the near-Earth environment and, more generally, the plasma conditions of the solar wind throughout the heliosphere, require detailed knowledge of these fluxes in the near-Sun corona. Here we present a short case for a space-based imaging-spectrometer coronagraph, which will have the ability to provide synoptic information on the coronal environment and provide strict constraints on the mass, energy, and momentum flux through the corona. The instrument would ideally achieve cadences of $\sim10$~s, spatial resolution of 1" and observe the corona out to 2~$R_{\sun}$. Such an instrument will enable significant progress in our understanding of MHD waves throughout complex plasmas, as well as potentially providing routine data products to aid Space Weather forecasting.

The propagation of magnetoacoustic waves in the neighbourhood of a 2D null point is investigated for both $\beta=0$ and $\beta \neq 0$ plasmas. Previous work has shown that the Alfv\'en speed, here $v_A \propto r$, plays a vital role in such systems and so a natural choice is to switch to polar coordinates. For the $\beta=0$ plasma, we derive an analytical solution for the behaviour of the fast magnetoacoustic wave in terms of the Klein-Gordon equation. We also solve the system with a semi-analytical WKB approximation which shows that the $\beta=0$ wave focuses on the null and contracts around it but, due to exponential decay, never reaches the null in a finite time. For the $\beta \neq 0$ plasma, we solve the system numerically and find the behaviour to be similar to that of the $\beta=0$ system at large radii, but completely different close to the null. We show that for an initially cylindrically-symmetric fast magnetoacoustic wave perturbation, there is a decrease in wave speed along the separatrices and so the perturbation starts to take on a quasi-diamond shape; with the corners located along the separatrices. This is due to the growth in pressure gradients that reach a maximum along the separatrices, which in turn reduces the acceleration of the fast wave along the separatrices leading to a deformation of the wave morphology.

We study the propagation of a fast magnetoacoustic wave in a 3D magnetic field created from two magnetic dipoles. The magnetic topology contains an X-line. We aim to contribute to the overall understanding of MHD wave propagation within inhomogeneous media, specifically around X-lines. We investigate the linearised, 3D MHD equations under the assumptions of ideal and cold plasma. We utilise the WKB approximation and Charpit's method during our investigation. It is found that the behaviour of the fast magnetoacoustic wave is entirely dictated by the local, inhomogeneous, equilibrium Alfv\'en speed profile. All parts of the wave experience refraction during propagation, where the magnitude of the refraction effect depends on the location of an individual wave element within the inhomogeneous magnetic field. The X-line, along which the Alfv\'en speed is identically zero, acts as a focus for the refraction effect. There are two main types of wave behaviour: part of the wave is either trapped by the X-line or escapes the system, and there exists a critical starting region around the X-line that divides these two types of behaviour. For the set-up investigated, it is found that $15.5\%$ of the fast wave energy is trapped by the X-line. We conclude that linear, $\beta=0$ fast magnetoacoustic waves can accumulate along X-lines and thus these will be specific locations of fast wave energy deposition and thus preferential heating. The work here highlights the importance of understanding the magnetic topology of a system. We also demonstrate how the 3D WKB technique described in this paper can be applied to other magnetic configurations.

The propagation of linear Alfv\'en wave pulses in an inhomogeneous plasma near a 2D coronal null point is investigated. When a uniform plasma density is considered, it is seen that an initially planar Alfv\'en wavefront remains planar, despite the varying equilibrium Alfv\'en speed, and that all the wave collects at the separatrices. Thus, in the non-ideal case, these Alfv\'enic disturbances preferentially dissipate their energy at these locations. For a non-uniform equilibrium density, it is found that the Alfv\'en wavefront is significantly distorted away from the initially planar geometry, inviting the possibility of dissipation due to phase mixing. Despite this however, we conclude that for the Alfv\'en wave, current density accumulation and preferential heating still primarily occur at the separatrices, even when an extremely non-uniform density profile is considered.

There is intense interest in determining the precise contribution of Alfv\'enic waves propagating along solar structures to the problems of coronal heating and solar wind acceleration. Since the launch of SDO/AIA, it has been possible to resolve transverse oscillations in off-limb solar polar plumes and recently McIntosh et al. (2011, Nature, 475, 477) concluded that such waves are energetic enough to play a role in heating the corona and accelerating the fast solar wind. However, this result is based on comparisons to Monte Carlo simulations and confirmation via direct measurements is still outstanding. Thus, this letter reports on the first direct measurements of transverse wave motions in solar polar plumes. Over a 4 hour period, we measure the transverse displacements, periods and velocity amplitudes of 596 distinct oscillations observed in the 171 \AA channel of SDO/AIA. We find a broad range of non-uniformly distributed parameter values which are well described by log-normal distributions with peaks at $234$ km, $121$ s and $8$ km s$^{-1}$, and mean and standard deviations of $407\pm297$ km, $173\pm118$ s and $14\pm10$ km s$^{-1}$. Within standard deviations, our direct measurements are broadly consistent with previous results. However, accounting for the whole of our observed non-uniform parameter distribution we calculate an energy flux of $9-24$ W m$^{-2}$, which is $4-10$ times below the energy requirement for solar wind acceleration. Hence, our results indicate that transverse MHD waves as resolved by SDO/AIA cannot be the dominant energy source for fast solar wind acceleration in the open-field corona.

The \textit{High resolution Coronal Imager (Hi-C)} has provided the sharpest view of the EUV corona to date. In this paper we exploit its impressive resolving power to provide the first analysis of the fine-scale structure of moss in an active region. The data reveal that the moss is made up of a collection of fine threads, that have widths with a mean and standard deviation of $440\pm190$~km (Full Width Half Maximum). {The brightest moss emission is located at the visible head of the fine-scale structure and the fine structure appears to extend into the lower solar atmosphere.} The emission decreases along the features implying the lower sections are most likely dominated by cooler transition region plasma. These threads appear to be the cool, lower legs of the hot loops. In addition, the increased resolution allows for the first direct observation {of physical displacements of the moss fine-structure in a direction transverse to its central axis. Some of these transverse displacements demonstrate periodic behaviour, which we interpret as a signature of kink (Alfv\'enic) waves. Measurements of the properties of the transverse motions are made and the wave motions have} means and standard deviations of $55\pm37$~km for the transverse displacement amplitude, $77\pm33$~s for the period and $4.7\pm2.5$~km/s for the velocity amplitude. The presence of waves in the transition region of hot loops could have important implications for the heating of active regions.

Context: MHD waves and magnetic null points are both prevalent in many astrophysical plasmas, including the solar atmosphere. Interaction between waves and null points has been implicated as a possible mechanism for localised heating events. Aims: Here we investigate the transient behaviour of the Alfven wave about fully 3D proper and improper 3D magnetic null points. Previously, the behaviour of fast magnetoacoustic waves at null points in 3D, cold MHD was considered by Thurgood & McLaughlin (Astronomy & Astrophysics, 2012, 545, A9). Methods: We introduce an Alfven wave into the vicinity of both proper and improper null points by numerically solving the ideal, $\beta=0$ MHD equations using the LARE3D code. A magnetic fieldline and flux-based coordinate system permits the isolation of resulting wave-modes and the analysis of their interaction. Results: We find that the Alfven wave propagates throughout the region and accumulates near the fan-plane, causing current build up. For different values of null point eccentricity, the qualitative behaviour changes only by the imposition of anisotropic pulse dilation, due to the differing rates at which fieldlines diverge from the spine. For all eccentricities, we find that the fast and Alfven waves are linearly decoupled. During the driving phase, an independently propagating fast wave is nonlinearly generated due to the ponderomotive force. Subsequently, no further excitation of fast waves occurs. Conclusions: We find that the key aspects of the theory of Alfv\'en waves about 2D null points extends intuitively to the fully 3D case; i.e. the wave propagates along fieldlines and thus accumulates at predictable parts of the topology. We also highlight that unlike in the 2D case, in 3D Alfven-wave pulses are always toroidal and thus any aspects of 2D Alfven-wave-null models that are pulse-geometry specific must be reconsidered in 3D.

Aims. Linear magnetohydrostatic (MHS) models of solar magnetic fields balance plasma pressure gradients, gravity and Lorentz forces where the current density is composed of a linear force-free component and a cross-field component that depends on gravitational stratification. In this paper, we investigate an efficient numerical procedure for calculating such equilibria. Methods. The MHS equations are reduced to two scalar elliptic equations - one on the lower boundary and the other within the interior of the computational domain. The normal component of the magnetic field is prescribed on the lower boundary and a multigrid method is applied on both this boundary and within the domain to find the poloidal scalar potential. Once solved to a desired accuracy, the magnetic field, plasma pressure and density are found using a finite difference method. Results. We investigate the effects of the cross-field currents on the linear MHS equilibria. Force-free and non-force-free examples are given to demonstrate the numerical scheme and an analysis of speed-up due to parallelization on a graphics processing unit (GPU) is presented. It is shown that speed-ups of x30 are readily achievable.

The recent launch of the High resolution Coronal imager (Hi-C) provided a unique opportunity to study the EUV corona with unprecedented spatial resolution. We utilize these observations to investigate the properties of low frequency (50-200 s) active region transverse waves, whose omnipresence had been suggested previously. The five-fold improvement in spatial resolution over SDO/AIA reveals coronal loops with widths 150-310 km and that these loops support transverse waves with displacement amplitudes <50 km. However, the results suggest that wave activity in the coronal loops is of low energy, with typical velocity amplitudes <3 km/s. An extended time-series of SDO data suggest low-energy wave behaviour is typical of the coronal structures both before and after the Hi-C observations.

Context : In the linear, {\beta}=0 MHD regime, the transient properties of MHD waves in the vicinity of 2D null points are well known. The waves are decoupled and accumulate at predictable parts of the magnetic topology: fast waves accumulate at the null point; whereas Alfv\'en waves cannot cross the separatricies. However, in nonlinear MHD mode conversion can occur at regions of inhomogeneous Alfv\'en speed, suggesting that the decoupled nature of waves may not extend to the nonlinear regime. Aims: We investigate the behaviour of low-amplitude Alfv\'en waves about a 2D magnetic null point in nonlinear, {\beta}= 0 MHD. Methods: We numerically simulate the introduction of low-amplitude Alfv\'en waves into the vicinity of a magnetic null point using the nonlinear LARE2D code. Results: Unlike in the linear regime, we find that the Alfv\'en wave sustains cospatial daughter disturbances, manifest in the transverse and longitudinal fluid velocity, owing to the action of nonlinear magnetic pressure gradients (viz. the ponderomotive force). These disturbances are dependent on the Alfv\'en wave and do not interact with the medium to excite magnetoacoustic waves, although the transverse daughter becomes focused at the null point. Additionally, an independently propagating fast magnetoacoustic wave is generated during the early stages, which transports some of the initial Alfv\'en wave energy towards the null point. Subsequently, despite undergoing dispersion and phase-mixing due to gradients in the Alfv\'en-speed profile (Grad c_A \neq 0) there is no further nonlinear generation of fast waves. Conclusions: We find that Alfv\'en waves at 2D cold null points behave largely as in the linear regime, however they sustain transverse and longitudinal disturbances - effects absent in the linear regime - due to nonlinear magnetic pressure gradients.

Where spatial gradients in the amplitude of an Alfven wave are non-zero, a nonlinear magnetic-pressure gradient acts upon the medium (the ponderomotive force). We investigate the nature of such a force in inhomogeneous 2.5D MHD plasmas by analysing source terms in the nonlinear wave equations for the general case of inhomogeneous B and {\rho}, and consider supporting nonlinear numerical simulations. Our equations indicate there are two distinct classes of ponderomotive effect induced by Alfven waves in general 2.5D MHD, each with both a longitudinal and transverse manifestation; i) Geometric Effects: Gradients in the pulse geometry relative to the background magnetic field cause the wave to sustain cospatial disturbances, the longitudinal and transverse daughter disturbances - where we report on the transverse disturbance for the first time. ii) Grad (c_A) Effects: Where a pulse propagates through an inhomogeneous region (where the non-zero gradients in the Alfven-speed profile are non-zero), the nonlinear magnetic-pressure gradient acts to accelerate the plasma. Transverse gradients (phase mixing regions) excite independently propagating fast magnetoacoustic waves (generalising the result of Nakariakov et al. 1997, Solar Physics, 175, 93) and longitudinal gradients perturb along the field (thus creating static disturbances in {\beta} = 0, and slow waves in {\beta} \= 0). We additionally demonstrate that mode conversion due the nonlinear Lorentz force is a one-way process, and does not act as a mechanism to nonlinearly generate Alfven waves due to propagating magnetoacoustic waves. We conclude that these ponderomotive effects are induced by an Alfven wave propagating in any MHD medium, and have the potential to have significant consequences on the dynamics of energy transport and aspects of dissipation provided the system is sufficiently nonlinear and inhomogeneous.

Oscillatory reconnection is a time-dependent magnetic reconnection mechanism that naturally produces periodic outputs from aperiodic drivers. This paper aims to quantify and measure the periodic nature of oscillatory reconnection for the first time. We solve the compressible, resistive, nonlinear MHD equations using 2.5D numerical simulations. We identify two distinct periodic regimes: the impulsive and stationary phases. In the impulsive phase, we find the greater the amplitude of the initial velocity driver, the longer the resultant current sheet and the earlier its formation. In the stationary phase, we find that the oscillations are exponentially decaying and for driving amplitudes 6.3 - 126.2 km/s, we measure stationary-phase periods in the range 56.3 - 78.9 s, i.e. these are high frequency (0.01 - 0.02 Hz) oscillations. In both phases, we find that the greater the amplitude of the initial velocity driver, the shorter the resultant period, but note that different physical processes and periods are associated with both phases. We conclude that the oscillatory reconnection mechanism behaves akin to a damped harmonic oscillator.

We investigate the excitation of magnetoacoustic-gravity waves generated from localized pulses in the gas pressure as well as in vertical component of velocity. These pulses are initially launched at the top of the solar photosphere that is permeated by a weak magnetic field. We investigate three different configurations of the background magnetic field lines: horizontal, vertical and oblique to the gravitational force. We numerically model magnetoacoustic-gravity waves by implementing a realistic (VAL-C) model of solar temperature. We solve two-dimensional ideal magnetohydrodynamic equations numerically with the use of the FLASH code to simulate the dynamics of the lower solar atmosphere. The initial pulses result in shocks at higher altitudes. Our numerical simulations reveal that a small-amplitude initial pulse can produce magnetoacoustic-gravity waves, which are later reflected from the transition region due to the large temperature gradient. The atmospheric cavities in the lower solar atmosphere are found to be the ideal places that may act as a resonator for various oscillations, including their trapping and leakage into the higher atmosphere. Our numerical simulations successfully model the excitation of such wave modes, their reflection and trapping, as well as the associated plasma dynamics.

Context: Coronal magnetic null points have been implicated as possible locations for localised heating events in 2D models. We investigate this possibility about fully 3D null points. Aims: We investigate the nature of the fast magnetoacoustic wave about a fully 3D magnetic null point, with a specific interest in its propagation, and we look for evidence of MHD mode coupling and/or conversion to the Alfv\'en mode. Methods: A special fieldline and flux-based coordinate system was constructed to permit the introduction of a pure fast magnetoacoustic wave in the vicinity of proper and improper 3D null points. We considered the ideal, {\beta} = 0, MHD equations, which are solved using the LARE3D numerical code. The constituent modes of the resulting wave were isolated and identified using the special coordinate system. Numerical results were supported by analytical work derived from perturbation theory and a linear implementation of the WKB method. Results: An initially pure fast wave is found to be permanently decoupled from the Alfv\'en mode both linearly and nonlinearly for both proper and improper 3D null points. The pure fast mode also generates and sustains a nonlinear disturbance aligned along the equilibrium magnetic field. The resulting pure fast magnetoacoustic pulse has transient behaviour, which is found to be governed by the (equilibrium) Alfv\'en-speed profile, and a refraction effect focuses all the wave energy towards the null point. Conclusions: Thus, the main results from previous 2D work do indeed carry over to the fully 3D magnetic null points and

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.

We investigate the long-term evolution of an initially buoyant magnetic flux tube emerging into a gravitationally-stratified coronal hole environment and report on the resulting oscillations and outflows. We perform 2.5D nonlinear numerical simulations, generalizing the models of McLaughlin et al. (2009) and Murray et al. (2009). We find that the physical mechanism of oscillatory reconnection naturally generates quasi-periodic vertical outflows, with a transverse/swaying aspect. The vertical outflows consist of both a periodic aspect and evidence of a positively-directed flow. The speed of the vertical outflow (20-60 km/s) is comparable to those reported in the observational literature. We also perform a parametric study varying the magnetic strength of the buoyant flux tube and find a range of associated periodicities: 1.75-3.5 min. Thus, the mechanism of oscillatory reconnection may provide a physical explanation to some of the high-speed, quasi-periodic, transverse outflows/jets recently reported by a multitude of authors and instruments.

A thin dark thread is observed in a UV/EUV solar jet in the 171A, 193A and 211A and partially in 304A. The dark thread appears to originate in the chromosphere but its temperature does not appear to lie within the passbands of the Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory (SDO). We therefore implement solar magnetoseismology to estimate the plasma parameters of the dark thread. A propagating fast kink (transverse) wave is observed to travel along the dark thread. The wave is tracked over a range of ~7000km by placing multiple slits along the axis of the dark thread. The phase speed and amplitude of the wave are estimated and magnetoseismological theory is employed to determine the plasma parameters. We are able to estimate the plasma temperature, density gradient, magnetic field gradient and sub-resolution expansion of the dark thread. The dark thread is found to be cool, T~3x10^4$ with both strong density and magnetic field gradients. The expansion of the flux tube along its length is ~300-400km.

We investigate the behaviour of nonlinear, nonideal Alfv\'en wave propagation within an inhomogeneous magnetic environment. The governing MHD equations are solved in 1D and 2D using both analytical techniques and numerical simulations. We find clear evidence for the ponderomotive effect and visco-resistive heating. The ponderomotive effect generates a longitudinal component to the transverse Alfv\'en wave, with a frequency twice that of the driving frequency. Analytical work shows the addition of resistive heating. This leads to a substantial increase in the local temperature and thus gas pressure of the plasma, resulting in material being pushed along the magnetic field. In 2D, our system exhibits phase mixing and we observe an evolution in the location of the maximum heating, i.e. we find a drifting of the heating layer. Considering Alfv\'en wave propagation in 2D with an inhomogeneous density gradient, we find that the equilibrium density profile is significantly modified by both the flow of density due to visco-resistive heating and the nonlinear response to the localised heating through phase mixing.

This paper extends the models of Craig & McClymont (1991) and McLaughlin & Hood (2004) to include finite $\beta$ and nonlinear effects. We investigate the nature of nonlinear fast magnetoacoustic waves about a 2D magnetic X-point. We solve the compressible and resistive MHD equations using a Lagrangian remap, shock capturing code (Arber et al. 2001) and consider an initial condition in $ {\bf{v}}\times{\bf{B}} \cdot {\hat{\bf{z}}}$ (a natural variable of the system). We observe the formation of both fast and slow oblique magnetic shocks. The nonlinear wave deforms the X-point into a 'cusp-like' point which in turn collapses to a current sheet. The system then evolves through a series of horizontal and vertical current sheets, with associated changes in connectivity, i.e. the system exhibits oscillatory reconnection. Our final state is non-potential (but in force balance) due to asymmetric heating from the shocks. Larger amplitudes in our initial condition correspond to larger values of the final current density left in the system. The inclusion of nonlinear terms introduces several new features to the system that were absent from the linear regime.

At this time there does not exist a robust set of rules connecting low and high $\beta$ waves across the $\beta \approx 1$ layer. The work here contributes specifically to what happens when a low $\beta$ fast wave crosses the $\beta \approx 1$ layer and transforms into high $\beta$ fast and slow waves. The nature of fast and slow magnetoacoustic waves is investigated in a finite $\beta$ plasma in the neighbourhood of a two-dimensional null point. The linearised equations are solved in both polar and cartesian forms with a two-step Lax-Wendroff numerical scheme. Analytical work (e.g. small $\beta$ expansion and WKB approximation) also complement the work. It is found that when a finite gas pressure is included in magnetic equilibrium containing an X-type null point, a fast wave is attracted towards the null by a refraction effect and that a slow wave is generated as the wave crosses the $\beta \approx 1$ layer. Current accumulation occurs close to the null and along nearby separatrices. The fast wave can now \emph{pass through the origin} due to the non-zero sound speed, an effect not previously seen in related papers but clear seen for larger values of $\beta$. Some of the energy can now leave the region of the null point and there is again generation of a slow wave component (we find that the fraction of the incident wave converted to a slow wave is proportional to $\beta$). We conclude that there are two competing phenomena; the refraction effect (due to the variable Alfv\'en speed) and the contribution from the non-zero sound speed. These experiments illustrate the importance of the magnetic topology and of the location of the $\beta \approx 1$ layer in the system.

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.

The nature of fast magnetoacoustic and Alfv\'en waves is investigated in a zero $\beta$ plasma in the neighbourhood of a pair of two-dimensional null points. This gives an indication of wave propagation in the low $\beta$ solar corona, for a more complicated magnetic configuration than that looked at by McLaughlin & Hood (2004). It is found that the fast wave is attracted to the null points and that the front of the wave slows down as it approaches the null point pair, with the wave splitting and part of the wave accumulating at one null and the rest at the other. Current density will then accumulate at these points and ohmic dissipation will then extract the energy in the wave at these points. This suggests locations where wave heating will occur in the corona. The Alfv\'en wave behaves in a different manner in that the wave accumulates along the separatrices. Hence, the current density will accumulate at this part of the topology and this is where wave heating will occur. However, the phenomenon of wave accumulation at a specific place is a feature of both wave types, and illustrates the importance of studying the topology of the corona when considering MHD wave propagation.

The nature of fast magnetoacoustic and Alfv\'en waves is investigated in a zero $\beta$ plasma. This gives an indication of wave propagation in the low $\beta$ solar corona. It is found that for a two-dimensional null point, the fast wave is attracted to that point and the front of the wave slows down as it approaches the null point, causing the current density to accumulate there and rise rapidly. Ohmic dissipation will extract the energy in the wave at this point. This illustrates that null points play an important role in the rapid dissipation of fast magnetoacoustic waves and suggests the location where wave heating will occur in the corona. The Alfv\'en wave behaves in a different manner in that the wave energy is dissipated along the separatrices. For Alfv\'en waves that are decoupled from fast waves, the value of the plasma $\beta$ is unimportant. However, the phenomenon of dissipating the majority of the wave energy at a specific place is a feature of both wave types.

This paper is the third in a series of investigations by the authors. The nature of fast magnetoacoustic and Alfv\'en waves is investigated in a 2D $\beta=0$ plasma in the neighbourhood of two dipoles. We use both numerical simulations (two-step Lax-Wendroff scheme) and analytical techniques (WKB approximation). It is found that the propagation of the linear fast wave is dictated by the Alfv\'en speed profile and that close to the null, the wave is attracted to the neutral point. However, it is also found that in this magnetic configuration some of the wave can escape the refraction effect; this had not been seen in previous investigations by the authors. The wave split occurs near the regions of very high Alfv\'en speed (found near the loci of the two dipoles). Also, for the set-up investigated it was found that 40% of the wave energy accumulates at the null. Ohmic dissipation will then extract the wave energy at this point. The Alfv\'en wave behaves in a different manner in that part of the wave accumulates along the separatrices and part escapes. Hence, the current density will accumulate at this part of the topology and this is where wave heating will occur. The phenomenon of wave accumulation at a specific place is a feature of both wave types, as is the result that a fraction of the wave can now escape the numerical box when propagating in this magnetic configuration.