Showing 1–50 of 165 results for author: Kontar, E

AtLAST Science Overview Report

Authors: Mark Booth, Pamela Klaassen, Claudia Cicone, Tony Mroczkowski, Martin A. Cordiner, Luca Di Mascolo, Doug Johnstone, Eelco van Kampen, Minju M. Lee, Daizhong Liu, John Orlowski-Scherer, Amélie Saintonge, Matthew W. L. Smith, Alexander Thelen, Sven Wedemeyer, Kazunori Akiyama, Stefano Andreon, Doris Arzoumanian, Tom J. L. C. Bakx, Caroline Bot, Geoffrey Bower, Roman Brajša, Chian-Chou Chen, Elisabete da Cunha, David Eden , et al. (59 additional authors not shown)

Abstract: Submillimeter and millimeter wavelengths provide a unique view of the Universe, from the gas and dust that fills and surrounds galaxies to the chromosphere of our own Sun. Current single-dish facilities have presented a tantalising view of the brightest (sub-)mm sources, and interferometers have provided the exquisite resolution necessary to analyse the details in small fields, but there are still… ▽ More Submillimeter and millimeter wavelengths provide a unique view of the Universe, from the gas and dust that fills and surrounds galaxies to the chromosphere of our own Sun. Current single-dish facilities have presented a tantalising view of the brightest (sub-)mm sources, and interferometers have provided the exquisite resolution necessary to analyse the details in small fields, but there are still many open questions that cannot be answered with current facilities. In this report we summarise the science that is guiding the design of the Atacama Large Aperture Submillimeter Telescope (AtLAST). We demonstrate how tranformational advances in topics including star formation in high redshift galaxies, the diffuse circumgalactic medium, Galactic ecology, cometary compositions and solar flares motivate the need for a 50m, single-dish telescope with a 1-2 degree field of view and a new generation of highly multiplexed continuum and spectral cameras. AtLAST will have the resolution to drastically lower the confusion limit compared to current single-dish facilities, whilst also being able to rapidly map large areas of the sky and detect extended, diffuse structures. Its high sensitivity and large field of view will open up the field of submillimeter transient science by increasing the probability of serendipitous detections. Finally, the science cases listed here motivate the need for a highly flexible operations model capable of short observations of individual targets, large surveys, monitoring programmes, target of opportunity observations and coordinated observations with other observatories. AtLAST aims to be a sustainable, upgradeable, multipurpose facility that will deliver orders of magnitude increases in sensitivity and mapping speeds over current and planned submillimeter observatories. △ Less

Submitted 1 July, 2024; originally announced July 2024.

Comments: 47 pages, 12 figures. For further details on AtLAST see https://atlast.uio.no

First determination of the angular dependence of rise and decay times of solar radio bursts using multi-spacecraft observations

Authors: Nicolina Chrysaphi, Milan Maksimovic, Eduard P. Kontar, Antonio Vecchio, Xingyao Chen, Aikaterini Pesini

Abstract: Radio photons interact with anisotropic density fluctuations in the heliosphere, which can alter their trajectory and influence properties deduced from observations. This is particularly evident in solar radio observations, where anisotropic scattering leads to highly-directional radio emissions. Consequently, observers at varying locations will measure different properties, including different so… ▽ More Radio photons interact with anisotropic density fluctuations in the heliosphere, which can alter their trajectory and influence properties deduced from observations. This is particularly evident in solar radio observations, where anisotropic scattering leads to highly-directional radio emissions. Consequently, observers at varying locations will measure different properties, including different source sizes, source positions, and intensities. However, it is not known if measurements of the decay time of solar radio bursts are also affected by the observer's position. Decay times are dominated by scattering effects, and so are frequently used as proxies of the level of density fluctuations in the heliosphere, making the identification of any location-related dependence crucial. We combine multi-vantage observations of interplanetary Type III bursts from four non-collinear, angularly-separated spacecraft with simulations, to investigate the dependence of both the decay- and rise-time measurements on the separation of the observer from the source. We propose a function to characterise the entire time profile of radio signals, allowing for the simultaneous estimation of the peak flux, decay time, and rise time, while demonstrating that the rise phase of radio bursts has a non-constant, non-exponential growth rate. We determine that the decay and rise times are independent of the observer's position, identifying them as the only properties to remain unaffected, thus not requiring corrections for the observer's location. Moreover, we examine the ratio between the rise and decay times, finding that it does not depend on the frequency. Therefore, we provide the first evidence that the rise phase is also significantly impacted by scattering effects, adding to our understanding of the plasma emission process. △ Less

Submitted 1 April, 2024; originally announced April 2024.

Comments: Accepted for publication in Astronomy & Astrophysics; 15 pages, 11 figures

Plasma motions and compressive wave energetics in the solar corona and solar wind from radio wave scattering observations

Authors: Francesco Azzollini, A. Gordon Emslie, Daniel L. Clarkson, Nicolina Chrysaphi, Eduard P. Kontar

Abstract: Radio signals propagating via the solar corona and solar wind are significantly affected by compressive waves, impacting properties of solar bursts as well as sources viewed through the turbulent solar atmosphere. While static fluctuations scatter radio waves elastically, moving, turbulent or oscillating density irregularities act to broaden the frequency of the scattered waves. Using a new anisot… ▽ More Radio signals propagating via the solar corona and solar wind are significantly affected by compressive waves, impacting properties of solar bursts as well as sources viewed through the turbulent solar atmosphere. While static fluctuations scatter radio waves elastically, moving, turbulent or oscillating density irregularities act to broaden the frequency of the scattered waves. Using a new anisotropic density fluctuation model from the kinetic scattering theory for solar radio bursts, we deduce the plasma velocities required to explain observations of spacecraft signal frequency broadening. The inferred velocities are consistent with motions that are dominated by the solar wind at distances $\gtrsim 10$ $R_\odot$, but the levels of frequency broadening for $\lesssim 10$ $R_\odot$ require additional radial speeds $\sim (100-300)$ km s$^{-1}$ and/or transverse speeds $\sim (20-70)$ km s$^{-1}$. The inferred radial velocities also appear consistent with the sound or proton thermal speeds, while the speeds perpendicular to the radial direction are consistent with non-thermal motions measured via coronal Doppler-line broadening, interpreted as Alfvénic fluctuations. Landau damping of parallel propagating ion-sound (slow MHD) waves allow an estimate of the proton heating rate. The energy deposition rates due to ion-sound wave damping peak at a heliocentric distance of $\sim(1-3)$ $R_\odot$ are comparable to the rates available from a turbulent cascade of Alfvénic waves at large scales, suggesting a coherent picture of energy transfer, via the cascade or/and parametric decay of Alfvén waves to the small scales where heating takes place. △ Less

Submitted 22 April, 2024; v1 submitted 19 March, 2024; originally announced March 2024.

Comments: 27 pages, 5 figures, submitted to The Astrophysical Journal

A Modelling Investigation for Solar Flare X-ray Stereoscopy with Solar Orbiter/STIX and Earth Orbiting Missions

Authors: Natasha L. S. Jeffrey, Säm Krucker, Morgan Stores, Eduard P. Kontar, Pascal Saint-Hilaire, Andrea F. Battaglia, Laura Hayes, Hannah Collier, Astrid Veronig, Yang Su, Srikar Paavan Tadepalli, Fanxiaoyu Xia

Abstract: The Spectrometer/Telescope for Imaging X-rays (STIX) on board Solar Orbiter (SolO) provides a unique opportunity to systematically perform stereoscopic X-ray observations of solar flares with current and upcoming X-ray missions at Earth. These observations will produce the first reliable measurements of hard X-ray (HXR) directivity in decades, providing a new diagnostic of the flare-accelerated el… ▽ More The Spectrometer/Telescope for Imaging X-rays (STIX) on board Solar Orbiter (SolO) provides a unique opportunity to systematically perform stereoscopic X-ray observations of solar flares with current and upcoming X-ray missions at Earth. These observations will produce the first reliable measurements of hard X-ray (HXR) directivity in decades, providing a new diagnostic of the flare-accelerated electron angular distribution and helping to constrain the processes that accelerate electrons in flares. However, such observations must be compared to modelling, taking into account electron and X-ray transport effects and realistic plasma conditions, all of which can change the properties of the measured HXR directivity. Here, we show how HXR directivity, defined as the ratio of X-ray spectra at different spacecraft viewing angles, varies with different electron and flare properties (e.g., electron angular distribution, highest energy electrons, and magnetic configuration), and how modelling can be used to extract these typically unknown properties from the data. Lastly, we present a preliminary HXR directivity analysis of two flares, observed by the Fermi Gamma-ray Burst Monitor (GBM) and SolO/STIX, demonstrating the feasibility and challenges associated with such observations, and how HXR directivity can be extracted by comparison with the modelling presented here. △ Less

Submitted 29 January, 2024; originally announced January 2024.

Comments: Accepted for publication in ApJ (January 2024)

Probing Turbulent Scattering Effects on Suprathermal Electrons in the Solar Wind: Modeling, Observations and Implications

Authors: Arnaud Zaslavsky, Justin C. Kasper, Eduard P. Kontar, Davin E. Larson, Milan Maksimovic, José M. D. C. Marques, Georgios Nicolaou, Christopher J. Owen, Orlando Romeo, Phyllis L. Whittlesey

Abstract: This study explores the impact of a turbulent scattering mechanism, akin to those influencing solar and galactic cosmic rays propagating in the interplanetary medium, on the population of suprathermal electrons in the solar wind. We employ a Fokker-Planck equation to model the radial evolution of electron pitch angle distributions under the action of magnetic focusing, which moves the electrons aw… ▽ More This study explores the impact of a turbulent scattering mechanism, akin to those influencing solar and galactic cosmic rays propagating in the interplanetary medium, on the population of suprathermal electrons in the solar wind. We employ a Fokker-Planck equation to model the radial evolution of electron pitch angle distributions under the action of magnetic focusing, which moves the electrons away from isotropy, and of a diffusion process that tends to bring them back to it. We compare the steady-state solutions of this Fokker-Planck equation with data obtained from the Solar Orbiter and Parker Solar Probe missions and find a remarkable agreement, varying the turbulent mean free path as the sole free parameter in our model. The obtained mean free paths are of the order of the astronomical unit, and display weak dependence on electron energy within the $100$ eV to $1$ keV range. This value is notably lower than Coulomb collision estimates but aligns well with observed mean free paths of low-rigidity solar energetic particles events. The strong agreement between our model and observations leads us to conclude that the hypothesis of turbulent scattering at work on electrons at all heliospheric distances is justified. We discuss several implications, notably the existence of a low Knudsen number region at large distances from the Sun, which offers a natural explanation for the presence of an isotropic ``halo'' component at all distances from the Sun -- electrons being isotropized in this distant region before travelling back into the inner part of the interplanetary medium. △ Less

Submitted 8 January, 2024; originally announced January 2024.

Comments: 23 pages, 13 figures, 1 table

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.

Comments: 27 pages, 12 figures

Spectroscopic Imaging of the Sun with MeerKAT: Opening a New Frontier in Solar Physics

Authors: Devojyoti Kansabanik, Surajit Mondal, Divya Oberoi, James O. Chibueze, N. E. Engelbrecht, R. D. Strauss, Eduard P. Kontar, Gert J. J. Botha, P. J. Steyn, Amore E. Nel

Abstract: Solar radio emissions provide several unique diagnostics to estimate different physical parameters of the solar corona, which are otherwise simply inaccessible. However, imaging the highly dynamic solar coronal emissions spanning a large range of angular scales at radio wavelengths is extremely challenging. At GHz frequencies, MeerKAT radio telescope is possibly globally the best-suited instrument… ▽ More Solar radio emissions provide several unique diagnostics to estimate different physical parameters of the solar corona, which are otherwise simply inaccessible. However, imaging the highly dynamic solar coronal emissions spanning a large range of angular scales at radio wavelengths is extremely challenging. At GHz frequencies, MeerKAT radio telescope is possibly globally the best-suited instrument at present for providing high-fidelity spectroscopic snapshot solar images. Here, we present the first published spectroscopic images of the Sun made using the observations with MeerKAT in the 880-1670 MHz band. This work demonstrates the high fidelity of spectroscopic snapshot MeerKAT solar images through a comparison with simulated radio images at MeerKAT frequencies. The observed images show extremely good morphological similarities with the simulated images. Our analysis shows that below ~900 MHz MeerKAT images can recover essentially the entire flux density from the large angular scale solar disc. Not surprisingly, at higher frequencies, the missing flux density can be as large as ~50%. However, it can potentially be estimated and corrected for. We believe once solar observation with MeerKAT is commissioned, it will enable a host of novel studies, open the door to a large unexplored phase space with significant discovery potential, and also pave the way for solar science with the upcoming Square Kilometre Array-Mid telescope, for which MeerKAT is a precursor. △ Less

Submitted 17 January, 2024; v1 submitted 4 July, 2023; originally announced July 2023.

Comments: Published at the Astrophysical Journal, 15 pages, 9 figures

Journal ref: ApJ 961 96 (2024)

Source positions of an interplanetary type III radio burst and anisotropic radio-wave scattering

Authors: Xingyao Chen, Eduard P. Kontar, Nicolina Chrysaphi, Peijin Zhang, Vratislav Krupar, Sophie Musset, Milan Maksimovic, Natasha L. S. Jeffrey, Francesco Azzollini, Antonio Vecchio

Abstract: Interplanetary solar radio type III bursts provide the means for remotely studying and tracking energetic electrons propagating in the interplanetary medium. Due to the lack of direct radio source imaging, several methods have been developed to determine the source positions from space-based observations. Moreover, none of the methods consider the propagation effects of anisotropic radio-wave scat… ▽ More Interplanetary solar radio type III bursts provide the means for remotely studying and tracking energetic electrons propagating in the interplanetary medium. Due to the lack of direct radio source imaging, several methods have been developed to determine the source positions from space-based observations. Moreover, none of the methods consider the propagation effects of anisotropic radio-wave scattering, which would strongly distort the trajectory of radio waves, delay their arrival times, and affect their apparent characteristics. We investigate the source positions and directivity of an interplanetary type III burst simultaneously observed by Parker Solar Probe, Solar Orbiter, STEREO, and Wind and compare the results of applying the intensity fit and timing methods with ray-tracing simulations of radio-wave propagation with anisotropic density fluctuations. The simulation calculates the trajectories of the rays, their time profiles at different viewing sites, and the apparent characteristics for various density fluctuation parameters. The results indicate that the observed source positions are displaced away from the locations where emission is produced, and their deduced radial distances are larger than expected from density models. This suggests that the apparent position is affected by anisotropic radio-wave scattering, which leads to an apparent position at a larger heliocentric distance from the Sun. The methods to determine the source positions may underestimate the apparent positions if they do not consider the path of radio-wave propagation and incomplete scattering at a viewing site close to the intrinsic source position. △ Less

Submitted 14 October, 2023; v1 submitted 15 June, 2023; originally announced June 2023.

Comments: 7 pages, 4 figures

Journal ref: A&A 680, A1 (2023)

The need for focused, hard X-ray investigations of the Sun

Authors: Lindsay Glesener, Albert Y. Shih, Amir Caspi, Ryan Milligan, Hugh Hudson, Mitsuo Oka, Juan Camilo Buitrago-Casas, Fan Guo, Dan Ryan, Eduard Kontar, Astrid Veronig, Laura A. Hayes, Andrew Inglis, Leon Golub, Nicole Vilmer, Dale Gary, Hamish Reid, Iain Hannah, Graham S. Kerr, Katharine K. Reeves, Joel Allred, Silvina Guidoni, Sijie Yu, Steven Christe, Sophie Musset , et al. (24 additional authors not shown)

Abstract: Understanding the nature of energetic particles in the solar atmosphere is one of the most important outstanding problems in heliophysics. Flare-accelerated particles compose a huge fraction of the flare energy budget; they have large influences on how events develop; they are an important source of high-energy particles found in the heliosphere; and they are the single most important corollary to… ▽ More Understanding the nature of energetic particles in the solar atmosphere is one of the most important outstanding problems in heliophysics. Flare-accelerated particles compose a huge fraction of the flare energy budget; they have large influences on how events develop; they are an important source of high-energy particles found in the heliosphere; and they are the single most important corollary to other areas of high-energy astrophysics. Despite the importance of this area of study, this topic has in the past decade received only a small fraction of the resources necessary for a full investigation. For example, NASA has selected no new Explorer-class instrument in the past two decades that is capable of examining this topic. The advances that are currently being made in understanding flare-accelerated electrons are largely undertaken with data from EOVSA (NSF), STIX (ESA), and NuSTAR (NASA Astrophysics). This is despite the inclusion in the previous Heliophysics decadal survey of the FOXSI concept as part of the SEE2020 mission, and also despite NASA's having invested heavily in readying the technology for such an instrument via four flights of the FOXSI sounding rocket experiment. Due to that investment, the instrumentation stands ready to implement a hard X-ray mission to investigate flare-accelerated electrons. This white paper describes the scientific motivation for why this venture should be undertaken soon. △ Less

Submitted 8 June, 2023; originally announced June 2023.

Comments: White paper submitted to the Decadal Survey for Solar and Space Physics (Heliophysics) 2024-2033; 15 pages, 5 figures

Journal ref: Bulletin of the AAS, Vol. 55, Issue 3, Whitepaper #129 (14pp); 2023 July 31

The efficiency of electron acceleration during the impulsive phase of a solar flare

Authors: Eduard P. Kontar, A. Gordon Emslie, Galina G. Motorina, Brian R. Dennis

Abstract: Solar flares are known to be prolific electron accelerators, yet identifying the mechanism(s) for such efficient electron acceleration in solar flare (and similar astrophysical settings) presents a major challenge. This is due in part to a lack of observational constraints related to conditions in the primary acceleration region itself. Accelerated electrons with energies above $\sim$20~keV are re… ▽ More Solar flares are known to be prolific electron accelerators, yet identifying the mechanism(s) for such efficient electron acceleration in solar flare (and similar astrophysical settings) presents a major challenge. This is due in part to a lack of observational constraints related to conditions in the primary acceleration region itself. Accelerated electrons with energies above $\sim$20~keV are revealed by hard X-ray (HXR) bremsstrahlung emission, while accelerated electrons with even higher energies manifest themselves through radio gyrosynchrotron emission. Here we show, for a well-observed flare on 2017~September~10, that a combination of \emph{RHESSI} hard X-ray and and SDO/AIA EUV observations provides a robust estimate of the fraction of the ambient electron population that is accelerated at a given time, with an upper limit of $\sim 10^{-2}$ on the number density of nonthermal ($\ge 20$~keV) electrons, expressed as a fraction of the number density of ambient protons in the same volume. This upper limit is about two orders of magnitude lower than previously inferred from microwave observations of the same event. Our results strongly indicate that the fraction of accelerated electrons in the coronal region at any given time is relatively small, but also that the overall duration of the HXR emission requires a steady resupply of electrons to the acceleration site. Simultaneous measurements of the instantaneous accelerated electron number density and the associated specific electron acceleration rate provide key constraints for a quantitative study of the mechanisms leading to electron acceleration in magnetic reconnection events. △ Less

Submitted 9 April, 2023; v1 submitted 3 April, 2023; originally announced April 2023.

Comments: 5 figures, 10 pages

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.

Comments: Accepted for publication in the Astrophysical Journal (on 16 Feb 2023), 17 pages, 15 figures

Journal ref: ApJ, 946, 33 (2023)

The frequency ratio and time delay of solar radio emissions with fundamental and harmonic components

Authors: Xingyao Chen, Eduard P. Kontar, Daniel L. Clarkson, Nicolina Chrysaphi

Abstract: Solar radio bursts generated through the plasma emission mechanism produce radiation near the local plasma frequency (fundamental emission) and double the plasma frequency (harmonic). While the theoretical ratio of these two frequencies is close to 2, simultaneous observations give ratios ranging from 1.6 to 2, suggesting either a ratio different from 2, a delay of the fundamental emission, or bot… ▽ More Solar radio bursts generated through the plasma emission mechanism produce radiation near the local plasma frequency (fundamental emission) and double the plasma frequency (harmonic). While the theoretical ratio of these two frequencies is close to 2, simultaneous observations give ratios ranging from 1.6 to 2, suggesting either a ratio different from 2, a delay of the fundamental emission, or both. To address this long-standing question, we conducted high frequency, high time resolution imaging spectroscopy of type III and type J bursts with fine structures for both the fundamental and harmonic components with LOFAR between 30 and 80 MHz. The short-lived and narrow frequency-band fine structures observed simultaneously at fundamental and harmonic frequencies give a frequency ratio of 1.66 and 1.73, similar to previous observations. However, frequency-time cross-correlations suggest a frequency ratio of 1.99 and 1.95 with a time delay between the F and H emissions of 1.00 and 1.67 s, respectively for each event. Hence, simultaneous frequency ratio measurements different from 2 are caused by the delay of the fundamental emission. Among the processes causing fundamental emission delays, anisotropic radio-wave scattering is dominant. Moreover, the levels of anisotropy and density fluctuations reproducing the delay of fundamental emissions are consistent with those required to simulate the source size and duration of fundamental emissions. Using these simulations we are able to, for the first time, provide quantitative estimates of the delay time of the fundamental emissions caused by radio-wave propagation effects at multiple frequencies, which can be used in future studies. △ Less

Submitted 26 January, 2023; originally announced January 2023.

Comments: 13 pages, 8 figures

Robust construction of differential emission measure profiles using a regularized maximum likelihood method

Authors: Paolo Massa, A. Gordon Emslie, Iain G. Hannah, Eduard P. Kontar

Abstract: Extreme-ultraviolet (EUV) observations provide considerable insight into evolving physical conditions in the active solar atmosphere. For a prescribed density and temperature structure, it is straightforward to construct the corresponding differential emission measure profile $ξ(T)$, such that $ξ(T) \, dT$ is proportional to the emissivity from plasma in the temperature range $[T, T + dT]$. Here w… ▽ More Extreme-ultraviolet (EUV) observations provide considerable insight into evolving physical conditions in the active solar atmosphere. For a prescribed density and temperature structure, it is straightforward to construct the corresponding differential emission measure profile $ξ(T)$, such that $ξ(T) \, dT$ is proportional to the emissivity from plasma in the temperature range $[T, T + dT]$. Here we study the inverse problem of obtaining a valid $ξ(T)$ profile from a set of EUV spectral line intensities observed at a pixel within a solar image. Our goal is to introduce and develop a regularized maximum likelihood (RML) algorithm designed to address the mathematically ill-posed problem of constructing differential emission measure profiles from a discrete set of EUV intensities in specified wavelength bands, specifically those observed by the Atmospheric Imaging Assembly (AIA) on the NASA Solar Dynamics Observatory. The RML method combines features of maximum likelihood and regularized approaches used by other authors. It is also guaranteed to produce a positive definite differential emission measure profile. We evaluate and compare the effectiveness of the method against other published algorithms, using both simulated data generated from parametric differential emission profile forms, and AIA data from a solar eruptive event on 2010 November 3. Similarities and differences between the differential emission measure profiles and maps reconstructed by the various algorithms are discussed. The RML inversion method is mathematically rigorous, computationally efficient, and produces acceptable measures of performance in the following three key areas: fidelity to the data, accuracy in the reconstruction, and robustness in the presence of data noise. As such, it shows considerable promise for computing differential emission measure profiles from datasets of discrete spectral lines. △ Less

Submitted 7 March, 2023; v1 submitted 11 January, 2023; originally announced January 2023.

Data-Constrained Solar Modeling with GX Simulator

Authors: Gelu M. Nita, Gregory D. Fleishman, Alexey A. Kuznetsov, Sergey A. Anfinogentov, Alexey G. Stupishin, Eduard P. Kontar, Samuel J. Schonfeld, James A. Klimchuk, Dale E. Gary

Abstract: To facilitate the study of solar active regions and flaring loops, we have created a modeling framework, the freely distributed GX Simulator IDL package, that combines 3D magnetic and plasma structures with thermal and non-thermal models of the chromosphere, transition region, and corona. The package has integrated tools to visualize the model data cubes, compute multi-wavelength emission maps fro… ▽ More To facilitate the study of solar active regions and flaring loops, we have created a modeling framework, the freely distributed GX Simulator IDL package, that combines 3D magnetic and plasma structures with thermal and non-thermal models of the chromosphere, transition region, and corona. The package has integrated tools to visualize the model data cubes, compute multi-wavelength emission maps from them, and quantitatively compare the resulting maps with observations. Its object-based modular architecture, which runs on Windows, Mac, and Unix/Linux platforms, offers capabilities that include the ability to either import 3D density and temperature distribution models, or to assign numerically defined coronal or chromospheric temperatures and densities, or their distributions to each individual voxel. The application integrates FORTRAN and C++ libraries for fast calculation of radio emission (free-free, gyroresonance, and gyrosynchrotron emission) along with soft and hard X-ray and EUV codes developed in IDL. To facilitate the creation of models, we have developed a fully automatic model production pipeline that downloads the required SDO/HMI vector magnetic field data and (optionally) the contextual SDO/AIA images, performs potential or nonlinear force free field extrapolations, populates the magnetic field skeleton with parameterized heated plasma coronal models that assume either steady-state or impulsive plasma heating, and generates non-LTE density and temperature distribution models of the chromosphere that are constrained by photospheric measurements. The standardized models produced by this pipeline may be further customized through a set of interactive tools provided by the graphical user interface. Here we describe the GX Simulator framework and its applications. △ Less

Submitted 2 January, 2023; originally announced January 2023.

Microwave Imaging of Quasi-periodic Pulsations at Flare Current Sheet

Authors: Yuankun Kou, Xin Cheng, Yulei Wang, Sijie Yu, Bin Chen, Eduard P. Kontar, Mingde Ding

Abstract: Quasi-periodic pulsations (QPPs) are frequently detected in solar and stellar flares, but the underlying physical mechanisms are still to be ascertained. Here, we show microwave QPPs during a solar flare originating from quasi-periodic magnetic reconnection at the flare current sheet. They appear as two vertically detached but closely related sources with the brighter ones located at flare loops a… ▽ More Quasi-periodic pulsations (QPPs) are frequently detected in solar and stellar flares, but the underlying physical mechanisms are still to be ascertained. Here, we show microwave QPPs during a solar flare originating from quasi-periodic magnetic reconnection at the flare current sheet. They appear as two vertically detached but closely related sources with the brighter ones located at flare loops and the weaker ones along the stretched current sheet. Although the brightness temperatures of the two microwave sources differ greatly, they vary in phase with periods of about 10--20 s and 30--60 s. The gyrosynchrotron-dominated microwave spectra also present a quasi-periodic soft-hard-soft evolution. These results suggest that relevant high-energy electrons are accelerated by quasi-periodic reconnection, likely arising from the modulation of magnetic islands within the current sheet as validated by a 2.5-dimensional magnetohydrodynamic simulation. △ Less

Submitted 16 December, 2022; originally announced December 2022.

Journal ref: Nature Communications (2022) 13:7680

The radial variation of the solar wind turbulence spectra near the kinetic break scale from Parker Solar Probe measurements

Authors: S. Lotz, A. E. Nel, R. T. Wicks, O. W. Roberts, N. E. Engelbrecht, R. D. Strauss, G. J. J. Botha, E. P. Kontar, A. Pitna, S. D. Bale

Abstract: In this study we examine the radial dependence of the inertial and dissipation range indices, as well as the spectral break separating the inertial and dissipation range in power density spectra of interplanetary magnetic field fluctuations using {\it Parker Solar Probe} data from the fifth solar encounter between $\sim$0.1 and $\sim$0.7 au. The derived break wavenumber compares reasonably well wi… ▽ More In this study we examine the radial dependence of the inertial and dissipation range indices, as well as the spectral break separating the inertial and dissipation range in power density spectra of interplanetary magnetic field fluctuations using {\it Parker Solar Probe} data from the fifth solar encounter between $\sim$0.1 and $\sim$0.7 au. The derived break wavenumber compares reasonably well with previous estimates at larger radial distances and is consistent with gyro-resonant damping of Alfvénic fluctuations by thermal protons. We find that the inertial scale power law index varies between approximately -1.65 and -1.45. This is consistent with either the Kolmogorov (-5/3) or Iroshnikov-Kraichnan (-3/2) values, has a very weak radial dependence with a possible hint that the spectrum becomes steeper closer to the Sun. The dissipation range power law index, however, has a clear dependence on radial distance (and turbulence age), decreasing from -3 near 0.7 au (4 days) to -4 [$\pm$0.3] at 0.1 au (0.75 days) closer to the Sun. △ Less

Submitted 5 December, 2022; originally announced December 2022.

Comments: Accepted for publication in ApJ

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.

Comments: Accepted for publication in the Astrophysical Journal

Spectral Analysis of Solar Radio Type III Bursts from 20 kHz to 410 MHz

Authors: K. Sasikumar Raja, Milan Maksimovic, Eduard P. Kontar, Xavier Bonnin, Philippe Zarka, Laurent Lamy, Hamish Reid, Nicole Vilmer, Alain Lecacheux, Vratislav Krupar, Baptiste Cecconi, Lahmiti Nora, Laurent Denis

Abstract: We present the statistical analysis of the spectral response of solar radio type III bursts over the wide frequency range between 20 kHz and 410 MHz. For this purpose, we have used observations that were carried out using both spaced-based (Wind/Waves) and ground-based (Nançay Decameter Array and Nançay Radioheliograph) facilities. In order to compare the flux densities observed by the different i… ▽ More We present the statistical analysis of the spectral response of solar radio type III bursts over the wide frequency range between 20 kHz and 410 MHz. For this purpose, we have used observations that were carried out using both spaced-based (Wind/Waves) and ground-based (Nançay Decameter Array and Nançay Radioheliograph) facilities. In order to compare the flux densities observed by the different instruments, we have carefully calibrated the data and displayed them in Solar Flux Units. The main result of our study is that type III bursts, in the metric to hectometric wavelength range, statistically exhibit a clear maximum of their median radio flux density around 2 MHz. Although this result was already reported by inspecting the spectral profiles of type III bursts in the frequency range 20 kHz - 20 MHz, our study extends such analysis for the first time to metric radio frequencies (i.e., from 20 kHz to 410 MHz) and confirms the maximum spectral response around 2 MHz. In addition, using a simple empirical model we show that the median radio flux $S$ of the studied dataset obeys the polynomial form $Y = 0.04 X^3 - 1.63 X^2 + 16.30 X -41.24$, with $X=\ln{(F_\text{MHz})}$ and with $Y=\ln{(S_\text{SFU})}$. Using the Sittler and Guhathakurtha model for coronal streamers \citep{Sit1999}, we have found that maximum of radio power falls therefore in the range 4 to 10 $R_{\odot}$, depending on whether the type III emissions are assumed to be at the fundamental or the harmonic. △ Less

Submitted 21 October, 2021; originally announced October 2021.

Comments: 21 Pages, 6 Figures, 1 Table, Accepted for publication in the Astrophysical Journal

The Spatial and Temporal Variations of Turbulence in a Solar Flare

Authors: Morgan Stores, Natasha L. S. Jeffrey, Eduard P. Kontar

Abstract: Magnetohydrodynamic (MHD) plasma turbulence is believed to play a vital role in the production of energetic electrons during solar flares and the non-thermal broadening of spectral lines is a key sign of this turbulence. Here, we determine how flare turbulence evolves in time and space using spectral profiles of Fe xxiv, Fe xxiii and Fe xvi, observed by Hinode/EIS. Maps of non-thermal velocity are… ▽ More Magnetohydrodynamic (MHD) plasma turbulence is believed to play a vital role in the production of energetic electrons during solar flares and the non-thermal broadening of spectral lines is a key sign of this turbulence. Here, we determine how flare turbulence evolves in time and space using spectral profiles of Fe xxiv, Fe xxiii and Fe xvi, observed by Hinode/EIS. Maps of non-thermal velocity are created for times covering the X-ray rise, peak, and decay. For the first time, the creation of kinetic energy density maps reveal where energy is available for energization, suggesting that similar levels of energy may be available to heat and/or accelerate electrons in large regions of the flare. We find that turbulence is distributed throughout the entire flare; often greatest in the coronal loop tops, and decaying at different rates at different locations. For hotter ions (Fe xxiv and Fe xxiii), the non-thermal velocity decreases as the flare evolves and during/after the X-ray peak shows a clear spatial variation decreasing linearly from the loop apex towards the ribbon. For the cooler ion (Fe xvi), the non-thermal velocity remains relativity constant throughout the flare, but steeply increases in one region corresponding to the southern ribbon, peaking just prior to the peak in hard X-rays before declining. The results suggest turbulence has a more complex temporal and spatial structure than previously assumed, while newly introduced turbulent kinetic energy maps show the availability of the energy and identify important spatial inhomogeneities in the macroscopic plasma motions leading to turbulence. △ Less

Submitted 4 October, 2021; originally announced October 2021.

Comments: Accepted for publication to Astrophysical Journal (October 1st 2021). Some figures may appear blurry in some pdf viewers, such as Preview. Please view pdf using Adobe

Simulations of radio-wave anisotropic scattering to interpret type III radio bursts measurements by Solar Orbiter, Parker Solar Probe, STEREO and Wind

Authors: S. Musset, M. Maksimovic, E. Kontar, V. Krupar, N. Chrysaphi, X. Bonnin, A. Vecchio, B. Cecconi, A. Zaslavsky, K. Issautier, S. D. Bale, M. Pulupa

Abstract: We use multi-spacecraft observations of invididual type III radio bursts in order to calculate the directivity of the radio emission, to be compared to the results of ray-tracing simulations of the radio-wave propagation and probe the plasma properties of the inner heliosphere. Ray-tracing simulations of radio-wave propagation with anisotropic scattering on density inhomogeneities are used to stud… ▽ More We use multi-spacecraft observations of invididual type III radio bursts in order to calculate the directivity of the radio emission, to be compared to the results of ray-tracing simulations of the radio-wave propagation and probe the plasma properties of the inner heliosphere. Ray-tracing simulations of radio-wave propagation with anisotropic scattering on density inhomogeneities are used to study the directivity of radio emissions. Simultaneous observations of type III radio bursts by four widely-separated spacecraft are used to calculate the directivity and position of the radio sources. The shape of the directivity pattern deduced for individual events is compared to the directivity pattern resulting from the ray-tracing simulations. We show that simultaneous observations of type radio III bursts by 4 different probes provide the opportunity to estimate the radio source positions and the directivity of the radio emission. The shape of the directivity varies from one event to another, and is consistent with anisotropic scattering of the radio-waves. △ Less

Submitted 29 September, 2021; v1 submitted 28 September, 2021; originally announced September 2021.

Journal ref: A&A 656, A34 (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.

Comments: Accepted for publication in the Astrophysical Journal Letters (on 4 Aug 2021), 8 pages, 5 figures

Journal ref: ApJL, 917, L32 (2021)

Energy budget of plasma motions, heating, and electron acceleration in a three-loop solar flare

Authors: Gregory D. Fleishman, Lucia Kleint, Galina G. Motorina, Gelu M. Nita, Eduard P. Kontar

Abstract: Non-potential magnetic energy promptly released in solar flares is converted to other forms of energy. This may include nonthermal energy of flare-accelerated particles, thermal energy of heated flaring plasma, and kinetic energy of eruptions, jets, up/down flows, and stochastic (turbulent) plasma motions. The processes or parameters governing partitioning of the released energy between these comp… ▽ More Non-potential magnetic energy promptly released in solar flares is converted to other forms of energy. This may include nonthermal energy of flare-accelerated particles, thermal energy of heated flaring plasma, and kinetic energy of eruptions, jets, up/down flows, and stochastic (turbulent) plasma motions. The processes or parameters governing partitioning of the released energy between these components is an open question. How these components are distributed between distinct flaring loops and what controls these spatial distributions is also unclear. Here, based on multi-wavelength data and 3D modeling, we quantify the energy partitioning and spatial distribution in the well observed SOL2014-02-16T064620 solar flare of class C1.5. Nonthermal emissions of this flare displayed a simple impulsive single-spike light curves lasting about 20\,s. In contrast, the thermal emission demonstrated at least three distinct heating episodes, only one of which was associated with the nonthermal component. The flare was accompanied by up and down flows and substantial turbulent velocities. The results of our analysis suggest that (i) the flare occurs in a multi-loop system that included at least three distinct flux tubes; (ii) the released magnetic energy is divided unevenly between the thermal and nonthermal components in these loops; (iii) only one of these three flaring loops contains an energetically important amount of nonthermal electrons, while two other loops remain thermal; (iv) the amounts of direct plasma heating and that due to nonthermal electron loss are comparable; (v) the kinetic energy in the flare footpoints constitute only a minor fraction compared with the thermal and nonthermal energies. △ Less

Submitted 1 April, 2021; originally announced April 2021.

Comments: ApJ accepted. 19 pages, 16 figures

Fine structure of type III solar radio bursts from Langmuir wave motion in turbulent plasma

Authors: Hamish A. S. Reid, Eduard P. Kontar

Abstract: The Sun frequently accelerates near-relativistic electron beams that travel out through the solar corona and interplanetary space. Interacting with their plasma environment, these beams produce type III radio bursts, the brightest astrophysical radio sources seen from the Earth. The formation and motion of type III fine frequency structures is a puzzle but is commonly believed to be related to pla… ▽ More The Sun frequently accelerates near-relativistic electron beams that travel out through the solar corona and interplanetary space. Interacting with their plasma environment, these beams produce type III radio bursts, the brightest astrophysical radio sources seen from the Earth. The formation and motion of type III fine frequency structures is a puzzle but is commonly believed to be related to plasma turbulence in the solar corona and solar wind. Combining a theoretical framework with kinetic simulations and high-resolution radio type III observations using the Low Frequency Array, we quantitatively show that the fine structures are caused by the moving intense clumps of Langmuir waves in a turbulent medium. Our results show how type III fine structure can be used to remotely analyse the intensity and spectrum of compressive density fluctuations, and can infer ambient temperatures in astrophysical plasma, both significantly expanding the current diagnostic potential of solar radio emission. △ Less

Submitted 15 March, 2021; originally announced March 2021.

Comments: 22 pages, 9 figures

Constraints on the acceleration region of type III radio bursts from decimetric radio spikes and faint X-ray bursts

Authors: Sophie Musset, Eduard Kontar, Lindsay Glesener, Nicole Vilmer, Abdallah Hamini

Abstract: We study the release of energy during the gradual phase of a flare, characterized by faint bursts of non-thermal hard X-ray (HXR) emission associated with decimetric radio spikes and type III radio bursts starting at high frequencies and extending to the heliosphere. We characterize the site of electron acceleration in the corona and study the radial evolution of radio source sizes in the high cor… ▽ More We study the release of energy during the gradual phase of a flare, characterized by faint bursts of non-thermal hard X-ray (HXR) emission associated with decimetric radio spikes and type III radio bursts starting at high frequencies and extending to the heliosphere. We characterize the site of electron acceleration in the corona and study the radial evolution of radio source sizes in the high corona. Imaging and spectroscopy of the HXR emission with Fermi and RHESSI provide a diagnostic of the accelerated electrons in the corona as well as a lower limit on the height of the acceleration region. Radio observations in the decimetric range with the ORFEES spectrograph provide radio diagnostics close to the acceleration region. Radio spectro-imaging with LOFAR in the meter range provide the evolution of the radio source sizes with their distance from the Sun, in the high corona. Non-thermal HXR bursts and radio spikes are well correlated on short timescales. The spectral index of non-thermal HXR emitting electrons is -4 and their number is about $2\times 10^{33}$ electrons/s. The density of the acceleration region is constrained between $1-5 \times 10^9$ cm$^{-3}$. Electrons accelerated upward rapidly become unstable to Langmuir wave production, leading to high starting frequencies of the type III radio bursts, and the elongation of the radio beam at its source is between 0.5 and 11.4 Mm. The radio source sizes and their gradient observed with LOFAR are larger than the expected size and gradient of the size of the electron beam, assuming it follows the expansion of the magnetic flux tubes. These observations support the idea that the fragmentation of the radio emission into spikes is linked to the fragmentation of the acceleration process itself. The combination of HXR and radio diagnostics in the corona provides strong constrains on the site of electron acceleration. △ Less

Submitted 19 January, 2021; originally announced January 2021.

Parametric simulation studies on the wave propagation of solar radio emission: the source size, duration, and position

Authors: PeiJin Zhang, ChuanBing Wang, Eduard P. Kontar

Abstract: The observed features of the radio sources indicate complex propagation effects embedded in the waves of solar radio bursts. In this work, we perform ray-tracing simulations on radio wave transport in the corona and interplanetary region with anisotropic electron density fluctuations. For the first time, the variation of the apparent source size, burst duration, and source position of both fundame… ▽ More The observed features of the radio sources indicate complex propagation effects embedded in the waves of solar radio bursts. In this work, we perform ray-tracing simulations on radio wave transport in the corona and interplanetary region with anisotropic electron density fluctuations. For the first time, the variation of the apparent source size, burst duration, and source position of both fundamental emission and harmonic emission at frequency 35 MHz are simulated as the function of the anisotropic parameter $α$ and the angular scattering rate coefficient $η=ε^2/h_0$, where $ε^2={\langle δn^2\rangle}/{n^2}$ is the density fluctuation level and $h_0$ is its correlation length near the wave exciting site. It is found that isotropic fluctuations produce a much larger decay time than a highly anisotropic fluctuation for fundamental emission. By comparing the observed duration and source size with the simulation results in the parameter space, we can estimate the scattering coefficient and the anisotropy parameter $η= 8.9\times 10^{-5}\, \mathrm{km^{-1}}$ and $α= 0.719$ with point pulse source assumption. Position offsets due to wave scattering and refraction can produce the co-spatial of fundamental and harmonic waves in observations of some type III radio bursts. The visual speed due to the wave propagation effect can reach 1.5\,$c$ for $η= 2.4\times 10^{-4}\, \mathrm{km^{-1}}$ and $α=0.2$ for fundamental emission in the sky plane, accompanying with large expansion rate of the source size. The visual speed direction is mostly identical to the offset direction, thus, for the observation aiming at obtaining the source position, the source centroid at the starting point is closer to the wave excitation point. △ Less

Submitted 4 January, 2021; originally announced January 2021.

Comments: 19 pages

Electron Acceleration during Macroscale Magnetic Reconnection

Authors: Harry Arnold, James Drake, Marc Swisdak, Fan Guo, Joel Dahlin, Bin Chen, Gregory Fleishman, Lindsay Glesener, Eduard Kontar, Tai Phan, Chengcai Shen

Abstract: The first self-consistent simulations of electron acceleration during magnetic reconnection in a macroscale system are presented. Consistent with solar flare observations the spectra of energetic electrons take the form of power-laws that extend more than two decades in energy. The drive mechanism for these nonthermal electrons is Fermi reflection in growing and merging magnetic flux ropes. A stro… ▽ More The first self-consistent simulations of electron acceleration during magnetic reconnection in a macroscale system are presented. Consistent with solar flare observations the spectra of energetic electrons take the form of power-laws that extend more than two decades in energy. The drive mechanism for these nonthermal electrons is Fermi reflection in growing and merging magnetic flux ropes. A strong guide field is found to suppress the production of nonthermal electrons by weakening the Fermi drive mechanism. For a weak guide field the total energy content of nonthermal electrons dominates that of the hot thermal electrons even though their number density remains small. Our results are benchmarked with the hard x-ray, radio and extreme ultra-violet (EUV) observations of the X8.2-class solar flare on September 10, 2017. △ Less

Submitted 22 January, 2021; v1 submitted 2 November, 2020; originally announced November 2020.

Journal ref: Phys. Rev. Lett. 126, 135101 (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.

Comments: 17 pages, 6 figures, accepted

The Solar Orbiter Science Activity Plan: translating solar and heliospheric physics questions into action

Authors: I. Zouganelis, A. De Groof, A. P. Walsh, D. R. Williams, D. Mueller, O. C. St Cyr, F. Auchere, D. Berghmans, A. Fludra, T. S. Horbury, R. A. Howard, S. Krucker, M. Maksimovic, C. J. Owen, J. Rodriiguez-Pacheco, M. Romoli, S. K. Solanki, C. Watson, L. Sanchez, J. Lefort, P. Osuna, H. R. Gilbert, T. Nieves-Chinchilla, L. Abbo, O. Alexandrova , et al. (160 additional authors not shown)

Abstract: Solar Orbiter is the first space mission observing the solar plasma both in situ and remotely, from a close distance, in and out of the ecliptic. The ultimate goal is to understand how the Sun produces and controls the heliosphere, filling the Solar System and driving the planetary environments. With six remote-sensing and four in-situ instrument suites, the coordination and planning of the operat… ▽ More Solar Orbiter is the first space mission observing the solar plasma both in situ and remotely, from a close distance, in and out of the ecliptic. The ultimate goal is to understand how the Sun produces and controls the heliosphere, filling the Solar System and driving the planetary environments. With six remote-sensing and four in-situ instrument suites, the coordination and planning of the operations are essential to address the following four top-level science questions: (1) What drives the solar wind and where does the coronal magnetic field originate? (2) How do solar transients drive heliospheric variability? (3) How do solar eruptions produce energetic particle radiation that fills the heliosphere? (4) How does the solar dynamo work and drive connections between the Sun and the heliosphere? Maximising the mission's science return requires considering the characteristics of each orbit, including the relative position of the spacecraft to Earth (affecting downlink rates), trajectory events (such as gravitational assist manoeuvres), and the phase of the solar activity cycle. Furthermore, since each orbit's science telemetry will be downloaded over the course of the following orbit, science operations must be planned at mission level, rather than at the level of individual orbits. It is important to explore the way in which those science questions are translated into an actual plan of observations that fits into the mission, thus ensuring that no opportunities are missed. First, the overarching goals are broken down into specific, answerable questions along with the required observations and the so-called Science Activity Plan (SAP) is developed to achieve this. The SAP groups objectives that require similar observations into Solar Orbiter Observing Plans (SOOPs), resulting in a strategic, top-level view of the optimal opportunities for science observations during the mission lifetime. △ Less

Submitted 22 September, 2020; originally announced September 2020.

Comments: 20 pages, 1 figure, accepted by Astronomy & Astrophysics

Journal ref: A&A 642, A3 (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.

Comments: Accepted for publication in the Astrophysical Journal (on 19 Aug 2020)

Probing solar flare accelerated electron distributions with prospective X-ray polarimetry missions

Authors: Natasha L. S. Jeffrey, Pascal Saint-Hilaire, Eduard P. Kontar

Abstract: Solar flare electron acceleration is an extremely efficient process, but the method of acceleration is not well constrained. Two of the essential diagnostics: electron anisotropy (velocity angle to the guiding magnetic field) and the high energy cutoff (highest energy electrons produced by the acceleration conditions: mechanism, spatial extent, time), are important quantities that can help to cons… ▽ More Solar flare electron acceleration is an extremely efficient process, but the method of acceleration is not well constrained. Two of the essential diagnostics: electron anisotropy (velocity angle to the guiding magnetic field) and the high energy cutoff (highest energy electrons produced by the acceleration conditions: mechanism, spatial extent, time), are important quantities that can help to constrain electron acceleration at the Sun but both are poorly determined. Here, using electron and X-ray transport simulations that account for both collisional and non-collisional transport processes such as turbulent scattering, and X-ray albedo, we show that X-ray polarization can be used to constrain the anisotropy of the accelerated electron distribution and the most energetic accelerated electrons together. Moreover, we show that prospective missions, e.g. CubeSat missions without imaging information, can be used alongside such simulations to determine these parameters. We conclude that a fuller understanding of flare acceleration processes will come from missions capable of both X-ray flux and polarization spectral measurements together. Although imaging polarimetry is highly desired, we demonstrate that spectro-polarimeters without imaging can also provide strong constraints on electron anisotropy and the high energy cutoff. △ Less

Submitted 18 December, 2020; v1 submitted 18 August, 2020; originally announced August 2020.

Comments: Accepted for publication in Astronomy & Astrophysics August 2020

Journal ref: A&A 642, A79 (2020)

Radio Echo in the Turbulent Corona and Simulations of Solar Drift-Pair Radio Bursts

Authors: Alexey A. Kuznetsov, Nicolina Chrysaphi, Eduard P. Kontar, Galina Motorina

Abstract: Drift-pair bursts are an unusual type of solar low-frequency radio emission, which appear in the dynamic spectra as two parallel drifting bright stripes separated in time. Recent imaging spectroscopy observations allowed for the quantitative characterization of the drifting pairs in terms of source size, position, and evolution. Here, the drift-pair parameters are qualitatively analyzed and compar… ▽ More Drift-pair bursts are an unusual type of solar low-frequency radio emission, which appear in the dynamic spectra as two parallel drifting bright stripes separated in time. Recent imaging spectroscopy observations allowed for the quantitative characterization of the drifting pairs in terms of source size, position, and evolution. Here, the drift-pair parameters are qualitatively analyzed and compared with the newly-developed Monte Carlo ray-tracing technique simulating radio-wave propagation in the inhomogeneous anisotropic turbulent solar corona. The results suggest that the drift-pair bursts can be formed due to a combination of the refraction and scattering processes, with the trailing component being the result of turbulent reflection (turbulent radio echo). The formation of drift-pair bursts requires an anisotropic scattering with the level of plasma density fluctuations comparable to that in type III bursts, but with a stronger anisotropy at the inner turbulence scale. The anisotropic radio-wave scattering model can quantitatively reproduce the key properties of drift-pair bursts: the apparent source size and its increase with time at a given frequency, the parallel motion of the source centroid positions, and the delay between the burst components. The trailing component is found to be virtually co-spatial and following the main component. The simulations suggest that the drift-pair bursts are likely to be observed closer to the disk center and below 100 MHz due to the effects of free-free absorption and scattering. The exciter of drift-pairs is consistent with propagating packets of whistlers, allowing for a fascinating way to diagnose the plasma turbulence and the radio emission mechanism. △ Less

Submitted 29 July, 2020; originally announced July 2020.

Comments: 13 pages, 7 figures

First Observation of a Type II Solar Radio Burst Transitioning Between a Stationary and Drifting State

Authors: Nicolina Chrysaphi, Hamish A. S. Reid, Eduard P. Kontar

Abstract: Standing shocks are believed to be responsible for stationary Type II solar radio bursts, whereas drifting Type II bursts are excited by moving shocks often related to coronal mass ejections (CMEs). Observations of either stationary or drifting Type II bursts are common, but a transition between the two states has not yet been reported. Here, we present a Type II burst which shows a clear, continu… ▽ More Standing shocks are believed to be responsible for stationary Type II solar radio bursts, whereas drifting Type II bursts are excited by moving shocks often related to coronal mass ejections (CMEs). Observations of either stationary or drifting Type II bursts are common, but a transition between the two states has not yet been reported. Here, we present a Type II burst which shows a clear, continuous transition from a stationary to a drifting state, the first observation of its kind. Moreover, band splitting is observed in the stationary parts of the burst, as well as intriguing negative and positive frequency-drift fine structures within the stationary emissions. The relation of the radio emissions to an observed jet and a narrow CME was investigated across multiple wavelengths, and the mechanisms leading to the transitioning Type II burst were determined. We find that a jet eruption generates a streamer-puff CME and that the interplay between the CME-driven shock and the streamer is likely to be responsible for the observed radio emissions. △ Less

Submitted 24 April, 2020; v1 submitted 24 March, 2020; originally announced March 2020.

Comments: Published in ApJ, 11 pages, 10 figures

Journal ref: ApJ, 893, 115 (2020)

Particle acceleration with anomalous pitch angle scattering in 3D separator reconnection

Authors: Alexei Borissov, Thomas Neukirch, Eduard Kontar, James Threlfall, Clare Parnell

Abstract: Understanding how the release of stored magnetic energy contributes to the generation of non-thermal high energy particles during solar flares is an important open problem in solar physics. Magnetic reconnection plays a fundamental role in the energy release and conversion processes taking place during flares. A common approach for investigating particle acceleration is to use test particles in fi… ▽ More Understanding how the release of stored magnetic energy contributes to the generation of non-thermal high energy particles during solar flares is an important open problem in solar physics. Magnetic reconnection plays a fundamental role in the energy release and conversion processes taking place during flares. A common approach for investigating particle acceleration is to use test particles in fields derived from magnetohydrodynamic (MHD) simulations of reconnection. These MHD simulations use anomalous resistivities that are much larger than the Spitzer resistivity based on Coulomb collisions. The processes leading to enhanced resistivity should also affect the test particles. We explore the link between resistivity and particle orbits building on a previous study using a 2D MHD simulation of magnetic reconnection. This paper extends the previous investigation to a 3D magnetic reconnection configuration and to study the effect on test particle orbits. We carried out orbit calculations using a 3D MHD simulation of separator reconnection. We use the relativistic guiding centre approximation including stochastic pitch angle scattering. The effects of varying the resistivity and the models for pitch angle scattering on particle orbit trajectories, final positions, energy spectra, final pitch angle distribution, and orbit duration are all studied in detail. Pitch angle scattering widens collimated beams of orbit trajectories, allowing orbits to access previously unaccessible field lines; this causes final positions to spread to topological structures that were previously inaccessible. Scattered orbit energy spectra are found to be predominantly affected by the level of anomalous resistivity, with the pitch angle scattering model only playing a role in isolated cases. Scattering is found to play a crucial role in determining the pitch angle and orbit duration distributions. △ Less

Submitted 21 January, 2020; originally announced January 2020.

Comments: 18 pages, 10 figures, accepted to Astronomy & Astrophysics

Journal ref: A&A 635, A63 (2020)

Density Fluctuations in the Solar Wind Based on Type III Radio Bursts Observed by Parker Solar Probe

Authors: Vratislav Krupar, Adam Szabo, Milan Maksimovic, Oksana Kruparova, Eduard P. Kontar, Laura A. Balmaceda, Xavier Bonnin, Stuart D. Bale, Marc Pulupa, David M. Malaspina, John W. Bonnell, Peter R. Harvey, Keith Goetz, Thierry Dudok de Wit, Robert J. MacDowall, Justin C. Kasper, Anthony W. Case, Kelly E. Korreck, Davin E. Larson, Roberto Livi, Michael L. Stevens, Phyllis L. Whittlesey, Alexander M. Hegedus

Abstract: Radio waves are strongly scattered in the solar wind, so that their apparent sources seem to be considerably larger and shifted than the actual ones. Since the scattering depends on the spectrum of density turbulence, better understanding of the radio wave propagation provides indirect information on the relative density fluctuations $ε=\langleδn\rangle/\langle n\rangle$ at the effective turbulenc… ▽ More Radio waves are strongly scattered in the solar wind, so that their apparent sources seem to be considerably larger and shifted than the actual ones. Since the scattering depends on the spectrum of density turbulence, better understanding of the radio wave propagation provides indirect information on the relative density fluctuations $ε=\langleδn\rangle/\langle n\rangle$ at the effective turbulence scale length. Here, we have analyzed 30 type III bursts detected by Parker Solar Probe (PSP). For the first time, we have retrieved type III burst decay times $τ_{\rm{d}}$ between 1 MHz and 10 MHz thanks to an unparalleled temporal resolution of PSP. We observed a significant deviation in a power-law slope for frequencies above 1 MHz when compared to previous measurements below 1 MHz by the twin-spacecraft Solar TErrestrial RElations Observatory (STEREO) mission. We note that altitudes of radio bursts generated at 1 MHz roughly coincide with an expected location of the Alfvén point, where the solar wind becomes super-Alfvénic. By comparing PSP observations and Monte Carlo simulations, we predict relative density fluctuations $ε$ at the effective turbulence scale length at radial distances between 2.5$R_\odot$ and 14$R_\odot$ to range from $0.22$ and $0.09$. Finally, we calculated relative density fluctuations $ε$ measured in situ by PSP at a radial distance from the Sun of $35.7$~$R_\odot$ during the perihelion \#1, and the perihelion \#2 to be $0.07$ and $0.06$, respectively. It is in a very good agreement with previous STEREO predictions ($ε=0.06-0.07$) obtained by remote measurements of radio sources generated at this radial distance. △ Less

Submitted 10 January, 2020; originally announced January 2020.

Comments: 12 pages, 10 figures, accepted for publication in ApJS

Spatio-temporal energy partitioning in a non-thermally dominated two-loop solar flare

Authors: Galina G. Motorina, Gregory D. Fleishman, Eduard P. Kontar

Abstract: Solar flares show remarkable variety of the energy partitioning between thermal and nonthermal components. Those with a prominent nonthermal component but only a modest thermal one are particularly well suited to study the direct effect of the nonthermal electrons on plasma heating. Here, we analyze such a well observed, impulsive single-spike nonthermal event, a SOL2013-11-05T035054 solar flare,… ▽ More Solar flares show remarkable variety of the energy partitioning between thermal and nonthermal components. Those with a prominent nonthermal component but only a modest thermal one are particularly well suited to study the direct effect of the nonthermal electrons on plasma heating. Here, we analyze such a well observed, impulsive single-spike nonthermal event, a SOL2013-11-05T035054 solar flare, where the plasma heating can be entirely attributed to the energy losses of these impulsively accelerated electrons. Evolution of the energy budget of thermal and nonthermal components during the flare is analysed using X-ray, microwave, and EUV observations and three-dimensional modeling. The results suggest that (i) the flare geometry is consistent with a two-loop morphology and the magnetic energy is likely released due to interaction between these two loops; (ii) the released magnetic energy converted to the nonthermal energy of accelerated electrons only, which is subsequently converted to the thermal energy of the plasma; (iii) the energy is partitioned in these two flaring loops in comparable amounts; (iv) one of these flaring loops remained relatively tenuous but rather hot, while the other remained relatively cool but denser than the first one. Therefore, this solar flare demonstrates an extreme efficiency of conversion of the free magnetic energy to the nonthermal energy of particle acceleration and the energy flow into two loops from the non-thermal to thermal one with a negligible direct heating. △ Less

Submitted 7 January, 2020; originally announced January 2020.

Comments: 23 pages, 10 figures, 1 table, ApJ, accepted

First imaging spectroscopy observations of solar drift pair bursts

Authors: Alexey Kuznetsov, Eduard Kontar

Abstract: Drift pairs are an unusual type of fine structure sometimes observed in dynamic spectra of solar radio emission. They appear as two identical short narrowband drifting stripes separated in time; both positive and negative frequency drifts are observed. Using the Low Frequency Array (LOFAR), we report unique observations of a cluster of drift pair bursts in the frequency range of 30-70 MHz made on… ▽ More Drift pairs are an unusual type of fine structure sometimes observed in dynamic spectra of solar radio emission. They appear as two identical short narrowband drifting stripes separated in time; both positive and negative frequency drifts are observed. Using the Low Frequency Array (LOFAR), we report unique observations of a cluster of drift pair bursts in the frequency range of 30-70 MHz made on 12 July 2017. Spectral imaging capabilities of the instrument have allowed us for the first time to resolve the temporal and frequency evolution of the source locations and sizes at a fixed frequency and along the drifting pair components. Sources of two components of a drift pair have been imaged and found to propagate in the same direction along nearly the same trajectories. Motion of the second component source is delayed in time with respect to that of the first one. The source trajectories can be complicated and non-radial; positive and negative frequency drifts correspond to opposite propagation directions. The drift pair bursts with positive and negative frequency drifts, as well as the associated broadband type-III-like bursts, are produced in the same regions. The visible source velocities are variable from zero to a few $10^4$ (up to ${\sim 10^5}$) km/s, which often exceeds the velocities inferred from the drift rate ($\sim 10^4$ km/s). The visible source sizes are of about $10'-18'$; they are more compact than typical type III sources at the same frequencies. The existing models of drift pair bursts cannot adequately explain the observed features. We discuss the key issues that need to be addressed, and in particular the anisotropic scattering of the radio waves. The broadband bursts observed simultaneously with the drift pairs differ in some aspects from common type III bursts and may represent a separate type of emission. △ Less

Submitted 22 October, 2019; originally announced October 2019.

Comments: Accepted for publication in Astronomy & Astrophysics Letters

On the Source Position and Duration of a Solar Type III Radio Burst Observed by LOFAR

Authors: PeiJin Zhang, SiJie Yu, Eduard Kontar, ChuanBing Wang

Abstract: Solar type III radio bursts are excited by electron beams propagating outward from the Sun. The flux of type III radio burst has a time profile of rising and decay phase at a given frequency, which has been actively studied since 1970s. Several factors that may influence the duration of a type III radio burst has been proposed. However, the major cause of the duration is still an open question. In… ▽ More Solar type III radio bursts are excited by electron beams propagating outward from the Sun. The flux of type III radio burst has a time profile of rising and decay phase at a given frequency, which has been actively studied since 1970s. Several factors that may influence the duration of a type III radio burst has been proposed. However, the major cause of the duration is still an open question. In this work, to study the dominant cause of the duration, we investigate the source positions of the front edge, the peak, and the tail edge in the dynamic spectrum of a single and clear type III radio burst. The duration of this type III burst at a given frequency is about 3 second for decameter wave. The beam-formed observations by the LOw-Frequency ARray (LOFAR) are used, which can provide the radio source positions and the dynamic spectra at the same time. We find that, for this burst, the source positions of the front edge, the peak, and the tail edge split with each other spatially. The radial speed of the front edge, the peak, and the tail edge is 0.42 c, 0.25 c, and 0.16 c, respectively. We estimate the influences of the corona density fluctuation and the electron-velocity dispersion on the duration, and the scattering effect by comparison with a few short-duration bursts from the same region. The analysis yields that, in the frequency range of 30 - 41 MHz, the electron-velocity dispersion is the dominant factor that determines the time duration of type III radio bursts with long duration, while scattering may play important role in the duration of short bursts. △ Less

Submitted 18 September, 2019; originally announced September 2019.

Comments: 12 pages

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.

Comments: 21 pages, 11 figures, accepted for publication in The Astrophysical Journal

Global Energetics of Solar Flares: VIII. The Low-Energy Cutoff

Authors: Markus Aschwanden, Eduard P. Kontar, Natasha L. S. Jeffrey

Abstract: One of the key problems in solar flare physics is the determination of the low-energy cut-off; the value that determines the energy of nonthermal electrons and hence flare energetics. We discuss different approaches to determine the low-energy cut-off in the spectrum of accelerated electrons: (i) the total electron number model, (ii) the time-of-flight model (based on the equivalence of the time-o… ▽ More One of the key problems in solar flare physics is the determination of the low-energy cut-off; the value that determines the energy of nonthermal electrons and hence flare energetics. We discuss different approaches to determine the low-energy cut-off in the spectrum of accelerated electrons: (i) the total electron number model, (ii) the time-of-flight model (based on the equivalence of the time-of-flight and the collisional deflection time); (iii) the warm target model of Kontar et al.~(2015), and (iv) the model of the spectral cross-over between thermal and nonthermal components. We find that the first three models are consistent with a low-energy cutoff with a mean value of $\approx 10$ keV, while the cross-over model provides an upper limit for the low-energy cutoff with a mean value of $ \approx 21$ keV. Combining the first three models we find that the ratio of the nonthermal energy to the dissipated magnetic energy in solar flares has a mean value of $q_E=0.57\pm0.08$, which is consistent with an earlier study based on the simplified approximation of the warm target model alone ($q_E=0.51\pm0.17$). This study corroborates the self-consistency between three different low-energy cutoff models in the calculation of nonthermal flare energies. △ Less

Submitted 13 June, 2019; originally announced June 2019.

Comments: 9 Figures

Coronal Loop Scaling Laws for Various Forms of Parallel Heat Conduction

Authors: Stephen J. Bradshaw, A. Gordon Emslie, Nicolas H. Bian, Eduard P. Kontar

Abstract: The solar atmosphere is dominated by loops of magnetic flux which connect the multi-million-degree corona to the much cooler chromosphere. The temperature and density structure of quasi-static loops is determined by the continuous flow of energy from the hot corona to the lower solar atmosphere. Loop scaling laws provide relationships between global properties of the loop (such as peak temperature… ▽ More The solar atmosphere is dominated by loops of magnetic flux which connect the multi-million-degree corona to the much cooler chromosphere. The temperature and density structure of quasi-static loops is determined by the continuous flow of energy from the hot corona to the lower solar atmosphere. Loop scaling laws provide relationships between global properties of the loop (such as peak temperature, pressure, and length); they follow from the physical variable dependencies of various terms in the energy equation, and hence the form of the loop scaling law provides insight into the key physics that controls the loop structure. Traditionally, scaling laws have been derived under the assumption of collision-dominated thermal conduction. Here we examine the impact of different regimes of thermal conduction -- collision-dominated, turbulence-dominated, and free-streaming -- on the form of the scaling laws relating the loop temperature and heating rate to its pressure and half-length. We show that the scaling laws for turbulence-dominated conduction are fundamentally different than those for collision-dominated and free-streaming conduction, inasmuch as the form of the scaling laws now depend primarily on conditions at the low-temperature, rather than high-temperature, part of the loop. We also establish regimes in temperature and density space in which each of the applicable scaling laws prevail. △ Less

Submitted 7 June, 2019; originally announced June 2019.

The Role of Energy Diffusion in the Deposition of Energetic Electron Energy in Solar and Stellar Flares

Authors: Natasha L. S. Jeffrey, Eduard P. Kontar, Lyndsay Fletcher

Abstract: During solar flares, a large fraction of the released magnetic energy is carried by energetic electrons that transfer and deposit energy in the Sun's atmosphere. Electron transport is often approximated by a cold thick-target model (CTTM), assuming that electron energy is much larger than the temperature of the ambient plasma, and electron energy evolution is modeled as a systematic loss. Using ki… ▽ More During solar flares, a large fraction of the released magnetic energy is carried by energetic electrons that transfer and deposit energy in the Sun's atmosphere. Electron transport is often approximated by a cold thick-target model (CTTM), assuming that electron energy is much larger than the temperature of the ambient plasma, and electron energy evolution is modeled as a systematic loss. Using kinetic modeling of electrons, we re-evaluate the transport and deposition of flare energy. Using a full collisional warm-target model (WTM), we account for electron thermalization and for the properties of the ambient coronal plasma such as its number density, temperature and spatial extent. We show that the deposition of non-thermal electron energy in the lower atmosphere is highly dependent on the properties of the flaring coronal plasma. In general, thermalization and a reduced WTM energy loss rate leads to an increase of non-thermal energy transferred to the chromosphere, and the deposition of non-thermal energy at greater depths. The simulations show that energy is deposited in the lower atmosphere initially by high energy non-thermal electrons, and later by lower energy non-thermal electrons that partially or fully thermalize in the corona, over timescales of seconds, unaccounted for in previous studies. This delayed heating may act as a diagnostic of both the injected non-thermal electron distribution and the coronal plasma, vital for constraining flare energetics. △ Less

Submitted 5 June, 2019; originally announced June 2019.

Comments: Accepted by the Astrophysical Journal, June 2019

Nonlinear wave interactions as a model for naturally enhanced ion-acoustic lines in the ionosphere

Authors: E. P. Kontar, H. L. Pecseli

Abstract: Incoherent radar scatter from the ionosphere will, for equilibrium conditions, show two symmetric ion-acoustic lines, one for each direction of wave propagation. Many observation, from the EISCAT Svalbard Radar (ESR) for instance, demonstrate that the symmetry of this ion line can be broken, accompanied by an enhanced, nonthermal, level of fluctuations, i.e., Naturally Enhanced Ion-Acoustic Lines… ▽ More Incoherent radar scatter from the ionosphere will, for equilibrium conditions, show two symmetric ion-acoustic lines, one for each direction of wave propagation. Many observation, from the EISCAT Svalbard Radar (ESR) for instance, demonstrate that the symmetry of this ion line can be broken, accompanied by an enhanced, nonthermal, level of fluctuations, i.e., Naturally Enhanced Ion-Acoustic Lines (NEIALs). Several models have been proposed for explaining these naturally enhanced lines. Here, we consider one of these, suggesting that decay of electron beam excited Langmuir waves gives rise to enhanced asymmetric ion lines in the ionosphere. We use a weak-turbulence approximation, and identify crucial parameters for Langmuir decay processes to be effective in generating the observed signals. △ Less

Submitted 3 May, 2019; originally announced May 2019.

Comments: 4 figures

Dynamics of electron beams in the solar corona plasma with density fluctuations

Authors: Eduard P. Kontar

Abstract: The problem of beam propagation in a plasma with small scale and low intensity inhomogeneities is investigated. It is shown that the electron beam propagates in a plasma as a beam-plasma structure and is a source of Langmuir waves. The plasma inhomogeneity changes the spatial distribution of the waves. The spatial distribution of the waves is fully determined by the distribution of plasma inhomoge… ▽ More The problem of beam propagation in a plasma with small scale and low intensity inhomogeneities is investigated. It is shown that the electron beam propagates in a plasma as a beam-plasma structure and is a source of Langmuir waves. The plasma inhomogeneity changes the spatial distribution of the waves. The spatial distribution of the waves is fully determined by the distribution of plasma inhomogeneities. The possible applications to the theory of radio emission associated with electron beams are discussed. △ Less

Submitted 11 April, 2019; originally announced April 2019.

Comments: 8 pages, 7 figures

Dynamics of electron beams in the inhomogeneous solar corona plasma

Authors: Eduard P. Kontar

Abstract: Dynamics of an spatially limited electron beam in the inhomogeneous solar corona plasma is considered in the framework of weak turbulence theory when the temperature of the beam significantly exceeds that of surrounding plasma. The numerical solution of kinetic equations manifests that generally the beam accompanied by Langmuir waves propagates as a beam-plasma structure with a decreasing velocity… ▽ More Dynamics of an spatially limited electron beam in the inhomogeneous solar corona plasma is considered in the framework of weak turbulence theory when the temperature of the beam significantly exceeds that of surrounding plasma. The numerical solution of kinetic equations manifests that generally the beam accompanied by Langmuir waves propagates as a beam-plasma structure with a decreasing velocity. Unlike the uniform plasma case the structure propagates with the energy losses in the form of Langmuir waves. The results obtained are compared with the results of observations of type III bursts. It is shown that the deceleration of type III sources can be explained by the corona inhomogeneity. The frequency drift rates of the type III sources are found in a good agreement with the numerical results of beam dynamics. △ Less

Submitted 21 March, 2019; originally announced March 2019.

Comments: 10 figures, 19 pages

Numerical consideration of quasilinear electron cloud dynamics in plasma

Authors: Eduard P. Kontar

Abstract: The dynamics of a hot electron cloud in the solar corona-like plasma based on the numerical solution of kinetic equations of weak turbulence theory is considered. Different finite difference schemes are examined to fit the exact analytical solutions of quasilinear equations in hydrodynamic limit (gas-dynamic solution). It is shown that the scheme suggested demonstrates correct asymptotic behavior… ▽ More The dynamics of a hot electron cloud in the solar corona-like plasma based on the numerical solution of kinetic equations of weak turbulence theory is considered. Different finite difference schemes are examined to fit the exact analytical solutions of quasilinear equations in hydrodynamic limit (gas-dynamic solution). It is shown that the scheme suggested demonstrates correct asymptotic behavior and can be employed to solve initial value problems for an arbitrary initial electron distribution function. △ Less

Submitted 20 March, 2019; originally announced March 2019.

Comments: 17 pages, 6 figures

Nonlinear development of electron-beam-driven weak turbulence in an inhomogeneous plasma

Authors: E. P. Kontar, H. L. Pecseli

Abstract: The self-consistent description of Langmuir wave and ion-sound wave turbulence in the presence of an electron beam is presented for inhomogeneous non-isothermal plasmas. Full numerical solutions of the complete set of kinetic equations for electrons, Langmuir waves, and ion-sound waves are obtained for a inhomogeneous unmagnetized plasma. The result show that the presence of inhomogeneity signific… ▽ More The self-consistent description of Langmuir wave and ion-sound wave turbulence in the presence of an electron beam is presented for inhomogeneous non-isothermal plasmas. Full numerical solutions of the complete set of kinetic equations for electrons, Langmuir waves, and ion-sound waves are obtained for a inhomogeneous unmagnetized plasma. The result show that the presence of inhomogeneity significantly changes the overall evolution of the system. The inhomogeneity is effective in shifting the wavenumbers of the Langmuir waves, and can thus switch between different process governing the weakly turbulent state. The results can be applied to a variety of plasma conditions, where we choose solar coronal parameters as an illustration, when performing the numerical analysis. △ Less

Submitted 20 March, 2019; originally announced March 2019.

Comments: 11 pages, 8 figures

A Fokker-Planck Framework for Studying the Diffusion of Radio Burst Waves in the Solar Corona

Authors: N. H. Bian, A. G. Emslie, E. P. Kontar

Abstract: Electromagnetic wave scattering off density inhomogeneities in the solar corona is an important process which determines both the apparent source size and the time profile of radio bursts observed at 1 AU. Here we model the scattering process using a Fokker-Planck equation and apply this formalism to several regimes of interest. In the first regime the density fluctuations are considered quasi-sta… ▽ More Electromagnetic wave scattering off density inhomogeneities in the solar corona is an important process which determines both the apparent source size and the time profile of radio bursts observed at 1 AU. Here we model the scattering process using a Fokker-Planck equation and apply this formalism to several regimes of interest. In the first regime the density fluctuations are considered quasi-static and diffusion in wavevector space is dominated by angular diffusion on the surface of a constant energy sphere. In the small-angle ("pencil beam") approximation, this diffusion further occurs over a small solid angle in wavevector space. The second regime corresponds to a much later time, by which scattering has rendered the photon distribution near-isotropic resulting in a spatial diffusion of the radiation. The third regime involves time-dependent fluctuations and, therefore, Fermi acceleration of photons. Combined, these results provide a comprehensive theoretical framework within which to understand several important features of propagation of radio burst waves in the solar corona: emitted photons are accelerated in a relatively small inner region and then diffuse outwards to larger distances. En route, angular diffusion results both in source sizes which are substantially larger than the intrinsic source, and in observed intensity-versus-time profiles that are asymmetric, with a sharp rise and an exponential decay. Both of these features are consistent with observations of solar radio bursts. △ Less

Submitted 12 February, 2019; v1 submitted 1 February, 2019; originally announced February 2019.

Comments: 28 pages , 1 figure, submitted to ApJ

Electron distribution and energy release in magnetic reconnection outflow regions during the pre-impulsive phase of a solar flare

Authors: Marina Battaglia, Eduard P. Kontar, Galina Motorina

Abstract: We present observations of electron energization in magnetic reconnection outflows during the pre-impulsive phase of solar flare SOL2012-07-19T05:58. During a time-interval of about 20 minutes, starting 40 minutes before the onset of the impulsive phase, two X-ray sources were observed in the corona, one above the presumed reconnection region and one below. For both of these sources, the mean elec… ▽ More We present observations of electron energization in magnetic reconnection outflows during the pre-impulsive phase of solar flare SOL2012-07-19T05:58. During a time-interval of about 20 minutes, starting 40 minutes before the onset of the impulsive phase, two X-ray sources were observed in the corona, one above the presumed reconnection region and one below. For both of these sources, the mean electron distribution function as a function of time is determined over an energy range from 0.1~keV up to several tens of keV, for the first time. This is done by simultaneous forward fitting of X-ray and EUV data. Imaging spectroscopy with RHESSI provides information on the high-energy tail of the electron distribution in these sources while EUV images from SDO/AIA are used to constrain the low specific electron energies. The measured electron distribution spectrum in the magnetic reconnection outflows is consistent with a time-evolving kappa-distribution with $κ=3.5-5.5$. The spectral evolution suggests that electrons are accelerated to progressively higher energies in the source above the reconnection region, while in the source below, the spectral shape does not change but an overall increase of the emission measure is observed, suggesting density increase due to evaporation. The main mechanisms by which energy is transported away from the source regions are conduction and free-streaming electrons. The latter dominates by more than one order of magnitude and is comparable to typical non-thermal energies during the hard X-ray peak of solar flares, suggesting efficient acceleration even during this early phase of the event. △ Less

Submitted 23 January, 2019; originally announced January 2019.

Comments: 9 pages, 5 figures, accepted for publication in ApJ

Determination of the total accelerated electron rate and power using solar flare hard X-ray spectra

Authors: Eduard P. Kontar, Natasha L. S. Jeffrey, A. Gordon Emslie

Abstract: Solar flare hard X-ray spectroscopy serves as a key diagnostic of the accelerated electron spectrum. However, the standard approach using the collisional cold thick-target model poorly constrains the lower-energy part of the accelerated electron spectrum, and hence the overall energetics of the accelerated electrons are typically constrained only to within one or two orders of magnitude. Here we d… ▽ More Solar flare hard X-ray spectroscopy serves as a key diagnostic of the accelerated electron spectrum. However, the standard approach using the collisional cold thick-target model poorly constrains the lower-energy part of the accelerated electron spectrum, and hence the overall energetics of the accelerated electrons are typically constrained only to within one or two orders of magnitude. Here we develop and apply a physically self-consistent warm-target approach which involves the use of both hard X-ray spectroscopy and imaging data. The approach allows an accurate determination of the electron distribution low-energy cutoff, and hence the electron acceleration rate and the contribution of accelerated electrons to the total energy released, by constraining the coronal plasma parameters. Using a solar flare observed in X-rays by the {\em RHESSI} spacecraft, we demonstrate that using the standard cold-target methodology, the low-energy cutoff (and hence the energy content in electrons) is essentially undetermined. However, the warm-target methodology can determine the low-energy electron cutoff with $\sim$7\% uncertainty at the $3σ$ level and hence permits an accurate quantitative study of the importance of accelerated electrons in solar flare energetics. △ Less

Submitted 4 January, 2019; v1 submitted 22 December, 2018; originally announced December 2018.

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

CME-driven Shock and Type II Solar Radio Burst Band Splitting

Authors: Nicolina Chrysaphi, Eduard P. Kontar, Gordon D. Holman, Manuela Temmer

Abstract: Coronal Mass Ejections (CMEs) are believed to be effective in producing shocks in the solar corona and the interplanetary space. One of the important signatures of shocks and shock acceleration are Type II solar radio bursts that drift with the shock speed and produce bands of fundamental and higher harmonic plasma radio emission. An intriguing aspect of Type II radio bursts is the occasional spli… ▽ More Coronal Mass Ejections (CMEs) are believed to be effective in producing shocks in the solar corona and the interplanetary space. One of the important signatures of shocks and shock acceleration are Type II solar radio bursts that drift with the shock speed and produce bands of fundamental and higher harmonic plasma radio emission. An intriguing aspect of Type II radio bursts is the occasional split of a harmonic band into thinner lanes, known as band-splitting. Here, we report a detailed imaging and spectroscopic observation of a CME-driven shock producing band-splitting in a Type II burst. Using the Low Frequency Array (LOFAR), we examine the spatial and temporal relation of the Type II burst to the associated CME event, use source imaging to calculate the apparent coronal density, and demonstrate how source imaging can be used to estimate projection effects. We consider two widely accepted band-splitting models that make opposing predictions regarding the locations of the true emission sources with respect to the shock front. Our observations suggest that the locations of the upper and lower sub-band sources are spatially separated by $\sim 0.2 \pm 0.05 \, \mathrm{R_\odot}$. However, we quantitatively show, for the first time, that such separation is consistent with radio-wave scattering of plasma radio emission from a single region, implying that the split-band Type II sources could originate from nearly co-spatial locations. Considering the effects of scattering, the observations provide supporting evidence for the model that interprets the band-splitting as emission originating in the upstream and downstream regions of the shock front, two virtually co-spatial areas. △ Less

Submitted 27 November, 2018; v1 submitted 18 October, 2018; originally announced October 2018.

Comments: Published by ApJ, 13 pages, 7 figures