The authors of this white paper met on 16-17 January 2020 at the New Jersey Institute of Technology, Newark, NJ, for a 2-day workshop that brought together a group of heliophysicists, data providers, expert modelers, and computer/data scientists. Their objective was to discuss critical developments and prospects of the application of machine and/or deep learning techniques for data analysis, modeling and forecasting in Heliophysics, and to shape a strategy for further developments in the field. The workshop combined a set of plenary sessions featuring invited introductory talks interleaved with a set of open discussion sessions. The outcome of the discussion is encapsulated in this white paper that also features a top-level list of recommendations agreed by participants.

In the standard model of solar flares, a large-scale reconnection current sheet is postulated as the central engine for powering the flare energy release and accelerating particles. However, where and how the energy release and particle acceleration occur remain unclear due to the lack of measurements for the magnetic properties of the current sheet. Here we report the measurement of spatially-resolved magnetic field and flare-accelerated relativistic electrons along a current-sheet feature in a solar flare. The measured magnetic field profile shows a local maximum where the reconnecting field lines of opposite polarities closely approach each other, known as the reconnection X point. The measurements also reveal a local minimum near the bottom of the current sheet above the flare loop-top, referred to as a "magnetic bottle". This spatial structure agrees with theoretical predictions and numerical modeling results. A strong reconnection electric field of ~4000 V/m is inferred near the X point. This location, however, shows a local depletion of microwave-emitting relativistic electrons. These electrons concentrate instead at or near the magnetic bottle structure, where more than 99% of them reside at each instant. Our observations suggest that the loop-top magnetic bottle is likely the primary site for accelerating and/or confining the relativistic electrons.

Nonthermal electrons accelerated in solar flares produce electromagnetic emission in two distinct, highly complementary domains - hard X-rays (HXRs) and microwaves (MWs). This paper reports MW imaging spectroscopy observations from the Expanded Owens Valley Solar Array of an M1.2 flare that occurred on 2017 September 9, from which we deduce evolving coronal parameter maps. We analyze these data jointly with the complementary Reuven Ramaty High-Energy Solar Spectroscopic Imager HXR data to reveal the spatially-resolved evolution of the nonthermal electrons in the flaring volume. We find that the high-energy portion of the nonthermal electron distribution, responsible for the MW emission, displays a much more prominent evolution (in the form of strong spectral hardening) than the low-energy portion, responsible for the HXR emission. We show that the revealed trends are consistent with a single electron population evolving according to a simplified trap-plus-precipitation model with sustained injection/acceleration of nonthermal electrons, which produces a double-powerlaw with steadily increasing break energy.

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.

We succeeded in observing two large spicules simultaneously with the Atacama Large Millimeter/submillimeter Array (ALMA), the Interface Region Imaging Spectrograph (IRIS), and the Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory. One is a spicule seen in the IRIS Mg II slit-jaw images and AIA 304\AA\ images (MgII/304A spicule). The other one is a spicule seen in the 100GHz images obtained with ALMA (100GHz spicule). Although the 100GHz spicule overlapped with the MgII/304A spicule in the early phase, it did not show any corresponding structures in the IRIS Mg II and AIA 304A images after the early phase. It suggests that the spicules are individual events and do not have a physical relationship. To obtain the physical parameters of the 100GHz spicule, we estimate the optical depths as a function of temperature and density using two different methods. One is using the observed brightness temperature by assuming a filling factor, and the other is using an emission model for the optical depth. As a result of comparing them, the kinetic temperature of the plasma and the number density of ionized hydrogens in the 100GHz spicule are ~6800 K and 2.2 x 10^10 cm^-3. The estimated values can explain the absorbing structure in the 193A image, which appear as a counterpart of the 100GHz spicule. These results suggest that the 100GHz spicule presented in this paper is classified to a macrospicule without a hot sheath in former terminology.

Particle acceleration in solar flares remains an outstanding problem in solar physics. It is yet unclear which of the acceleration mechanisms dominates and how exactly is the excessive magnetic energy transferred to the nonthermal and other forms of energy. We emphasize, that the ultimate acceleration mechanism must be capable of efficiently working in the most extreme conditions, such as the shortest detected time scales and the highest acceleration efficiency. Here we focus on detailed multiwavelength analysis of a very initial phase of the SOL2011-08-04 flare, which demonstrated prominent short subpeaks of nonthermal emission during filament eruption associated with the flare. We demonstrate that the three-dimensional configuration of the flare, combined with timing and spectral behavior of the rapidly varying component, put very stringent constraints on the acceleration regime. Specifically, the rapid subpeaks are generated by short injections of nonthermal electrons with a reasonably hard, single power-law spectrum and a relatively narrow spread of pitch-angles along the mean magnetic field. The acceleration site is a compact volume located near the top of extended coronal loop(s). The electrons are accelerated up to several hundreds of keV promptly, with the characteristic acceleration time shorter than 50 ms. We show, that these properties are difficult to reconcile with widely adopted stochastic acceleration models, while the data inescapably require acceleration by a super-Dreicer electric field, whether regular or random.

At the Sun, the strongest magnetic fields are routinely detected at dark sunspots. The magnitude of the field is typically about 3000 G, with only a few exceptions that reported the magnetic field in excess of 5000 G. Given that the magnetic field decreases with height in the solar atmosphere, no coronal magnetic field above ~2000 G was ever reported. Here, we present imaging microwave observations of an anomalously strong magnetic field of about 4000 G at the base of the corona in solar active region NOAA 12673 on 06 September 2017. Combining the photospheric vector measurements of the magnetic field and the coronal probing, we created and validated a nonlinear force-free field coronal model, with which we quantify the record-breaking coronal magnetic field at various coronal heights.

We report hard X-ray and gamma-ray observations of the impulsive phase of the SOL2017-09-06T11:55 X9.3 solar flare. We focus on a high-energy part of the spectrum, >100 keV, and perform time resolved spectral analysis for a portion of the impulsive phase, recorded by the Konus-Wind experiment, that displayed prominent gamma-ray emission. Given a variety of possible emission components contributing to the gamma-ray emission, we employ a Bayesian inference to build the most probable fitting model. The analysis confidently revealed contributions from nuclear deexcitation lines, electron-positron annihilation line at 511 keV, and a neutron capture line at 2.223 MeV along with two components of the bremsstrahlung continuum. The revealed time evolution of the spectral components is particularly interesting. The low-energy bremsstrahlung continuum shows a soft-hard-soft pattern typical for impulsive flares, while the high-energy one shows a persistent hardening at the course of the flare. The neutron capture line emission shows an unusually short time delay relative to the nuclear deexcitation line component, which implies that the production of neutrons was significantly reduced soon after the event onset. This in turn may imply a prominent softening of the accelerated proton spectrum at the course of the flare, similar to the observed softening of the low-energy component of the accelerated electrons responsible for the low-energy bremsstrahlung continuum. We discuss possible physical scenarios, which might result in the obtained relationships between these gamma-ray components.

The Sun is an active star that can have a direct impact on the Earth, its magnetosphere, and the technological infrastructure on which modern society depends. Among the phenomena that drive "space weather" are fast solar wind streams and co-rotating interaction regions, solar flares, coronal mass ejections, the shocks they produce, and the energetic particles they accelerate. Radio emission from these and associated phenomena offer unique diagnostic possibilities that complement those available at other wavelengths. Here, the relevant space weather drivers are briefly described, the potential role of radio observations is outlined, and the requirements of an instrument to provide them are provided: specifically, ultrabroadband imaging spectropolarimetry. The insights provided by radio observations of space weather drivers will not only inform the science of space weather, they will pave the way for new tools for forecasting and "nowcasting" space weather. They will also serve as an important touchstone against which local environment of exoplanets and the impact of "exo-space weather" can be evaluated.

Identification of the mechanisms responsible for heating the solar chromosphere and corona remains an outstanding problem, one of great relevance to late-type stars as well. There has been tremendous progress in the past decade, largely driven by new instruments, new observations, and sophisticated modeling efforts. Despite this progress, gaps remain. We briefly discuss the need for radio coverage of the 3D solar atmosphere and discuss the requirements.

Magnetic reconnection is a fundamental physical process in many laboratory, space, and astrophysical plasma contexts. Solar flares serve as an outstanding laboratory to study the magnetic reconnection and the associated energy release and conversion processes under plasma conditions difficult to reproduce in the laboratory, and with considerable spatiotemporal details not possible elsewhere in astrophysics. Here we emphasize the unique power of remote-sensing observations of solar flares at radio wavelengths. In particular, we discuss the transformative technique of broadband radio dynamic imaging spectroscopy in making significant contributions to addressing several outstanding challenges in magnetic reconnection, including the capability of pinpointing magnetic reconnection sites, measuring the time-evolving reconnecting magnetic fields, and deriving the spatially and temporally resolved distribution function of flare-accelerated electrons.

3D picture of the coronal magnetic field remains an outstanding problem in solar physics, particularly, in active regions. Nonlinear force-free field reconstructions that employ routinely available full-disk photospheric vector magnetograms represent state-of-the-art coronal magnetic field modeling. Such reconstructions, however, suffer from an inconsistency between a force-free coronal magnetic field and non-force-free photospheric boundary condition, from which the coronal reconstruction is performed. In this study we focus on integrating the additional chromospheric and / or coronal magnetic field data with the vector photospheric magnetograms with the goal of improving the reliability of the magnetic field reconstructions. We develop a corresponding modification of the available optimization codes described in Fleishman et al. (2017) and test their performance using a full-fledged MHD model obtained from the Bifrost code by performing a `voxel-by-voxel' comparison between the reconstructed and the model magnetic fields. We demonstrate that adding even an incomplete set of chromospheric magnetic field data can measurably improve the reconstruction of the coronal magnetic field, greatly improve reconstructions of the magnetic connectivity and of the coronal electric current.

The authors of this report met on 28-30 March 2018 at the New Jersey Institute of Technology, Newark, New Jersey, for a 3-day workshop that brought together a group of data providers, expert modelers, and computer and data scientists, in the solar discipline. Their objective was to identify challenges in the path towards building an effective framework to achieve transformative advances in the understanding and forecasting of the Sun-Earth system from the upper convection zone of the Sun to the Earth's magnetosphere. The workshop aimed to develop a research roadmap that targets the scientific challenge of coupling observations and modeling with emerging data-science research to extract knowledge from the large volumes of data (observed and simulated) while stimulating computer science with new research applications. The desire among the attendees was to promote future trans-disciplinary collaborations and identify areas of convergence across disciplines. The workshop combined a set of plenary sessions featuring invited introductory talks and workshop progress reports, interleaved with a set of breakout sessions focused on specific topics of interest. Each breakout group generated short documents, listing the challenges identified during their discussions in addition to possible ways of attacking them collectively. These documents were combined into this report-wherein a list of prioritized activities have been collated, shared and endorsed.

Sudden jets of collimated plasma arise from many locations on the Sun, including active regions. The magnetic field along which a jet emerges is often open to interplanetary space, offering a clear "escape route" for any flare-accelerated electrons and making jets lucrative targets for studying particle acceleration and the solar sources of transient heliospheric events. Bremsstrahlung hard X-rays (HXRs) could, in principle, trace the accelerated electrons that escape along the paths of the jets, but measurements of the escaping electron beams are customarily difficult due to the low densities of the corona. In this work, we augment HXR observations with gyrosynchrotron emission observed in microwaves, as well as extreme ultraviolet (EUV) emission and modeling to investigate flare-accelerated electrons in a coronal jet. HXR and microwave data from RHESSI and OVSA, respectively, give complementary insight into electron spectra and locations, including the presence of accelerated electrons in the jet itself. High-time-resolution HXR data from the Konus-Wind instrument suggest electron acceleration timescales on the order of 1 second or shorter. We model the energetic electron distributions in the GX Simulator framework using SoHO/MDI, RHESSI, TRACE, and OVSA data as constraints. The result is a modeled distribution, informed and constrained by measurements, of accelerated electrons as they escape the Sun. Combining the detection of microwave gyrosynchrotron emission from an open, rather than closed, magnetic configuration, with realistic 3D modeling constrained by magnetograms, EUV, and X-ray emission, we obtain the most stringent constraints to date on the accelerated electrons within a solar jet.

Solar flares are due to the catastrophic release of magnetic energy in the Sun's corona, resulting in plasma heating, mass motions, particle acceleration, and radiation emitted from radio to $\gamma$-ray wavelengths. They are associated with global coronal eruptions of plasma into the interplanetary medium---coronal mass ejections---that can result in a variety of "space weather" phenomena. Flares release energy over a vast range of energies, from $\sim\!10^{23}$ ergs (nanoflares) to more than $10^{32}$ ergs. Solar flares are a phenomenon of general astrophysical interest, allowing detailed study of magnetic energy release, eruptive processes, shock formation and propagation, particle acceleration and transport, and radiative processes. Observations at radio wavelengths offer unique diagnostics of the physics of flares. To fully exploit these diagnostics requires the means of performing time-resolved imaging spectropolarimetry. Recent observations with the Jansky Very Large Array (JVLA) and the Expanded Owens Valley Solar Array (EOVSA), supported by extensive development in forward modeling, have demonstrated the power of the approach. The ngVLA has the potential to bring our understanding of flare processes to a new level through its combination of high spatial resolution, broad frequency range, and imaging dynamic range---especially when used in concert with multi-wavelength observations and data at hard X-ray energies.

Energy stored in the magnetic field in the solar atmosphere above active regions is a key driver of all solar activity (e.g., solar flares and coronal mass ejections), some of which can affect life on Earth. Radio observations provide a unique diagnostic of the coronal magnetic fields that make them a critical tool for the study of these phenomena, using the technique of broadband radio imaging spectropolarimetry. Observations with the ngVLA will provide unique observations of coronal magnetic fields and their evolution, key inputs and constraints for MHD numerical models of the solar atmosphere and eruptive processes, and a key link between lower layers of the solar atmosphere and the heliosphere. In doing so they will also provide practical "research to operations" guidance for space weather forecasting.

Observations of solar flares at sub-THz frequencies (mm and sub-mm wavelengths) over the last two decades often show a spectral component rising with frequency. Unlike a typical gyrosynchrotron spectrum decreasing with frequency, or a weak thermal component from hot coronal plasma, the observations can demonstrate a high flux level (up to ~10^4 s.f.u. at 0.4 THz) and fast variability on sub-second time scales. Although, many models has been put forward to explain the puzzling observations, none of them have clear observational support. Here we propose a scenario to explain the intriguing sub-THz observations. We show that the model, based on free-free emission from the plasma of flare ribbons at temperatures 10^4-10^6 K, is consistent with all existing observations of frequency-rising sub-THz flare emission. The model provides a temperature diagnostic of the flaring chromosphere and suggests fast heating and cooling of the dense transition region plasma.

Quantifying coronal magnetic field remains a central problem in solar physics. Nowadays the coronal magnetic field is often modelled using nonlinear force-free field (NLFFF) reconstructions, whose accuracy has not yet been comprehensively assessed. Here we perform a detailed casting of the NLFFF reconstruction tools, such as pi-disambiguation, photospheric field preprocessing, and volume reconstruction methods using a 3D snapshot of the publicly available full-fledged radiative MHD model. Specifically, from the MHD model we know the magnetic field vector in the entire 3D domain, which enables us to perform "voxel-by-voxel" comparison of the restored and the true magnetic field in the 3D model volume. Our tests show that the available pi-disambiguation methods often fail at the quiet sun areas dominated by small-scale magnetic elements, while they work well at the AR photosphere and (even better) chromosphere. The preprocessing of the photospheric magnetic field, although does produce a more force-free boundary condition, also results in some effective `elevation' of the magnetic field components. This `elevation' height is different for the longitudinal and transverse components, which results in a systematic error in absolute heights in the reconstructed magnetic data cube. The extrapolations performed starting from actual AR photospheric magnetogram are free from this systematic error, while have other metrics comparable with those for extrapolations from the preprocessed magnetograms. This finding favors the use of extrapolations from the original photospheric magnetogram without preprocessing. Our tests further suggest that extrapolations from a force-free chromospheric boundary produce measurably better results, than those from the photospheric boundary.

The spectrum of gyrosynchrotron emission from solar flares generally peaks in the microwave range. Its optically-thin, high-frequency component, above the spectral peak, is often used for diagnostics of the nonthermal electrons and the magnetic field in the radio source. Under favorable conditions, its low-frequency counterpart brings additional, complementary information about these parameters as well as thermal plasma diagnostics, either through gyrosynchrotron self-absorption, free-free absorption by the thermal plasma, or the suppression of emission through the so-called Razin effect. However, their effects on the low-frequency spectrum are often masked by spatial nonuniformity. To disentangle the various contributions to low-frequency gyrosynchrotron emission, a combination of spectral and imaging data is needed. To this end, we have investigated Owens Valley Solar Array (OVSA) multi-frequency images for 26 solar bursts observed jointly with Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) during the first half of 2002. For each, we examined dynamic spectra, time- and frequency-synthesis maps, RHESSI images with overlaid OVSA contours, and a few representative single-frequency snapshot OVSA images. We focus on the frequency dependence of microwave source sizes derived from the OVSA images and their effect on the low-frequency microwave spectral slope. We succeed in categorizing 18 analyzed events into several groups. Four events demonstrate clear evidence of being dominated by gyrosynchrotron self-absorption, with an inferred brightness temperature of $\geq10^8$~K. The low-frequency spectra in the remaining events are affected to varying degree by Razin suppression. We find that many radio sources are rather large at low frequencies, which can have important implications for solar energetic particle production and escape.

We report the first science results from the newly completed Expanded Owens Valley Solar Array (EOVSA), which obtained excellent microwave imaging spectroscopy observations of SOL2017-09-10, a classic partially-occulted solar limb flare associated with an erupting flux rope. This event is also well-covered by the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) in hard X-rays (HXRs). We present an overview of this event focusing on microwave and HXR data, both associated with high-energy nonthermal electrons, and discuss them within the context of the flare geometry and evolution revealed by extreme ultraviolet (EUV) observations from the Atmospheric Imaging Assembly aboard the Solar Dynamics Observatory (SDO/AIA). The EOVSA and RHESSI data reveal the evolving spatial and energy distribution of high-energy electrons throughout the entire flaring region. The results suggest that the microwave and HXR sources largely arise from a common nonthermal electron population, although the microwave imaging spectroscopy provides information over a much larger volume of the corona.

Understanding nonthermal particle generation, transport, and escape in solar flares requires detailed quantification of the particle evolution in the realistic 3D domain where the flare takes place. Rather surprisingly, apart of standard flare scenario and integral characteristics of the nonthermal electrons, not much is known about actual evolution of nonthermal electrons in the 3D spatial domain. This paper attempts to begin to remedy this situation by creating sets of evolving 3D models, the synthesized emission from which matches the evolving observed emission. Here we investigate two contrasting flares: a dense, "coronal-thick-target" flare SOL2002-04-12T17:42, that contained a single flare loop observed in both microwave and X-ray, and a more complex flare, SOL2015-06-22T17:50, that contained at least four distinct flaring loops needed to consistently reproduce the microwave and X-ray emission. Our analysis reveals differing evolution pattern of the nonthermal electrons in the dense and tenuous loops; however, both of which imply the central role of resonant wave-particle interaction with turbulence. These results offer new constraints for theory and models of the particle acceleration and transport in solar flares.

Why the tenuous solar outer atmosphere, or corona, is much hotter than the underlying layers remains one of the greatest challenge for solar modeling. Detailed diagnostics of the coronal thermal structure come from extreme ultraviolet (EUV) emission. The EUV emission is produced by heavy ions in various ionization states and depends on the amount of these ions and on plasma temperature and density. Any nonuniformity of the elemental distribution in space or variability in time affects thermal diagnostics of the corona. Here we theoretically predict ionized chemical element concentrations in some areas of the solar atmosphere, where the electric current is directed upward. We then detect these areas observationally, by comparing the electric current density with the EUV brightness in an active region. We found a significant excess in EUV brightness in the areas with positive current density rather than negative. Therefore, we report the observational discovery of substantial concentrations of heavy ions in current-carrying magnetic flux tubes, which might have important implications for the elemental fractionation in the solar corona known as the first ionization potential (FIP) effect. We call such areas of heavy ion concentration the "ion traps." These traps hold enhanced ion levels until they are disrupted by a flare whether large or small.

Solar flares often happen after a preflare / preheating phase, which is almost or entirely thermal. In contrast, there are the so-called early impulsive flares that do not show a (significant) preflare heating but instead often show the Neupert effect--a relationship where the impulsive phase is followed by a gradual, cumulative-like, thermal response. This has been interpreted as a dominance of nonthermal energy release at the impulsive phase, even though a similar phenomenology is expected if the thermal and nonthermal energies are released in comparable amounts at the impulsive phase. Nevertheless, some flares do show a good quantitative correspondence between the nonthermal electron energy input and plasma heating, in such cases the thermal response was weak, which results in calling them "cold" flares. We undertook a systematic search of such events among early impulsive flares registered by Konus-Wind instrument in the triggered mode from 11/1994 to 04/2017 and selected 27 cold flares based on relationships between HXR (Konus-Wind) and SXR (GOES) emission. For these events we put together all available microwave data from different instruments. We obtained temporal and spectral parameters of HXR and microwave emissions of the events and examined correlations between them. We found that, compared with a `mean' flare, the cold flares: (i) are weaker, shorter, and harder in the X-ray domain, (ii) are harder and shorter, but not weaker in the microwaves, (iii) have a significantly higher spectral peak frequencies in the microwaves. We discuss the possible physical reasons for these distinctions and implication of the finding.

The well-established notion of a "common population" of the accelerated electrons simultaneously producing the hard X-ray (HXR) and the microwave (MW) emission during the flare impulsive phase has been challenged by some studies reporting the discrepancies between the HXR-inferred and the MW-inferred electron energy spectra. The traditional methods of their spectral inversion have some problems that can be mainly attributed to the unrealistic and the oversimplified treatment of the flare emission. To properly address this problem, we use a Non-linear Force Free Field (NLFFF) model extrapolated from an observed photospheric magnetogram as input to the three-dimensional, multi-wavelength modeling platform GX Simulator, and create a unified electron population model that can simultaneously reproduce the observed HXR and MW observations. We model the end of the impulsive phase of the 2015-06-22 M6.5 flare, and constrain the modeled electron spatial and energy parameters using observations made by the highest-resolving instruments currently available in two wavelengths, the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) for HXR and the Expanded Owens Valley Solar Array (EOVSA) for MW. Our results suggest that the HXR-emitting electron population model fits the standard flare model with a broken power-law spectrum (E_break ~ 200 keV) that simultaneously produces the HXR footpoint emission and the MW high frequency emission. The model also includes an "HXR invisible" population of nonthermal electrons that are trapped in a large volume of magnetic field above the HXR-emitting loops, which is observable by its gyrosynchrotron (GS) radiation emitting mainly in MW low frequency range.

The new generations of multiwavelength radioheliographs with high spatial resolution will employ microwave imaging spectropolarimetry to recover flare topology and plasma parameters in the flare sources and along the wave propagation paths. The recorded polarization depends on the emission mechanism and emission regime (optically thick or thin), the emitting particle properties, and propagation effects. Here, we report an unusual flare, SOL2012-07-06T01:37, whose optically thin gyrosynchrotron emission of the main source displays an apparently ordinary mode sense of polarization in contrast to the classical theory that favors the extraordinary mode. This flare produced copious nonthermal emission in hard X-rays and in high-frequency microwaves up to 80 GHz. It is found that the main flare source corresponds to an interaction site of two loops with greatly different sizes. We have investigated the three possible known reasons of the circular polarization sense reversal - mode coupling, positron contribution, and the effect of beamed angular distribution. We excluded polarization reversal due to contribution of positrons because there was no relevant response in the X-ray emission. We find that a beam-like electron distribution can produce the observed polarization behavior, but the source thermal density must be much higher than the estimate from to the X-ray data. We conclude that the apparent ordinary wave emission in the optically thin mode is due to radio wave propagation across the quasi-transverse (QT) layer. The abnormally high transition frequency (above 35 GHz) can be achieved reasonably low in the corona where the magnetic field value is high and transverse to the line of sight. This places the microwave source below this QT layer, i.e. very low in the corona.

Ever-increasing multi-frequency imaging of solar observations suggests that solar flares often involve more than one magnetic fluxtube. Some of the fluxtubes are closed, while others can contain open field. The relative proportion of nonthermal electrons among those distinct loops is highly important for understanding the energy release, particle acceleration, and transport. The access of nonthermal electrons to the open field is further important as the open field facilitates the solar energetic particle (SEP) escape from the flaring site, and thus controls the SEP fluxes in the solar system, both directly and as seed particles for further acceleration. The large-scale fluxtubes are often filled with a tenuous plasma, which is difficult to detect in either EUV or X-ray wavelengths; however, they can dominate at low radio frequencies, where a modest component of nonthermal electrons can render the source optically thick and, thus, bright enough to be observed. Here we report detection of a large-scale `plume' at the impulsive phase of an X-class solar flare, SOL2001-08-25T16:23, using multi-frequency radio data from Owens Valley Solar Array. To quantify the flare spatial structure, we employ 3D modeling utilizing force-free-field extrapolations from the line-of-sight SOHO/MDI magnetograms with our modeling tool GX Simulator. We found that a significant fraction of the nonthermal electrons accelerated at the flare site low in the corona escapes to the plume, which contains both closed and open field. We propose that the proportion between the closed and open field at the plume is what determines the SEP population escaping into interplanetary space.

Observations of the Sun at millimeter and submillimeter wavelengths offer a unique probe into the structure, dynamics, and heating of the chromosphere; the structure of sunspots; the formation and eruption of prominences and filaments; and energetic phenomena such as jets and flares. High-resolution observations of the Sun at millimeter and submillimeter wavelengths are challenging due to the intense, extended, low- contrast, and dynamic nature of emission from the quiet Sun, and the extremely intense and variable nature of emissions associated with energetic phenomena. The Atacama Large Millimeter/submillimeter Array (ALMA) was designed with solar observations in mind. The requirements for solar observations are significantly different from observations of sidereal sources and special measures are necessary to successfully carry out this type of observations. We describe the commissioning efforts that enable the use of two frequency bands, the 3 mm band (Band 3) and the 1.25 mm band (Band 6), for continuum interferometric-imaging observations of the Sun with ALMA. Examples of high-resolution synthesized images obtained using the newly commissioned modes during the solar commissioning campaign held in December 2015 are presented. Although only 30 of the eventual 66 ALMA antennas were used for the campaign, the solar images synthesized from the ALMA commissioning data reveal new features of the solar atmosphere that demonstrate the potential power of ALMA solar observations. The ongoing expansion of ALMA and solar-commissioning efforts will continue to enable new and unique solar observing capabilities.

Solar flares are generally believed to be powered by free magnetic energy stored in the corona, but the build up of coronal energy alone may be insufficient for the imminent flare occurrence. The flare onset mechanism is a critical but less understood problem, insights into which could be gained from small-scale energy releases known as precursors, which are observed as small pre-flare brightenings in various wavelengths, and also from certain small-scale magnetic configurations such as the opposite polarity fluxes, where magnetic orientation of small bipoles is opposite to that of the ambient main polarities. However, high-resolution observations of flare precursors together with the associated photospheric magnetic field dynamics are lacking. Here we study precursors of a flare using unprecedented spatiotemporal resolution of the 1.6 m New Solar Telescope, complemented by novel microwave data. Two episodes of precursor brightenings are initiated at a small-scale magnetic channel (a form of opposite polarity fluxes) with multiple polarity inversions and enhanced magnetic fluxes and currents, lying near the footpoints of sheared magnetic loops. The low-atmospheric origin of these precursor emissions is corroborated by microwave spectra. We propose that the emerging magnetic channel field interacts with the sheared arcades to cause precursor brightenings at the main flare core region. These high-resolution results provide evidence of low-atmospheric small-scale energy release and possible relationship to the onset of the main flare.

We use state-of-the-art, three-dimensional non-local thermodynamic equilibrium (non-LTE) radiative magnetohydrodynamic simulations of the quiet solar atmosphere to carry out detailed tests of chromospheric magnetic field diagnostics from free-free radiation at millimeter and submillimeter wavelengths (mm/submm). The vertical component of the magnetic field was deduced from the mm/submm brightness spectra and the degree of circular polarization synthesized at millimeter frequencies. We used the frequency bands observed by the Atacama Large Millimeter/Submillimeter Array (ALMA) as a convenient reference. The magnetic field maps obtained describe the longitudinal magnetic field at the effective formation heights of the relevant wavelengths in the solar chromosphere. The comparison of the deduced and model chromospheric magnetic fields at the spatial resolution of both the model and current observations demonstrates a good correlation, but has a tendency to underestimate the model field. The systematic discrepancy of about 10 percent is probably due to averaging of the restored field over the heights contributing to the radiation, weighted by the strength of the contribution. On the whole, the method of probing the longitudinal component of the magnetic field with free-free emission at mm/submm wavelengths is found to be applicable to measurements of the weak quiet-Sun magnetic fields. However, successful exploitation of this technique requires very accurate measurements of the polarization properties (primary beam and receiver polarization response) of the antennas, which will be the principal factor that determines the level to which chromospheric magnetic fields can be measured. Consequently, high-resolution and high-precision observations of circularly polarized radiation at millimeter wavelengths can be a powerful tool for producing chromospheric longitudinal magnetograms.

Recently, in a few case studies we demonstrated that gyrosynchrotron microwave emission can be detected directly from the acceleration region when the trapped electron component is insignificant. For the statistical study reported here, we have identified events with steep (narrowband) microwave spectra that do not show a significant trapped component and at the same time show evidence of source uniformity, which simplifies the data analysis greatly. Initially, we identified a subset of more than 20 radio bursts with such narrow spectra, having low- and high-frequency spectral indices larger than 3 in absolute value. A steep low-frequency spectrum implies that the emission is nonthermal (for optically-thick thermal emission, the spectral index cannot be steeper than 2), and the source is reasonably dense and uniform. A steep high-frequency spectrum implies that no significant electron trapping occurs; otherwise a progressive spectral flattening would be observed. Roughly half of these radio bursts have RHESSI data, which allows for detailed, joint diagnostics of the source parameters and evolution. Based on an analysis of radio-to-X-ray spatial relationships, timing, and spectral fits, we conclude that the microwave emission in these narrowband bursts originates directly from the acceleration regions, which have relatively strong magnetic field, high density, and low temperature. In contrast, the thermal X-ray emission comes from a distinct loop with smaller magnetic field, lower density, but higher temperature. Therefore, these flares occurred likely due to interaction between two (or more) magnetic loops.

This document was created by the Solar Simulations for the Atacama Large Millimeter Observatory Network (SSALMON) in preparation of the first regular observations of the Sun with the Atacama Large Millimeter/submillimeter Array (ALMA), which are anticipated to start in ALMA Cycle 4 in October 2016. The science cases presented here demonstrate that a large number of scientifically highly interesting observations could be made already with the still limited solar observing modes foreseen for Cycle 4 and that ALMA has the potential to make important contributions to answering long-standing scientific questions in solar physics. With the proposal deadline for ALMA Cycle 4 in April 2016 and the Commissioning and Science Verification campaign in December 2015 in sight, several of the SSALMON Expert Teams composed strategic documents in which they outlined potential solar observations that could be feasible given the anticipated technical capabilities in Cycle 4. These documents have been combined and supplemented with an analysis, resulting in recommendations for solar observing with ALMA in Cycle 4. In addition, the detailed science cases also demonstrate the scientific priorities of the solar physics community and which capabilities are wanted for the next observing cycles. The work on this White Paper effort was coordinated in close cooperation with the two international solar ALMA development studies led by T. Bastian (NRAO, USA) and R. Brajsa, (ESO). This document will be further updated until the beginning of Cycle 4 in October 2016. In particular, we plan to adjust the technical capabilities of the solar observing modes once finally decided and to further demonstrate the feasibility and scientific potential of the included science cases by means of numerical simulations of the solar atmosphere and corresponding simulated ALMA observations.

Recently, a number of peculiar flares have been reported, which demonstrate significant non-thermal particle signatures with a low, if any, thermal emission, that implies close association of the observed emission with the primary energy release/electron acceleration region. This paper presents a flare that appears a "cold" one at the impulsive phase, while displaying a delayed heating later on. Using HXR data from \kw, microwave observations by SSRT, RSTN, NoRH and NoRP, context observations, and 3D modeling, we study the energy release, particle acceleration and transport, and the relationships between the nonthermal and thermal signatures. The flaring process is found to involve interaction between a small and a big loop and the accelerated particles divided in roughly equal numbers between them. Precipitation of the electrons from the small loop produced only weak thermal response because the loop volume was small, while the electrons trapped in the big loop lost most of their energy in the coronal part of the loop, which resulted in the coronal plasma heating but no or only weak chromospheric evaporation, and thus unusually weak soft X-ray emission. Energy losses of fast electrons in the big tenuous loop were slow resulting in the observed delay of the plasma heating. We determined that the impulsively accelerated electron population had a beamed angular distribution in the direction of electric force along the magnetic field of the small loop. The accelerated particle transport in big loop was primarily mediated by turbulent waves like in the other reported cold flares.

The Atacama Large Millimeter/submillimeter Array (ALMA) is a new powerful tool for observing the Sun at high spatial, temporal, and spectral resolution. These capabilities can address a broad range of fundamental scientific questions in solar physics. The radiation observed by ALMA originates mostly from the chromosphere - a complex and dynamic region between the photosphere and corona, which plays a crucial role in the transport of energy and matter and, ultimately, the heating of the outer layers of the solar atmosphere. Based on first solar test observations, strategies for regular solar campaigns are currently being developed. State-of-the-art numerical simulations of the solar atmosphere and modeling of instrumental effects can help constrain and optimize future observing modes for ALMA. Here we present a short technical description of ALMA and an overview of past efforts and future possibilities for solar observations at submillimeter and millimeter wavelengths. In addition, selected numerical simulations and observations at other wavelengths demonstrate ALMA's scientific potential for studying the Sun for a large range of science cases.

We present detailed 3D modeling of a dense, coronal thick target X-ray flare using the GX Simulator tool, photospheric magnetic measurements, and microwave imaging and spectroscopy data. The developed model offers a remarkable agreement between the synthesized and observed spectra and images in both X-ray and microwave domains, which validates the entire model. The flaring loop parameters are chosen to reproduce the emission measure, temperature, and the nonthermal electron distribution at low energies derived from the X-ray spectral fit, while the remaining parameters, unconstrained by the X-ray data, are selected such as to match the microwave images and total power spectra. The modeling suggests that the accelerated electrons are trapped in the coronal part of the flaring loop, but away from where the magnetic field is minimal, and, thus, demonstrates that the data are clearly inconsistent with electron magnetic trapping in the weak diffusion regime mediated by the Coulomb collisions. Thus, the modeling supports the interpretation of the coronal thick-target sources as sites of electron acceleration in flares and supplies us with a realistic 3D model with physical parameters of the acceleration region and flaring loop.

We developed an efficient algorithm integrated in our 3D modeling tool, GX Simulator (Nita et al. 2015), allowing quick computation of the synthetic intensity and polarization maps of solar active regions (AR) in the ALMA spectral range. The algorithm analyzes the photospheric input (white light and magnetogram) to classify a given photospheric pixel to belong to a given photospheric structure. Then, a 1D chromospheric model (Fontenla et al. 2009) is added on top of each pixel, which forms a chromospheric model of the AR. Next step is computation of the mm and sub-mm emission produced from this chromosphere model. A huge advantage of this approach is that emission from any given AR can be synthesized very fast, on the order of a few minutes after the AR selection. Using the GX Simulator tool it is also possible to produce synthetic maps of the microwave (gyroresonance) and EUV emission from the same AR model and compare them with the ALMA synthetic maps and with the corresponding observed microwave and/or EUV data.

The Solar Simulations for the Atacama Large Millimeter Observatory Network (SSALMON) was initiated in 2014 in connection with two ALMA development studies. The Atacama Large Millimeter/submillimeter Array (ALMA) is a powerful new tool, which can also observe the Sun at high spatial, temporal, and spectral resolution. The international SSALMONetwork aims at coordinating the further development of solar observing modes for ALMA and at promoting scientific opportunities for solar physics with particular focus on numerical simulations, which can provide important constraints for the observing modes and can aid the interpretation of future observations. The radiation detected by ALMA originates mostly in the solar chromosphere - a complex and dynamic layer between the photosphere and corona, which plays an important role in the transport of energy and matter and the heating of the outer layers of the solar atmosphere. Potential targets include active regions, prominences, quiet Sun regions, flares. Here, we give a brief overview over the network and potential science cases for future solar observations with ALMA.

We have simulated the Expanded Owens Valley Solar Array (EOVSA) radio images generated at multiple frequencies from a model solar active region, embedded in a realistic solar disk model, and explored the resulting datacube for different spectral analysis schemes to evaluate the potential for realizing one of EOVSA's most important scientific goals--coronal magnetography. In this paper, we focus on modeling the gyroresonance and free-free emission from an on-disk solar active region model with realistic complexities in electron density, temperature and magnetic field distribution. We compare the magnetic field parameters extrapolated from the image datacube along each line of sight after folding through the EOVSA instrumental profile with the original (unfolded) parameters used in the model. We find that even the most easily automated, image-based analysis approach (Level 0) provides reasonable quantitative results, although they are affected by systematic effects due to finite sampling in the Fourier (uv) plane. Finally, we note the potential for errors due to misidentified harmonics of the gyrofrequency, and discuss the prospects for applying a more sophisticated spectrally-based analysis scheme (Level 1) to resolve the issue in cases where improved uv coverage and spatial resolution are available.

This paper presents a solely thermal flare, which we detected in the microwave range from the thermal gyro- and free-free emission it produced. An advantage of analyzing thermal gyro emission is its unique ability to precisely yield the magnetic field in the radiating volume. When combined with observationally-deduced plasma density and temperature, these magnetic field measurements offer a straightforward way of tracking evolution of the magnetic and thermal energies in the flare. For the event described here, the magnetic energy density in the radio-emitting volume declines over the flare rise phase, then stays roughly constant during the extended peak phase, but recovers to the original level over the decay phase. At the stage where the magnetic energy density decreases, the thermal energy density increases; however, this increase is insufficient, by roughly an order of magnitude, to compensate for the magnetic energy decrease. When the magnetic energy release is over, the source parameters come back to nearly their original values. We discuss possible scenarios to explain this behavior.

We have undertaken a major enhancement of our IDL-based simulation tools developed earlier for modeling microwave and X-ray emission. The object-based architecture provides an interactive graphical user interface that allows the user to import photospheric magnetic field maps and perform magnetic field extrapolations to almost instantly generate 3D magnetic field models, to investigate the magnetic topology of these models by interactively creating magnetic field lines and associated magnetic flux tubes, to populate the flux tubes with user-defined nonuniform thermal plasma and anisotropic, nonuniform, nonthermal electron distributions; to investigate the spatial and spectral properties of radio and X-ray emission calculated from the model, and to compare the model-derived images and spectra with observational data. The application integrates shared-object libraries containing fast gyrosynchrotron emission codes developed in FORTRAN and C++, soft and hard X-ray codes developed in IDL, a FORTRAN-based potential-field extrapolation routine and an IDL-based linear force free field extrapolation routine. The interactive interface allows users to add any user-defined radiation code that adheres to our interface standards, as well as user-defined magnetic field extrapolation routines. Here we use this tool to analyze a simple single-loop flare and use the model to constrain the 3D structure of the magnetic flaring loop and 3D spatial distribution of the fast electrons inside this loop. We iteratively compute multi-frequency microwave and multi-energy X-ray images from realistic magnetic fluxtubes obtained from an extrapolation of a magnetogram taken prior to the flare, and compare them with imaging data obtained by SDO, NoRH, and RHESSI instruments. We use this event to illustrate use of the tool for general interpretation of solar flares to address disparate problems in solar physics.

Spectra derived from fast Fourier transform (FFT) analysis of time-domain data intrinsically contain statistical fluctuations whose distribution depends on the number of accumulated spectra contributing to a measurement. The tail of this distribution, which is essential for separation of the true signal from the statistical fluctuations, deviates noticeably from the normal distribution for a finite number of the accumulations. In this paper we develop a theory to properly account for the statistical fluctuations when fitting a model to a given accumulated spectrum. The method is implemented in software for the purpose of automatically fitting a large body of such FFT-derived spectra. We apply this tool to analyze a portion of a dense cluster of spikes recorded by our FST instrument during a record-breaking event that occurred on 06 Dec 2006. The outcome of this analysis is briefly discussed.

Currently there is a concern about ability of the classical thermal (Maxwellian) distribution to describe quasi-steady-state plasma in solar atmosphere including active regions. In particular, other distributions have been proposed to better fit observations, for example, kappa- and $n$-distributions. If present, these distributions will generate radio emissions with different observable properties compared with the classical gyroresonance (GR) or free-free emission, which implies a way of remote detecting these non-Maxwellian distributions in the radio observations. Here we present analytically derived GR and free-free emissivities and absorption coefficients for the kappa- and $n$-distributions and discuss their properties, which are in fact remarkably different from each other and from the classical Maxwellian plasma. In particular, the radio brightness temperature from a gyrolayer increases with the optical depth $\tau$ for kappa-distribution, but decreases with $\tau$ for $n$-distribution. This property has a remarkable consequence allowing a straightforward observational test: the gyroresonance radio emission from the non-Maxwellian distributions is supposed to be noticeably polarized even in the optically thick case, where the emission would have strictly zero polarization in the case of Maxwellian plasma. This offers a way of remote probing the plasma distribution in astrophysical sources including solar active regions as a vivid example.

Relativistic antiparticles can be created in high-energy nuclear interactions; thus, detection of antiparticles in an astrophysical source can tell us something remarkable about the underlying high-energy processes and nuclear interactions. However, once created, the antiparticles remain a minor fraction of their conjugant normal particles, so the detection of the antiparticles represents a big science challenge. To address this challenge we employ imaging and polarimetry of microwave radiation produced as the positrons gyrate in the ambient magnetic field. The key property of the radiation used in this method is that the oppositely charged particles, electrons and positrons, produce radiation with opposite helicity, easily distinguishable by currently operating radio facilities. Analysis of available spatially resolved microwave data augmented by independent magnetic field measurements allows us to remotely detect the relativistic positron component in several solar flares.

The thermal-to-nonthermal partition was found to vary greatly from one flare to another resulting in a broad variety of cases from 'heating without acceleration' to 'acceleration without heating'. Recent analysis of microwave data of these differing cases suggests that a similar acceleration mechanism, forming a power-law nonthermal tail up to a few MeV or even higher, operates in all the cases. However, the level of this nonthermal spectrum compared to the original thermal distribution differs significantly from one case to another, implying a highly different thermal-to-nonthermal energy partition in various cases. This further requires a specific mechanism capable of extracting the charged particles from the thermal pool and supplying them to a bulk acceleration process to operate in flares \textit{in addition} to the bulk acceleration process itself, which, in contrast, efficiently accelerates the seed particles, while cannot accelerate the thermal particles. Within this 'microwave' view on the flare energy partition and particle acceleration I present a few contrasting examples of acceleration regions detected with microwave data and compare them with the most popular acceleration mechanisms---in DC fields, in collapsing traps, and stochastic acceleration by a turbulence spectrum - to identify the key elements needed to conform with observations. In particular, I point out that the turbulence needed to drive the particle acceleration is generated in nonpotential magnetic structures, which results in nonzero helicity of the turbulence. This helicity, in its turn, produces a nonzero mean DC electric field on top of stochastic turbulent fields driving the main stochastic acceleration; thus, acceleration by helical turbulence combines properties of the standard stochastic acceleration with some features of acceleration in DC electric fields, exactly what is demanded by observation.

Based on detailed analysis of radio and X-ray observations of a flare on 2002 April 11 augmented by realistic 3D modeling, we have identified a radio emission component produced directly at the flare acceleration region. This acceleration region radio component has distinctly different (i) spectrum, (ii) light curves, (iii) spatial location, and, thus, (iv) physical parameters from those of the separately identified, trapped or precipitating electron components. To derive evolution of physical parameters of the radio sources we apply forward fitting of the radio spectrum time sequence with the gyrosynchrotron source function with 5 to 6 free parameters. At the stage when the contribution from the acceleration region dominates the radio spectrum, the X-ray- and radio- derived electron energy spectral indices agree well with each other. During this time the maximum energy of the accelerated electron spectrum displays a monotonic increase with time from ~ 300 keV to ~ 2 MeV over roughly one minute duration indicative of an acceleration process in the form of growth of the power-law tail; the fast electron residence time in the acceleration region is about 2 - 4 s, which is much longer than the time of flight and so requires a strong diffusion mode there to inhibit free-streaming propagation. The acceleration region has a relatively strong magnetic field, B ~ 120 G, and a low thermal density, less than 2x10^9 cm^-3. These acceleration region properties are consistent with a stochastic acceleration mechanism.

We have developed a general framework for modeling gyrosynchrotron and free-free emission from solar flaring loops and used it to test the premise that 2D maps of source parameters, particularly magnetic field, can be deduced from spatially resolved microwave spectropolarimetry data. In this paper we show quantitative results for a flaring loop with a realistic magnetic geometry, derived from a magnetic field extrapolation, and containing an electron distribution with typical thermal and nonthermal parameters, after folding through the instrumental profile of a realistic interferometric array. We compare the parameters generated from forward fitting a homogeneous source model to each line of sight through the folded image data cube with both the original parameters used in the model and with parameters generated from forward fitting a homogeneous source model to the original (unfolded) image data cube. We find excellent agreement in general, but with systematic effects that can be understood as due to finite resolution in the folded images and the variation of parameters along the line of sight, which are ignored in the homogeneous source model. We discuss the use of such 2D parameter maps within a larger framework of 3D modeling, and the prospects for applying these methods to data from a new generation of multifrequency radio arrays now or soon to be available.

Flaring release of magnetic energy in solar corona is only possible if the magnetic field deviates from a potential one. We show that the linear MHD modes excited on top of the non-potential magnetic field possess a nonzero kinetic helicity. Accordingly, this necessarily results in a noticeable kinetic helicity of the turbulence, composed of these linear modes with various scales and random phases, generated at the flare site by the primary energy release, which may be important for many applications. In particular, a nonzero turbulence helicity has a potentially strong effect on the particle acceleration because the helical component of the turbulence induces a mean regular large-scale (DC) electric field capable of directly accelerating the charged particles in addition to the commonly considered stochastic turbulent electric field. In this paper, we derive the kinetic helicity density of the linear MHD modes excited on top of a twisted large-scale magnetic field, estimate the corresponding turbulence helicity and take its effect on stochastic particle acceleration by the turbulence into consideration; in particular, we compare this induced mean electric field with the electron and estimated effective ion Dreicer fields. We have discovered that this, so far missing but highly important, ingredient of the turbulence at the flare site can be responsible for the thermal-to-nonthermal energy partition in flares by controlling the process of particle extraction from the thermal pool and formation of the seed particle population to be then stochastically accelerated to higher energies. In addition, it is naturally consistent with such puzzling flare manifestations as spatial separation of electron and proton emission sites, electron beam formation, and enrichment of the accelerated particle population by $^3He$ and other rare ions.

In some flares the thermal component appears much earlier than the nonthermal component in X-ray range. Using sensitive microwave observations we revisit this finding made by Battaglia et al. (2009) based on RHESSI data analysis. We have found that nonthermal microwave emission produced by accelerated electrons with energy of at least several hundred keV, appears as early as the thermal soft X-ray emission indicative that the electron acceleration takes place at the very early flare phase. The non-detection of the hard X-rays at that early stage of the flares is, thus, an artifact of a limited RHESSI sensitivity. In all considered events, the microwave emission intensity increases at the early flare phase. We found that either thermal or nonthermal gyrosynchrotron emission can dominate the low-frequency part of the microwave spectrum below the spectral peak occurring at 3-10 GHz. In contrast, the high-frequency optically thin part of the spectrum is always formed by the nonthermal, accelerated electron component, whose power-law energy spectrum can extend up to a few MeV at this early flare stage. This means that even though the total number of accelerated electrons is small at this stage, their nonthermal spectrum is fully developed. This implies that an acceleration process of available seed particles is fully operational. While, creation of this seed population (the process commonly called `injection' of the particles from the thermal pool into acceleration) has a rather low efficiency at this stage, although, the plasma heating efficiency is high. This imbalance between the heating and acceleration (in favor of the heating) is difficult to reconcile within most of available flare energization models. Being reminiscent of the tradeoff between the Joule heating and runaway electron acceleration, it puts additional constraints on the electron injection into the acceleration process.

A quantitative study of the observable radio signatures of the sausage, kink, and torsional MHD oscillation modes in flaring coronal loops is performed. Considering first non-zero order effect of these various MHD oscillation modes on the radio source parameters such as magnetic field, line of sight, plasma density and temperature, electron distribution function, and the source dimensions, we compute time dependent radio emission (spectra and light curves). The radio light curves (of both flux density and degree of polarization) at all considered radio frequencies are than quantified in both time domain (via computation of the full modulation amplitude as a function of frequency) and in Fourier domain (oscillation spectra, phases, and partial modulation amplitude) to form the signatures specific to a particular oscillation mode and/or source parameter regime. We found that the parameter regime and the involved MHD mode can indeed be distinguished using the quantitative measures derived in the modeling. We apply the developed approach to analyze radio burst recorded by Owens Valley Solar Array and report possible detection of the sausage mode oscillation in one (partly occulted) flare and kink or torsional oscillations in another flare.

Microwave emission of solar flares is formed primarily by incoherent gyrosynchrotron radiation generated by accelerated electrons in coronal magnetic loops. The resulting emission depends on many factors, including pitch-angle distribution of the emitting electrons and the source geometry. In this work, we perform systematic simulations of solar microwave emission using recently developed tools (GS Simulator and fast gyrosynchrotron codes) capable of simulating maps of radio brightness and polarization as well as spatially resolved emission spectra. A 3D model of a symmetric dipole magnetic loop is used. We compare the emission from isotropic and anisotropic (of loss-cone type) electron distributions. We also investigate effects caused by inhomogeneous distribution of the emitting particles along the loop. It is found that effect of the adopted moderate electron anisotropy is the most pronounced near the footpoints and it also depends strongly on the loop orientation. Concentration of the emitting particles at the loop top results in a corresponding spatial shift of the radio brightness peak, thus reducing effects of the anisotropy. The high-frequency (around 50 GHz) emission spectral index is specified mainly by the energy spectrum of the emitting electrons; however, at intermediate frequencies (around 10-20 GHz), the spectrum shape is strongly dependent on the electron anisotropy, spatial distribution, and magnetic field nonuniformity. The implications of the obtained results for the diagnostics of the energetic electrons in solar flares are discussed.

We use rotation stereoscopy to estimate the height of a steady-state solar feature relative to the photosphere, based on its apparent motion in the image plane recorded over several days of observation. The stereoscopy algorithm is adapted to work with either one- or two-dimensional data (i.e. from images or from observations that record the projected position of the source along an arbitrary axis). The accuracy of the algorithm is tested on simulated data, and then the algorithm is used to estimate the coronal radio source heights associated with the active region NOAA 10956, based on multifrequency imaging data over 7 days from the Siberian Solar Radio Telescope near 5.7 GHz, the Nobeyama Radio Heliograph at 17 GHz, as well as one-dimensional scans at multiple frequencies spanning the 5.98--15.95 GHz frequency range from the RATAN-600 instrument. The gyroresonance emission mechanism, which is sensitive to the coronal magnetic field strength, is applied to convert the estimated radio source heights at various frequencies, h(f), to information about magnetic field vs. height B(h), and the results are compared to a magnetic field extrapolation derived from photospheric magnetic field observations obtained by Hinode and MDI. We found that the gyroresonant emission comes from the heights exceeding location of the third gyrolayer irrespectively on the magnetic extrapolation method; implications of this finding for the coronal magnetography and coronal plasma physics are discussed.

We report the observation of an unusual cold, tenuous solar flare, which reveals itself via numerous and prominent non-thermal manifestations, while lacking any noticeable thermal emission signature. RHESSI hard X-rays and 0.1-18 GHz radio data from OVSA and Phoenix-2 show copious electron acceleration (10^35 electrons per second above 10 keV) typical for GOES M-class flares with electrons energies up to 100 keV, but GOES temperatures not exceeding 6.1 MK. The imaging, temporal, and spectral characteristics of the flare have led us to a firm conclusion that the bulk of the microwave continuum emission from this flare was produced directly in the acceleration region. The implications of this finding for the flaring energy release and particle acceleration are discussed.

Gyrosynchrotron (GS) emission of charged particles spiraling in magnetic fields plays an exceptionally important role in astrophysics. In particular, this mechanism makes a dominant contribution to the continuum solar and stellar radio emissions. However, the available exact equations describing the emission process are extremely slow computationally, thus limiting the diagnostic capabilities of radio observations. In this work, we present approximate GS codes capable of fast calculating the emission from anisotropic electron distributions. The computation time is reduced by several orders of magnitude compared with the exact formulae, while the computation error remains within a few percent. The codes are implemented as the executable modules callable from IDL; they are made available for users via web sites.

The coming decade will see the routine use of solar data of unprecedented spatial and spectral resolution, time cadence, and completeness. To capitalize on the new (or soon to be available) facilities such as SDO, ATST and FASR, and the challenges they present in the visualization and synthesis of multi-wavelength datasets, we propose that realistic, sophisticated, 3D active region and flare modeling is timely and critical, and will be a forefront of coronal studies over the coming decade. To make such modeling a reality, a broad, concerted effort is needed to capture the wealth of information resulting from the data, develop a synergistic modeling effort, and generate the necessary visualization, interpretation and model-data comparison tools to accurately extract the key physics.

The coming years will see routine use of solar data of unprecedented spatial and spectral resolution, time cadence, and completeness in the wavelength domain. To capitalize on the soon to be available radio facilities such as the expanded OVSA, SSRT and FASR, and the challenges they present in the visualization and synthesis of the multi-frequency datasets, we propose that realistic, sophisticated 3D active region and flare modeling is timely now and will be a forefront of coronal studies over the coming years. Here we summarize our 3D modeling efforts, aimed at forward fitting of imaging spectroscopy data, and describe currently available 3D modeling tools. We also discuss plans for future generalization of our modeling tools.

Fireball model of the gamma-ray bursts (GRBs) predicts generation of numerous internal shocks, which efficiently accelerate charged particles and generate relatively small-scale stochastic magnetic and electric fields. The accelerated particles diffuse in space due to interaction with the random waves and so emit so called Diffusive Synchrotron Radiation (DSR) in contrast to standard synchrotron radiation they would produce in a large-scale regular magnetic fields. In this contribution we present key results of detailed modeling of the GRB spectral parameters, which demonstrate that the non-perturbative DSR emission mechanism in a strong random magnetic field is consistent with observed distributions of the Band parameters and also with cross-correlations between them.

Fireball model of the GRBs predicts generation of numerous internal shocks, which then efficiently accelerate charged particles and generate magnetic and electric fields. These fields are produced in the form of relatively small-scale stochastic ensembles of waves, thus, the accelerated particles diffuse in space due to interaction with the random waves and so emit so called Diffusive Synchrotron Radiation (DSR) in contrast to standard synchrotron radiation they would produce in a large-scale regular magnetic fields. In this paper we present first results of comprehensive modeling of the GRB spectral parameters within the fireball/internal shock concept. We have found that the non-perturbative DSR emission mechanism in a strong random magnetic field is consistent with observed distributions of the Band parameters and also with cross-correlations between them; this analysis allowed to restrict GRB physical parameters from the requirement of consistency between the model and observed distributions.

Observations in the sub-THz range of large solar flares have revealed a mysterious spectral component increasing with frequency and hence distinct from the microwave component commonly accepted to be produced by gyrosynchrotron (GS) emission from accelerated electrons. Evidently, having a distinct sub-THz component requires either a distinct emission mechanism (compared to the GS one), or different properties of electrons and location, or both. We find, however, that the list of possible emission mechanisms is incomplete. This Letter proposes a more complete list of emission mechanisms, capable of producing a sub-THz component, both well-known and new in this context and calculates a representative set of their spectra produced by a) free-free emission, b) gyrosynchrotron emission, c) synchrotron emission from relativistic positrons/electrons, d) diffusive radiation, and e) Cherenkov emission. We discuss the possible role of the mechanisms in forming the sub-THz emission and emphasize their diagnostics potential for flares.

The Letter takes up a question of what radio emission is produced by electrons at the very acceleration site of a solar flare. Specifically, we calculate incoherent radio emission produced within two competing acceleration models--stochastic acceleration by cascading MHD turbulence and regular acceleration in collapsing magnetic traps. Our analysis clearly demonstrates that the radio emission from the acceleration sites: (i) has sufficiently strong intensity to be observed by currently available radio instruments and (ii) has spectra and light curves which are distinctly different in these two competing models, which makes them observationally distinguishable. In particular, we suggest that some of the narrowband microwave and decimeter continuum bursts may be a signature of the stochastic acceleration in solar flares.

We develop a practical forward fitting method based on the SIMPLEX algorithm with shaking, which allows the derivation of the magnetic field and other parameters along a solar flaring loop using microwave imaging spectroscopy of gyrosynchrotron emission. We illustrate the method using a model loop with spatially varying magnetic field, filled with uniform ambient density and an evenly distributed fast electron population with an isotropic, power-law energy distribution.

We propose a nonlinear self-consistent model of the turbulent non-resonant particle acceleration in solar flares. We simulate temporal evolution of the spectra of charged particles accelerated by strong long-wavelength MHD turbulence taking into account back reaction of the accelerated particles on the turbulence. The main finding is that the nonlinear coupling of accelerated particles and MHD turbulence result in prominent evolution of the spectra of accelerated particles, which can be either soft-hard-soft or soft-hard-harder depending on the particle injection efficiency. Such evolution patterns are widely observed in hard X-ray and gamma-ray emission from solar flares.

We describe microwave and hard X-ray observations of strong quasiperiodic pulsations from the GOES X1.3 solar flare on 15 June 2003. The radio observations were made jointly by the Owens Valley Solar Array (OVSA), the Nobeyama Polarimeter (NoRP), and the Nobeyama Radioheliograph (NoRH). Hard X-ray observations were made by the Ramaty High Energy Solar Spectroscopic Imager (RHESSI). Using Fourier analysis, we study the frequency- and energy-dependent oscillation periods, differential phase, and modulation amplitudes of the radio and X-ray pulsations. Focusing on the more complete radio observations, we also examine the modulation of the degree of circular polarization and of the radio spectral index. The observed properties of the oscillations are compared with those derived from two simple models for the radio emission. In particular, we explicitly fit the observed modulation amplitude data to the two competing models. The first model considers the effects of MHD oscillations on the radio emission. The second model considers the quasi-periodic injection of fast electrons. We demonstrate that quasiperiodic acceleration and injection of fast electrons is the more likely cause of the quasiperiodic oscillations observed in the radio and hard X-ray emission, which has important implications for particle acceleration and transport in the flaring sources.

We analyze a dense cluster of solar radio spikes registered at ~ 4.5 -- 6 GHz by the Purple Mountain Observatory spectrometer (Nanjing, China) operating in the 4.5 -- 7.5 GHz range with the 5 ms temporal resolution. To handle with the data from the spectrometer we developed a new technique utilizing a nonlinear multi-Gaussian spectral fit based on chi-squared criteria to extract individual spikes from the originally recorded spectra. Applying this method to the experimental raw data we eventually identified about 3000 spikes for this event, which allows for a detailed statistical analysis. Various statistical characteristics of the spikes have been evaluated, including intensity distributions, spectral bandwidth distributions, and distribution of the spike mean frequencies. The most striking finding of this analysis is distributions of the spike bandwidth, which are remarkably asymmetric. To reveal the underlaying microphysics we explore the local trap model with the renormalized theory of spectral profile of the electron cyclotron maser (ECM) emission peak in a source with random magnetic irregularities. The distribution of the solar spikes relative bandwidth calculated within the local trap model represents an excellent fit to the experimental data. Accordingly, the developed technique may offer a new tool of studying very low levels of the magnetic turbulence in the spike sources, when the ECM mechanism of the spike cluster is confirmed.

Theoretical and experimental study of fast electron beams attracts a lot of attention in the astrophysics and laboratory. In the case of solar flares the problem of reliable beam detection and diagnostics is of exceptional importance. This paper explores the fact that the electron beams moving oblique to the magnetic field or along the field with some angular scatter around the beam propagation direction can generate microwave continuum bursts via gyrosynchrotron mechanism. The characteristics of the microwave bursts produced by beams differ from those in case of isotropic or loss-cone distributions, which suggests a new tool for quantitative diagnostics of the beams in the solar corona. To demonstrate the potentiality of this tool, we analyze here a radio burst occurred during an impulsive flare 1B/M6.7 on 10 March 2001 (AR 9368, N27W42). Based on detailed analysis of the spectral, temporal, and spatial relationships, we obtained firm evidence that the microwave continuum burst is produced by electron beams. For the first time we developed and applied a new forward fitting algorithm based on exact gyrosynchrotron formulae and employing both the total power and polarization measurements to solve the inverse problem of the beam diagnostics. We found that the burst is generated by a oblique beam in a region of reasonably strong magnetic field ($\sim 200-300$ G) and the burst is observed at a quasi-transverse viewing angle. We found that the life time of the emitting electrons in the radio source is relatively short, $\tau_l \approx 0.5$ s, consistent with a single reflection of the electrons from a magnetic mirror at the foot point with the stronger magnetic field. We discuss the implications of these findings for the electron acceleration in flares and for beam diagnostics.

Modern observations and models of various astrophysical objects suggest that many of their physical parameters fluctuate substantially at different spatial scales. The rich variety of the emission processes, including Transition Radiation but not limited to it, arising in such turbulent media constitutes the scope of Stochastic Theory of Radiation. We review general approaches applied in the stochastic theory of radiation and specific methods used to calculate the transition radiation produced by fast particles in the magnetized randomly inhomogeneous plasma. The importance of the theory of transition radiation for astrophysics is illustrated by one example of its detailed application to a solar radio burst, including specially designed algorithms of the spectral forward fitting.

Anisotropic distributions of charged particles including two-stream distributions give rise to generation of either stochastic electric fields (in the form of Langmuir waves, Buneman instability) or random quasi-static magnetic fields (Weibel and filamentation instabilities) or both. These two-stream instabilities are known to play a key role in collisionless shock formation, shock-shock interactions, and shock-induced electromagnetic emission. This paper applies the general non-perturbative stochastic theory of radiation to study electromagnetic emission produced by relativistic particles, which random walk in the stochastic electric fields of the Langmuir waves. This analysis takes into account the cumulative effect of uncorrelated Langmuir waves on the radiating particle trajectory giving rise to angular diffusion of the particle, which eventually modifies the corresponding radiation spectra. We demonstrate that the radiative process considered is probably relevant for emission produced in various kinds of astrophysical jets, in particular, prompt gamma-ray burst spectra, including X-ray excesses and prompt optical flashes.

We calculate spectra of radiation produced by a relativistic particle in the presence of one-dimensional Langmuir turbulence which might be generated by a streaming instability in the plasma, in particular, in the shock front or at the shock-shock interactions. The shape of the radiation spectra is shown to depend sensitively on the angle between the particle velocity and electric field direction. The radiation spectrum in the case of exactly transverse particle motion is degenerate and similar to that of spatially uniform Langmuir oscillations. In case of oblique propagation, the spectrum is more complex, it consists of a number of power-law regions and may contain a distinct high-frequency spectral peak. %at $\omega=2\omega\pe \gamma^2$. The emission process considered is relevant to various laboratory plasma settings and for astrophysical objects as gamma-ray bursts and collimated jets.

We present radio and X-ray observations of an impulsive solar flare that was moderately intense in microwaves, yet showed very meager EUV and X-ray emission. The flare occurred on 2001 Oct 24 and was well-observed at radio wavelengths by the Nobeyama Radioheliograph (NoRH), the Nobeyama Radio Polarimeters (NoRP), and by the Owens Valley Solar Array (OVSA). It was also observed in EUV and X-ray wavelength bands by the TRACE, GOES, and Yohkoh satellites. We find that the impulsive onset of the radio emission is progressively delayed with increasing frequency relative to the onset of hard X-ray emission. In contrast, the time of flux density maximum is progressively delayed with decreasing frequency. The decay phase is independent of radio frequency. The simple source morphology and the excellent spectral coverage at radio wavelengths allowed us to employ a nonlinear chi-squared minimization scheme to fit the time series of radio spectra to a source model that accounts for the observed radio emission in terms of gyrosynchrotron radiation from MeV-energy electrons in a relatively dense thermal plasma. We discuss plasma heating and electron acceleration in view of the parametric trends implied by the model fitting. We suggest that stochastic acceleration likely plays a role in accelerating the radio-emitting electrons.

Diffusive Synchrotron Radiation (DSR) is produced by charged particles as they random walk in a stochastic magnetic field. The spectrum of the radiation produced by particles in such fields differs substantially from those of standard synchrotron emission because the corresponding particle trajectories deviate significantly from gyration in a regular field. The Larmor radius, therefore, is no longer a good measure of the particle trajectory. In this paper we analyze a special DSR regime which arises as highly relativistic electrons move through magnetic fields which have only random structure on a wide range of spatial scales. Such stochastic fields arise in turbulent processes, and are likely present in pulsar wind nebulae (PWNe). We show that DSR generated by a single population of electrons can reproduce the observed broad-band spectra of PWNe from the radio to the X-ray, in particular producing relatively flat spectrum radio emission as is usually observed in PWNe. DSR can explain the existence of several break frequencies in the broad-band emission spectrum without recourse to breaks in the energy spectrum of the relativistic particles. The shape of the radiation spectrum depends on the spatial spectrum of the stochastic magnetic field. The implications of the presented DSR regime for PWN physics are discussed.

The spectrum of electromagnetic emission generated by relativistic electrons scattered on small-scale random magnetic fields, implied by current models of the magnetic field generation in the gamma-ray burst sources, is considered. The theory developed includes both perturbative and non-perturbative versions and, therefore, suggests a general treatment of the radiation in arbitrary small-scale random field. It is shown that a general treatment of the random nature of the small-scale magnetic field, as well as angular diffusion of the electrons due to multiple scattering by magnetic inhomogeneities (i.e., non-perturbative effects), give rise to a radiation spectrum that differs significantly from so-called "jitter" spectrum. The spectrum of diffusive synchrotron radiation seems to be consistent with the low energy spectral index distribution of the gamma-ray bursts.

Flattenings of nonthermal radiation spectra observed from knots and interknot locations of the jets of 3C273 and M87 in UV and X-ray bands are discussed within modern models of magnetic field generation in the relativistic jets. Specifically, we explicitly take into account the effect of the small-scale random magnetic field, probably present in such jets, which gives rise to emission of Diffusive Synchrotron Radiation, whose spectrum deviates substantially from the standard synchrotron spectrum, especially at high frequencies. The calculated spectra agree well with the observed ones if the energy densities contained in small-scale and large-scale magnetic fields are comparable. The implications of this finding for magnetic field generation, particle acceleration, and jet composition are discussed.

The review presents basic elements of the theory of electromagnetic emission by nonthermal particles propagating in stochastic media. In particular, general theory and a few special regimes of Diffusive Synchrotron Radiation (DSR)arising in the presence of random magnetic/electric fields are analyzed in detail. Possible applications of the DSR theory to astrophysical sources as well as other emission processes typical for the random media are discussed briefly.

Microturbulence, i.e. enhanced fluctuations of plasma density, electric and magnetic fields, is of great interest in astrophysical plasmas, but occurs on spatial scales far too small to resolve by remote sensing, e.g., at ~ 1-100 cm in the solar corona. This paper reports spatially resolved observations that offer strong support for the presence in solar flares of a suspected radio emission mechanism, resonant transition radiation, which is tightly coupled to the level of microturbulence and provides direct diagnostics of the existence and level of fluctuations on decimeter spatial scales. Although the level of the microturbulence derived from the radio data is not particularly high, <\Delta n^2 >/n^2 ~ 10^{-5}$, it is large enough to affect the charged particle diffusion and give rise to effective stochastic acceleration. This finding has exceptionally broad astrophysical implications since modern sophisticated numerical models predict generation of much stronger turbulence in relativistic objects, e.g., in gamma-ray burst sources.

The effect of pitch-angle anisotropy of fast electrons on generation of nonthermal radio emission is studied. Incoherent gyrosynchrotron radiation is shown to depend strongly on the anisotropy. In particular, the spectral index of gyrosynchrotron radiation increases up to a factor of 3-4 compared with the isotropic case. The degree of polarization (X-mode) increases for loss-cone distributions, while decreases for beam-like distributions. Ample evidence of the pitch-angle anisotropy effect on (1) spatial distribution of the radio brightness, (2) spatially resolved light curves of the intensity and polarization, (3) spectral hardness of microwave bursts, is found by exploring observations performed with Nobeyama Radioheliograph and Owens Valley Solar Array. Above some threshold in the angular gradient, the electron cyclotron maser instability (coherent gyrosynchrotron emission) can develop, provided that standard gyrosynchrotron emission is accompanied by a lower-frequency intense coherent emission. The coherent emission is studied in detail for a realistic distribution of fast electrons over the energy and pitch-angle. In agreement with observations of narrowband radio spikes, it is found that (1) hard (power-law) distributions over energy are preferable to produce the coherent emission, (2) the threshold of the instability corresponds to quite an anisotropic electron distribution, thus, the pitch-angle anisotropy derived from the properties of the continuum gyrosynchrotron radiation will not necessarily give rise to coherent emission (either enhanced isotropization implied by quasilnear saturation of the instability).