Showing 1–39 of 39 results for author: Pontin, D I

Exploring the generation and annihilation of three dimensional nulls through MHD simulations in initially chaotic magnetic field devoid of nulls

Authors: Yogesh Kumar Maurya, Ramit Bhattacharyya, David I. Pontin, Sanjay Kumar

Abstract: Three-dimensional (3D) magnetic nulls are abundant in the solar atmosphere, as been firmly established through contemporary observations. They are established to be important magnetic structures in, for example, jets and circular ribbon flares. While simulations and extrapolations support this, the mechanisms behind 3D null generation remain an open question. Recent magnetohydrodynamics (MHD) simu… ▽ More Three-dimensional (3D) magnetic nulls are abundant in the solar atmosphere, as been firmly established through contemporary observations. They are established to be important magnetic structures in, for example, jets and circular ribbon flares. While simulations and extrapolations support this, the mechanisms behind 3D null generation remain an open question. Recent magnetohydrodynamics (MHD) simulations propose that magnetic reconnection is responsible for both generating and annihilating 3D nulls, a novel concept. However, these simulations began with initial magnetic fields already supporting pre-existing nulls, raising the question of whether magnetic reconnection can create nulls in fields initially devoid of them. Previously, this question was briefly explored in a simulation with an initial chaotic magnetic field. However, the study failed to precisely identify locations, topological degrees, and natures (spiral or radial) of nulls, and it approximated magnetic reconnection without fully tracking field line in time. In this paper these findings are revisited in light of recent advancements and tools used to locate and trace nulls, along with the tracing of field lines, through which the concept of generation/annihilation of 3D nulls from chaotic fields is established in a precise manner. △ Less

Submitted 10 June, 2024; originally announced June 2024.

Generation and annihilation of three dimensional magnetic nulls in extrapolated solar coronal magnetic field: Data-based Implicit Large Eddy Simulation

Authors: Yogesh Kumar Maurya, Ramit Bhattacharyya, David I. Pontin

Abstract: Three-dimensional magnetic nulls are the points where magnetic field vanishes and are preferential sites for magnetic reconnection: a process which converts magnetic energy into heat and accelerates charged particles along with a rearrangement of magnetic field lines. In the solar corona, the reconnections manifest as coronal transients including solar flares, coronal mass ejections and coronal je… ▽ More Three-dimensional magnetic nulls are the points where magnetic field vanishes and are preferential sites for magnetic reconnection: a process which converts magnetic energy into heat and accelerates charged particles along with a rearrangement of magnetic field lines. In the solar corona, the reconnections manifest as coronal transients including solar flares, coronal mass ejections and coronal jets. The nulls are generally found to be collocated with complex active regions on the solar photosphere. Extrapolation of magnetic field from corresponding photospheric magnetogram indicate an abundance of these nulls in solar atmosphere. Nevertheless, their generation is still not well understood. Recently, Maurya et al. (2023) have demonstrated magnetic reconnection to be a cause for generation and annihilation of magnetic nulls through magnetohydrodynamics simulation, where the initial magnetic field is idealized to have a single radial null. This article further extends the study in a more realistic scenario where the initial magnetic field is constructed by extrapolating photospheric magnetogram data and hence, incorporates field line complexities inherent to a complex active region. For the purpose, the active region NOAA 11977 hosting a C6.6 class flare is selected. The simulation is initiated using non-force-free extrapolated magnetic field from the photospheric vector magnetogram at around 02:48:00 UT on 17 February 2014, 16 minutes before the flare peak. The generation, annihilation and dynamics of nulls are explored by a complimentary usage of trilinear null detection technique and tracing of magnetic field line dynamics. It is found that the nulls can spontaneously generate/annihilate in pairs while preserving the topological degree and can have observational implications like footpoint brightenings. Magnetic reconnection is found to be the cause of such generation and annihilation. △ Less

Submitted 18 April, 2024; originally announced April 2024.

Coronal Heating and Solar Wind Generation by Flux Cancellation Reconnection

Authors: D. I. Pontin, E. R. Priest, L. P. Chitta, V. S. Titov

Abstract: In this paper we propose that flux cancellation on small granular scales ($\lesssim 1000$ km) ubiquitously drives reconnection at a multitude of sites in the low solar atmosphere, contributing to chromospheric/coronal heating and the generation of the solar wind. We analyse the energy conversion in these small-scale flux cancellation events using both analytical models and three-dimensional, resis… ▽ More In this paper we propose that flux cancellation on small granular scales ($\lesssim 1000$ km) ubiquitously drives reconnection at a multitude of sites in the low solar atmosphere, contributing to chromospheric/coronal heating and the generation of the solar wind. We analyse the energy conversion in these small-scale flux cancellation events using both analytical models and three-dimensional, resistive magnetohydrodynamic (MHD) simulations. The analytical models -- in combination with the latest estimates of flux cancellation rates -- allow us to estimate the energy release rates due to cancellation events, which are found to be of the order $10^6-10^7$ erg cm$^{-2}$sec$^{-1}$, sufficient to heat the chromosphere and corona of the Quiet Sun and active regions, and to power the solar wind. The MHD simulations confirm the conversion of energy in reconnecting current sheets, in a geometry representing a small-scale bipole being advected towards an intergranular lane. A ribbon-like jet of heated plasma is accelerated upwards, that could also escape the Sun as the solar wind in an open-field configuration. We conclude that through two phases of atmospheric energy release -- pre-cancellation and cancellation -- the cancellation of photospheric magnetic flux fragments and the associated magnetic reconnection may provide a substantial energy and mass flux contribution to coronal heating and solar wind generation. △ Less

Submitted 21 January, 2024; originally announced January 2024.

Comments: 19 pages, 11 figures, published in ApJ

Journal ref: The Astrophysical Journal, 960:51 (19pp), 2024 January 1

2.5-D MHD Simulation of the Formation and Evolution of Plasmoids in Coronal Current Sheets

Authors: Sripan Mondal, Abhishek K Srivastava, David I. Pontin, Ding Yuan, Eric R. Priest

Abstract: In the present paper, using MPI-AMRVAC, we perform a 2.5-D numerical MHD simulation of the dynamics and associated thermodynamical evolution of an initially force-free Harris current sheet subjected to an external velocity perturbation under the condition of uniform resistivity. The amplitude of the magnetic field is taken to be 10 Gauss, typical of the solar corona. We impose a Gaussian velocity… ▽ More In the present paper, using MPI-AMRVAC, we perform a 2.5-D numerical MHD simulation of the dynamics and associated thermodynamical evolution of an initially force-free Harris current sheet subjected to an external velocity perturbation under the condition of uniform resistivity. The amplitude of the magnetic field is taken to be 10 Gauss, typical of the solar corona. We impose a Gaussian velocity pulse across this current sheet mimicking the interaction of fast magnetoacoustic waves with a current sheet in corona. This leads to a variety of dynamics and plasma processes in the current sheet, which is initially quasi-static. The initial pulse interacts with the current sheet and splits into a pair of counter-propagating wavefronts, which forms a rarefied region and leads to inflow and a thinning of the current sheet. The thinning results in Petschek-type magnetic reconnection followed by tearing instability and plasmoid formation. The reconnection outflows containing outward-moving plasmoids have accelerated motions with velocities ranging from 105-303 km/s. The average temperature and density of the plasmoids are found to be 8 MK and twice the background density of the solar corona, respectively. These estimates of velocity, temperature and density of plasmoids are similar to values reported from various solar coronal observations. Therefore, we infer that the external triggering of a quasi-static current sheet by a single velocity pulse is capable of initiating magnetic reconnection and plasmoid formation in the absence of a localized enhancement of resistivity in the solar corona. △ Less

Submitted 13 January, 2024; originally announced January 2024.

Comments: 20 pages, 10 figures, Accepted for publication in The Astrophysical Journal

Reconnection generated plasma flows in the quasi-separatrix layer in localised solar corona

Authors: Sripan Mondal, A. K. Srivastava, Sudheer K. Mishra, K. Sangal, Pradeep Kayshap, Yang Guo, David I. Pontin, Vadim M. Uritsky, Leon Ofman, T. -J. Wang, Ding Yuan

Abstract: Multiwavelength observations of the propagating disturbances (PDs), discovered by Atmospheric Imaging Assembly (AIA) onboard Solar Dynamics Observatory (SDO), are analyzed to determine its driving mechanism and physical nature. Two magnetic strands in the localised corona are observed to approach and merge with each other followed by the generation of brightening, which further propagates in a cus… ▽ More Multiwavelength observations of the propagating disturbances (PDs), discovered by Atmospheric Imaging Assembly (AIA) onboard Solar Dynamics Observatory (SDO), are analyzed to determine its driving mechanism and physical nature. Two magnetic strands in the localised corona are observed to approach and merge with each other followed by the generation of brightening, which further propagates in a cusp-shaped magnetic channel. Differential emission measure analysis shows an occurrence of heating in this region-of-interest (ROI). We extrapolate potential magnetic field lines at coronal heights from observed Helioseismic and Magnetic Imager (HMI) vector magnetogram via Green's function method using MPI-AMRVAC. We analyze the field to locate magnetic nulls and quasi-separatrix layers (QSLs) which are preferential locations for magnetic reconnection. Dominant QSLs including a magnetic null are found to exist and match the geometry followed by PDs, therefore, it provides conclusive evidence of magnetic reconnection. In addition, spectroscopic analysis of Interface Region Imaging Spectrograph (IRIS) Si IV 1393.77 Å line profiles show a rise of line-width in the same time range depicting presence of mass motion in the observed cusp-shaped region. PDs are observed to exhibit periodicities of around four minutes. The speeds of PDs measured by Surfing Transform Technique are almost close to each other in four different SDO/AIA bandpasses, i.e., 304, 171, 193 and 131 Å excluding the interpretation of PDs in terms of slow magnetoacoustic waves. We describe comprehensively the observed PDs as quasi-periodic plasma flows generated due to periodic reconnection in vicinity of a coronal magnetic null. △ Less

Submitted 3 May, 2023; originally announced May 2023.

Comments: 21 Pages; 10 Figures; Accepted for the Publication in The Astrophysical Journal

Parallel plasma loops and the energization of the solar corona

Authors: Hardi Peter, Lakshmi Pradeep Chitta, Feng Chen, David I. Pontin, Amy R. Winebarger, Leon Golub, Sabrina L. Savage, Laurel A. Rachmeler, Ken Kobayashi, David H. Brooks, Jonathan W. Cirtain, Bart De Pontieu, David E. McKenzie, Richard J. Morton, Paola Testa, Sanjiv K. Tiwari, Robert W. Walsh, Harry P. Warren

Abstract: The outer atmosphere of the Sun is composed of plasma heated to temperatures well in excess of the visible surface. We investigate short cool and warm (<1 MK) loops seen in the core of an active region to address the role of field-line braiding in energising these structures. We report observations from the High-resolution Coronal imager (Hi-C) that have been acquired in a coordinated campaign wit… ▽ More The outer atmosphere of the Sun is composed of plasma heated to temperatures well in excess of the visible surface. We investigate short cool and warm (<1 MK) loops seen in the core of an active region to address the role of field-line braiding in energising these structures. We report observations from the High-resolution Coronal imager (Hi-C) that have been acquired in a coordinated campaign with the Interface Region Imaging Spectrograph (IRIS). In the core of the active region, the 172 A band of Hi-C and the 1400 A channel of IRIS show plasma loops at different temperatures that run in parallel. There is a small but detectable spatial offset of less than 1 arcsec between the loops seen in the two bands. Most importantly, we do not see observational signatures that these loops might be twisted around each other. Considering the scenario of magnetic braiding, our observations of parallel loops imply that the stresses put into the magnetic field have to relax while the braiding is applied: the magnetic field never reaches a highly braided state on these length-scales comparable to the separation of the loops. This supports recent numerical 3D models of loop braiding in which the effective dissipation is sufficiently large that it keeps the magnetic field from getting highly twisted within a loop. △ Less

Submitted 31 May, 2022; originally announced May 2022.

Comments: Accepted for publication in the Astrophysical Journal, 24 pages, 18 figures

The Dynamic Coupling of Streamers and Pseudostreamers to the Heliosphere

Authors: V. Aslanyan, D. I. Pontin, A. K. Higginson, P. F. Wyper, R. B. Scott, S. K. Antiochos

Abstract: The slow solar wind is generally believed to result from the interaction of open and closed coronal magnetic flux at streamers and pseudostreamers. We use 3-dimensional magnetohydrodynamic simulations to determine the detailed structure and dynamics of open-closed interactions that are driven by photospheric convective flows. The photospheric magnetic field model includes a global dipole giving ri… ▽ More The slow solar wind is generally believed to result from the interaction of open and closed coronal magnetic flux at streamers and pseudostreamers. We use 3-dimensional magnetohydrodynamic simulations to determine the detailed structure and dynamics of open-closed interactions that are driven by photospheric convective flows. The photospheric magnetic field model includes a global dipole giving rise to a streamer together with a large parasitic polarity region giving rise to a pseudostreamer that separates a satellite coronal hole from the main polar hole. Our numerical domain extends out to 30 solar radii and includes an isothermal solar wind, so that the coupling between the corona and heliosphere can be calculated rigorously. This system is driven by imposing a large set of quasi-random surface flows that capture the driving of coronal flux in the vicinity of streamer and pseudostreamer boundaries by the supergranular motions. We describe the resulting structures and dynamics. Interchange reconnection dominates the evolution at both streamer and pseudostreamer boundaries, but the details of the resulting structures are clearly different from one another. Additionally, we calculate in situ signatures of the reconnection and determine the dynamic mapping from the inner heliosphere back to the Sun for a test spacecraft orbit. We discuss the implications of our results for interpreting observations from inner heliospheric missions, such as Parker Solar Probe and Solar Orbiter, and for space weather modeling of the slow solar wind. △ Less

Submitted 27 April, 2022; v1 submitted 7 January, 2022; originally announced January 2022.

Comments: This article was submitted under duress from an unreasonable reviewer for a separate article

Is flare-ribbon fine structure related to tearing in the flare current sheet?

Authors: P. F. Wyper, D. I. Pontin

Abstract: Observations of solar flare ribbons show significant fine structure in the form of breaking wave-like perturbations and spirals. The origin of this structure is not well understood, but one possibility is that it is related to the tearing instability in the flare current sheet. Here we study this connection by constructing an analytical three-dimensional magnetic field representative of an eruptin… ▽ More Observations of solar flare ribbons show significant fine structure in the form of breaking wave-like perturbations and spirals. The origin of this structure is not well understood, but one possibility is that it is related to the tearing instability in the flare current sheet. Here we study this connection by constructing an analytical three-dimensional magnetic field representative of an erupting flux rope with a flare current sheet below it. We introduce small-scale flux ropes representative of those formed during a tearing instability in the current layer, and use the squashing factor on the solar surface to identify the shape of the presumed flare ribbons and fine structure. Our analysis suggests there is a direct link between flare-ribbon fine structure and flare current sheet tearing, with the majority of the ribbon fine structure related to oblique tearing modes. Depending upon the size, location and twist of the small-scale flux ropes, breaking wave-like and spiral features within the hooks and straight sections of the flare ribbon can be formed that are qualitatively similar to observations. We also show that the handedness of the spirals/waves must be the same as the handedness of the hooks of the main ribbon. We conclude that tearing in the flare current layer is a likely explanation for spirals and wave-like features in flare ribbons. △ Less

Submitted 24 August, 2021; originally announced August 2021.

Comments: ApJ, in press, 16 pages, 14 figures

Non-thermal line broadening due to braiding-induced turbulence in solar coronal loops

Authors: D. I. Pontin, H. Peter, L. P. Chitta

Abstract: Aims: Emission line profiles from solar coronal loops exhibit properties that are unexplained by current models. We investigate the non-thermal broadening associated with plasma heating in coronal loops that is induced by magnetic field line braiding. Methods: We describe the coronal loop by a 3D magnetohydrodynamic model of the turbulent decay of an initially-braided magnetic field. From this we… ▽ More Aims: Emission line profiles from solar coronal loops exhibit properties that are unexplained by current models. We investigate the non-thermal broadening associated with plasma heating in coronal loops that is induced by magnetic field line braiding. Methods: We describe the coronal loop by a 3D magnetohydrodynamic model of the turbulent decay of an initially-braided magnetic field. From this we synthesize the Fe XII line at 193 Angstroms that forms around 1.5 MK. Results: Key features of current observations of extreme UV lines from the corona are reproduced in the synthesised spectra: (i) Typical non-thermal widths range from 15 to 20 km/s. (ii) The widths are approximately independent of the size of the field of view. (iii) There is a correlation between line intensity and the non-thermal broadening. (iv) Spectra are found to be non-Gaussian, with enhanced power in the wings of order 10-20%. Conclusions: Our model provides an explanation that self-consistently connects the heating process to the observed non-thermal line broadening. The non-Gaussian nature of the spectra is a consequence of the non-Gaussian nature of the underlying velocity fluctuations, interpreted as a signature of intermittency in the turbulence. △ Less

Submitted 27 August, 2020; originally announced August 2020.

Comments: Published in A&A

Journal ref: A&A 639, A21 (2020)

On the periodicity of linear and nonlinear oscillatory reconnection

Authors: J. O. Thurgood, D. I. Pontin, J. A. McLaughlin

Abstract: (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 be… ▽ More (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. △ Less

Submitted 21 November, 2018; originally announced November 2018.

Comments: Accepted in A&A, 5 Figures (1 animated)

Journal ref: A&A 621, A106 (2019)

Resistively-limited current sheet implosions in planar anti-parallel (1D) and null-point containing (2D) magnetic field geometries

Authors: J. O. Thurgood, D. I. Pontin, J. A. McLaughlin

Abstract: 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,… ▽ More 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. △ Less

Submitted 11 July, 2018; v1 submitted 21 June, 2018; originally announced June 2018.

Comments: Published open access in Physics of Plasmas - go to main journal for animations https://doi.org/10.1063/1.5035489

Journal ref: Physics of Plasmas, 25, 072105 (2018)

Magnetic Structures at the Boundary of the Closed Corona: Interpretation of S-web Arcs

Authors: Roger B. Scott, David I. Pontin, Anthony R. Yeates, Peter F. Wyper

Abstract: The topology of magnetic fields near the open-closed flux boundary in the Sun's corona is an important influencing factor in the process of interchange reconnection, whereby plasma is exchanged between open and closed flux domains. Maps of the magnetic squashing factor at the radial outer boundary in coronal field models reveal the presence of the so-called `S-web', and suggest that interchange re… ▽ More The topology of magnetic fields near the open-closed flux boundary in the Sun's corona is an important influencing factor in the process of interchange reconnection, whereby plasma is exchanged between open and closed flux domains. Maps of the magnetic squashing factor at the radial outer boundary in coronal field models reveal the presence of the so-called `S-web', and suggest that interchange reconnection could potentially deposit closed coronal material into high-latitude regions far from the heliospheric current sheet. Here we demonstrate that certain features of the S-web reveal the underlying topological structure of the magnetic field. Specifically, in order for the arcing bands of highly squashed magnetic flux of the S-web to terminate or intersect away from the helmet streamer apex, there must be a null spine line that maps a finite segment of the photospheric open-closed boundary up to a singular point in the open flux domain. We propose that this association between null spine lines and arc termination points may be used to identify locations in the heliosphere that are preferential for the appearance of solar energetic particles or slow solar wind plasma with certain characteristics. △ Less

Submitted 11 May, 2018; originally announced May 2018.

Implosive collapse about magnetic null points: A quantitative comparison between 2D and 3D nulls

Authors: J. O. Thurgood, D. I. Pontin, J. A. McLaughlin

Abstract: 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.… ▽ More 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. △ Less

Submitted 20 February, 2018; originally announced February 2018.

Comments: Accepted in ApJ, will be published gold open access, refer to main journal

Three-dimensional oscillatory magnetic reconnection

Authors: J. O. Thurgood, D. I. Pontin, J. A. McLaughlin

Abstract: 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… ▽ More 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 △ Less

Submitted 29 June, 2017; originally announced June 2017.

Comments: Accepted in ApJ, will be published gold open access, refer to main journal for animated content

Effects of Fieldline Topology on Energy Propagation in the Corona

Authors: Simon Candelaresi, David Ian Pontin, Gunnar Hornig

Abstract: We study the effect of photospheric footpoint motions on magnetic field structures containing magnetic nulls. The footpoint motions are prescribed on the photospheric boundary as a velocity field which entangles the magnetic field. We investigate the propagation of the injected energy, the conversion of energy, emergence of current layers and other consequences of the non-trivial magnetic field to… ▽ More We study the effect of photospheric footpoint motions on magnetic field structures containing magnetic nulls. The footpoint motions are prescribed on the photospheric boundary as a velocity field which entangles the magnetic field. We investigate the propagation of the injected energy, the conversion of energy, emergence of current layers and other consequences of the non-trivial magnetic field topology in this situation. These boundary motions lead initially to an increase in magnetic and kinetic energy. Following this, the energy input from the photosphere is partially dissipated and partially transported out of the domain through the Poynting flux. The presence of separatrix layers and magnetic null-points fundamentally alters the propagation behavior of disturbances from the photosphere into the corona. Depending on the field line topology close to the photosphere, the energy is either trapped or free to propagate into the corona. △ Less

Submitted 18 October, 2017; v1 submitted 10 November, 2016; originally announced November 2016.

Comments: 14 pages, 15 figures

Journal ref: (2016) ApJ, 832, 150

Why are flare ribbons associated with the spines of magnetic null points generically elongated?

Authors: D. I. Pontin, K. Galsgaard, P. Demoulin

Abstract: Coronal magnetic null points exist in abundance as demonstrated by extrapolations of the coronal field, and have been inferred to be important for a broad range of energetic events. These null points and their associated separatrix and spine field lines represent discontinuities of the field line mapping, making them preferential locations for reconnection. This field line mapping also exhibits st… ▽ More Coronal magnetic null points exist in abundance as demonstrated by extrapolations of the coronal field, and have been inferred to be important for a broad range of energetic events. These null points and their associated separatrix and spine field lines represent discontinuities of the field line mapping, making them preferential locations for reconnection. This field line mapping also exhibits strong gradients adjacent to the separatrix (fan) and spine field lines, that can be analysed using the `squashing factor', $Q$. In this paper we make a detailed analysis of the distribution of $Q$ in the presence of magnetic nulls. While $Q$ is formally infinite on both the spine and fan of the null, the decay of $Q$ away from these structures is shown in general to depend strongly on the null-point structure. For the generic case of a non-radially-symmetric null, $Q$ decays most slowly away from the spine/fan in the direction in which $|{\bf B}|$ increases most slowly. In particular, this demonstrates that the extended, elliptical high-$Q$ halo around the spine footpoints observed by Masson et al. (Astrophys. J., 700, 559, 2009) is a generic feature. This extension of the $Q$ halos around the spine/fan footpoints is important for diagnosing the regions of the photosphere that are magnetically connected to any current layer that forms at the null. In light of this, we discuss how our results can be used to interpret the geometry of observed flare ribbons in `circular ribbon flares', in which typically a coronal null is implicated. We conclude that both the physics in the vicinity of the null and how this is related to the extension of $Q$ away from the spine/fan can be used in tandem to understand observational signatures of reconnection at coronal null points. △ Less

Submitted 18 May, 2016; originally announced May 2016.

Comments: Pre-print version of article accepted for publication in Solar Physics

Braided magnetic fields: equilibria, relaxation and heating

Authors: D. I. Pontin, S. Candelaresi, A. J. B. Russell, G. Hornig

Abstract: 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, the… ▽ More 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. △ Less

Submitted 13 April, 2016; v1 submitted 18 December, 2015; originally announced December 2015.

Comments: To appear in Plasma Physics and Controlled Fusion

The effect of reconnection on the structure of the Sun's open-closed-flux boundary

Authors: D. I. Pontin, P. F. Wyper

Abstract: Global magnetic field extrapolations are now revealing the huge complexity of the Sun's corona, and in particular the structure of the boundary between open and closed magnetic flux. Moreover, recent developments indicate that magnetic reconnection in the corona likely occurs in highly fragmented current layers, and that this typically leads to a dramatic increase in the topological complexity bey… ▽ More Global magnetic field extrapolations are now revealing the huge complexity of the Sun's corona, and in particular the structure of the boundary between open and closed magnetic flux. Moreover, recent developments indicate that magnetic reconnection in the corona likely occurs in highly fragmented current layers, and that this typically leads to a dramatic increase in the topological complexity beyond that of the equilibrium field. In this paper we investigate the consequences of reconnection at the open-closed flux boundary ("interchange reconnection") in a fragmented current layer. We demonstrate that it leads to a situation in which magnetic flux (and therefore plasma) from open and closed field regions is efficiently mixed together. This corresponds to an increase in the length and complexity of the open-closed boundary. Thus, whenever reconnection occurs at a null point or separator of the open-closed boundary, the associated separatrix arc of the so-called "S-web" in the high corona becomes not a single line but a band of finite thickness within which the open-closed flux boundary is highly structured. This has significant implications for the acceleration of the slow solar wind, for which the interaction of open and closed field is thought to be important, and may also explain the coronal origins of certain solar energetic particles. The topological structures examined contain magnetic null points, separatrices and separators, and include a model for a pseudo-streamer. The potential for understanding both the large scale morphology and fine structure observed in flare ribbons associated with coronal nulls is also discussed. △ Less

Submitted 25 May, 2015; v1 submitted 4 February, 2015; originally announced February 2015.

Journal ref: The Astrophysical Journal, 805:39 (2015)

The structure of current layers and degree of field line braiding in coronal loops

Authors: D. I. Pontin, G. Hornig

Abstract: One proposed resolution to the long-standing problem of solar coronal heating involves the buildup of magnetic energy in the corona due to turbulent motions at the photosphere that braid the coronal field, and the subsequent release of this energy via magnetic reconnection. In this paper the ideal relaxation of braided magnetic fields modelling solar coronal loops is followed. A sequence of loops… ▽ More One proposed resolution to the long-standing problem of solar coronal heating involves the buildup of magnetic energy in the corona due to turbulent motions at the photosphere that braid the coronal field, and the subsequent release of this energy via magnetic reconnection. In this paper the ideal relaxation of braided magnetic fields modelling solar coronal loops is followed. A sequence of loops with increasing braid complexity is considered, with the aim of understanding how this complexity influences the development of small scales in the magnetic field, and thus the energy available for heating. It is demonstrated that the ideally accessible force-free equilibrium for these braided fields contains current layers of finite thickness. It is further shown that for any such braided field, if a force-free equilibrium exists then it should contain current layers whose thickness is determined by length scales in the field line mapping. The thickness and intensity of the current layers follow scaling laws, and this allows us to extrapolate beyond the numerically accessible parameter regime and to place an upper bound on the braid complexity possible at coronal plasma parameters. At this threshold level the braided loop contains $10^{26}$--$10^{28}{\rm ergs}$ of available free magnetic energy, more than sufficient for a large nanoflare. △ Less

Submitted 31 March, 2015; v1 submitted 11 November, 2014; originally announced November 2014.

Comments: To appear in ApJ. 20 pages, 10 figures

Dynamic Topology and Flux Rope Evolution During Non-linear Tearing of 3D Null Point Current Sheets

Authors: P. F. Wyper, D. I Pontin

Abstract: In this work the dynamic magnetic field within a tearing-unstable three-dimensional (3D) current sheet about a magnetic null point is described in detail. We focus on the evolution of the magnetic null points and flux ropes that are formed during the tearing process. Generally, we find that both magnetic structures are created prolifically within the layer and are non-trivially related. We examine… ▽ More In this work the dynamic magnetic field within a tearing-unstable three-dimensional (3D) current sheet about a magnetic null point is described in detail. We focus on the evolution of the magnetic null points and flux ropes that are formed during the tearing process. Generally, we find that both magnetic structures are created prolifically within the layer and are non-trivially related. We examine how nulls are created and annihilated during bifurcation processes, and describe how they evolve within the current layer. The type of null bifurcation first observed is associated with the formation of pairs of flux ropes within the current layer. We also find that new nulls form within these flux ropes, both following internal reconnection and as adjacent flux ropes interact. The flux ropes exhibit a complex evolution, driven by a combination of ideal kinking and their interaction with the outflow jets from the main layer. The finite size of the unstable layer also allows us to consider the wider effects of flux rope generation. We find that the unstable current layer acts as a source of torsional MHD waves and dynamic braiding of magnetic fields. The implications of these results to several areas of heliophysics are discussed. △ Less

Submitted 23 June, 2014; originally announced June 2014.

Comments: Submitted to Physics of Plasmas (12 pages, 12 figures)

Current singularities in line-tied three-dimensional magnetic fields

Authors: I. J. D. Craig, D. I. Pontin

Abstract: This paper considers the current distributions that derive from finite amplitude perturbations of line-tied magnetic fields comprising hyperbolic field structures. The initial equilibrium on which we principally focus is a planar magnetic X-point threaded by a uniform axial field. This field is line-tied on all surfaces but subject to three-dimensional (3D) disturbances that alter the initial topo… ▽ More This paper considers the current distributions that derive from finite amplitude perturbations of line-tied magnetic fields comprising hyperbolic field structures. The initial equilibrium on which we principally focus is a planar magnetic X-point threaded by a uniform axial field. This field is line-tied on all surfaces but subject to three-dimensional (3D) disturbances that alter the initial topology. Results of ideal relaxation simulations are presented which illustrate how intense current structures form that can be related, through the influence of line-tying, to the quasi-separatrix layers (QSLs) of the initial configuration. It is demonstrated that the location within the QSL that attracts the current, and its scaling properties, are strongly dependent on the relative dimensions of the QSL with respect to the line-tied boundaries. These results are contrasted with the behavior of a line-tied 3D field containing an isolated null point. In this case, it is found that the dominant current always forms at the null, but that the collapse is inhibited when the null is closer to a line-tied boundary. △ Less

Submitted 5 June, 2014; originally announced June 2014.

Comments: To appear in ApJ

On the nature of reconnection at a solar coronal null point above a separatrix dome

Authors: D. I. Pontin, E. R. Priest, K. Galsgaard

Abstract: Three-dimensional magnetic null points are ubiquitous in the solar corona, and in any generic mixed-polarity magnetic field. We consider magnetic reconnection at an isolated coronal null point, whose fan field lines form a dome structure. We demonstrate using analytical and computational models several features of spine-fan reconnection at such a null, including the fact that substantial magnetic… ▽ More Three-dimensional magnetic null points are ubiquitous in the solar corona, and in any generic mixed-polarity magnetic field. We consider magnetic reconnection at an isolated coronal null point, whose fan field lines form a dome structure. We demonstrate using analytical and computational models several features of spine-fan reconnection at such a null, including the fact that substantial magnetic flux transfer from one region of field line connectivity to another can occur. The flux transfer occurs across the current sheet that forms around the null point during spine-fan reconnection, and there is no separator present. Also, flipping of magnetic field lines takes place in a manner similar to that observed in quasi-separatrix layer or slip-running reconnection. △ Less

Submitted 25 July, 2013; originally announced July 2013.

Comments: Accepted for publication in the Astrophysical Journal

Journal ref: ApJ, 774, 154 (2013)

Kelvin-Helmholtz instability in a current-vortex sheet at a 3D magnetic null

Authors: P. F. Wyper, D. I. Pontin

Abstract: We report here, for the first time, an observed instability of a Kelvin-Helmholtz (KH) nature occurring in a fully three-dimensional (3D) current-vortex sheet at the fan plane of a 3D magnetic null point. The current-vortex layer forms self-consistently in response to foot point driving around the spine lines of the null. The layer first becomes unstable at an intermediate distance from the null p… ▽ More We report here, for the first time, an observed instability of a Kelvin-Helmholtz (KH) nature occurring in a fully three-dimensional (3D) current-vortex sheet at the fan plane of a 3D magnetic null point. The current-vortex layer forms self-consistently in response to foot point driving around the spine lines of the null. The layer first becomes unstable at an intermediate distance from the null point, with the instability being characterized by a rippling of the fan surface and a filamentation of the current density and vorticity in the shear layer. Owing to the 3D geometry of the shear layer, a branching of the current filaments and vortices is observed. The instability results in a mixing of plasma between the two topologically distinct regions of magnetic flux on either side of the fan separatrix surface, as flux is reconnected across this surface. We make a preliminary investigation of the scaling of the system with the dissipation parameters. Our results indicate that the fan plane separatrix surface is an ideal candidate for the formation of current-vortex sheets in complex magnetic fields and, therefore, the enhanced heating and connectivity change associated with the instabilities of such layers. △ Less

Submitted 25 March, 2013; originally announced March 2013.

Comments: 12 pages, 14 figures, to appear in Physics of Plasmas

Journal ref: Phys. Plasmas 20, 032117 (2013)

On the formation of current sheets in response to the compression or expansion of a potential magnetic field

Authors: D. I. Pontin, Y. -M. Huang

Abstract: The compression or expansion of a magnetic field that is initially potential is considered. It was recently suggested by Janse & Low [2009, ApJ, 690, 1089] that, following the volumetric deformation, the relevant lowest energy state for the magnetic field is another potential magnetic field that in general contains tangential discontinuities (current sheets). Here we examine this scenario directly… ▽ More The compression or expansion of a magnetic field that is initially potential is considered. It was recently suggested by Janse & Low [2009, ApJ, 690, 1089] that, following the volumetric deformation, the relevant lowest energy state for the magnetic field is another potential magnetic field that in general contains tangential discontinuities (current sheets). Here we examine this scenario directly using a numerical relaxation method that exactly preserves the topology of the magnetic field. It is found that of the magnetic fields discussed by Janse & Low, only those containing magnetic null points develop current singularities during an ideal relaxation, while the magnetic fields without null points relax toward smooth force-free equilibria with finite non-zero current. △ Less

Submitted 4 July, 2012; originally announced July 2012.

Comments: Accepted for publication in ApJ

Theory of magnetic reconnection in solar and astrophysical plasmas

Authors: D. I. Pontin

Abstract: Magnetic reconnection is a fundamental process in a plasma that facilitates the release of energy stored in the magnetic field by permitting a change in the magnetic topology. In this article we present a review of the current state of understanding of magnetic reconnection. We discuss theoretical results regarding the formation of current sheets in complex 3D magnetic fields, and describe the fun… ▽ More Magnetic reconnection is a fundamental process in a plasma that facilitates the release of energy stored in the magnetic field by permitting a change in the magnetic topology. In this article we present a review of the current state of understanding of magnetic reconnection. We discuss theoretical results regarding the formation of current sheets in complex 3D magnetic fields, and describe the fundamental differences between reconnection in two and three dimensions. We go on to outline recent developments in modelling of reconnection with kinetic theory, as well as in the MHD framework where a number of new 3D reconnection regimes have been identified. We discuss evidence from observations and simulations of solar system plasmas that support this theory, and summarise some prominent locations in which this new reconnection theory is relevant in astrophysical plasmas. △ Less

Submitted 3 July, 2012; v1 submitted 17 February, 2012; originally announced February 2012.

Journal ref: Phil. Trans R. Soc. A, 370, 3169-3192 (2012)

Heating of braided coronal loops

Authors: A. L. Wilmot-Smith, D. I. Pontin, A. R. Yeates, G. Hornig

Abstract: We investigate the relaxation of braided magnetic loops in order to find out how the type of braiding via footpoint motions affects resultant heating of the loop. Two magnetic loops, braided in different ways, are used as initial conditions in resistive MHD simulations and their subsequent evolution is studied. The fields both undergo a resistive relaxation in which current sheets form and fragmen… ▽ More We investigate the relaxation of braided magnetic loops in order to find out how the type of braiding via footpoint motions affects resultant heating of the loop. Two magnetic loops, braided in different ways, are used as initial conditions in resistive MHD simulations and their subsequent evolution is studied. The fields both undergo a resistive relaxation in which current sheets form and fragment and the system evolves towards a state of lower energy. In one case this relaxation is very efficient with current sheets filling the volume and homogeneous heating of the loop occurring. In the other case fewer current sheets develop, less magnetic energy is released in the process and a patchy heating of the loop results. The two cases, although very similar in their setup, can be distinguished by the mixing properties of the photospheric driver. The mixing can be measured by the topological entropy of the plasma flow, an observable quantity. △ Less

Submitted 4 November, 2011; originally announced November 2011.

Comments: 16 pages, 13 Figures. To be published in Astronomy and Astrophysics

Generalised models for torsional spine and fan magnetic reconnection

Authors: D. I. Pontin, A. K. Al-Hachami, K. Galsgaard

Abstract: Three-dimensional null points are present in abundance in the solar corona, and the same is likely to be true in other astrophysical environments. Recent studies suggest that reconnection at such 3D nulls may play an important role in the coronal dynamics. In this paper the properties of the torsional spine and torsional fan modes of magnetic reconnection at 3D nulls are investigated. New analytic… ▽ More Three-dimensional null points are present in abundance in the solar corona, and the same is likely to be true in other astrophysical environments. Recent studies suggest that reconnection at such 3D nulls may play an important role in the coronal dynamics. In this paper the properties of the torsional spine and torsional fan modes of magnetic reconnection at 3D nulls are investigated. New analytical models are developed, which for the first time include a current layer that is spatially localised around the null, extending along either the spine or the fan of the null. These are complemented with numerical simulations. The principal aim is to investigate the effect of varying the degree of asymmetry of the null point magnetic field on the resulting reconnection process - where previous studies always considered a non-generic radially symmetric null. The geometry of the current layers within which torsional spine and torsional fan reconnection occur is found to be strongly dependent on the symmetry of the magnetic field. Torsional spine reconnection still occurs in a narrow tube around the spine, but with elliptical cross-section when the fan eigenvalues are different, and with the short axis of the ellipse being along the strong field direction. The spatiotemporal peak current, and the peak reconnection rate attained, are found not to depend strongly on the degree of asymmetry. For torsional fan reconnection, the reconnection occurs in a planar disk in the fan surface, which is again elliptical when the symmetry of the magnetic field is broken. The short axis of the ellipse is along the weak field direction, with the current being peaked in these weak field regions. The peak current and peak reconnection rate in this case are clearly dependent on the asymmetry, with the peak current increasing but the reconnection rate decreasing as the degree of asymmetry is increased. △ Less

Submitted 27 September, 2011; v1 submitted 13 May, 2011; originally announced May 2011.

Journal ref: Astronomy & Astrophysics, vol. 533, article A78, 2011

Steady state reconnection at a single 3D magnetic null point

Authors: K. Galsgaard, D. I. Pontin

Abstract: To systematically stress a rotationally symmetric 3D magnetic null point by advecting the opposite footpoints of the spine axis in opposite directions. This stress eventually concentrates in the vicinity of the null point forming a local current sheet through which magnetic reconnection takes place. The aim is to look for a steady state evolution of the current sheet dynamics which may provide sca… ▽ More To systematically stress a rotationally symmetric 3D magnetic null point by advecting the opposite footpoints of the spine axis in opposite directions. This stress eventually concentrates in the vicinity of the null point forming a local current sheet through which magnetic reconnection takes place. The aim is to look for a steady state evolution of the current sheet dynamics which may provide scaling relations for various characteristic parameters of the system. The evolution is followed by solving numerically the non-ideal MHD equations in a Cartesian domain. The null point is embedded in an initially constant density and temperature plasma. It is shown that a quasi-steady reconnection process can be set up at a 3D null by continuous shear driving. It appears that a true steady state in unlikely to be realised as the current layer tends to grow until restricted by the geometry of the computational domain and imposed driving profile. However, ratios between characteristic quantities clearly settle after some time to stable values -- so that the evolution is quasi-steady. The experiments show a number of scaling relations, but they do not provide a clear consensus for extending to lower magnetic resistivity or faster driving velocities. More investigations are needed to fully clarify the properties of current sheets at magnetic null points. △ Less

Submitted 11 February, 2011; originally announced February 2011.

Comments: 9 pages, 11 figures, accepted for publication in A&A

Three-dimensional magnetic reconnection regimes: A review

Authors: D. I. Pontin

Abstract: The magnetic field in many astrophysical plasmas -- such as the Solar corona and Earth's magnetosphere -- has been shown to have a highly complex, three-dimensional structure. Recent advances in theory and computational simulations have shown that reconnection in these fields also has a three-dimensional nature, in contrast to the widely used two-dimensional (or 2.5-dimensional) models. Here we di… ▽ More The magnetic field in many astrophysical plasmas -- such as the Solar corona and Earth's magnetosphere -- has been shown to have a highly complex, three-dimensional structure. Recent advances in theory and computational simulations have shown that reconnection in these fields also has a three-dimensional nature, in contrast to the widely used two-dimensional (or 2.5-dimensional) models. Here we discuss the underlying theory of three-dimensional magnetic reconnection. We also review a selection of new models that illustrate the current state of the art, as well as highlighting the complexity of energy release processes mediated by reconnection in complicated three-dimensional magnetic fields. △ Less

Submitted 5 January, 2011; originally announced January 2011.

Comments: Accepted for publication in Advances in Space Research. (Based on paper presented at 38th COSPAR meeting, Bremen, 2010.)

Dynamics of braided coronal loops II: Cascade to multiple small-scale reconnection events

Authors: D. I. Pontin, A. L. Wilmot-Smith, G. Hornig, K. Galsgaard

Abstract: Aims: Our aim is to investigate the resistive relaxation of a magnetic loop that contains braided magnetic flux but no net current or helicity. The loop is subject to line-tied boundary conditions. We investigate the dynamical processes that occur during this relaxation, in particular the magnetic reconnection that occurs, and discuss the nature of the final equilibrium. Methods: The three-dimensi… ▽ More Aims: Our aim is to investigate the resistive relaxation of a magnetic loop that contains braided magnetic flux but no net current or helicity. The loop is subject to line-tied boundary conditions. We investigate the dynamical processes that occur during this relaxation, in particular the magnetic reconnection that occurs, and discuss the nature of the final equilibrium. Methods: The three-dimensional evolution of a braided magnetic field is followed in a series of resistive MHD simulations. Results: It is found that, following an instability within the loop, a myriad of thin current layers forms, via a cascade-like process. This cascade becomes more developed and continues for a longer period of time for higher magnetic Reynolds number. During the cascade, magnetic flux is reconnected multiple times, with the level of this `multiple reconnection' positively correlated with the magnetic Reynolds number. Eventually the system evolves into a state with no more small-scale current layers. This final state is found to approximate a non-linear force-free field consisting of two flux tubes of oppositely-signed twist embedded in a uniform background field. △ Less

Submitted 21 October, 2010; v1 submitted 30 March, 2010; originally announced March 2010.

Comments: To be published in A&A

Dynamics of braided coronal loops - I. Onset of magnetic reconnection

Authors: A. L. Wilmot-Smith, D. I. Pontin, G. Hornig

Abstract: The response of the solar coronal magnetic field to small-scale photospheric boundary motions including the possible formation of current sheets via the Parker scenario is one of open questions of solar physics. Here we address the problem via a numerical simulation. The three-dimensional evolution of a braided magnetic field which is initially close to a force-free state is followed using a res… ▽ More The response of the solar coronal magnetic field to small-scale photospheric boundary motions including the possible formation of current sheets via the Parker scenario is one of open questions of solar physics. Here we address the problem via a numerical simulation. The three-dimensional evolution of a braided magnetic field which is initially close to a force-free state is followed using a resistive MHD code. A long-wavelength instability takes place and leads to the formation of two thin current layers. Magnetic reconnection occurs across the current sheets with three-dimensional features shown, including an elliptic magnetic field structure about the reconnection site, and results in an untwisting of the global field structure. △ Less

Submitted 18 March, 2010; v1 submitted 11 January, 2010; originally announced January 2010.

Comments: 18 pages, 9 figures

3D Null Point Reconnection Regimes

Authors: E. R. Priest, D. I. Pontin

Abstract: Recent advances in theory and computational experiments have shown the need to refine the previous categorisation of magnetic reconnection at three-dimensional null points -- points at which the magnetic field vanishes. We propose here a division into three different types, depending on the nature of the flow near the spine and fan of the null. The spine is an isolated field line which approache… ▽ More Recent advances in theory and computational experiments have shown the need to refine the previous categorisation of magnetic reconnection at three-dimensional null points -- points at which the magnetic field vanishes. We propose here a division into three different types, depending on the nature of the flow near the spine and fan of the null. The spine is an isolated field line which approaches the null (or recedes from it), while the fan is a surface of field lines which recede from it (or approach it). So-called "torsional spine reconnection" occurs when field lines in the vicinity of the fan rotate, with current becoming concentrated along the spine, so that nearby field lines undergo rotational slippage. In "torsional fan reconnection" field lines near the spine rotate and create a current that is concentrated in the fan with a rotational flux mismatch and rotational slippage. In both of these regimes, the spine and fan are perpendicular and there is no flux transfer across spine or fan. The third regime, called "spine-fan reconnection", is the most common in practice and combines elements of the previous spine and fan models. In this case, in response to a generic shearing motion, the null point collapses to form a current sheet that is focused at the null itself, in a sheet that locally spans both the spine and fan. In this regime the spine and fan are no longer perpendicular and there is flux transfer across both of them. △ Less

Submitted 16 October, 2009; originally announced October 2009.

Comments: 16 pages, 10 figures. Accepted for publication in Physics of Plasmas

Journal ref: Phys. Plasmas 16, 122101 (2009)

Magnetic reconnection at 3D null points: effect of magnetic field asymmetry

Authors: A. K. Al-Hachami, D. I. Pontin

Abstract: The aim of this paper is to investigate the properties of magnetic reconnection at a 3D null point, with respect to their dependence on the symmetry of the magnetic field around the null. In particular we examine the rate of flux transport across the null point with symmetric/asymmetric diffusion regions, as well as how the current sheet forms in time, and its properties. Mathematical modelling… ▽ More The aim of this paper is to investigate the properties of magnetic reconnection at a 3D null point, with respect to their dependence on the symmetry of the magnetic field around the null. In particular we examine the rate of flux transport across the null point with symmetric/asymmetric diffusion regions, as well as how the current sheet forms in time, and its properties. Mathematical modelling and finite difference resistive MHD simulations are used. It is found that the basic structure of the mode of magnetic reconnection considered is unaffected by varying the magnetic field symmetry, that is, the plasma flow is found cross both the spine and fan of the null. However, the peak intensity and dimensions of the current sheet are dependent on the symmetry/ asymmetry of the field lines. As a result, the reconnection rate is also found to be strongly dependent on the field asymmetry. △ Less

Submitted 1 February, 2010; v1 submitted 31 August, 2009; originally announced August 2009.

Comments: 23 pages, 9 figures

Magnetic Braiding and Quasi-Separatrix Layers

Authors: A. L. Wilmot-Smith, G. Hornig, D. I. Pontin

Abstract: The squashing factor Q (Titov et al. 2002), a property of the magnetic field line mapping, has been suggested as an indicator for the formation of current sheets, and subsequently magnetic reconnection, in astrophysical plasmas. Here we test this hypothesis for a particular class of braided magnetic fields which serve as a model for solar coronal loops. We explore the relationship between Quasi-… ▽ More The squashing factor Q (Titov et al. 2002), a property of the magnetic field line mapping, has been suggested as an indicator for the formation of current sheets, and subsequently magnetic reconnection, in astrophysical plasmas. Here we test this hypothesis for a particular class of braided magnetic fields which serve as a model for solar coronal loops. We explore the relationship between Quasi-Separatrix Layers (QSLs), that is, layer-like structures with high Q value, electric currents and integrated parallel currents, the latter being a quantity closely related to the reconnection rate. It is found that as the degree of braiding of the magnetic field is increased the maximum values of Q increase exponentially. At the same time, the distribution of Q becomes increasingly filamentary, with the width of the high-Q layers exponentially decreasing. This is accompanied by an increase in the number of layers so that as the field is increasingly braided the volume becomes occupied by a myriad of thin QSLs. QSLs are not found to be good predictors of current features in this class of braided fields. Indeed, despite the presence of multiple QSLs, the current associated with the field remains smooth and large-scale under ideal relaxation; the field dynamically adjusts to a smooth equilibrium. Regions of high Q are found to be better related to regions of high integrated parallel current than to actual current sheets. △ Less

Submitted 22 July, 2009; originally announced July 2009.

Comments: 17 pages, 8 figures

Lagrangian relaxation schemes for calculating force-free magnetic fields, and their limitations

Authors: D. I. Pontin, G. Hornig, A. L. Wilmot-Smith, I. J. D. Craig

Abstract: Force-free magnetic fields are important in many astrophysical settings. Determining the properties of such force-free fields -- especially smoothness and stability properties -- is crucial to understanding many key phenomena in astrophysical plasmas, for example energy release processes that heat the plasma and lead to dynamic or explosive events. Here we report on a serious limitation on the c… ▽ More Force-free magnetic fields are important in many astrophysical settings. Determining the properties of such force-free fields -- especially smoothness and stability properties -- is crucial to understanding many key phenomena in astrophysical plasmas, for example energy release processes that heat the plasma and lead to dynamic or explosive events. Here we report on a serious limitation on the computation of force-free fields that has the potential to invalidate the results produced by numerical force-free field solvers even for cases in which they appear to converge (at fixed grid resolution) to an equilibrium magnetic field. In the present work we discuss this problem within the context of a Lagrangian relaxation scheme that conserves magnetic flux and div(B) identically. Error estimates are introduced to assess the quality of the calculated equilibrium. We go on to present an algorithm, based on re-writing the curl operation via Stokes' theorem, for calculating the current which holds great promise for improving dramatically the accuracy of the Lagrangian relaxation procedure. △ Less

Submitted 6 March, 2009; originally announced March 2009.

Comments: Submitted to ApJ

Journal ref: Astrophys.J.700:1449-1455,2009

Magnetic Braiding and Parallel Electric Fields

Authors: A. L. Wilmot-Smith, G. Hornig, D. I. Pontin

Abstract: The braiding of the solar coronal magnetic field via photospheric motions - with subsequent relaxation and magnetic reconnection -- is one of the most widely debated ideas of solar physics. We readdress the theory in the light of developments in three-dimensional magnetic reconnection theory. It is known that the integrated parallel electric field along field lines is the key quantity determinin… ▽ More The braiding of the solar coronal magnetic field via photospheric motions - with subsequent relaxation and magnetic reconnection -- is one of the most widely debated ideas of solar physics. We readdress the theory in the light of developments in three-dimensional magnetic reconnection theory. It is known that the integrated parallel electric field along field lines is the key quantity determining the rate of reconnection, in contrast with the two-dimensional case where the electric field itself is the important quantity. We demonstrate that this difference becomes crucial for sufficiently complex magnetic field structures. A numerical method is used to relax a braided magnetic field to an ideal force-free equilibrium; that equilibrium is found to be smooth, with only large- scale current structures. However, the equilibrium is shown to have a highly filamentary integrated parallel current structure with extremely short length- scales. An analytical model is developed to show that, in a coronal situation, the length scales associated with the integrated parallel current structures will rapidly decrease with increasing complexity, or degree of braiding, of the magnetic field. Analysis shows the decrease in these length scales will, for any finite resistivity, eventually become inconsistent with the stability of a force- free field. Thus the inevitable consequence of the magnetic braiding process is shown to be a loss of equilibrium of the coronal field, probably via magnetic reconnection events. △ Less

Submitted 8 October, 2008; originally announced October 2008.

Comments: 19 pages with 7 figures

Journal ref: Astrophys.J.696:1339-1347,2009

Current sheets at three-dimensional magnetic nulls: Effect of compressibility

Authors: D. I. Pontin, A. Bhattacharjee, K. Galsgaard

Abstract: The nature of current sheet formation in the vicinity of three-dimensional (3D) magnetic null points is investigated. The particular focus is upon the effect of the compressibility of the plasma on the qualitative and quantitative properties of the current sheet. An initially potential 3D null is subjected to shearing perturbations, as in a previous paper [Pontin et al., Phys. Plasmas, in press… ▽ More The nature of current sheet formation in the vicinity of three-dimensional (3D) magnetic null points is investigated. The particular focus is upon the effect of the compressibility of the plasma on the qualitative and quantitative properties of the current sheet. An initially potential 3D null is subjected to shearing perturbations, as in a previous paper [Pontin et al., Phys. Plasmas, in press (2007)]. It is found that as the incompressible limit is approached, the collapse of the null point is suppressed, and an approximately planar current sheet aligned to the fan plane is present instead. This is the case regardless of whether the spine or fan of the null is sheared. Both the peak current and peak reconnection rate are reduced. The results have a bearing on previous analytical solutions for steady-state reconnection in incompressible plasmas, implying that fan current sheet solutions are dynamically accessible, while spine current sheet solutions are not. △ Less

Submitted 26 April, 2007; v1 submitted 17 January, 2007; originally announced January 2007.

Comments: to appear in Physics of Plasmas. This version contains updated figures and references, additional discussion, and typos are fixed. This is the second in a series of papers - the first of which (by the same authors) is located at astro-ph/0701462. A version with higher quality figures can be found at http://www.maths.dundee.ac.uk/~dpontin/

Current amplification and magnetic reconnection at a 3D null point. I - Physical characteristics

Authors: D. I. Pontin, K. Galsgaard

Abstract: The behaviour of magnetic perturbations of an initially potential three-dimensional equilibrium magnetic null point are investigated. The basic components which constitute a typical disturbance are taken to be rotations and shears, in line with previous work. The spine and fan of the null point (the field lines which asymptotically approach or recede from the null) are subjected to such rotation… ▽ More The behaviour of magnetic perturbations of an initially potential three-dimensional equilibrium magnetic null point are investigated. The basic components which constitute a typical disturbance are taken to be rotations and shears, in line with previous work. The spine and fan of the null point (the field lines which asymptotically approach or recede from the null) are subjected to such rotational and shear perturbations, using three-dimensional magnetohydrodynamic simulations. It is found that rotations of the fan plane and about the spine lead to current sheets which are spatially diffuse in at least one direction, and form in the locations of the spine and fan. However, shearing perturbations lead to 3D-localised current sheets focused at the null point itself. In addition, rotations are associated with a growth of current parallel to the spine, driving rotational flows and a type of rotational reconnection. Shears, on the other hand, cause a current through the null which is parallel to the fan plane, and are associated with stagnation-type flows and field line reconnection across both the spine and fan. The importance of the parallel electric field, and its meaning as a reconnection rate, are discussed. △ Less

Submitted 19 January, 2007; originally announced January 2007.

Comments: Accepted for publication in JGR (Space Physics)

Journal ref: J.Geophys.Res.112:A03103,2007

Current sheet formation and non-ideal behaviour at three-dimensional magnetic null points

Authors: D. I. Pontin, A. Bhattacharjee, K. Galsgaard

Abstract: The nature of the evolution of the magnetic field, and of current sheet formation, at three-dimensional (3D) magnetic null points is investigated. A kinematic example is presented which demonstrates that for certain evolutions of a 3D null (specifically those for which the ratios of the null point eigenvalues are time-dependent) there is no possible choice of boundary conditions which renders th… ▽ More The nature of the evolution of the magnetic field, and of current sheet formation, at three-dimensional (3D) magnetic null points is investigated. A kinematic example is presented which demonstrates that for certain evolutions of a 3D null (specifically those for which the ratios of the null point eigenvalues are time-dependent) there is no possible choice of boundary conditions which renders the evolution of the field at the null ideal. Resistive MHD simulations are described which demonstrate that such evolutions are generic. A 3D null is subjected to boundary driving by shearing motions, and it is shown that a current sheet localised at the null is formed. The qualitative and quantitative properties of the current sheet are discussed. Accompanying the sheet development is the growth of a localised parallel electric field, one of the signatures of magnetic reconnection. Finally, the relevance of the results to a recent theory of turbulent reconnection is discussed. △ Less

Submitted 22 March, 2007; v1 submitted 16 January, 2007; originally announced January 2007.

Comments: to appear in Phys. Plasmas. A version with higher quality figures can be found at http://www.maths.dundee.ac.uk/~dpontin/ In this replacement version, typos have been corrected, and in addition references and some further discussion added