The SunPy Project developed a 13-question survey to understand the software and hardware usage of the solar physics community. 364 members of the solar physics community, across 35 countries, responded to our survey. We found that 99$\pm$0.5% of respondents use software in their research and 66% use the Python scientific software stack. Students are twice as likely as faculty, staff scientists, and researchers to use Python rather than Interactive Data Language (IDL). In this respect, the astrophysics and solar physics communities differ widely: 78% of solar physics faculty, staff scientists, and researchers in our sample uses IDL, compared with 44% of astrophysics faculty and scientists sampled by Momcheva and Tollerud (2015). 63$\pm$4% of respondents have not taken any computer-science courses at an undergraduate or graduate level. We also found that most respondents utilize consumer hardware to run software for solar-physics research. Although 82% of respondents work with data from space-based or ground-based missions, some of which (e.g. the Solar Dynamics Observatory and Daniel K. Inouye Solar Telescope) produce terabytes of data a day, 14% use a regional or national cluster, 5% use a commercial cloud provider, and 29% use exclusively a laptop or desktop. Finally, we found that 73$\pm$4% of respondents cite scientific software in their research, although only 42$\pm$3% do so routinely.

Extreme ultraviolet (EUV) waves are large-scale propagating disturbances observed in the solar corona, frequently associated with coronal mass ejections and flares. They appear as faint, extended structures propagating from a source region across the structured solar corona. Since their discovery, over two hundred papers discussing their properties, causes and physical nature have been published. However, despite this their fundamental properties and the physics of their interactions with other solar phenomena are still not understood. To further the understanding of EUV waves, we have constructed the Automated Wave Analysis and REduction (AWARE) algorithm for the measurement of EUV waves. AWARE is implemented in two stages. In the first stage, we use a new type of running difference image, the running difference persistence image, which enables the efficient isolation of propagating, brightening wavefronts as they propagate across the corona. In the second stage, AWARE detects the presence of a wavefront, and measures the distance, velocity and acceleration of that wavefront across the Sun. The fit of propagation models to the wave progress isolated in the first stage is achieved using the Random Sample and Consensus (RANSAC) algorithm. AWARE is tested against simulations of EUV wave propagation, and is applied to measure EUV waves in observational data from the Atmospheric Imaging Assembly (AIA). We also comment on unavoidable systematic errors that bias the estimation of wavefront velocity and acceleration. In addition, the full AWARE software suite comes with a package that creates simulations of waves propagating across the disk from arbitrary starting points.

Torsional motions are ubiquitous in the solar atmosphere. In this work, we perform 3D numerical simulations which mimic a vortex-type photospheric driver with a Gaussian spatial profile. This driver is implemented to excite MHD waves in an axially symmetric, 3D magnetic flux tube embedded in a realistic solar atmosphere. The Gaussian width of the driver is varied and the resulting perturbations are compared. Velocity vectors were decomposed into parallel, perpendicular and azimuthal components with respect to pre-defined magnetic flux surfaces. These components correspond broadly to the fast, slow and Alfven modes, respectively. From these velocities the corresponding wave energy fluxes are calculated, allowing us to estimate the contribution of each mode to the energy flux. For the narrowest driver ($0.15$ Mm) the parallel component accounts for $\sim 55 - 65\%$ of the flux. This contribution increases smoothly with driver width up to nearly $90\%$ for the widest driver ($0.35$ Mm). The relative importance of the perpendicular and azimuthal components decrease at similar rates. The azimuthal energy flux varied between $\sim 35\%$ for the narrowest driver and $< 10\%$ for the widest one. Similarly, the perpendicular flux was $\sim 25 - 10\%$. We also demonstrate that the fast mode corresponds to the sausage wave in our simulations. Our results therefore show that the fast sausage wave is easily excited by this driver and that it carries the majority of the energy transported. For this vortex-type driver the Alfven wave does not contribute a significant amount of energy.

The Astropy project supports and fosters the development of open-source and openly-developed Python packages that provide commonly-needed functionality to the astronomical community. A key element of the Astropy project is the core package Astropy, which serves as the foundation for more specialized projects and packages. In this article, we provide an overview of the organization of the Astropy project and summarize key features in the core package as of the recent major release, version 2.0. We then describe the project infrastructure designed to facilitate and support development for a broader ecosystem of inter-operable packages. We conclude with a future outlook of planned new features and directions for the broader Astropy project.

The QCD Axion is a particle postulated to exist since the 1970s to explain the Strong-CP problem in particle physics. It could also account for all of the observed Dark Matter in the universe. The Axion Resonant InterAction DetectioN Experiment (ARIADNE) experiment intends to detect the QCD axion by sensing the fictitious "magnetic field" created by its coupling to spin. The experiment must be sensitive to magnetic fields below the $10^{-19}$ T level to achieve its design sensitivity, necessitating tight control of the experiment's magnetic environment. We describe a method for controlling three aspects of that environment which would otherwise limit the experimental sensitivity. Firstly, a system of superconducting magnetic shielding is described to screen ordinary magnetic noise from the sample volume at the $10^8$ level. Secondly, a method for reducing magnetic field gradients within the sample up to $10^2$ times is described, using a simple and cost-effective design geometry. Thirdly, a novel coil design is introduced which allows the generation of fields similar to those produced by Helmholtz coils in regions directly abutting superconducting boundaries. The methods may be generally useful for magnetic field control near superconducting boundaries in other experiments where similar considerations apply.

The Axion Resonant InterAction Detection Experiment (ARIADNE) is a collaborative effort to search for the QCD axion using techniques based on nuclear magnetic resonance. In the experiment, axions or axion-like particles would mediate short-range spin-dependent interactions between a laser-polarized 3He gas and a rotating (unpolarized) tungsten source mass, acting as a tiny, fictitious "magnetic field". The experiment has the potential to probe deep within the theoretically interesting regime for the QCD axion in the mass range of 0.1-10 meV, independently of cosmological assumptions. The experiment relies on a stable rotary mechanism and superconducting magnetic shielding, required to screen the 3He sample from ordinary magnetic noise. Progress on testing the stability of the rotary mechanism is reported, and the design for the superconducting shielding is discussed.

This paper presents SunPy (version 0.5), a community-developed Python package for solar physics. Python, a free, cross-platform, general-purpose, high-level programming language, has seen widespread adoption among the scientific community, resulting in the availability of a large number of software packages, from numerical computation (NumPy, SciPy) and machine learning (scikit-learn) to visualisation and plotting (matplotlib). SunPy is a data-analysis environment specialising in providing the software necessary to analyse solar and heliospheric data in Python. SunPy is open-source software (BSD licence) and has an open and transparent development workflow that anyone can contribute to. SunPy provides access to solar data through integration with the Virtual Solar Observatory (VSO), the Heliophysics Event Knowledgebase (HEK), and the HELiophysics Integrated Observatory (HELIO) webservices. It currently supports image data from major solar missions (e.g., SDO, SOHO, STEREO, and IRIS), time-series data from missions such as GOES, SDO/EVE, and PROBA2/LYRA, and radio spectra from e-Callisto and STEREO/SWAVES. We describe SunPy's functionality, provide examples of solar data analysis in SunPy, and show how Python-based solar data-analysis can leverage the many existing tools already available in Python. We discuss the future goals of the project and encourage interested users to become involved in the planning and development of SunPy.

High-resolution observations of the solar photosphere have identified a wide variety of spiralling motions in the solar plasma. These spirals vary in properties, but are observed to be abundant at the solar surface. In this work these spirals are studied for their potential as magnetohydrodynamic (MHD) wave generation mechanisms. The inter-granular lanes, where these spirals are commonly observed, are also regions where the magnetic field strength is higher than average. This combination of magnetic field and spiralling plasma is a recipe for the generation of Alfv\'en waves and other MHD waves. This work employs numerical simulations of a self-similar magnetic flux tube embedded in a realistic, gravitationally stratified, solar atmosphere to study the effects of a single magnetic flux tube perturbed by a logarithmic velocity spiral driver. The expansion factor of the logarithmic spiral driver is varied and multiple simulations are run for a range of values of the expansion factor centred around observational constraints. The simulations are analysed using `flux surfaces' constructed from the magnetic field lines so that the vectors perpendicular, parallel and azimuthal to the local magnetic field vector can be calculated. The results of this analysis show that the Alfv\'en wave is the dominant wave for lower values of the expansion factor, whereas, for the higher values the parallel component is dominant. This transition occurs within the range of the observational constraints, meaning that spiral drivers, as observed in the solar photosphere, have the potential to generate a variety of MHD wave modes.

Recent ground- and space-based observations reveal the presence of small-scale motions between convection cells in the solar photosphere. In these regions small-scale magnetic flux tubes are generated due to the interaction of granulation motion and background magnetic field. This paper studies the effects of these motions, on magnetohydrodynamic wave excitation from broadband photospheric drivers. Numerical experiments of linear magnetohydrodynamic wave propagation in a magnetic flux tube embedded in a realistic gravitationally stratified solar atmosphere between the photosphere and the low choromosphere (above $\beta=1$) are performed. Horizontal and vertical velocity field drivers mimic granular buffeting and solar global oscillations. A uniform torsional driver as well as Archimedean and logarithmic spiral drivers mimic observed torsional motions in the solar photosphere. The results are analysed using a novel method for extracting the parallel, perpendicular and azimuthal components of the perturbations, which caters for both the linear and non-linear case. Employing this method yields the identification of the wave modes excited in the numerical simulations and enables a comparison of excited modes via velocity perturbations and wave energy flux. The wave energy flux distribution is calculated, to enable the quantification of the relative strengths of excited modes. The torsional drivers excite primarily Alfv\'en modes (~60% of the total flux) with small contributions from the slow kink mode, and, for the logarithmic spiral driver, small amounts of slow sausage mode. The horizontal and vertical drivers excite primarily slow kink or fast sausage modes respectively, with small variations dependent upon flux surface radius.

There have been ubiquitous observations of wave-like motions in the solar atmosphere for decades. Recent improvements to space- and ground-based observatories have allowed the focus to shift to smaller magnetic structures on the solar surface. In this paper, high-resolution ground-based data taken using the Swedish 1 m Solar Telescope is combined with co-spatial and co-temporal data from the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory (SDO) satellite to analyze running penumbral waves (RPWs). RPWs have always been thought to be radial wave propagation that occurs within sunspots. Recent research has suggested that they are in fact upwardly propagating field-aligned waves (UPWs). Here, RPWs within a solar pore are observed for the first time and are interpreted as UPWs due to the lack of a penumbra that is required to support RPWs. These UPWs are also observed co-spatially and co-temporally within several SDO/AIA elemental lines that sample the transition region and low corona. The observed UPWs are traveling at a horizontal velocity of around 17 +- 0.5 km s-1 and a minimum vertical velocity of 42 +- 21 km s-1. The estimated energy of the waves is around 150 W m-2, which is on the lower bound required to heat the quiet-Sun corona. This is a new, yet unconsidered source of wave energy within the solar chromosphere and low corona.

A single open magnetic flux tube spanning the solar photosphere (solar radius R) and the lower corona (R + 10 Mm) is modelled in magnetohydrostatic equilibrium within a realistic stratified atmosphere subject to solar gravity. Such flux tubes are observed to remain relatively stable for up to a day or more, and it is our aim to apply the model as the background condition for numerical studies of energy transport mechanisms from the surface to the corona. We solve analytically an axially symmetric 3D structure for the model, with magnetic field strength, plasma density, pressure and temperature all consistent with observational and theoretical estimates. The self similar construction ensures the magnetic field is divergence free. The equation of pressure balance for this particular set of flux tubes can be integrated analytically to find the pressure and density corrections required to preserve the magnetohydrostatic equilibrium. The model includes a number of free parameters, which makes the solution applicable to a variety of other physical problems and it may therefore be of more general interest.

The properties of Ellerman bombs (EBs), small-scale brightenings in the H-alpha line wings, have proved difficult to establish due to their size being close to the spatial resolution of even the most advanced telescopes. Here, we aim to infer the size and lifetime of EBs using high-resolution data of an emerging active region collected using the Interferometric BIdimensional Spectrometer (IBIS) and Rapid Oscillations of the Solar Atmosphere (ROSA) instruments as well as the Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO). We develop an algorithm to track EBs through their evolution, finding that EBs can often be much smaller (around 0.3") and shorter lived (less than 1 minute) than previous estimates. A correlation between G-band magnetic bright points and EBs is also found. Combining SDO/HMI and G-band data gives a good proxy of the polarity for the vertical magnetic field. It is found that EBs often occur both over regions of opposite polarity flux and strong unipolar fields, possibly hinting at magnetic reconnection as a driver of these events.The energetics of EB events is found to follow a power-law distribution in the range of "nano-flare" (10^{22-25} ergs).