Issue 
A&A
Volume 609, January 2018



Article Number  A41  
Number of page(s)  9  
Section  The Sun  
DOI  https://doi.org/10.1051/00046361/201730546  
Published online  05 January 2018 
Radio signatures of shockaccelerated electron beams in the solar corona
^{1} LeibnizInstitut für Astrophysik Potsdam (AIP), An der Sternwarte 16, 14482 Potsdam, Germany
email: GMann@aip.de
^{2} National Academy of Sciences of Ukraine, Institute of Radio Astronomy of NAS, 4 Chervonopraporna Str., 61002 Kharkov, Ukraine
^{3} Commission for Astronomy, Austrian Academy of Sciences, Schmiedlstrasse 6, 8042 Graz, Austria
^{4} Poltava Gravimetrical Observatory of Institute of Geophysics of NASU, Myasoedova, 27/29, 36014 Poltava, Ukraine
Received: 2 February 2017
Accepted: 11 November 2017
Context. The Sun’s activity can appear in terms of radio bursts. In the frequency range 8−33 MHz the radio telescope URAN2 observed special fine structures appearing as a chain of stripes of enhanced radio emission in the dynamic radio spectrum. The chain drifts slowly from 26 to 23 MHz within 4 min. The individual structures consist of a “head” at the highfrequency edge and a “tail” rapidly drifting from the “head” to lower frequencies over an extent of ≈10 MHz within 8 s. Since they resemble the wellknown “herring bones” in type II radio bursts, they are interpreted as shock accelerated electron beams.
Aims. The electron beams generating these fine structures are considered to be produced by shock drift acceleration (SDA). The beam electrons excite Langmuir waves which are converted into radio waves by nonlinear waveplasma processes. That is called plasma emission. The aim of this paper is to link the radio spectral data of these fine structures to the theoretical results in order to gain a better understanding of the generation of energetic electrons by shocks in the solar corona.
Methods. Adopting SDA for generating energetic electrons, the accelerated electrons establish a beamlike velocity distribution. Plasma emission requires the excitation of Langmuir waves, which is efficient if the velocity of the beam electrons exceeds a few times thermal electron speed. That is the case if the angle between the shock normal and the upstream magnetic field is nearly perpendicular. Hence, the RankineHugoniot relationships, which describe the shock transition in the framework of magnetohydrodynamics, are evaluated for the special case of nearly perpendicular shocks under coronal circumstances.
Results. The radio data deduced from the dynamic radio spectrum can be related in the best way to the theoretical results, if the electron beams, which generate these fine structures, are generated via SDA at an almost perpendicular shock, which is traveling nearly horizontally to the surface of the Sun.
Key words: Sun: flares / Sun: radio radiation / Sun: particle emission / shock waves / acceleration of particles / plasmas
© ESO, 2017
1. Introduction
The Sun’s activity appears not only in the wellknown 11yr sunspot cycle but also in eruptive events such as flares and coronal mass ejections, as well as radio bursts. All these events are accompanied by an enhanced radio emission. Therefore the study of the Sun’s radio radiation provides important information on plasma processes associated with the Sun’s activity.
On April 7, 2011, special fine structures (see Fig. 1) were recorded in the solar radio radiation in the frequency range 8–33 MHz with the radio spectrometer at the radio telescope URAN2 of Poltava Gravimetrical Observatory of Institute of Geophysics of National Academy of Sciences of Ukraine (Brazhenko et al. 2005). Figure 1 shows the dynamic radio spectrum of these fine structures (Brazhenko et al. 2012).
Fig. 1
Dynamic radio spectrum in the frequency range 8–33 MHz as recorded during the period 11:43–11:50 UT on 7 April 2011 with the grayscaled intensity (top), which is given in solar flux units (1 s.f.u. = 10^{22} W m^{2} Hz^{1}), and degree of circular polarization (bottom). 
They appear as a chain of stripes of enhanced radio emission with a pronounced circular polarization. The chain drifts slowly from high towards low frequencies. The individual structures consist of two components namely the “head” at the highfrequency edge and the “tail” extending to lower frequencies (see also Fig. 2). In other events in the decameter range, sometimes opposite cases are observed, that is, the heads appear at the lowfrequency edge, whereas the tails drift rapidly towards higher frequencies (V. N. Melnik, priv. comm.). In the case discussed here, the high frequency edges of these heads drift slowly from 26 MHz towards 23 MHz in the period 11:44−11:48 UT (see Fig. 1) leading to a drift rate D_{f,sd} = −12.5 kHz/s. The tails extend over a frequency range of ≈10 MHz with drift rates of about D_{f,tail} = −1.25 MHz/s (see Fig. 2).
These structures resemble the “herring bones” that appear as fine structures in solar typeII radio bursts. In dynamic radio spectra, solar typeII radio bursts manifest themselves as stripes of enhanced radio emission slowly drifting from high to low frequencies and are considered as the radio signature of shock waves traveling through the corona (Smerd et al. 1962; Wild & Smerd 1972; see as reviews e.g. Nelson & Melrose 1985; and Mann 1995). Approximately, 21% of all typeII radio bursts show these socalled herring bones (Cane & White 1989). They are rapidly drifting emission stripes shooting away from the slowly drifting “backbone” towards lower and higher frequencies (see e.g., Nelson & Melrose 1985; Cairns & Robinson 1987). Herring bones are regarded as the radio signature of electron beams generated by the shock wave associated with the backbone (see e.g., Nelson & Melrose 1985; Mann 1995).
The fine structures shown in Fig. 1 can be considered as individual herring bones due to the following reasons.

Cairns & Robinson (1987) described herring bones (HB) indynamic radio spectra: “... while the duration of a typeIII burstinvariably increases towards lower frequencies, resulting in afanlike shape, a HB burst has a longer duration nearer the startingfrequency, thereby resembling a thorn. This observation allowsus to safely distinguish normal HB bursts from the typeIIIlikebursts ...” (see Cairns & Robinson 1987, p. 367). The finestructures discussed here show the same spectral morphology,that is, a thornlike shape (or in other words a headtail structure).

HBs exhibit a pronounced circular polarization (Stewart & Magun 1980; Suzuki et al. 1982; Cairns & Robinson 1987; Dorovskyy et al. 2015). The fundamental radiation of HBs is usually more circularly polarized than the harmonic one (Cairns & Robinson 1987). The fine structures discussed here show also a significantly circular polarization as seen at the bottom of Fig. 1.

The slow drift of –12.5 kHz/s is typical for typeII bursts in the decameter range (see e.g., Mann 1995; Melnik et al. 2004; Dorovskyy et al. 2015).

The fast drift rate of the tails of –1.25 MHz/s is typical of HBs in the decameter range (Melnik et al. 2004; Dorovskyy et al. 2015). We note that the drift rate of typeIII bursts is roughly about –2.36 MHz/s at 20 MHz according to Alvarez & Haddock (1973). Typically, HBs have a lower drift rate than typeIII bursts in the same frequency range (Mann & Klassen 2002).
Since the fine structures discussed here show all typical features of HBs in the decameter range, it is justified to also consider them as the radio signature of shockaccelerated electron beams (Nelson & Melrose 1985; Cairns & Robinson 1987; Zlobec et al. 1993; Mann 1995, 2006; Aurass 1996; Cairns et al. 2003; Nindos et al. 2008; Dorovsky et al. 2015). Hence, the slow drift of the highfrequency edge of the chain of heads and the rapid drift of the tails represent the radio signature of the associated shock wave propagating outwards in the corona and the shock accelerated electron beams shooting away from the shock, respectively.
The structures of discussion are superimposed on a longlasting typeIV continuum accompanied by a behindlimb coronal mass ejection (CME). It is not unusual for typeII bursts to occur with a typeIV continuum as reported by Melnik et al. (2008) and Dorovskyy et al. (2015).
Fig. 2
Sketch of the shape of the fine structures discussed here. 
In the solar corona, shock waves can be generated either by the blast wave due to the pressure pulse of a flare (Uchida et al. 1973; Vršnak et al. 1995; Magdalenic et al. 2012) or driven by a CME (Stewart et al. 1974a, 1974b; Melnik et al. 2008; Dorovskyy et al. 2015); they can be observed as typeII radio bursts in the metric and decameter wave ranges (Wild & McCready 1950; Uchida 1960; Melnik et al. 2004; Dorovskyy et al. 2015).
The mechanism by which electrons are accelerated up to high energies is currently under debate, since energetic electrons are the source of the nonthermal radio and Xray radiation. The aim of this paper is to examine what can be learnt from these radio observations for the generation mechanism of energetic electrons at shocks in space. Here, we assume that the electrons are accelerated by shock drift acceleration (SDA; see e.g., Holman & Pesses 1983; Ball & Melrose 2001; Mann & Klassen 2005).
The radio telescope URAN2 is briefly described in Sect. 2, where the observational data deduced from the dynamic radio spectrum in Fig. 1 are also presented. Data analysis is presented in Sect. 3. The necessary conditions of radio emission via plasma emission under coronal circumstances are discussed in Sect. 4. In Sect. 5, the SDA is briefly introduced and the properties of shock accelerated electron beams are derived. The link between the observational data and the theoretical results is discussed in Sect. 6. A Summary is presented in Sect. 7. A summary of the RankineHugoniot relationships describing the relation between the physical quantities in the up and downstream regions of a shock is given in Appendix A.
2. Observations
The observations were performed with the radio telescope URAN2 (Brazhenko et al. 2005) of the Poltava Gravimetrical Observatory of the Geophysics Institute in the National Academy of Sciences of Ukraine. The telescope was built near Poltava (Ukraine) in 1992 and is part of the allUkrainian interferometer URAN (Ukrainian Radio Interferometer of National Academy of Sciences). The radio telescope consists of 512 broadband dipoles. The telescope area comprises 238 m × 118 m (corresponding to 28 000 m^{2}). The beam of the radio telescope is 7°× 3° at 25 MHz and is operated electronically. The working frequency band is from 8 to 33 MHz. The radio telescope URAN2 has the capacity to measure polarization of radio emission. It is used both as a part of the interferometer for observations of extragalactic objects such as quasars and radio galaxies, and as an autonomous instrument for observations of the Sun and Jupiter. Observations of solar radio emission are carried out daily at intervals of approximately 3 h and more centered to the meridian. The tracking of the Sun is automatic with a switch frequency of once every four minutes. DSP (Digital Spectral Polarimeter; Ryabov et al. 2010) with time resolution 1–100 ms and frequency resolution 4 kHz was used in solar observations during 2011. The dynamic range of DSP is 90 dB. The sensitivity of the radio telescope URAN2 is about 500 Jy during solar observations.
With this instrument, special fine structures are observed on April 7, 2011, as presented in Fig. 1. In dynamic spectra, the morphology of such individual structures is illustrated in Fig. 2. These consist of a “head” at its high frequency edge and a “tail” rapidly drifting towards lower frequencies. The high frequency edges (heads) of these structures drift from 26 MHz to 23 MHz in the period 11:44–11:48 UT, that is, within 4 min, resulting in a drift rate of –12.5 kHz/s. The structure occurring on 11:46:40 UT in the dynamic radio spectrum (see Fig. 1) is selected to determine their typical spectral features. The temporal leading edge of its “tail” drifts from 25 MHz towards 15 MHz within 8 s leading to a drift rate of –1.25 MHz/s. In Fig. 1, tail structures are also seen in the frequency range 28−32 MHz indicating the emission in the 25−15 MHz range to be the fundamental one. This assumption is supported by the argument that they show a pronounced circular polarization (see Fig. 1) as is typical for the fundamental HB emission (Cairns & Robinson 1987).
3. Data analysis
In the solar corona, the nonthermal radio radiation in the MHzrange is regarded to be plasma emission, that is, the emission of the frequency f takes place near the local electron plasma frequency ω_{pe} = (4πe^{2}N_{e}/m_{e})^{1/2} (e, elementary charge; N_{e}, electron number density; m_{e}, electron mass) and/or its harmonics (Ginzburg & Zheleznyakov 1958; Melrose 1985), that is, f ≈ nf_{pe} (f_{pe} = ω_{pe}/ 2π) with n = 1 and n = 2 for fundamental and harmonic emission, respectively. Since the electron plasma frequency depends solely on the electron number density N_{e} and because of the gravitational density stratification of the corona, a radio source causes a negative or positive frequency drift in the dynamic radio spectrum, if it is moving up or downwards in the corona, respectively. In order to derive a radial source velocity from the frequency drift, a density model of the corona has to be employed. By means of such a density model, a radial distance r from the center of the Sun is related to a definite electron number density and, hence, electron plasma frequency.
The 25 MHzlevel corresponds to an electron number density N_{e} = 7.75 × 10^{6} cm^{3} assuming fundamental (n = 1) emission. We note that the assumption of fundamental emission is justified as already discussed in Sect. 2.
Since the 25 MHz level is probably located in the outer corona (Mann et al. 1999; Warmuth & Mann 2005; Zucca et al. 2014) the onefold Newkirk model (1961), (1)with N_{0} = 4.2 × 10^{4} cm^{3} seems to be appropriate in this region. R_{⊙} denotes the radius of the Sun. The Newkirk (1961) model results from white light scattering measurements of the corona and describes the radial behavior of the electron number density very well in the middle of the corona (Koutchmy 1994). It corresponds to a hydrostatic (or barometric) density model with a temperature T = 1.4 × 10^{6} K (Mann et al. 1999). Employing this model, the 26 and 23 MHzlevels are located at radial distances r = 1.878 R_{⊙}, and 1.969 R_{⊙} from the center of the Sun, respectively. Thus, a distance of 0.091R_{⊙} is covered by the slowly drifting radio source within a period of 240 s resulting in a radial velocity of V_{r,sd} = 264 km s^{1}. The tails are drifting from 25 MHz towards 15 MHz within 8 s. Since the 25 MHz and 15 MHzlevels are located at radial distances r = 1.906 R_{⊙} and 2.370 R_{⊙} from the center of the Sun according to Eq. (1), respectively, the radial source velocity of the tail is found to be V_{r,tail} = 40 400 km s^{1}.
At the 25 MHzlevel, that is, at 1.906 R_{⊙}, a magnetic field B of 0.58 G is derived by means of (2)with B_{0} = 0.5 G (Dulk & McLean 1978). This formula approximately describes the radial behavior of the magnetic field strength above active regions. Hence, one gets an Alfvén speed (for definition see Appendix A) of 423 km s^{1} and a ratio between the electron plasma frequency ω_{pe} and the electron cyclotron frequency ω_{ce} of ω_{pe}/ω_{ce} = 15 at the 25 MHz level (ω_{ce} = eB/m_{e}c; c, velocity of light). For comparison, a sound speed (see Appendix A) and a thermal electron velocity v_{th,e} = (k_{B}T/m_{e})^{1/2} (k_{B}, Boltzmann’s constant; m_{p}, proton mass) of 179 km s^{1} and 4600 km s^{1} are obtained for a coronal temperature of T = 1.4 × 10^{6} K, γ = 5/3, and (Priest 1982), respectively. With such values for the sound and Alfvén speed, a plasma beta of 0.2 is obtained according to Eq. (A11).
4. Plasma emission
As already mentioned in the previous section, the radio emission of the fine structures is assumed to be plasma emission. In this case, energetic electrons excite electrostatic highfrequency plasma waves such as, for example, Langmuir waves which are converted into escaping radio waves by nonlinear plasma processes (Melrose 1985). Hence, the appearance of energetic electrons is a necessary condition for nonthermal radio emission.
Langmuir waves with the frequency ω_{L} and wave number k are governed by the dispersion relation (3)(Baumjohann & Treumann 1997). These waves can be excited by electronbeam plasma instability. The resonance condition, (4)(Treumann & Baumjohann 1997) leads to (5)Here, ω_{pe}/ω_{ce} ≫ 1 (see Sect. 3) has been taken into account. V_{b} denotes the velocity of the electron beam along the ambient magnetic field. The inspection of Eq. (5) shows, that V_{b}> 3^{1/2}·v_{th,e} is a necessary condition for the excitation of Langmuir waves by an electron beam.
In the case discussed here, the electron velocity distribution function is composed of background electrons with a particle number density N_{0} and temperature T_{0} and the beam electrons with a particle number density N_{b} and the temperature T_{b}; it can be written for the complete electron velocity distribution function (6)with v_{th,0} = (k_{B}T_{0}/m_{e})^{1/2} and v_{th,b} = (k_{B}T_{b}/m_{e})^{1/2}. We note that the velocity distribution function is normalized to unity here. The requirement of a currentfree plasma leads to V_{0} = −V_{b}(N_{b}/N_{0}). Now, all velocities are normalized to the thermal speed of the background electrons v_{th,0}, that is, U = V/v_{th,0}, U_{0} = V_{0}/v_{th,0}, and U_{b} = V_{b}/v_{th,0}. Taking into account dV = v_{th,0}dU, Eq. (6) can be rewritten as (7)with ν = N_{b}/N_{0} and w = v_{th,b}/v_{th,0}. Thus, w denotes the width of the beam in velocity space. The occurrence of an instability for exciting Langmuir waves requires that there is a region where the velocity distribution function fulfills dF/ dU> 0 (Treumann & Baumjohann 1997) leading to (8)by means of Eq. (7). Equation (8) shows that instability occurs only for U<U_{b}. The lefthand side of Eq. (8) has a maximum at U = U_{b}−w. Inserting U = U_{b}−w into Eq. (8), the inequality, (9)results in a sufficient condition for exciting Langmuir waves.
5. Shock drift acceleration
Here, we assume that the electron beams related to the discussed fine structures (HBs) are produced by SDA. A fast magnetosonic shock wave is accompanied by a magnetic field compression. Thus, it is a moving magnetic mirror, at which charged particles can be reflected and accelerated. The acceleration happens by the electric field, which is induced in the shock transition region.
For a brief introduction to SDA, we refer to the papers by Holman & Pesses (1983), Ball & Melrose (2001), and Mann & Klassen (2005) for a much more detailed description.
Since the tail’s associated electron beams have typical velocities ≈40 000 km s^{1} (or energies ≈5 keV), they can be treated in a nonrelativistic manner.
The SDA is usually described in the de HoffmannTeller frame, in which the shock wave is at rest and the motional electric field is removed. Therefore, the reflection process can be treated by conservation of the kinetic energy and the magnetic moment in this special frame. As a result, the velocity gain due to SDA is given by (10)with v_{s} being the shock speed in the laboratory frame (Holman & Pesses 1983). Here, V_{i, ∥} and V_{r, ∥} denote the particle velocity parallel to the ambient magnetic field before and after the reflection process, respectively. θ is the angle between the upstream magnetic field B_{1} and the shock normal n_{s}. Because of the conservation of the magnetic moment, the particle velocity perpendicular to the magnetic field stays unchanged during the reflection process. (11)Assuming a Maxwellian distribution (it is normalized to unity here), (12)which is the initial state in the upstream region, after the reflection process, the velocity distribution function of the accelerated electrons has the form (see Eq. (15) in Mann & Klassen 2005) (13)with (14)where V_{s} = v_{s}secθ and V_{e} = (2eφ_{HT}/m_{e})^{1/2}. Here, Θ is the wellknown step function, that is, Θ(x) = 1 for x ≥ 0 and Θ(x) = 0 for x< 0. We note that V_{s} represents the shock speed in the de HoffmannTeller frame. The losscone angle α_{lc} is given by α_{lc} = arcsin [(B_{1}/B_{2})^{1/2}] with B_{1} and B_{2} as the magnetic field strength in the up and downstream regions, respectively.
Due to the different inertia of electrons and protons, a crossshock potential is established at the shock transition region (Goddrich & Scudder 1984; Kunic et al. 2001). In the de HoffmannTeller frame it is given by (15)where γ is the ratio of the specific heats (see Eq. (10) in Mann & Klassen 2005), and T_{1} and T_{2} denote the temperatures in the up and downstream regions, respectively.
The distribution function of Eq. (13) represents a socalled shifted losscone distribution (Leroy & Mangeney 1984; Wu 1984). It is nonzero in the region defined by V_{∥} ≥ V_{s} and V_{⊥} ≥ V_{⊥ ,lc} in the V_{⊥}−V_{∥} plane as illustrated in Fig. 3.
Fig. 3
The velocity distribution function of the accelerated electrons is nonzero in the hatched region in the V_{⊥}−V_{∥} −plane for V_{S} = 4v_{th,e}, V_{e} = 5v_{th,e} and α_{lc} = 45°. The dashed line represents the line V_{⊥} = (V_{∥}−V_{s})·tanα_{lc}. The full line represents the function V_{⊥} = V_{⊥ .lc}(V_{∥}) for the above given parameters. 
The reduced distribution function is defined by (16)with f as the usual distribution function in the whole velocity space. For the accelerated electrons it is found to be (17)with (18)(see Eq. (16) in Mann & Klassen 2005). The reduced distribution function depends solely on the velocity component V_{∥}. For illustration, the shape of the reduced distribution (see Eq. (17)) is presented in Fig. 4;
Fig. 4
Shape of the reduced distribution function for V_{S} = 4v_{th,e}, V_{e} = 5v_{th,e} and α_{lc} = 45°. Its beamlike shape is evident. 
it represents a distribution function of an electron beam with the velocity V_{b} = V_{s}(1 + cos^{2}α_{lc}) and the width 2^{1/2}v_{th,e}cosα_{lc}. F_{acc} is not vanishing for V_{∥} ≥ V_{s}. At V_{∥} = V_{s}, it takes the value (19)Once again, the distribution function of the accelerated electrons has the shape of a shifted losscone distribution. There are regions with ∂f_{acc}/∂V_{∥}> 0 and ∂f_{acc}/∂V_{⊥}> 0, which gives rise to an instability as a mixture of both a losscone and beam instability. As already discussed in Sect. 3, such instabilities lead to exciting Langmuir waves (Treumann & Baumjohann 1997), which are nonlinearly converted into radio waves. Hence, the electron beams produced by SDA are able to emit radio waves.
6. Discussion
As already mentioned in Sect. 1, the chain of fine structures in Fig. 1 is interpreted as the radio signature of a chain of electron beams generated by a shock wave traveling outwards in the corona. Hence, V_{b} = V_{s}(1 + cos^{2}α_{lc}) and v_{s} can be interpreted as the velocity of the electron beam and the associated shock, respectively.
Because the shockaccelerated electron beams must be able to excite Langmuir waves in order to emit radio waves, their velocity V_{b} must necessarily fulfill V_{b}> 3^{1/2}v_{th,e} (see Eq. (5)) leading to the requirement (20)Here, v_{th,e} = 4600 km s^{1} and v_{A} = 423 km s^{1} (see Sect. 3) have been used for the thermal electron speed v_{th,e} and Alfvén speed v_{A}, respectively. As a result, the appearance of radio emission from shockaccelerated electron beams requires that the associated shock must be a high AlfvénMach number (i.e., M_{A} ≫ 1) in the de HoffmannTeller frame and, hence, a nearly perpendicular one (i.e., θ → 90°). That is the reason why the RankineHugoniot relationships are especially evaluated for the case of an almost perpendicular shock in Appendix A, that is, in the regime M_{A} ≫ 1 and θ → 90°. We note that the relationship between the quantities of the up and downstream regions of the shock can be studied by the magnetohydrodynamic RankineHugoniot relationships in a magnetized plasma such as the corona (see Appendix A).
Figure 5 presents the shock geometry in the corona.
Fig. 5
Shock geometry in the corona. 
The shock is oriented perpendicular to the image plane. The shock normal n_{s} takes an angle α to the radial direction n_{r} and an angle θ to the upstream magnetic field B_{1}.
According to Fig. 5, the radial velocities of the shock v_{r,s} and of the electron beam V_{r,b} are given by v_{r,s} = v_{s}·cosα and V_{r,b} = V_{b}·cos(θ−α), respectively. They are identified with the velocities of the slowly drifting chain, that is, V_{r,sd} = v_{r,s}, and of the tails, that is, V_{r,tail} = V_{r,b}, as derived from the dynamic radio spectrum (see Fig. 1).
Then, the real shock speed is given by v_{s} = v_{r,s}/ cosα resulting in (21)with as the AlfvénMach number in the laboratory frame. Here, v_{r,s} = V_{r,sd} = 264 km s^{1} and v_{A1} = 463 km s^{1} have been used. Thus, the angle α and the shock speed can be determined for a given value of . Furthermore, the relationship (22)can be derived for calculating the angle θ for a given value of . The angle α_{lc} is given by α_{lc} = arcsin [(B_{1}/B_{2})^{1/2}] = arcsin [X^{− 1/2}] for a nearly perpendicular shock (see Eq. (A.9)). We note that the density jump X across the shock is a function of according to Eq. (A.10). With the sofound angles α, θ, and α_{lc}, the velocity of the electron beam V_{b} = v_{s}·secθ·(1 + cos^{2}α_{lc}) can be determined. As a result, all quantities can be fixed for a given value of or X taking into account the results of the RankineHugoniot relationships (see Appendix). Varying X in the range 1.4−2.6, the best fit with respect to the observations, that is, V_{r,sd} = v_{r,s} = 264 km s^{1} and V_{r,tail} = V_{r,b} = 40 400 km s^{1}, is found for X = 2.3 resulting in: , v_{s} = 1010 km s^{1}, α = 74.8°, α_{lc} = 41.3°, θ = 87.8°, V_{s} = 26.500 km s^{1}, and V_{b} = 41 400 km s^{1}. This leads to an AlfvénMach number of the associated shock M_{A} = V_{s}/v_{A} = 63 in the de HoffmannTeller frame as should be expected (i.e., M_{A} ≫ 1; see the discussion with respect to Eq. (20)). Employing all these parameters, the radial velocity of the electron beam V_{r,b} = V_{b}·cos(θ−α) is 40 400 km s^{1} and agrees well with the radial speed of the tail’s radio source as derived from the dynamic radio spectra in Fig. 1. Thus, the cycle of argumentation is closed in the best way for this fit. Consequently, it is justified to consider the derived shock parameters as the real ones.
By means of SDA the electrons are accelerated up to velocities V_{∥} ≥ V_{s}. We note that both V_{s} and V_{b} are greater than 3^{1/2}v_{th,e} (i.e., V_{s} = 26 500 km s^{1}> 5.8v_{th,e} and V_{b} = 41 400 km s^{1}> 9.0v_{th,e} with v_{th,e} = 4600 km s^{1}) as required for exciting Langmuir waves (see the discussion in Sect. 4). The ratio between the number density of the beam electrons and background ones N_{b}/N_{0} can be determined to be by comparing the velocity distribution function of shockaccelerated electrons (see Eq. (17)) with that of an electron beam (see Eq. (6)). F_{b} is given by Eq. (18). With the above given parameters, one gets ν = N_{b}/N_{0} = 4.9 × 10^{11} with V_{e} = 6.54v_{th,e} according to Eq. (A.13). Then, the sufficient condition for exciting Langmuir waves (see Eq. (9)) is well fulfilled by the shockaccelerated electrons, that is, ν = 4.9 × 10^{11} > 1.3 × 10^{14}.
In the special case of V_{∥} ≈ V_{s}, the accelerated electrons are at rest in the de HoffmannTeller frame, or, in other words, they are attached at the shock, that is, they cannot leave the shock front. The inspection of Eqs. (13) and (17) as well as of Fig. 4 reveals that there is a substantial number of accelerated electrons at velocities of V_{∥} ≈ V_{s}. They show a positive slope in the velocity distribution (see Eq. (17) and Fig. 4), that is, dF_{acc}/ dV_{∥}> 0 , giving rise to Langmuir wave excitation (Treumann & Baumjohann 1997). Since they cannot leave the shock front, they move with the shock through the corona. These electrons could cause the radio radiation of the heads, whereas the tails are generated by the electrons forming the beam at velocities V_{∥} ≈ V_{b}. Thus, this approach is also able to explain the head and tail radiation as different components of the fine structures.
To emit radio waves via plasma emission, the excitation of Langmuir waves is necessary. As discussed in Sect. 4, the velocity of the electron beam generating Langmuir waves must exceed at least 3^{1/2}v_{th,e}, that is, (23)This condition is dependent on the shock speed v_{s} and the thermal electron velocity v_{th,e} as well as the angles θ and α_{lc}. The angle θ is the most crucial parameter. The condition (23) is well fulfilled for θ> 87° in the case discussed here. Once again, θ is the angle between the shock normal and the upstream magnetic field. This angle can vary during the motion of the shock through the corona. If the value of θ drops below this threshold, the excitation of Langmuir waves and, consequently, the associated radio emission are switched off. Thus, the appearance of the discussed fine structures is crucially dependent on the value of θ. This could explain the bursty occurrence of them leading to a slowly drifting chain of such structures in the dynamic radio spectrum.
Shock drift acceleration only produces electron beams in the upstream region. Such beams are necessary to excite Langmuir waves and, consequently, for the emission of radio waves, as already mentioned in Sect. 1. Herring bones drift away from the backbone towards low and high frequencies in dynamic radio spectra. This phenomenon was interpreted by Holman & Pesses (1983) as a curved shock wave traveling through the corona. It can generate electron beams at different parts of the shock as illustrated in Fig. 6.
Fig. 6
Sketch of a shock propagating laterally with respect to the solar surface in the corona. The shock is traveling into the direction n. At the flanks, shocks are formed with different angles θ between the shock normal n_{s} and the upstream magnetic field B. Electrons accelerated due to SDA can propagate along the upstream magnetic field towards both the higher and/or lower corona. 
That leads to a drift towards lower and/or higher frequencies of the associated HBs in dynamic radio spectra. But the emission of radio waves can only happen if Langmuir waves are excited, which requires the fulfilment of condition (23). If this condition is not fulfilled, for example in the case discussed where θ< 87°, then there is no radio emission. That may explain why only negatively drifting tails are seen in the observations discussed here.
Energetic electrons in the solar corona are also observed in terms of typeIII bursts in the solar radio radiation (Wild & McCready 1950; Reid & Ratcliffe 2014, for review). They appear as emission stripes rapidly drifting from high to low frequencies (sometimes down to a few tens of kHz) in dynamic radio spectra (Suzuki & Dulk 1985). These bursts are regarded as the radio signature of beams of energetic electrons traveling along magnetic field lines through the corona and sometimes up to the interplanetary space. These energetic electrons are produced as a result of flares (or magnetic reconnection). Initially, they have a power law distribution in velocity space. Due to their spatiotemporal evolution, they form a beamlike distribution, that is, with positive slope in the velocity distribution, after a distance from the initial acceleration region as discussed by Reid et al. (2011, 2014). Then, they are able to excite Langmuir waves which are converted into radio waves by nonlinear waveparticle interactions (Melrose 1985). In contradiction to the energetic electrons associated with typeIII bursts, SDA initially generates a beamlike electron distribution (see Eq. (17)). Thus, the excitation of Langmuir waves and the emission of radio waves in terms of HBs happens immediately in the vicinity of the shock. In the case discussed here (see Fig. 1), the tails appear in the frequency range 25−15 MHz, that is, the associated electron beam travels over a distance of 323 Mm before the radio emission is switched off. According to the discussion by Reid & Kontar (2013, 2015), such electron beams are influenced by Coulomb collisions, scattering at density fluctuation in the background plasma, and inhomogeneities of the magnetic field as well as emissions and absorptions of Langmuir waves. Because they have energies of about 5 keV, as in the present case, Coulomb collisions and scattering at density fluctuations could be the most important influences on the electron beam (see the discussion in Mann & Klassen 2005; and the results of Reid & Kontar 2013, 2015) resulting in a destruction of the beam followed by the interruption of Langmuir wave excitation, and, finally, in the switching off of the radio emission.
As a result, the HB structures observed in the decameter range are generated by electron beams produced by a nearly perpendicular shock, which laterally propagates in the corona. Such shocks can be driven by the expanding flanks of CMEs moving upwards through the corona. Such a scenario is confirmed by observations (SalaMatamoros et al. 2016), theoretical studies (Schmidt & Cairns 2012a, b) and numerical simulations (Schmidt et al. 2014).
7. Summary
The Ukrainian radio telescope URAN2 recorded special fine structures in the solar radio radiation on April 7, 2011. They appear as slowly drifting chains of stripes of enhanced radio emission in the dynamic radio spectrum in the frequency range 8−33 MHz (see Fig. 1). The chain drifts from 26 to 23 MHz within 240 s. The discussed fine structures consist of two components, namely, the head at the high frequency edge and the tail, which rapidly drifts from high to low frequencies, that is, over ≈10 MHz within ≈8 s. The slow drift rate of the chains corresponds to those of solar typeII radio bursts in this frequency range. Solar typeII radio bursts are regarded as the radio signature of shock waves traveling outwards through the corona. The tails resembling HBs appear as special fine structures in solar typeII radio bursts. Herringbones are considered to be the radio signatures of electron beams accelerated at the typeIIburstassociated shock wave. For these reasons, the tails are interpreted as shockaccelerated electron beams. According to this view, the slow drift is identified as a shock wave moving outwards in the corona. As mentioned in Sect. 1, there are also cases in which the tails drift towards higher frequencies. These can be explained in the following way: generally, in the case of SDA, the electrons are scattered back into the upstream region, where they appear as beams of energetic electrons. If the upstream magnetic field is directed outwards, that is, into regions with lower electron number densities, the tails show a negative drift in the dynamic radio spectrum. In the opposite case, in which the upstream magnetic field is directed towards the lower corona, that is, into regions with higher electrons number densities, the tails drift towards higher frequencies.
Since the radio emission in this frequency range, that is, 8−33 MHz, is thought to be plasma emission (as discussed in Sect. 4), these radio waves are emanated from the outer corona. The radial velocities of the shockaccelerated electron beams and the associated shock waves are derived to 40 400 km s^{1} and 264 km s^{1}, respectively.
The electrons are thought to be accelerated by SDA. In the framework of SDA, a shifted losscone distribution is established for the accelerated electrons. This results in a beamlike distribution with respect to the electron velocities parallel to the ambient magnetic field. Such a distribution is unstable for exciting Langmuir waves. This requires the velocity of the electron beam to exceed 3^{1/2} times the thermal electron velocity. If the electron beam is generated by SDA, then only a nearlyperpendicular fast mode shock fulfills this condition. This is why the RankineHugoniot relationships are evaluated under the assumption of an almost perpendicular shock (see Appendix A).
Taking all of the above into account, the most robust scenario linking the observational data, that is, the radial velocities of the shock and the electron beam, to the theoretical results, involves an almost perpendicular shock, traveling with a speed of 1010 km s^{1} nearly horizontal with respect to the surface of the Sun, generating electron beams with a velocity of 41 400 km s^{1} by means of SDA.
Acknowledgments
A. I. Brazhenko, A. A. Konovalenko, G. Mann, and V. N. Melnik express their thanks to H. O. Rucker and the staff of the Space Research Institute of the Austrian Academy of Sciences (Graz, Austria) for their kind hospitality; work was carried out at this institute. G.M. performed this work under the framework of grant No. 50QL0001 of the Deutsches Zentrum für Luft und Raumfahrt (DLR). V.N.M. and A.A.K. made this work in the frame of project “Spectr3” by National Academy of Sciences of Ukraine.
References
 Alvarez, H., & Haddock, F. T. 1973, Sol. Phys., 29, 197 [NASA ADS] [CrossRef] [Google Scholar]
 Aurass, H. 1996, in Coronal Physics from Radio and Space Observations, ed. G. Trottet (Heidelberg: Springer Verlag), Lect. Not. Phys., 135 [Google Scholar]
 Ball, L., & Melrose, D. B. 2001, PASA, 18, 361 [NASA ADS] [CrossRef] [Google Scholar]
 Baumjohann, W., & Treumann, R. A. 1997, Basic Space Plasma Physics, (London: Imperial College Press) 202 [Google Scholar]
 Brazhenko, A. I., Bulatsen, V. G., Vashchishin, R. V., et al. 2005, Kinematika i Fizika Nebesnykh Tel, Suppl., 5, 43 [Google Scholar]
 Brazhenko, A., Melnik, V., Konovalenko, A., et al. 2012, Proc. European Planetary Science Congress (September 23–28, 2012 in Madrid, Spain), 544 [Google Scholar]
 Cairns, I. H., & Robinson, R. D. 1987, Sol. Phys., 111, 365 [NASA ADS] [CrossRef] [Google Scholar]
 Cairns, I. H., Knock, S. A., Robinson, P. A., & Kunic, Z. 2003, Space Sci. Rev., 107, 27 [NASA ADS] [CrossRef] [Google Scholar]
 Cane, H. V., & White, S. M. 1989, Sol. Phys., 120, 137 [NASA ADS] [CrossRef] [Google Scholar]
 Dorovskyy, V. V., Melnik, V. N., Konovalenko, A. A., et al. 2015, Sol. Phys., 290, 2031 [NASA ADS] [CrossRef] [Google Scholar]
 Dulk, G. A., & McLean, D. J. 1978, A&A, 66, 315 [NASA ADS] [Google Scholar]
 Goodrich, C. C., & Scudder 1984, J. Geophys. Res., 89, 6654 [NASA ADS] [CrossRef] [Google Scholar]
 Holman, G. D., & Pesses, M. E. 1983, ApJ, 267, 837 [NASA ADS] [CrossRef] [Google Scholar]
 Koutchmy, S. 1994, Adv. Space Res., 14, 29 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Kunic, Z., Cairns, I. H., & Knock, S. A. 2001, J. Geophys. Res., 107, 101029 [Google Scholar]
 Leroy, M. M., & Mangeney, A. 1984, Ann. Geophys., 2, 449 [NASA ADS] [Google Scholar]
 Magdalenic, J., Marque, C., Zhukov, A. N., Vršnak, B., & Veronig, A. 2012, ApJ, 746, 152 [NASA ADS] [CrossRef] [Google Scholar]
 Mann, G. 1995, in Coronal Magnetic Energy Release, eds. A. O. Benz, & A. Krüger (Heidelberg: Springer Verlag), Lect. Not. Phys. 183 [Google Scholar]
 Mann, G. 2006, in Solar Eruptions and Energetic Particles (AGU, Washington (DC)), Geophys. Monograph Series, 165, 221 [Google Scholar]
 Mann, G., & Klassen, A. 2002, Proc. 10th European Solar Physics Meeting, Prague, Czech Rep., ESA SP, 506, 245 [Google Scholar]
 Mann, G., & Klassen, A. 2005, A&A, 441, 319 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mann, G., Jansen, F., MacDowall, R. J., Kaiser, M. L., & Stone, R. G. 1999, A&A, 348, 614 [NASA ADS] [Google Scholar]
 Melnik, V. N., Konovalenko, A. A., Rucker, H. O., et al. 2004, Sol. Phys., 222, 151 [NASA ADS] [CrossRef] [Google Scholar]
 Melnik, V. N., Rucker, H. O., Konovalenko, A. A., et al. 2008, in Solar Physics Research Trends, ed. P. Wang (New York: Nova Science Publishers), 287 [Google Scholar]
 Melrose, D. 1985, in Solar Radio Physics, eds. D. J. McLean, & N. R. Labrum (Cambridge: Cambridge Univ. Press), 177 [Google Scholar]
 Nelson, G. S., & Melrose, D. 1985, in Solar Radio Physics, eds. D. J. McLean, & N. R. Labrum (Cambridge: Cambridge Univ. Press), 333 [Google Scholar]
 Newkirk, G. A. 1961, ApJ, 133, 983 [NASA ADS] [CrossRef] [Google Scholar]
 Nindos, A., Aurass, H., Klein, K.L., & Trottet, G. 2008, Sol. Phys., 253, 3 [NASA ADS] [CrossRef] [Google Scholar]
 Priest, E. R. 1982, Solar Magnetohydrodynamics (Dordrecht: Kluwer Academic Publishers), 1999 [Google Scholar]
 Reid, H. A. S., & Kontar, E. P. 2013, Sol. Phys., 285, 217 [NASA ADS] [CrossRef] [Google Scholar]
 Reid, H. A. S., & Kontar, E. P. 2015, A&A, 577, A124 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Reid, H. A. S., & Ratcliffe, H. 2014, RA&A, 14, 773 [Google Scholar]
 Reid, H. A. S., Vilmer, N., & Kontar, E. P. 2011, A&A, 529, A66 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Reid, H. A. S., Vilmer, N., & Kontar, E. P. 2014, A&A, 567, A85 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Ryabov, V. B., Vavriv, D. M., Zarka, P., et al. 2010, A&A, 510, A16 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 SalasMatamoros, C., Klein, K.L., & Rouillard, A. P. 2016, A&A, 590, A135 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Schmidt, J. M., & Cairns, I. H. 2012a, J. Geophys. Res., 117, A04106 [NASA ADS] [Google Scholar]
 Schmidt, J. M., & Cairns, I. H. 2012b, J. Geophys. Res., 117, A11104 [NASA ADS] [Google Scholar]
 Schmidt, J. M., Cairns, I. H., & Lobzin, V. V. 2014, J. Geophys. Res., 119, 6042 [CrossRef] [Google Scholar]
 Smerd, S. F., Wild, J. P., & Sheridan, K. V. 1962, Austr. J. Phys., 15, 180 [NASA ADS] [CrossRef] [Google Scholar]
 Stewart, R. T., Howard, R. A., Hansen, F., Gergely, T., & Kundu, M. R. 1974a, Sol. Phys., 36, 219 [NASA ADS] [CrossRef] [Google Scholar]
 Stewart, R. T., McCabe, M., Koomen, M. J., Hansen, F., & Dulk, G. 1974b, Sol. Phys., 26, 203 [NASA ADS] [CrossRef] [Google Scholar]
 Stewart, R. T., & Magun, A. 1980, PASA, 3, 375 [Google Scholar]
 Suzuki, S., & Dulk, G. A. 1985, in Solar Radio Physics, eds. D. J. McLean, & N. R. Labrum (Cambridge: Cambridge Univ. Press), 289 [Google Scholar]
 Suzuki, S., Stewart, R. T., & Magun, A. 1982, in Radikophysics of the Sun, eds. M. R. Kundu, & T. R. Gergely, IAU Symp., 86 [Google Scholar]
 Treumann, R. A., & Baumjohann, W. 1997, Advance Space Plasma Physics (London: Imperial College Press), 69 [Google Scholar]
 Uchida, Y. 1960, PASJ, 12, 376 [NASA ADS] [Google Scholar]
 Uchida, Y., Altschuler, D., & Newkirk, G. 1973, Sol. Phys., 28, 495 [Google Scholar]
 Vršnak, B., Ruždjak, V., & Aurass, H. 1995, Sol. Phys., 158, 331 [NASA ADS] [Google Scholar]
 Warmuth, A., & Mann, G. 2005, A&A, 435, 1123 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Wild, J. P., & McCready, L. L. 1950, Austr. J. Sci. Res., Ser. A, 3, 387 [Google Scholar]
 Wild, J. P., & Smerd, S. F. 1972, ARA&A, 10, 159 [NASA ADS] [CrossRef] [Google Scholar]
 Wu, C. S. 1984, J. Geophys. Res., 89, 8857 [NASA ADS] [CrossRef] [Google Scholar]
 Zlobec, P., Messerotti, M., Karlicky, M., & Urbarz, H. 1993, Sol. Phys., 144, 373 [NASA ADS] [CrossRef] [Google Scholar]
 Zucca, P., Carley, E. P., Bloomfield, D. Sh., & Gallagher, P. T. 2014, A&A, 564, A47 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
Appendix A: RankineHugoniot relationships in magnetohydrodynamics
A shock wave is basically a discontinuity with a mass flux across it. It is additionally a dissipative structure. In a magnetized plasma such as the corona, shock waves can be described by conservation of mass, momentum, energy, and magnetic flux in terms of the magnetohydrodynamic (MHD) equations. This is usually done in the de HoffmannTeller frame, in which the shock is at rest and the motional electric field is removed (see e.g., Priest 1982).
The framework of the shock is given in the following way: The shock normal n_{s} is directed along the xaxis, that is, the shock itself is located in the y−z plane. The magnetic field B is put in the x−z plane and takes an angle θ to the xaxis in the upstream region. All quantities with the subscript 1 and 2 denote those in the up and downstream regions, respectively. Then, the conservation of mass, x and zcomponent of the momentum, energy and magnetic flux leads to the jump relationships in terms of X = ρ_{2}/ρ_{1}: (Priest 1982). Here, the Alfvén and sound speed are given by v_{A1} = B_{1}/ (4πρ_{1})^{1/2} and c_{s1} = (γp_{1}/ρ_{1})^{1/2} (γ, ratio of the specific heats) in the upstream region, respectively. ρ denotes the mass density, which is related to the full particle number density N by (m_{p}, proton mass; , mean molecular weight). Additionally, the equation of state p = Nk_{B}T has been used. It relates the pressure p to the full particle number N and the temperature T. Then, the sound speed is also given by . The AlfvénMach number M_{A} is defined by M_{A} = (v_{s}/v_{A1})·secθ; it is actually the AlfvénMach number in the de HoffmannTeller frame. Here, v_{s} denotes the shock speed in the laboratory frame. These jump relationships are supplemented by a polynominal of the third order in X: (A.6)(Priest 1982). In the limit X → 1, this provides the wellknown three MHD modes, namely the Alfvén wave (M_{A} = 1) as well as the fast (M_{A}> 1) and slow (M_{A}< 1) magnetoacoustic waves.
By means of Eqs. (A.1), (A.3), and (A.4), one gets (A.7)At a shock, there must be a positive jump of the density (i.e., X> 1) in order to guarantee an increase of the entropy (see Priest 1982), which is necessary for a dissipative structure. Then, according to Eq. (A.7), B_{1}<B_{2}, that is, magnetic field energy increases across the fastmode shocks (i.e., M_{A}>X).
In the special case of a nearly perpendicular shock, that is, M_{A} ≫ 1 and θ → 90°, Eqs. (A.5) and (A.7) reduce to (A.8)and (A.9)with . In this limit, the polynominal (A.6) provides (A.10)Here, β denotes the ratio between the thermal and the magnetic pressure, that is, (A.11)The inspection of Eq. (A.10) shows if X → (γ + 1)/(γ−1) (=4 for γ = 5/3). The dependence of on X according to Eq. (A.10) is illustrated for β = 0.2 in Fig. A.1, that is, for 423 km s^{1} and 179 km s^{1} as typical values of the Alfvén and sound speed at the 25 MHz level in the corona, respectively. For this case, .
Fig. A.1
Dependence of on X according to Eq. (A.10) for β = 0.2. 
With the expression of the jump of the temperature across the shock (see Eq. (A.8)), the cross shock potential (see Eq. (15)) is found to be (A.12)Then, one gets for the velocity V_{e} = (2eφ_{HT}/m_{e})^{1/2}(A.13)
All Figures
Fig. 1
Dynamic radio spectrum in the frequency range 8–33 MHz as recorded during the period 11:43–11:50 UT on 7 April 2011 with the grayscaled intensity (top), which is given in solar flux units (1 s.f.u. = 10^{22} W m^{2} Hz^{1}), and degree of circular polarization (bottom). 

In the text 
Fig. 2
Sketch of the shape of the fine structures discussed here. 

In the text 
Fig. 3
The velocity distribution function of the accelerated electrons is nonzero in the hatched region in the V_{⊥}−V_{∥} −plane for V_{S} = 4v_{th,e}, V_{e} = 5v_{th,e} and α_{lc} = 45°. The dashed line represents the line V_{⊥} = (V_{∥}−V_{s})·tanα_{lc}. The full line represents the function V_{⊥} = V_{⊥ .lc}(V_{∥}) for the above given parameters. 

In the text 
Fig. 4
Shape of the reduced distribution function for V_{S} = 4v_{th,e}, V_{e} = 5v_{th,e} and α_{lc} = 45°. Its beamlike shape is evident. 

In the text 
Fig. 5
Shock geometry in the corona. 

In the text 
Fig. 6
Sketch of a shock propagating laterally with respect to the solar surface in the corona. The shock is traveling into the direction n. At the flanks, shocks are formed with different angles θ between the shock normal n_{s} and the upstream magnetic field B. Electrons accelerated due to SDA can propagate along the upstream magnetic field towards both the higher and/or lower corona. 

In the text 
Fig. A.1
Dependence of on X according to Eq. (A.10) for β = 0.2. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.