High-energy particles enter the solar atmosphere from Galactic or solar coronal sources, producing an "albedo'' source from the quiet Sun, now observable across a wide range of photon energies. The interaction of high-energy particles in a stellar atmosphere depends essentially upon the joint variation of the magnetic field and the gas, which heretofore has been characterized parametrically as P ~ B^alpha with P the gas pressure and B the magnitude of the magnetic field. We re-examine that parametrization by using a self-consistent 3D MHD model (Bifrost) and show that this relationship tends to P ~ B^{2.9+-0.1} based on the visible portions of the sample of open-field flux tubes in such a model, but with large variations from point to point. This scatter corresponds to the strong meandering of the open-field flux tubes in the lower atmosphere, which will have a strong effect on the prediction of the emission anisotropy (limb brightening). The simulations show that much of the open flux in coronal holes originates in weak-field regions within the granular pattern of the convective motions seen in the simulations.

We use the Monte Carlo particle physics code FLUKA (Fluktuierende Kaskade) to calculate $\gamma$-ray spectra expected from solar flare energetic ion distributions. The FLUKA code includes robust physics-based models for electromagnetic, hadronic and nuclear interactions, sufficiently detailed for it to be a useful tool for calculating nuclear de-excitation, positron annihilation and neutron capture line fluxes and shapes, as well as $\approx \; {\rm GeV}$ continuum radiation from pion decay products. We show nuclear de-excitation $\gamma$-ray line model spectra from a range of assumed primary accelerated ion distributions and find them to be in good agreement with those found using the code of Murphy et al. (2009). We also show full $\gamma$-ray model spectra which exhibit all the typical structures of $\gamma$-ray spectra observed in solar flares. From these model spectra we build templates which are incorporated into the software package Objective Spectral Executive (OSPEX) and used to fit the combined Fermi Gamma-ray Burst Monitor (GBM)/Large Area Telescope (LAT) spectrum of the 2010 June 12 solar flare, providing a statistically acceptable result. To the best of our knowledge, the fit carried out with the FLUKA templates for the full $\gamma$-ray spectrum can be regarded as the first attempt to use a single code to implement a self-consistent treatment of the several spectral components in the photon energy range from $\approx 100$s ${\rm keV}$ to $\approx 100$s ${\rm MeV}$.

Solar observations in the infrared domain can bring important clues on the response of the low solar atmosphere to primary energy released during flares. At present the infrared continuum has been detected at 30 THz (10 $\mu$m) in only a few flares. In this work we present a detailed multi-frequency analysis of SOL2012-03-13, including observations at radio millimeter and sub-millimeter wavelengths, in hard X-rays (HXR), gamma-rays (GR), H-alpha, and white-light. HXR/GR spectral analysis shows that the event is a GR line flare and allows estimating the numbers of and energy contents in electrons, protons and alpha particles produced during the flare. The energy spectrum of the electrons producing the HXR/GR continuum is consistent with a broken power-law with an energy break at ~800 keV. It is shown that the high-energy part (above ~800 keV) of this distribution is responsible for the high-frequency radio emission (> 20 GHz) detected during the flare. By comparing the 30 THz emission expected from semi-empirical and time-independent models of the quiet and flare atmospheres, we find that most (~80%) of the observed 30 THz radiation can be attributed to thermal free-free emission of an optically-thin source. Using the F2 flare atmospheric model this thin source is found to be at temperatures T~8000 K and is located well above the minimum temperature region. We argue that the chromospheric heating, which results in 80% of the 30 THz excess radiation, can be due to energy deposition by non-thermal flare accelerated electrons, protons and alpha particles. The remaining 20% of the 30 THz excess emission is found to be radiated from an optically-thick atmospheric layer at T~5000 K, below the temperature minimum region, where direct heating by non-thermal particles is insufficient to account for the observed infrared radiation.

Context. The Sun is an active source of radio emission that is often associated with energetic phenomena ranging from nanoflares to coronal mass ejections (CMEs). At low radio frequencies (<100 MHz), numerous millisecond duration radio bursts have been reported, such as radio spikes or solar S bursts (where S stands for short). To date, these have neither been studied extensively nor imaged because of the instrumental limitations of previous radio telescopes. Aims. Here, Low Frequency Array (LOFAR) observations were used to study the spectral and spatial characteristics of a multitude of S bursts, as well as their origin and possible emission mechanisms. Methods. We used 170 simultaneous tied-array beams for spectroscopy and imaging of S bursts. Since S bursts have short timescales and fine frequency structures, high cadence (~50 ms) tied-array images were used instead of standard interferometric imaging, that is currently limited to one image per second. Results. On 9 July 2013, over 3000 S bursts were observed over a time period of ~8 hours. S bursts were found to appear as groups of short-lived (<1 s) and narrow-bandwidth (~2.5 MHz) features, the majority drifting at ~3.5 MHz/s and a wide range of circular polarisation degrees (2-8 times more polarised than the accompanying Type III bursts). Extrapolation of the photospheric magnetic field using the potential field source surface (PFSS) model suggests that S bursts are associated with a trans-equatorial loop system that connects an active region in the southern hemisphere to a bipolar region of plage in the northern hemisphere. Conclusions. We have identified polarised, short-lived solar radio bursts that have never been imaged before. They are observed at a height and frequency range where plasma emission is the dominant emission mechanism, however they possess some of the characteristics of electron-cyclotron maser emission.

The Sun is an active source of radio emission which is often associated with energetic phenomena such as solar flares and coronal mass ejections (CMEs). At low radio frequencies (<100 MHz), the Sun has not been imaged extensively because of the instrumental limitations of previous radio telescopes. Here, the combined high spatial, spectral and temporal resolution of the Low Frequency Array (LOFAR) was used to study solar Type III radio bursts at 30-90 MHz and their association with CMEs. The Sun was imaged with 126 simultaneous tied-array beams within 5 solar radii of the solar centre. This method offers benefits over standard interferometric imaging since each beam produces high temporal (83 ms) and spectral resolution (12.5 kHz) dynamic spectra at an array of spatial locations centred on the Sun. LOFAR's standard interferometric output is currently limited to one image per second. Over a period of 30 minutes, multiple Type III radio bursts were observed, a number of which were found to be located at high altitudes (4 solar radii from the solar center at 30 MHz) and to have non-radial trajectories. These bursts occurred at altitudes in excess of values predicted by 1D radial electron density models. The non-radial high altitude Type III bursts were found to be associated with the expanding flank of a CME. The CME may have compressed neighbouring streamer plasma producing larger electron densities at high altitudes, while the non-radial burst trajectories can be explained by the deflection of radial magnetic fields as the CME expanded in the low corona.

This paper reports on a search for flare emission via charge-exchange radiation in the wings of the Lyman-alpha line of He ii at 304 A, as originally suggested for hydrogen by Orrall and Zirker. Via this mechanism a primary alpha particle that penetrates into the neutral chromosphere can pick up an atomic electron and emit in the He ii bound-bound spectrum before it stops. The Extreme-ultraviolet Variability Experiment (EVE) on board the Solar Dynamics Observatory (SDO) gives us our first chance to search for this effect systematically. The Orrall-Zirker mechanism has great importance for flare physics because of the essential roles that particle acceleration plays; this mechanism is one of the few proposed that would allow remote sensing of primary accelerated particles below a few MeV/nucleon. We study ten events in total, including the gamma-ray events SOL2010-06-12 (M2.0) and SOL2011-02-24 (M3.5) (the latter a limb flare), seven X-class flares, and one prominent M-class event that produced solar energetic particles (SEPs). The absence of charge-exchange line wings may point to a need for more complete theoretical work. Some of the events do have broad-band signatures, which could correspond to continua from other origins, but these do not have the spectral signatures expected from the Orrall-Zirker mechanism.

We present the numerical simulations for an electron-beam-driven and loss-cone-driven electron-cyclotron maser (ECM) with different plasma parameters and different magnetic field strengths for a relatively small region and short time-scale in an attempt to interpret the recent discovered intense radio emission from ultracool dwarfs. We find that a large amount of electromagnetic field energy can be effectively released from the beam-driven ECM, which rapidly heats the surrounding plasma. A rapidly developed high-energy tail of electrons in velocity space (resulting from the heating process of the ECM) may produce the radio continuum depending on the initial strength of the external magnetic field and the electron beam current. Both significant linear polarization and circular polarization of electromagnetic waves can be obtained from the simulations. The spectral energy distributions of the simulated radio waves show that harmonics may appear from 10 to 70$\nu_{\rm pe}$ ($\nu_{\rm pe}$ is the electron plasma frequency) in the non-relativistic case and from 10 to 600$\nu_{\rm pe}$ in the relativistic case, which makes it difficult to find the fundamental cyclotron frequency in the observed radio frequencies. A wide frequency band should therefore be covered by future radio observations.

X-ray observations are a powerful diagnostic tool for transport, acceleration, and heating of electrons in solar flares. Height and size measurements of X-ray footpoints sources can be used to determine the chromospheric density and constrain the parameters of magnetic field convergence and electron pitch-angle evolution. We investigate the influence of the chromospheric density, magnetic mirroring and collisional pitch-angle scattering on the size of X-ray sources. The time-independent Fokker-Planck equation for electron transport is solved numerically and analytically to find the electron distribution as a function of height above the photosphere. From this distribution, the expected X-ray flux as a function of height, its peak height and full width at half maximum are calculated and compared with RHESSI observations. A purely instrumental explanation for the observed source size was ruled out by using simulated RHESSI images. We find that magnetic mirroring and collisional pitch-angle scattering tend to change the electron flux such that electrons are stopped higher in the atmosphere compared with the simple case with collisional energy loss only. However, the resulting X-ray flux is dominated by the density structure in the chromosphere and only marginal increases in source width are found. Very high loop densities (>10^{11} cm^{-3}) could explain the observed sizes at higher energies, but are unrealistic and would result in no footpoint emission below about 40 keV, contrary to observations. We conclude that within a monolithic density model the vertical sizes are given mostly by the density scale-height and are predicted smaller than the RHESSI results show.

A fundamental problem in astrophysics is the interaction between magnetic turbulence and charged particles. It is now possible to use \emph{Ramaty High Energy Solar Spectroscopic Imager (RHESSI)} observations of hard X-rays (HXR) emitted by electrons to identify the presence of turbulence and to estimate the magnitude of the magnetic field line diffusion coefficient at least in dense coronal flaring loops.} {We discuss the various possible regimes of cross-field transport of non-thermal electrons resulting from broadband magnetic turbulence in coronal loops. The importance of the Kubo number $K$ as a governing parameter is emphasized and results applicable in both the large and small Kubo number limits are collected.} {Generic models, based on concepts and insights developed in the statistical theory of transport, are applied to the coronal loops and to the interpretation of hard X-ray imaging data in solar flares. The role of trapping effects, which become important in the non-linear regime of transport, is taken into account in the interpretation of the data.} For this flaring solar loop, we constrain the ranges of parallel and perpendicular correlation lengths of turbulent magnetic fields and possible Kubo numbers. We show that a substantial amount of magnetic fluctuations with energy $\sim 1%$ (or more) of the background field can be inferred from the measurements of the magnetic diffusion coefficient inside thick-target coronal loops.

Gamma-rays and neutrons are the only sources of information on energetic ions present during solar flares and on properties of these ions when they interact in the solar atmosphere. The production of {\gamma}-rays and neutrons results from convolution of the nuclear cross-sections with the ion distribution functions in the atmosphere. The observed {\gamma}-ray and neutron fluxes thus provide useful diagnostics for the properties of energetic ions, yielding strong constraints on acceleration mechanisms as well as properties of the interaction sites. The problem of ion transport between the accelerating and interaction sites must also be addressed to infer as much information as possible on the properties of the primary ion accelerator. In the last couple of decades, both theoretical and observational developments have led to substantial progress in understanding the origin of solar {\gamma}-rays and neutrons. This chapter reviews recent developments in the study of solar {\gamma}-rays and of solar neutrons at the time of the RHESSI era. The unprecedented quality of the RHESSI data reveals {\gamma}-ray line shapes for the first time and provides {\gamma}-ray images. Our previous understanding of the properties of energetic ions based on measurements from the former solar cycles is also summarized. The new results-obtained owing both to the gain in spectral resolution (both with RHESSI and with the non solar-dedicated INTEGRAL/SPI instrument) and to the pioneering imaging technique in the {\gamma}-ray domain-are presented in the context of this previous knowledge. Still open questions are emphasized in the last section of the chapter and future perspectives on this field are briefly discussed.

Recently unanticipated magnetic activity in ultracool dwarfs (UCDs, spectral classes later than M7) have emerged from a number of radio observations. The highly (up to 100%) circularly polarized nature and high brightness temperature of the emission has been interpreted as an effective amplification mechanism of the high-frequency electromagnetic waves, the electron cyclotron maser instability (ECMI). In order to understand the magnetic topology and the properties of the radio emitting region and associated plasmas in these ultracool dwarfs and interpret the origin of radio pulses and their radiation mechanism, we built an active region model, based on the rotation of the UCD and the ECMI mechanism. ECMI mechanism is responsible for the radio bursts from the magnetic tubes and the rotation of the dwarf can modulate the integral of flux with respect to time. The high degree of variability in the brightness and the diverse profile of pulses can be interpreted in terms of a large-scale hot active region with extended magnetic structure existing in the magnetosphere of TVLM 513-46546. We suggest the time profile of the radio light curve is in the form of power law in the model. The radio emitting region consists of complicated substructure. With this model, we can determine the nature (e.g. size, temperature, density) of the radio emitting region and plasma. The magnetic topology can also be constrained. We compare our predicted X-ray flux with Chandra X-ray observation of TVLM 513-46546. Although the X-ray detection is only marginally significant, our predicted flux is significantly lower than the observed flux. We suggest more observations at multi-wavelength will help us understand the magnetic field structure and plasma behavior on the ultracool dwarf.

Most theoretical descriptions of the production of solar flare bremsstrahlung radiation assume the collision of dilute accelerated particles with a cold, dense target plasma, neglecting interactions of the fast particles with each other. This is inadequate for situations where collisions with this background plasma are not completely dominant, as may be the case in, for example, low-density coronal sources. We aim to formulate a model of a self-interacting, entirely fast electron population in the absence of a dense background plasma, to investigate its implications for observed bremsstrahlung spectra and the flare energy budget. We derive approximate expressions for the time-dependent distribution function of the fast electrons using a Fokker-Planck approach. We use these expressions to generate synthetic bremsstrahlung X-ray spectra as would be seen from a corresponding coronal source. We find that our model qualitatively reproduces the observed behaviour of some flares. As the flare progresses, the model's initial power-law spectrum is joined by a lower energy, thermal component. The power-law component diminishes, and the growing thermal component proceeds to dominate the total emission over timescales consistent with flare observations. The power-law exhibits progressive spectral hardening, as is seen in some flare coronal sources. We also find that our model requires a factor of 7 - 10 fewer accelerated electrons than the cold, thick target model to generate an equivalent hard X-ray flux. This model forms the basis of a treatment of self-interactions among flare fast electrons, a process which affords a more efficient means to produce bremsstrahlung photons and so may reduce the efficiency requirements placed on the particle acceleration mechanism. It also provides a useful description of the thermalisation of fast electrons in coronal sources.

In solar flares, inverse Compton scattering (ICS) of photospheric photons might give rise to detectable hard X-ray photon fluxes from the corona where ambient densities are too low for significant bremsstrahlung or recombination. Gamma-ray lines and continuum in some large flares imply the presence of the necessary ~100 MeV electrons and positrons, the latter as by-products of GeV energy ions. Recent observations of coronal hard X-ray sources in particular prompt us to reconsider here the possible contribution of ICS. We aim to evaluate the ICS X-ray fluxes to be expected from prescribed populations of relativistic electrons and positrons in the solar corona. The ultimate aim is to determine if ICS coronal X-ray sources might offer a new diagnostic window on relativistic electrons and ions in flares. We use the complete formalism of ICS to calculate X-ray fluxes from possible populations of flare primary electrons and secondary positrons, paying attention to the incident photon angular distribution near the solar surface and thus improving on the assumption of isotropy made in previous solar discussions. Both primary electrons and secondary positrons produce very hard ICS X-ray spectra. The anisotropic primary radiation field results in pronounced centre-to-limb variation in predicted fluxes and spectra, with the most intense spectra, extending to the highest photon energies, expected from limb flares. Acceptable numbers of electrons or positrons could account for RHESSI coronal X/gamma-ray sources. Some coronal X-ray sources at least might be interpreted in terms of ICS by relativistic electrons or positrons, particularly when sources appear at such low ambient densities that bremsstrahlung appears implausible.

The collisional thick target model (CTTM) of solar hard X-ray (HXR) bursts has become an almost 'Standard Model' of flare impulsive phase energy transport and radiation. However, it faces various problems in the light of recent data, particularly the high electron beam density and anisotropy it involves.} {We consider how photon yield per electron can be increased, and hence fast electron beam intensity requirements reduced, by local re-acceleration of fast electrons throughout the HXR source itself, after injection.} {We show parametrically that, if net re-acceleration rates due to e.g. waves or local current sheet electric (${\cal E}$) fields are a significant fraction of collisional loss rates, electron lifetimes, and hence the net radiative HXR output per electron can be substantially increased over the CTTM values. In this local re-acceleration thick target model (LRTTM) fast electron number requirements and anisotropy are thus reduced. One specific possible scenario involving such re-acceleration is discussed, viz, a current sheet cascade (CSC) in a randomly stressed magnetic loop.} {Combined MHD and test particle simulations show that local ${\cal E}$ fields in CSCs can efficiently accelerate electrons in the corona and and re-accelerate them after injection into the chromosphere. In this HXR source scenario, rapid synchronisation and variability of impulsive footpoint emissions can still occur since primary electron acceleration is in the high Alfv\'{e}n speed corona with fast re-acceleration in chromospheric CSCs. It is also consistent with the energy-dependent time-of-flight delays in HXR features.

The magnetic field of the solar corona has a large-scale dipole character, which maps into the bipolar field in the solar wind. Using standard representations of the coronal field, we show that high-energy ions can be trapped stably in these large-scale closed fields. The drift shells that describe the conservation of the third adiabatic invariant may have complicated geometries. Particles trapped in these zones would resemble the Van Allen Belts and could have detectable consequences. We discuss potential sources of trapped particles.

A novel method of using hard X-rays as a diagnostic for chromospheric density and magnetic structures is developed to infer sub-arcsecond vertical variation of magnetic flux tube size and neutral gas density.Using Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) X-ray data and the newly developed X-ray visibilities forward fitting technique we find the FWHM and centroid positions of hard X-ray sources with sub-arcsecond resolution ($\sim 0.2"$) for a solar limb flare. We show that the height variations of the chromospheric density and the magnetic flux densities can be found with unprecedented vertical resolution of $\sim$ 150 km by mapping 18-250 keV X-ray emission of energetic electrons propagating in the loop at chromospheric heights of 400-1500 km. Our observations suggest that the density of the neutral gas is in good agreement with hydrostatic models with a scale height of around $140\pm 30$ km. FWHM sizes of the X-ray sources decrease with energy suggesting the expansion (fanning out) of magnetic flux tube in the chromosphere with height. The magnetic scale height $B(z)(dB/dz)^{-1}$ is found to be of the order of 300 km and strong horizontal magnetic field is associated with noticeable flux tube expansion at a height of $\sim$ 900 km.

Particle acceleration consequences from fluctuating electric fields superposed on an X-type magnetic field in collisionless solar plasma are studied. Such a system is chosen to mimic generic features of dynamic reconnection, or the reconnective dissipation of a linear disturbance. We explore numerically the consequences for charged particle distributions of fluctuating electric fields superposed on an X-type magnetic field. Particle distributions are obtained by numerically integrating individual charged particle orbits when a time varying electric field is superimposed on a static X-type neutral point. This configuration represents the effects of the passage of a generic MHD disturbance through such a system. Different frequencies of the electric field are used, representing different possible types of wave. The electric field reduces with increasing distance from the X-type neutral point as in linear dynamic magnetic reconnection. The resulting particle distributions have properties that depend on the amplitude and frequency of the electric field. In many cases a bimodal form is found. Depending on the timescale for variation of the electric field, electrons and ions may be accelerated to different degrees and often have energy distributions of different forms. Protons are accelerated to $\gamma$-ray producing energies and electrons to and above hard X-ray producing energies in timescales of 1 second. The acceleration mechanism is possibly important for solar flares and solar noise storms but is also applicable to all collisionless plasmas.

The observed hard X-ray (HXR) flux spectrum $I(\epsilon)$ from solar flares is a combination of primary bremsstrahlung photons $I_P(\epsilon)$ with a spectrally modified component from photospheric Compton backscatter of downward primary emission. The latter can be significant, distorting or hiding the true features of the primary spectrum which are key diagnostics for acceleration and propagation of high energy electrons and of their energy budget. For the first time in solar physics, we use a Green's function approach to the backscatter spectral deconvolution problem, constructing a Green's matrix including photoelectric absorption. This approach allows spectrum-independent extraction of the primary spectrum for several HXR flares observed by the {\it Ramaty High Energy Solar Spectroscopic Imager} (RHESSI). We show that the observed and primary spectra differ very substantially for flares with hard spectra close to the disk centre. We show in particular that the energy dependent photon spectral index $\gamma (\epsilon)=-d \log I/d \log \epsilon$ is very different for $I_P(\epsilon)$ and for $I(\epsilon)$ and that inferred mean source electron spectra ${\bar F}(E)$ differ greatly. Even for a forward fitting of a parametric ${\bar F}(E)$ to the data, a clear low-energy cutoff required to fit $I(\epsilon)$ essentially disappears when the fit is to $I_P(\epsilon)$ - i.e. when albedo correction is included. The self-consistent correction for backscattered photons is thus shown to be crucial in determining the energy spectra of flare accelerated electrons, and hence their total number and energy.

The deduction from solar flare X-ray photon spectroscopic data of the energy dependent model-independent spectral index is considered as an inverse problem. Using the well developed regularization approach we analyze the energy dependency of spectral index for a high resolution energy spectrum provided by Ramaty High Energy Solar Spectroscopic Imager (RHESSI). The regularization technique produces much smoother derivatives while avoiding additional errors typical of finite differences. It is shown that observations imply a spectral index varying significantly with energy, in a way that also varies with time as the flare progresses. The implications of these findings are discussed in the solar flare context.

Finite thermal velocity modifications to electron slowing-down rates may be important for the deduction of solar flare total electron energy. Here we treat both slowing-down and velocity diffusion of electrons in the corona at flare temperatures, for the case of a simple, spatially homogeneous source. Including velocity diffusion yields a consistent treatment of both `accelerated' and `thermal' electrons. It also emphasises that one may not invoke finite thermal velocity target effects on electron lifetimes without simultaneously treating the contribution to the observed X-ray spectrum from thermal electrons. We present model calculations of the X-ray spectra resulting from injection of a power-law energy distribution of electrons into a source with finite temperature. Reducing the power-law distribution low-energy cutoff to lower and lower energies only increases the relative magnitude of the thermal component of the spectrum, because the lowest energy electrons simply join the background thermal distribution. Acceptable fits to RHESSI flare data are obtained using this model. These also demonstrate, however, that observed spectra may in consequence be acceptably consistent with rather a wide range of injected electron parameters.

Possible entropy constraints on particle acceleration spectra are discussed. Solar flare models invoke a variety of initial distributions of the primary energy release over the particles of the flare plasma -- ie., the partition of the energy between thermal and nonthermal components. It is suggested that, while this partition can take any value as far as energy is concerned, the entropy of a particle distribution may provide a useful measure of the likelihood of its being produced for a prescribed total energy. The Gibbs' entropy is calculated for several nonthermal isotropic distribution functions $f$, for a single particle species, and compared with that of a Maxwellian, all distributions having the same total number and energy of particles. Speculations are made on the relevance of some of the results to the cosmic ray power-law spectrum, on their relation to the observed frequency distribution of nonthermal flare hard X-ray spectrum parameters and on the additional energy release required to achieve lower entropy $f$s.