-
Helicity-conserving relaxation in unstable and merging magnetic flux ropes
Authors:
Philippa Browning,
Mykola Gordovskyy,
Alan Hood
Abstract:
Twisted magnetic flux ropes are reservoirs of free magnetic energy. In a highly-conducting plasma such as the solar corona, energy release through multiple magnetic reconnections can be modelled as a helicity-conserving relaxation to a minimum energy state. One possible trigger for this relaxation is the ideal kink instability in a twisted flux rope. We show that this provides a good description f…
▽ More
Twisted magnetic flux ropes are reservoirs of free magnetic energy. In a highly-conducting plasma such as the solar corona, energy release through multiple magnetic reconnections can be modelled as a helicity-conserving relaxation to a minimum energy state. One possible trigger for this relaxation is the ideal kink instability in a twisted flux rope. We show that this provides a good description for confined solar flares, and develop from idealised cylindrical models to realistic models of coronal loops. Using 3D magnetohydrodynamic simulations combined with test-particle simulations of non-thermal electrons and ions, we predict multiple observational signatures of such flares. We then show how interactions and mergers of flux ropes can release free magnetic energy, using relaxation theory to complement simulations of merging-compression formation in spherical tokamaks and heating avalanches in the solar corona.
△ Less
Submitted 16 August, 2023;
originally announced August 2023.
-
An Anisotropic Density Turbulence Model from the Sun to 1 au Derived From Radio Observations
Authors:
Eduard P. Kontar,
A. Gordon Emslie,
Daniel L. Clarkson,
Xingyao Chen,
Nicolina Chrysaphi,
Francesco Azzollini,
Natasha L. S. Jeffrey,
Mykola Gordovskyy
Abstract:
Solar radio bursts are strongly affected by radio-wave scattering on density inhomogeneities, changing their observed time characteristics, sizes, and positions. The same turbulence causes angular broadening and scintillation of galactic and extra-galactic compact radio sources observed through the solar atmosphere. Using large-scale simulations of radio-wave transport, the characteristics of anis…
▽ More
Solar radio bursts are strongly affected by radio-wave scattering on density inhomogeneities, changing their observed time characteristics, sizes, and positions. The same turbulence causes angular broadening and scintillation of galactic and extra-galactic compact radio sources observed through the solar atmosphere. Using large-scale simulations of radio-wave transport, the characteristics of anisotropic density turbulence from $0.1 \, R_\odot$ to $1$ au are explored. For the first time, a profile of heliospheric density fluctuations is deduced that accounts for the properties of extra-solar radio sources, solar radio bursts, and in-situ density fluctuation measurements in the solar wind at $1$ au. The radial profile of the spectrum-weighted mean wavenumber of density fluctuations (a quantity proportional to the scattering rate of radio-waves) is found to have a broad maximum at around $(4-7) \, R_\odot$, where the slow solar wind becomes supersonic. The level of density fluctuations at the inner scale (which is consistent with the proton resonance scale) decreases with heliocentric distance as $\langleδ{n_i}^2 \rangle (r) \simeq 2 \times 10^7 \, (r/R_\odot-1)^{-3.7}$ cm$^{-6}$. Due to scattering, the apparent positions of solar burst sources observed at frequencies between $0.1$ and $300$ MHz are computed to be essentially cospatial and to have comparable sizes, for both fundamental and harmonic emission. Anisotropic scattering is found to account for the shortest solar radio burst decay times observed, and the required wavenumber anisotropy is $q_\parallel/q_\perp =0.25-0.4$, depending on whether fundamental or harmonic emission is involved. The deduced radio-wave scattering rate paves the way to quantify intrinsic solar radio burst characteristics.
△ Less
Submitted 29 August, 2023; v1 submitted 10 August, 2023;
originally announced August 2023.
-
Particle acceleration and their escape into the heliosphere in solar flares with open magnetic field
Authors:
Mykola Gordovskyy,
Philippa K. Browning,
Kanya Kusano,
Satoshi Inoue,
Gregory E. Vekstein
Abstract:
Energetic particle populations in the solar corona and in the heliosphere appear to have different characteristics even when produced in the same solar flare. It is not clear what causes this difference: properties of the acceleration region, the large-scale magnetic field configuration in the flare, or particle transport effects, such as scattering. In this study we use a combination of magnetohy…
▽ More
Energetic particle populations in the solar corona and in the heliosphere appear to have different characteristics even when produced in the same solar flare. It is not clear what causes this difference: properties of the acceleration region, the large-scale magnetic field configuration in the flare, or particle transport effects, such as scattering. In this study we use a combination of magnetohydrodynamic and test-particle approaches to investigate magnetic reconnection, particle acceleration and transport in two solar flares: an M-class flare on June 19th, 2013, and an X-class flare on September 6th, 2011. We show that in both events , the same regions are responsible for the acceleration of particles remaining in the coronal and being ejected towards the heliosphere. However, the magnetic field structure around the acceleration region acts as a filter, resulting in different characteristics (such as energy spectra) acquired by these two populations. We argue that this effect is an intrinsic property of particle acceleration in the current layers created by the interchange reconnection and, therefore, may be ubiquitous, particularly, in non-eruptive solar flares with substantial particle emission into the heliosphere.
△ Less
Submitted 30 May, 2023;
originally announced May 2023.
-
Solar Radio Spikes and Type IIIb Striae Manifestations of Sub-second Electron Acceleration Triggered by a Coronal Mass Ejection
Authors:
Daniel L. Clarkson,
Eduard P. Kontar,
Nicole Vilmer,
Mykola Gordovskyy,
Xingyao Chen,
Nicolina Chrysaphi
Abstract:
Understanding electron acceleration associated with magnetic energy release at sub-second scales presents a major challenges in solar physics. Solar radio spikes observed as sub-second, narrow bandwidth bursts with $Δ{f}/f\sim10^{-3}-10^{-2}$ are indicative of sub-second evolution of the electron distribution. We present a statistical analysis of frequency, and time-resolved imaging of individual…
▽ More
Understanding electron acceleration associated with magnetic energy release at sub-second scales presents a major challenges in solar physics. Solar radio spikes observed as sub-second, narrow bandwidth bursts with $Δ{f}/f\sim10^{-3}-10^{-2}$ are indicative of sub-second evolution of the electron distribution. We present a statistical analysis of frequency, and time-resolved imaging of individual spikes and Type IIIb striae associated with a coronal mass ejection (CME). LOFAR imaging reveals that co-temporal ($<2$ s) spike and striae intensity contours almost completely overlap. On average, both burst types have similar source size with fast expansion at millisecond scales. The radio source centroid velocities are often superluminal, and independent of frequency over 30-45 MHz. The CME perturbs the field geometry, leading to increased spike emission likely due to frequent magnetic reconnection. As the field restores towards the prior configuration, the observed sky-plane emission locations drift to increased heights over tens of minutes. Combined with previous observations above 1 GHz, average decay time and source size estimates follow $\sim1/f$ dependency over three decades in frequency, similar to radio-wave scattering predictions. Both time and spatial characteristics of the bursts between 30-70 MHz are consistent with radio-wave scattering with strong anisotropy of the density fluctuation spectrum. Consequently, the site of radio-wave emission does not correspond to the observed burst locations and implies acceleration and emission near the CME flank. The bandwidths suggest intrinsic emission source sizes $<1$ arcsec at 30 MHz, and magnetic field strengths a factor of two larger than average in events that produce decameter spikes.
△ Less
Submitted 22 February, 2023;
originally announced February 2023.
-
Oscillatory reconnection and waves driven by merging magnetic flux ropes in solar flares
Authors:
J. Stewart,
P. K. Browning,
M. Gordovskyy
Abstract:
Oscillatory reconnection is a process that has been suggested to underlie several solar and stellar phenomena, and is likely to play an important role in transient events such as flares. Quasi-periodic pulsations (QPPs) in flare emissions may be a manifestation of oscillatory reconnection, but the underlying mechanisms remain uncertain. In this paper, we present 2D magnetohydrodynamic (MHD) simula…
▽ More
Oscillatory reconnection is a process that has been suggested to underlie several solar and stellar phenomena, and is likely to play an important role in transient events such as flares. Quasi-periodic pulsations (QPPs) in flare emissions may be a manifestation of oscillatory reconnection, but the underlying mechanisms remain uncertain. In this paper, we present 2D magnetohydrodynamic (MHD) simulations of two current-carrying magnetic flux ropes with an out-of-plane magnetic field undergoing oscillatory reconnection in which the two flux ropes merge into a single flux rope. We find that oscillatory reconnection can occur intrinsically without an external oscillatory driver during flux rope coalescence, which may occur both during large-scale coronal loop interactions and the merging of plasmoids in fragmented current sheets. Furthermore, we demonstrate that radially propagating non-linear waves are produced in the aftermath of flux rope coalescence, due to the post-reconnection oscillations of the merged flux rope. The behaviour of these waves is found to be almost independent of the initial out-of-plane magnetic field. It is estimated that the waves emitted through merging coronal loops and merging plasmoids in loop-top current sheets would have a typical phase speed of 90 km/s and 900 km/s respectively. It is possible that the properties of the waves emitted during flux rope coalescence could be used as a diagnostic tool to determine physical parameters within a coalescing region.
△ Less
Submitted 6 May, 2022;
originally announced May 2022.
-
Pulsations of microwave emission from a solar flare in a twisted loop caused by intrinsic MHD oscillations
Authors:
C. Smith,
M. Gordovskyy,
P. K. Browning
Abstract:
We present results revealing microwave pulsations produced in a model of a flaring twisted solar coronal loop, without any external oscillatory driver. Two types of oscillations are identified: slowly-decaying oscillations with a period of about 70-75s and amplitude of about 5-10% seen in loops both with and without energetic electrons, and oscillations with period of about 40s and amplitude of a…
▽ More
We present results revealing microwave pulsations produced in a model of a flaring twisted solar coronal loop, without any external oscillatory driver. Two types of oscillations are identified: slowly-decaying oscillations with a period of about 70-75s and amplitude of about 5-10% seen in loops both with and without energetic electrons, and oscillations with period of about 40s and amplitude of a few tens of percent observed only in loops with energetic electrons for about 100s after onset of fast energy release. We interpret the longer-period oscillations as the result of a standing kink mode modulating the average magnetic field strength in the loop, whilst the short-period intermittent oscillations associated with energetic electrons are likely to be produced by fast variations of the electric field which produces energetic electrons in this scenario. The slowly-decaying oscillations can explain the quasi-periodic pulsations often observed in the flaring corona.
△ Less
Submitted 20 January, 2022;
originally announced January 2022.
-
Sizes and shapes of sources in solar metric radio bursts
Authors:
M. Gordovskyy,
E. P. Kontar,
D. L. Clarkson,
N. Chrysaphi,
P. K. Browning
Abstract:
Metric and decametric radio-emissions from the Sun are the only direct source of information about the dynamics of non-thermal electrons in the upper corona. In addition, the combination of spectral and imaging (sizes, shapes, and positions) observations of low-frequency radio sources can be used as a unique diagnostic tool to probe plasma turbulence in the solar corona and inner heliosphere. The…
▽ More
Metric and decametric radio-emissions from the Sun are the only direct source of information about the dynamics of non-thermal electrons in the upper corona. In addition, the combination of spectral and imaging (sizes, shapes, and positions) observations of low-frequency radio sources can be used as a unique diagnostic tool to probe plasma turbulence in the solar corona and inner heliosphere. The geometry of the low-frequency sources and its variation with frequency are still not understood, primarily due to the relatively low spatial resolution available for solar observations. Here we report the first detailed multi-frequency analysis of the sizes of solar radio sources observed by the Low-Frequency Array (LOFAR). Furthermore, we investigate the source shapes by approximating the derived intensity distributions using 2D Gaussian profiles with elliptical half-maximum contours. These measurements have been made possible by a novel empirical method for evaluating the instrumental and ionospheric effects on radio maps based on known source observations. The obtained deconvolved sizes of the sources are found to be smaller than previous estimations, and often show higher ellipticity. The sizes and ellipticities of the sources inferred using 2D Gaussian approximation, and their variation with frequency are consistent with models of anisotropic radio-wave scattering in the solar corona.
△ Less
Submitted 15 November, 2021;
originally announced November 2021.
-
First Frequency-Time-Resolved Imaging Spectroscopy Observations of Solar Radio Spikes
Authors:
Daniel L. Clarkson,
Eduard P. Kontar,
Mykola Gordovskyy,
Nicolina Chrysaphi,
Nicole Vilmer
Abstract:
Solar radio spikes are short duration and narrow bandwidth fine structures in dynamic spectra observed from GHz to tens of MHz range. Their very short duration and narrow frequency bandwidth are indicative of sub-second small-scale energy release in the solar corona, yet their origin is not understood. Using the LOw Frequency ARray (LOFAR), we present spatially, frequency and time resolved observa…
▽ More
Solar radio spikes are short duration and narrow bandwidth fine structures in dynamic spectra observed from GHz to tens of MHz range. Their very short duration and narrow frequency bandwidth are indicative of sub-second small-scale energy release in the solar corona, yet their origin is not understood. Using the LOw Frequency ARray (LOFAR), we present spatially, frequency and time resolved observations of individual radio spikes associated with a coronal mass ejection (CME). Individual radio spike imaging demonstrates that the observed area is increasing in time and the centroid positions of the individual spikes move superluminally parallel to the solar limb. Comparison of spike characteristics with that of individual Type IIIb striae observed in the same event show similarities in duration, bandwidth, drift rate, polarization and observed area, as well the spike and striae motion in the image plane suggesting fundamental plasma emission with the spike emission region on the order of ${\sim}\:10^8$ cm, with brightness temperature as high as $10^{13}$ K. The observed spatial, spectral, and temporal properties of the individual spike bursts are also suggesting the radiation responsible for spikes escapes through anisotropic density turbulence in closed loop structures with scattering preferentially along the guiding magnetic field oriented parallel to the limb in the scattering region. The dominance of scattering on the observed time profile suggests the energy release time is likely to be shorter than what is often assumed. The observations also imply that the density turbulence anisotropy along closed magnetic field lines is higher than along open field lines.
△ Less
Submitted 16 August, 2021; v1 submitted 13 August, 2021;
originally announced August 2021.
-
Sub-second time evolution of Type III solar radio burst sources at fundamental and harmonic frequencies
Authors:
Xingyao Chen,
Eduard P. Kontar,
Nicolina Chrysaphi,
Natasha L. S. Jeffrey,
Mykola Gordovskyy,
Yihua Yan,
Baolin Tan
Abstract:
Recent developments in astronomical radio telescopes opened new opportunities in imaging and spectroscopy of solar radio bursts at sub-second timescales. Imaging in narrow frequency bands has revealed temporal variations in the positions and source sizes that do not fit into the standard picture of type III solar radio bursts, and require a better understanding of radio-wave transport. In this pap…
▽ More
Recent developments in astronomical radio telescopes opened new opportunities in imaging and spectroscopy of solar radio bursts at sub-second timescales. Imaging in narrow frequency bands has revealed temporal variations in the positions and source sizes that do not fit into the standard picture of type III solar radio bursts, and require a better understanding of radio-wave transport. In this paper, we utilise 3D Monte Carlo ray-tracing simulations that account for the anisotropic density turbulence in the inhomogeneous solar corona to quantitatively explain the image dynamics at the fundamental (near plasma frequency) and harmonic (double) plasma emissions observed at \sim 32~MHz. Comparing the simulations with observations, we find that anisotropic scattering from an instantaneous emission point source can account for the observed time profiles, centroid locations, and source sizes of the fundamental component of type III radio bursts (generated where f_{pe} \approx 32~MHz). The best agreement with observations is achieved when the ratio of the perpendicular to the parallel component of the wave vector of anisotropic density turbulence is around 0.25. Harmonic emission sources observed at the same frequency (\sim 32~MHz, but generated where f_{pe} \approx 16~MHz) have apparent sizes comparable to those produced by the fundamental emission, but demonstrate a much slower temporal evolution. The simulations of radio-wave propagation make it possible to quantitatively explain the variations of apparent source sizes and positions at sub-second time-scales both for the fundamental and harmonic emissions, and can be used as a diagnostic tool for the plasma turbulence in the upper corona.
△ Less
Submitted 18 November, 2020; v1 submitted 17 October, 2020;
originally announced October 2020.
-
Forward modelling of particle acceleration and transport in an individual solar flare
Authors:
Mykola Gordovskyy,
Philippa K. Browning,
Satoshi Inoue,
Eduard P. Kontar,
Kanya Kusano,
Grigory E. Vekstein
Abstract:
The aim of this study is to generate maps of the hard X-ray emission produced by energetic electrons in a solar flare and compare them with observations. The ultimate goal is to test the viability of the combined MHD/test-particle approach for data-driven modelling of active events in the solar corona and their impact on the heliosphere. Based on an MHD model of X-class solar flare observed on the…
▽ More
The aim of this study is to generate maps of the hard X-ray emission produced by energetic electrons in a solar flare and compare them with observations. The ultimate goal is to test the viability of the combined MHD/test-particle approach for data-driven modelling of active events in the solar corona and their impact on the heliosphere. Based on an MHD model of X-class solar flare observed on the 8th of September 2017, we calculate trajectories of a large number of electrons and protons using the relativistic guiding-centre approach. Using the obtained particle trajectories, we deduce the spatial and energy distributions of energetic electrons and protons, and calculate bremsstrahlung hard X-ray emission using the 'thin target' approximation. Our approach predicts some key characteristics of energetic particles in the considered flare, including the size and location of the acceleration region, energetic particle trajectories and energy spectra. Most importantly, the hard X-ray bremsstrahlung intensity maps predicted by the model are in a good agreement with those observed by RHESSI. Furthermore, the locations of proton and electron precipitation appear to be close to the sources of helioseismic response detected in this flare. Therefore, the adopted approach can be used for observationally-driven modelling of individual solar flares, including manifestations of energetic particles in the corona, as well as inner heliosphere.
△ Less
Submitted 21 September, 2020;
originally announced September 2020.
-
Predicting the time variation of radio emission from MHD simulations of a flaring T-Tauri star
Authors:
C. O. G. Waterfall,
P. K. Browning,
G. A. Fuller,
M. Gordovskyy,
S. Orlando,
F. Reale
Abstract:
We model the time dependent radio emission from a disk accretion event in a T-Tauri star using 3D, ideal magnetohydrodynamic simulations combined with a gyrosynchrotron emission and radiative transfer model. We predict for the first time, the multi-frequency (1$-$1000 GHz) intensity and circular polarisation from a flaring T-Tauri star. A flux tube, connecting the star with its circumstellar disk,…
▽ More
We model the time dependent radio emission from a disk accretion event in a T-Tauri star using 3D, ideal magnetohydrodynamic simulations combined with a gyrosynchrotron emission and radiative transfer model. We predict for the first time, the multi-frequency (1$-$1000 GHz) intensity and circular polarisation from a flaring T-Tauri star. A flux tube, connecting the star with its circumstellar disk, is populated with a distribution of non-thermal electrons which is allowed to decay exponentially after a heating event in the disk and the system is allowed to evolve. The energy distribution of the electrons, as well as the non-thermal power law index and loss rate, are varied to see their effect on the overall flux. Spectra are generated from different lines of sight, giving different views of the flux tube and disk. The peak flux typically occurs around 20$-$30 GHz and the radio luminosity is consistent with that observed from T-Tauri stars. For all simulations, the peak flux is found to decrease and move to lower frequencies with elapsing time. The frequency-dependent circular polarisation can reach 10$-$30$\%$ but has a complex structure which evolves as the flare evolves. Our models show that observations of the evolution of the spectrum and its polarisation can provide important constraints on physical properties of the flaring environment and associated accretion event.
△ Less
Submitted 9 June, 2020;
originally announced June 2020.
-
Using the Stokes V widths of Fe I lines for diagnostics of the intrinsic solar photospheric magnetic field
Authors:
M. Gordovskyy,
S. Shelyag,
P. K. Browning,
V. G. Lozitsky
Abstract:
The goal of this study is to explore a novel method for the solar photospheric magnetic field diagnostics using Stokes V widths of different magnetosensitive Fe~I spectral lines. We calculate Stokes I and V profiles of several Fe I lines based on a one-dimensional photospheric model VAL C using the NICOLE radiative transfer code. These profiles are used to produce calibration curves linking the in…
▽ More
The goal of this study is to explore a novel method for the solar photospheric magnetic field diagnostics using Stokes V widths of different magnetosensitive Fe~I spectral lines. We calculate Stokes I and V profiles of several Fe I lines based on a one-dimensional photospheric model VAL C using the NICOLE radiative transfer code. These profiles are used to produce calibration curves linking the intrinsic magnetic field values with the widths of blue peaks of Stokes V profiles. The obtained calibration curves are then tested using the Stokes profiles calculated for more realistic photospheric models based on MHD models of magneto-convection. It is shown that the developed Stokes V widths (SVW) method can be used with various optical and near-infrared lines. Out of six lines considered in this study, FeI 6301 line appears to be the most effective: it is sensitive to fields over ~200G and does not show any saturation up to ~2kG. Other lines considered can also be used for the photospheric field diagnostics with this method, however, only in narrower field value ranges, typically from about 100G to 700-1000G. The developed method can be a useful alternative to the classical magnetic line ratio method, particularly when the choice of lines is limited.
△ Less
Submitted 6 December, 2019;
originally announced December 2019.
-
Anisotropic Radio-Wave Scattering and the Interpretation of Solar Radio Emission Observations
Authors:
Eduard P. Kontar,
Xingyao Chen,
Nicolina Chrysaphi,
Natasha L. S. Jeffrey,
A. Gordon Emslie,
Vratislav Krupar,
Milan Maksimovic,
Mykola Gordovskyy,
Philippa K. Browning
Abstract:
The observed properties (i.e., source size, source position, time duration, decay time) of solar radio emission produced through plasma processes near the local plasma frequency, and hence the interpretation of solar radio bursts, are strongly influenced by propagation effects in the inhomogeneous turbulent solar corona. In this work, a 3D stochastic description of the propagation process is prese…
▽ More
The observed properties (i.e., source size, source position, time duration, decay time) of solar radio emission produced through plasma processes near the local plasma frequency, and hence the interpretation of solar radio bursts, are strongly influenced by propagation effects in the inhomogeneous turbulent solar corona. In this work, a 3D stochastic description of the propagation process is presented, based on the Fokker-Planck and Langevin equations of radio-wave transport in a medium containing anisotropic electron density fluctuations. Using a numerical treatment based on this model, we investigate the characteristic source sizes and burst decay times for Type III solar radio bursts. Comparison of the simulations with the observations of solar radio bursts shows that predominantly perpendicular density fluctuations in the solar corona are required, with an anisotropy factor $\sim 0.3$ for sources observed at around 30~MHz. The simulations also demonstrate that the photons are isotropized near the region of primary emission, but the waves are then focused by large-scale refraction, leading to plasma radio emission directivity that is characterized by a half-width-half-maximum of about 40~degrees near 30~MHz. The results are applicable to various solar radio bursts produced via plasma emission.
△ Less
Submitted 4 September, 2019; v1 submitted 1 September, 2019;
originally announced September 2019.
-
Forced magnetic reconnection and plasmoid coalescence: I - MHD Simulations
Authors:
Max Potter,
Philippa Browning,
Mykola Gordovskyy
Abstract:
Forced magnetic reconnection, a reconnection event triggered by external perturbation, should be ubiquitous in the solar corona. Energy released during such cases can be much greater than that which was introduced by the perturbation. It is unclear how the properties of the external perturbation and the initial current sheet affect the reconnection region properties, and thereby the reconnection d…
▽ More
Forced magnetic reconnection, a reconnection event triggered by external perturbation, should be ubiquitous in the solar corona. Energy released during such cases can be much greater than that which was introduced by the perturbation. It is unclear how the properties of the external perturbation and the initial current sheet affect the reconnection region properties, and thereby the reconnection dynamics and energy release profile. We investigate the effect of the form of the external perturbation and initial current sheet on the evolution of the reconnection region and the energy release process. Chiefly we explore the non-linear interactions between multiple, simultaneous perturbations, which represent more realistic scenarios. Future work will use these results in test particle simulations to investigate particle acceleration over multiple reconnection events. Simulations are performed using Lare2d, a 2.5D Lagrangian-remap solver for the visco-resistive MHD equations. The model of forced reconnection is extended to include superpositions of sinusoidal driving disturbances, including localised Gaussian perturbations. A transient perturbation is applied to the boundaries of a region containing a force-free current sheet. The simulation domain is sufficiently wide to allow multiple magnetic islands to form and coalesce. Island coalescence contributes significantly to energy release and involves rapid reconnection. Long wavelength modes in perturbations dominate the evolution, without the presence of which reconnection is either slow, as in the case of short wavelength modes, or the initial current sheet remains stable, as in the case of noise perturbations. Multiple perturbations combine in a highly non-linear manner: reconnection is typically faster than when either disturbance is applied individually, with multiple low-energy events contributing to the same total energy release.
△ Less
Submitted 8 January, 2019;
originally announced January 2019.
-
Combining MHD and kinetic modelling of solar flares
Authors:
M. Gordovskyy,
P. K. Browning,
R. F. Pinto
Abstract:
Solar flares are explosive events in the solar corona, representing fast conversion of magnetic energy into thermal and kinetic energy, and hence radiation, due to magnetic reconnection. Modelling is essential for understanding and predicting these events. However, self-consistent modelling is extremely difficult due to the vast spatial and temporal scale separation between processes involving the…
▽ More
Solar flares are explosive events in the solar corona, representing fast conversion of magnetic energy into thermal and kinetic energy, and hence radiation, due to magnetic reconnection. Modelling is essential for understanding and predicting these events. However, self-consistent modelling is extremely difficult due to the vast spatial and temporal scale separation between processes involving thermal plasma (normally considered using magnetohydrodynamic (MHD) approach) and non-thermal plasma (requiring a kinetic approach). In this mini-review we consider different approaches aimed at bridging the gap between fluid and kinetic modelling of solar flares. Two types of approaches are discussed: combined MHD/test-particle (MHDTP) models, which can be used for modelling the flaring corona with relatively small numbers of energetic particles, and hybrid fluid-kinetic methods, which can be used for modelling stronger events with higher numbers of energetic particles. Two specific examples are discussed in more detail: MHDTP models of magnetic reconnection and particle acceleration in kink-unstable twisted coronal loops, and a novel reduced-kinetic model of particle transport.
△ Less
Submitted 26 September, 2018; v1 submitted 15 September, 2018;
originally announced September 2018.
-
Analysis of unresolved photospheric magnetic field structure using Fe I 6301 and 6302 lines
Authors:
M. Gordovskyy.,
S. Shelyag,
P. K. Browning,
V. G. Lozitsky
Abstract:
Early magnetographic observations indicated that magnetic field in the solar photosphere has unresolved small-scale structure. Near-infrared and optical data with extremely high spatial resolution show that these structures have scales of few tens of kilometres, which are not resolved in the majority of solar observations. The goal of this study is to establish the effect of unresolved photospheri…
▽ More
Early magnetographic observations indicated that magnetic field in the solar photosphere has unresolved small-scale structure. Near-infrared and optical data with extremely high spatial resolution show that these structures have scales of few tens of kilometres, which are not resolved in the majority of solar observations. The goal of this study is to establish the effect of unresolved photospheric magnetic field structure on Stokes profiles observed with relatively low spatial resolution. Ultimately, we aim to develop methods for fast estimation of the photospheric magnetic filling factor and line-of-sight gradient of the photospheric magnetic field, which can be applied to large observational data sets. We exploit 3D MHD models of magneto-convection developed using MURAM code. Corresponding profiles of Fe I 6301.5 and 6302.5 $\mathrmÅ$ spectral lines are calculated using NICOLE radiative transfer code. The resulting I and V Stokes [x,y,$λ$] cubes with reduced spatial resolution of 150 km are used to calculate magnetic field values as they would be obtained in observations with Hinode/SOT or SDO/HMI. Three different methods of the magnetic filling factor estimation are considered: the magnetic line ratio method, Stokes V width method and a simple statistical method. We find that the statistical method and the Stokes V width method are sufficiently reliable for fast filling factor estimations. Furthermore, we find that Stokes $I\pm V$ bisector splitting gradient can be used for fast estimation of line-of-sight gradient of the photospheric magnetic field.
△ Less
Submitted 21 August, 2018;
originally announced August 2018.
-
Polarisation of microwave emission from reconnecting twisted coronal loops
Authors:
Mykola Gordovskyy,
Philippa Browning,
Eduard Kontar
Abstract:
Magnetic reconnection and particle acceleration due to the kink instability in twisted coronal loops can be a viable scenario for confined solar flares. Detailed investigation of this phenomenon requires reliable methods for observational detection of magnetic twist in solar flares, which may not be possible solely through extreme UV and soft X-ray thermal emission. Polarisation of microwave emiss…
▽ More
Magnetic reconnection and particle acceleration due to the kink instability in twisted coronal loops can be a viable scenario for confined solar flares. Detailed investigation of this phenomenon requires reliable methods for observational detection of magnetic twist in solar flares, which may not be possible solely through extreme UV and soft X-ray thermal emission. Polarisation of microwave emission in flaring loops can be used as one of the detection criteria. The aim of this study is to investigate the effect of magnetic twist in flaring coronal loops on the polarisation of gyro-synchrotron microwave (GSMW) emission, and determine whether it could provide a means for magnetic twist detection. We use time-dependent magnetohydrodynamic and test-particle models developed using LARE3D and GCA codes to investigate twisted coronal loops relaxing following the kink-instability. Synthetic GSMW emission maps (I and V Stokes components) are calculated using GX simulator. It is found that flaring twisted coronal loops produce GSMW radiation with a gradient of circular polarisation across the loop. However, these patterns may be visible only for a relatively short period of time due to fast magnetic reconfiguration after the instability. Their visibility also depends on the orientation and position of the loop on solar disk. Typically, it would be difficult to see these characteristic polarisation pattern in a twisted loop seen from the top (close to the centre of the solar disk), but easier in a twisted loop seen from the side (i.e. observed very close to the limb).
△ Less
Submitted 29 March, 2017; v1 submitted 7 November, 2016;
originally announced November 2016.
-
Reduced drift-kinetics with thermal velocity distribution across magnetic field
Authors:
Mykola Gordovskyy,
Philippa Browning
Abstract:
The goal of this study is to develop an approximate self-consistent description of particle motion in strongly magnetised solar corona. We derive a set of reduced drift-kinetic equations based on the assumption that the gyro-velocity distribution is Maxwellian. The equations are tested using simple 1D models.
The goal of this study is to develop an approximate self-consistent description of particle motion in strongly magnetised solar corona. We derive a set of reduced drift-kinetic equations based on the assumption that the gyro-velocity distribution is Maxwellian. The equations are tested using simple 1D models.
△ Less
Submitted 31 January, 2016;
originally announced February 2016.
-
Plasma motions and non-thermal line broadening in flaring twisted coronal loops
Authors:
Mykola Gordovskyy,
Eduard Kontar,
Philippa Browning
Abstract:
Observation of coronal EUV spectral lines sensitive to different temperatures offers an opportunity to evaluate the thermal structure and flows in flaring atmospheres. This can be used to estimate the partitioning between the thermal and kinetic energies. Our aim is to model large-scale (50-10000km) velocity distributions in order to interpret non-thermal broadening of different spectral EUV lines…
▽ More
Observation of coronal EUV spectral lines sensitive to different temperatures offers an opportunity to evaluate the thermal structure and flows in flaring atmospheres. This can be used to estimate the partitioning between the thermal and kinetic energies. Our aim is to model large-scale (50-10000km) velocity distributions in order to interpret non-thermal broadening of different spectral EUV lines. The developed models allow us to understand the origin of the observed spectral line shifts and broadening, and link these features to particular physical effects in flaring atmospheres. We use ideal MHD to derive twisted magnetic fluxtube configurations in a stratified atmosphere. The evolution of these fluxtubes is followed using resistive MHD, with anomalous resistivity depending on the local density and temperature. The model also takes into account the thermal conduction and radiative losses. The model allows us to evaluate average velocities and velocity dispersions, which would be interpreted as `non-thermal' velocities in observations, at different temperatures for different parts of the models. Our models show qualitative and quantitative agreement with observations. The line-of-sight (LOS) velocity dispersions demonstrate substantial correlation with the temperature, increasing from about 20-30 km/s around 1 MK to about 200-400 km/s near 10-20 MK. The average LOS velocities also correlate with velocity dispersions, although they demonstrate a very strong scattering, compared to observations. We also note that near foot-points the velocity dispersions across the magnetic field are systematically lower than those along the field. We conclude, that the correlation between the flow velocities, velocity dispersions and temperatures are likely to indicate that the same heating mechanism is responsible for heating the plasma, its turbulisation and expansion/evaporation.
△ Less
Submitted 1 February, 2016; v1 submitted 26 August, 2015;
originally announced August 2015.
-
Energy release in driven twisted coronal loops
Authors:
M. R. Bareford,
M. Gordovskyy,
P. K. Browning,
A. W. Hood
Abstract:
In the present study we investigate magnetic reconnection in twisted magnetic fluxtubes with different initial configurations. In all considered cases, energy release is triggered by the ideal kink instability, which is itself the result of applying footpoint rotation to an initially potential field. The main goal of this work is to establish the influence of the field topology and various thermod…
▽ More
In the present study we investigate magnetic reconnection in twisted magnetic fluxtubes with different initial configurations. In all considered cases, energy release is triggered by the ideal kink instability, which is itself the result of applying footpoint rotation to an initially potential field. The main goal of this work is to establish the influence of the field topology and various thermodynamic effects on the energy release process. Specifically, we investigate convergence of the magnetic field at the loop footpoints, atmospheric stratification, as well as thermal conduction. In all cases, the application of vortical driving at the footpoints of an initally potential field leads to an internal kink instability. With the exception of the curved loop with high footpoint convergence, the global geometry of the loop change little during the simulation. Footpoint convergence, curvature and atmospheric structure clearly influences the rapidity with which a loop achieves instability as well as the size of the subsequent energy release. Footpoint convergence has a stabilising influence and thus the loop requires more energy for instability, which means that the subsequent relaxation has a larger heating effect. Large-scale curvature has the opposite result: less energy is needed for instability and so the amount of energy released from the field is reduced. Introducing a stratified atmosphere gives rise to decaying wave phenomena during the driving phase, and also results in a loop that is less stable.
△ Less
Submitted 3 June, 2015;
originally announced June 2015.
-
Thermal and non-thermal emission from reconnecting twisted coronal loops
Authors:
R. F. Pinto,
M. Gordovskyy,
P. K. Browning,
N. Vilmer
Abstract:
Twisted magnetic fields should be ubiquitous in flare-producing active regions where the magnetic fields are strongly non-potential. It has been shown that reconnection in helical magnetic coronal loops results in plasma heating and particle acceleration distributed within a large volume, including the lower coronal and chromospheric sections of the loops. This scenario can be an alternative to th…
▽ More
Twisted magnetic fields should be ubiquitous in flare-producing active regions where the magnetic fields are strongly non-potential. It has been shown that reconnection in helical magnetic coronal loops results in plasma heating and particle acceleration distributed within a large volume, including the lower coronal and chromospheric sections of the loops. This scenario can be an alternative to the standard flare model, where particles are accelerated only in a small volume located in the upper corona. We use a combination of MHD simulations and test-particle methods, which describe the development of kink instability and magnetic reconnection in twisted coronal loops using resistive compressible MHD, and incorporate atmospheric stratification and large-scale loop curvature. The resulting distributions of hot plasma let us estimate thermal X-ray emission intensities. The electric and magnetic fields obtained are used to calculate electron trajectories using the guiding-centre approximation. These trajectories combined with the MHD plasma density distributions let us deduce synthetic HXR bremsstrahlung intensities. Our simulations emphasise that the geometry of the emission patterns produced by hot plasma in flaring twisted coronal loops can differ from the actual geometry of the underlying magnetic fields. The twist angles revealed by the emission threads (SXR) are consistently lower than the field-line twist present at the onset of the kink-instability. HXR emission due to the interaction of energetic electrons with the stratified background are concentrated at the loop foot-points in these simulations, even though the electrons are accelerated everywhere within the coronal volume of the loop. The maximum of HXR emission consistently precedes that of SXR emission, with the HXR light-curve being approximately proportional to the temporal derivative of the SXR light-curve.
△ Less
Submitted 11 December, 2015; v1 submitted 3 June, 2015;
originally announced June 2015.
-
Effect of collisions and magnetic convergence on electron acceleration and transport in reconnecting twisted solar flare loops
Authors:
M. Gordovskyy,
P. K. Browning,
E. P. Kontar,
N. H. Bian
Abstract:
We study a model of particle acceleration coupled with an MHD model of magnetic reconnection in unstable twisted coronal loops. The kink instability leads to the formation of helical currents with strong parallel electric fields resulting in electron acceleration. The motion of electrons in the electric and magnetic fields of the reconnecting loop is investigated using a test-particle approach tak…
▽ More
We study a model of particle acceleration coupled with an MHD model of magnetic reconnection in unstable twisted coronal loops. The kink instability leads to the formation of helical currents with strong parallel electric fields resulting in electron acceleration. The motion of electrons in the electric and magnetic fields of the reconnecting loop is investigated using a test-particle approach taking into account collisional scattering. We discuss the effects of Coulomb collisions and magnetic convergence near loop footpoints on the spatial distribution and energy spectra of high-energy electron populations and possible implications on the hard X-ray emission in solar flares.
△ Less
Submitted 26 January, 2015;
originally announced January 2015.
-
Particle acceleration and transport in reconnecting twisted loops in a stratified atmosphere
Authors:
Mykola Gordovskyy,
Philippa Browning,
Eduard Kontar,
Nicolas Bian
Abstract:
Twisted coronal loops should be ubiquitous in the solar corona. Twisted magnetic fields contain excess magnetic energy, which can be released during magnetic reconnection, causing solar flares. The aim of this work is to investigate magnetic reconnection, and particle acceleration and transport in kink-unstable twisted coronal loops, with a focus on the effects of resistivity, loop geometry and at…
▽ More
Twisted coronal loops should be ubiquitous in the solar corona. Twisted magnetic fields contain excess magnetic energy, which can be released during magnetic reconnection, causing solar flares. The aim of this work is to investigate magnetic reconnection, and particle acceleration and transport in kink-unstable twisted coronal loops, with a focus on the effects of resistivity, loop geometry and atmospheric stratification. Another aim is to perform forward-modelling of bremsstrahlung emission and determine the structure of hard X-ray sources. We use a combination of magnetohydrodynamic (MHD) and test-particle methods. First, the evolution of the kinking coronal loop is considered using resistive MHD model, incorporating atmospheric stratification and loop curvature. Then, the obtained electric and magnetic fields and density distributions are used to calculate electron and proton trajectories using a guiding-centre approximation, taking into account Coulomb collisions. It is shown that electric fields in twisted coronal loops can effectively accelerate protons and electrons to energies up to 10 MeV. High-energy particles have hard, nearly power-law energy spectra. The volume occupied by high-energy particles demonstrates radial expansion, which results in the expansion of the visible hard X-ray loop and a gradual increase in hard X-ray footpoint area. Synthesised hard X-ray emission reveals strong footpoint sources and the extended coronal source, whose intensity strongly depends on the coronal loop density.
△ Less
Submitted 26 January, 2015;
originally announced January 2015.
-
Observations of Unresolved Photospheric Magnetic Fields in Solar Flares Using Fe I and Cr I Lines
Authors:
Mykola Gordovskyy,
Vsevolod G. Lozitsky
Abstract:
The structure of the photospheric magnetic field during solar flares is examined using echelle spectropolarimetric observations. The study is based on several Fe I and Cr I lines observed at locations corresponding to brightest H$α$ emission during thermal phase of flares. The analysis is performed by comparing magnetic field values deduced from lines with different magnetic sensitivities, as well…
▽ More
The structure of the photospheric magnetic field during solar flares is examined using echelle spectropolarimetric observations. The study is based on several Fe I and Cr I lines observed at locations corresponding to brightest H$α$ emission during thermal phase of flares. The analysis is performed by comparing magnetic field values deduced from lines with different magnetic sensitivities, as well as by examining the fine structure of $I\pm V$ Stokes profiles splitting. It is shown that the field has at least two components, with stronger unresolved flux tubes embedded in weaker ambient field. Based on a two-component magnetic field model, we compare observed and synthetic line profiles and show that the field strength in small-scale flux tubes is about $2-3$ kG. Furthermore, we find that the small-scale flux tubes are associated with flare emission, which may have implications for flare phenomenology.
△ Less
Submitted 15 April, 2014;
originally announced April 2014.
-
Two-fluid simulations of driven reconnection in the Mega-Ampere Spherical Tokamak
Authors:
A. Stanier,
P. Browning,
M. Gordovskyy,
K. G. McClements,
M. P. Gryaznevich,
V. S. Lukin
Abstract:
In the merging-compression method of plasma start-up, two flux-ropes with parallel toroidal current are formed around in-vessel poloidal field coils, before merging to form a spherical tokamak plasma. This start-up method, used in the Mega-Ampere Spherical Tokamak (MAST), is studied as a high Lundquist number and low plasma-beta magnetic reconnection experiment.
In this paper, 2D fluid simulatio…
▽ More
In the merging-compression method of plasma start-up, two flux-ropes with parallel toroidal current are formed around in-vessel poloidal field coils, before merging to form a spherical tokamak plasma. This start-up method, used in the Mega-Ampere Spherical Tokamak (MAST), is studied as a high Lundquist number and low plasma-beta magnetic reconnection experiment.
In this paper, 2D fluid simulations are presented of this merging process in order to understand the underlying physics, and better interpret the experimental data. These simulations examine the individual and combined effects of tight-aspect ratio geometry and two-fluid physics on the merging. The ideal self-driven flux-rope dynamics are coupled to the diffusion layer physics, resulting in a large range of phenomena. For resistive MHD simulations, the flux-ropes enter the sloshing regime for normalised resistivity eta < 1E-5. In Hall-MHD three regimes are found for the qualitative behaviour of the current sheet, depending on the ratio of the current sheet width to the ion-sound radius. These are a stable collisional regime, an open X-point regime, and an intermediate regime that is highly unstable to tearing-type instabilities.
In toroidal axisymmetric geometry, the final state after merging is a MAST-like spherical tokamak with nested flux-surfaces. It is also shown that the evolution of simulated 1D radial density profiles closely resembles the Thomson scattering electron density measurements in MAST. An intuitive explanation for the origin of the measured density structures is proposed, based upon the results of the toroidal Hall-MHD simulations.
△ Less
Submitted 13 August, 2013;
originally announced August 2013.
