A&A 425, 281-285 (2004)
DOI: 10.1051/0004-6361:20040174

The jump of reverse drift rate of V-shaped structure and the large-scale inhomogeneity of plasma density in the solar atmosphere

F. Y. Xu1,2 - H. A. Wu1,2


1 - Purple Mountain Observatory, Chinese Academy of Sciences, Nanjing 210008, PR China
2 - National Astronomical Observatories, Chinese Academy of Sciences, Beijing 100012, PR China

Received 16 October 2002 / Accepted 31 March 2004

Abstract
A rare solar radio V-shaped structure was observed on August 25, 1999 by the 4.5-7.5 GHz spectrograph of the Purple Mountain Observatory, China. It consists of a reverse drift component followed by a normal drift one, which seems to be an unusual frequency drift event because of the jump of reverse drift rate from 3.22 to 1.53 GHz/s and the fastest drift rate of -3.53 GHz/s. According to their different drift rates, the V-shaped structure can be divided into three parts. By the consideration of plasma emission of an electron beam, a three-component solar atmospheric model with large and small scale lengths $\lambda $, $\lambda _{1}$ and $\lambda _{2}$ to describe the equilibrium solar atmosphere and local inhomogeneity as well as the energy loss of the electron beam due to electron-electron collision is proposed. The local scale lengths $\lambda _{1}=5.5\times 10^{9}$, $\lambda _{2}=7.7\times 10^{9}$ cm and $\overline{v}_{1}=0.41$ c, $\overline{v}_{2}=0.25$ c at Parts 1 and 2 of the reverse drift component are obtained. The variations of scale length of $\lambda
_{1}<\lambda _{2}$ and $\overline{v}_{1}>\overline{v}_{2}$, especially $\overline{v}_{1}>\overline{v}_{2}$, could be responsible for the jump of drift rate. The normal drift rate would be interpreted as the increase of average energy $E^{\prime
}$ of the electron beam at the mirror due to the decrease in electrons of lower energy via energy loss of collision. The three-component solar atmospheric model and the energy loss of the electron beam are discussed.

Key words: sun: flares - Sun: radio radiation

1 Introduction

Earlier observations of solar radio type III bursts were made mainly in meter/decimeter wavebands, and their properties were characterized by narrow bandwidth, frequency drift and harmonic structure etc. (Suzuki & Dulk 1985). The characters of narrow bandwidth and frequency drift suggest that they may originate from a coherent emission mechanism such as from local plasma emission at its characteristic frequency, $\omega
_{\rm p}$, excited by a source moving upward/downward along the field line. The exciting Langmuir wave would be converted into fundamental or harmonic escaping radiation at characteristic frequencies $\omega =\omega _{\rm p}$, $\omega =2\omega _{\rm p}$ by scattering or mode-mode coupling processes (Zheleznyakov 1970; Melorose 1980). The frequency of escaping radiation follows the variation of exciting local plasma frequency which decreases/increases with time as the electron beam moves upward/downward, and leads to a detectable normal/reverse drift type III burst. Such a picture is still accepted today for type III or type III-like bursts.

Since the 1980s, the surveys of solar radio type III bursts have been conducted in microwave bands, such as 0.1-3 GHz, 4-8.5 GHz and 6-8.5 GHz (Stähli & Benz 1987; Alloart et al. 1999; Bruggmann et al. 1990). The observed characters of microwave type III (III$_{\mu }$) bursts are similar to those of meter/decimeter type III bursts except that the drift rate of type III$_{\mu }$ bursts is almost two orders of magnitude faster than that of meter/decimeter type III bursts, and also the dominant drift direction is a reverse drift in the microwave band but normal one in the meter and decimeter waveband. What electron density is required for a type III$_{\mu }$ burst? For example, in the range of 4.5-7.5 GHz, the required electron density is about $6\times 10^{10}$ to $2\times 10^{11}/\rm cm^{3}$ if harmonic emission is assumed. This probably corresponds to the region between the lower corona and higher chromosphere. Because type III and type III-like structures may be the reflection of the electron beam trace or magnetic field configuration, the type III$_{\mu }$ and type III$_{\mu }$-like structures in the range of 4.5-7.5 GHz may play an important role in plasma diagnosis in the transition region.

In this paper a rare type III$_{\mu }$-like structure, a V-shaped structure, observed on August 25, 1999 will be presented. The variation of reverse drift rate will be analyzed by means of the larger-scale inhomogeneity of the plasma density distribution in the solar atmosphere. The different effects of drift rate on both reverse and normal drift components as well as the parameters associated with the source region and environment are discussed.

2 Observations and data analysis

2.1 Observations

The solar flare (1N/M 3.6) on August 25, 1999 occurred near the region of NOAA 8673 and 8674. The accompanying radio burst consists of three Parts I, II and III (Fig. 1), and contains rich fine structures. To search for fine structures, especially for faint fine structures, the dynamic spectrogram corresponding to Part I was enlarged interval by interval, and its chromatic scale was also adjusted until the contrast between emissions of faint fine structure and background was large enough. Two rare V-shaped structures located in the interval of 01:34:44.396-01:34:44.960 UT and frequency range of 5.57-6.54 GHz were detected, and their spectrogram together with the partial enlargement of the 2nd V-shaped spectrum as well as the time profile are given in Figs. 2a, b and c. As seen from Fig. 2a, the 1st V-shaped structure is a faint featureless structure and is due to an electron beam moving downward along a magnetic loop then mirroring upward. The 2nd V-shaped structure, however, is an unusual one. The 2nd one began at 01:34:44.480 UT with a starting frequency of 5.57 GHz, then drifted downward to the frequency $\sim $6.01 GHz with a reverse drift rate 3.22 Ghz/s. It then drifted downward to the turnover frequency 6.26 GHz. In this latter part, however, the reverse drift rate was about 1.53 GHz/s and was slower than that of the initial part. From the turn-over frequency, the normal drift component began with the fastest drift rate of -3.53 GHz/s which led to an increase of the instantaneous bandwidth after the change of frequency drift. Such a feature will be discussed in Sect. 3. The variation of drift rate is a rare phenomenon.

  \begin{figure}
\par\includegraphics[height=4.9cm, width=7.8cm,clip]{3194.f1.eps}\end{figure} Figure 1: Time profile of the solar radio event on August 25, 1999. It consists of three bursts, complex, strong impulse and weaker complex burst, denoted by I, II and III, respectively.
Open with DEXTER


  \begin{figure}
\par\includegraphics[height=14.61cm, width=8.5cm,clip]{3194.f2.eps}\end{figure} Figure 2: Dynamic spectrogram and corresponding time profile of V-shaped structures with time resolution of 5 ms and frequency resolution of 10 MHz. a) Spectrogram of two V-shaped structures. The 1st V-shaped structure located at frequency range of $\sim $6.01-6.6 GHz is a weak and featureless structure but the 2nd is an unusual one with a jump of drift rate on the reverse drift component and a faster normal drift rate (compared with the reverse drift rate). b) The partial enlargement of the 2nd V-shaped structure, which could show the jump of drift rate clearly, Note that there is a strong diffusion in the frequency range of $\sim $5.90 to 6.26 GHz (only the higher frequency of 6.18 GHz is shown here). c) Corresponding time profile of V-shaped structures at selected frequency. The 2nd V-shaped structure could be seen still. Also note that in the frequency range of $\sim $5.9 to 6.26 GHz the peak of drift structure emission is superposed to form a continual enhancement emission, which may be caused by diffusion due to deflected electrons of lower energy.
Open with DEXTER

Benz et al. (2001) stated that some artificial effects may produce a spurious symmetric pattern in the solar radio spectrogram. Therefore, a series of methods has been adopted to avoid the origin of spurious spectra during the design of the solar radio spectrometer at 4.5-7.5 GHz. More details about this technology have been reported by Yao et al. (2002). The V-shaped structure does not seem to be produced by an artificial effect due to simultaneous observations, because the simultaneous observations of the V-shaped structures were independently obtained in the same interval and frequency range by two spectrometers located at two observatories (PMO and NAO) separated by a distance of over one thousand km. Therefore, the V-shaped structure presented here has a solar origin and hence is reliable.

2.2 Analysis

Under the assumption of an isothermal atmospheric model, the electron density exponentially decreases with increasing height above photosphere i.e. $n(h)=n_{0}{\rm e}^{-h/\lambda }$, the frequency drift rate of type III or type III-like burst may be written as

\begin{displaymath}\frac{{\rm d}f}{{\rm d}t}=-\frac{fv_{0}\cos \theta }{2\lambda },
\end{displaymath} (1)

where ${\rm d}f/{\rm d}t$, f, the frequency drift rate and frequency, are the measurable quantities, v0, the initial velocity of the electron beam in the 0.1-0.3 c range, sometimes up to 0.6 c (c-light speed), $\lambda $, the scale length to describe the density gradient, usually assumed to be ${\sim} 10^{10}$ cm and $\theta $, the angle between the zenith at the source location and source motion direction, which is invariable except for the source location at the mirror or near the top of the magnetic loop. As seen from (1), Eq. (1) is valid only if v0 = constant. Therefore, ${\rm d}f/{\rm d}t$ depends on $\lambda $.

The electron density in the microwave source region, however, is two orders of magnitude higher than that in the meter or decimeter source region. Thus, energy loss of the electron beam and hence the corresponding velocity decrease via electron-electron collision must be taken into account. For an average electron of energy E, the energy loss rate then may be given by

\begin{displaymath}\frac{{\rm d}E}{{\rm d}t}=-\sigma _{\rm ee}Env,
\end{displaymath} (2)

where v is the velocity corresponding to E, n, the local electron density and $\sigma _{\rm ee}$, the collision cross-section for electron on electron, and $\sigma _{\rm ee}$ may be expressed as

\begin{displaymath}\sigma _{\rm ee}=\frac{K}{E^{2}},
\end{displaymath} (3)

where $K=55.7\pi e^{4}$ (Brown 1971), and e is the electronic charge. Then energy loss may be obtained by integration of Eq. (2)

\begin{displaymath}E=\left(E_{0}^{3/2}-\frac{3}{2}\sqrt{\frac{2}{m_{\rm e}}}n\Delta t\right)^{2/3},
\end{displaymath} (4)

where $m_{\rm e}$ is the electronic mass, E, the remaining energy after energy loss, which depends on both the initial energy E0 and n if $\Delta t$ is given. As shown by Eq. (4), under the assumption of an isothermal atmospheric model, the energy E(corresponding to v) changing with time is also a slowly exponentially varying function, and hence cannot yield a jump of drift rate.
  \begin{figure}
\par\includegraphics[height=6cm, width=7.8cm,clip]{3194.f3.eps}\end{figure} Figure 3: Model of a three-component atmosphere. The density distribution in the magnetic loop is characterized by the scale length $\lambda $, while the inhomogeneity is described by two small local scale lengths  $\lambda _{1}$ and  $\lambda _{2}$.
Open with DEXTER

The solar atmosphere is not in equilibrium during a flare and thus an isothermal atmospheric model seems to be invalid. Based on the two-component model suggested by Aschwanden & Benz (1986), a similar three-component model described by three exponential parts with differing scale lengths (see Fig. 3) will be assumed. The first component with a large scale length $\lambda $ represents an isothermal atmospheric model in a magnetic loop, the second and the third components with small local scale lengths $\lambda _{1}$and $\lambda _{2}$ describe an inhomogeneous region superimposed on the first component.

Table 1: Parameters and drift rates for Parts 1 and 2 calculated in comparison with the observations.

Such a model will apply to not only the interpretation of the change of frequency drift but also the diagnosis for plasma in the source region e.g. the diagnosis for the deviation of plasma density distribution from the equilibrium solar atmospheric model as well as for the origin of the brightened knot in the soft X-ray image etc. The density in the inhomogeneous region then may be written as

\begin{displaymath}n_{i}=n_{i+1}{\rm e}^{-h_{i}/\lambda _{i}}~~~~~~~(i=1,2\cdots),
\end{displaymath} (5)

where ni corresponds to the character frequency fi(assuming harmonic emission) of the V-shaped structure. For the sake of clearness the 2-nd V-shaped structure will also be divided into three parts according to their different drift rate denoted by numbers of 1, 2 and 3 (as shown schematically in Fig. 4). The parameters such as f, n and v etc. are denoted by their subscripts. For example, f1 denotes the starting frequency, f2, the frequency at which the reverse drift rate becomes slower, f3, the turnover frequency and f4, the lower-frequency of normal drift component. To estimate the energy loss at Parts 1 and 2 of V-shaped structure, Eq. (4) will be rewritten as

\begin{displaymath}E_{ei}=\left(E_{oi}^{3/2}-\frac{3}{2}K\sqrt{\frac{2}{m_{\rm e}}}n_{i}\Delta
t_{i}\right)^{2/3},
\end{displaymath} (6)

where E0i, Eei, ni and $\Delta$ti are the initial and remaining (i.e. at the end of part i) energy as well as the constant density and time interval in part i. E02 is equivalent to Ee1, while E01 can be derived from the requirement that the deflection time $t_{\rm D}$ of the electron beam must exceed the total duration $t_{\rm d}$ of the type III burst. Setting $t_{\rm D}=t_{\rm d}$, v0 (corresponding to E0) may be determined by $t_{\rm d}$ (Benz et al. 1992)

\begin{displaymath}v_{0}=\left(\frac{10^{20}nt_{\rm d}}{3.1}\right)^{1/3},
\end{displaymath} (7)

where $t_{\rm d}$ is a measurable quantity and n the density corresponding to the whole frequency range of the type III burst. For an exponential distribution of $n(h)=n_{0}{\rm e}^{-h/\lambda }$, however, provided $h/ \lambda$ is small enough, e.g. $h/\lambda\leq0.2$, then the local electron density may be approximately expressed by a linear function of $h/ \lambda$ i.e. $n\approx n_{0}(1-h/ \lambda)$, which is equivalent to its arithmetic mean ${\sim} n_{0}(1-h/2\lambda)$ . Therefore, $\overline{n}_{i}$ can be expressed by (ni+ni+1)/2, where $n_{i}=(f_{i}/(2\times 8.98\times10^{3}))^{2}$. Such an approximation is also applicable for our case because $h_{1}/ \lambda_{1}$ and $h_{2}/ \lambda_{2}$ calculated from Eq. (5) are 0.15 and 0.08. Therefore, ni and n in Eqs. (6) and (7) will be replaced by $\overline{n}_{i}$ and $\overline{n}$. For the sake of generality, an average density $\overline{n}\approx1\times10^{11}/\rm cm^{3}$ and an appropriate duration $t_{\rm d}=0.65$ s according to the typical duration range of a type III burst given by Isliker & Benz (1994) will be adopted in Eq. (7) for the evaluation of the initial energy E0. By means of Eq. (6) the energy loss and the corresponding velocity (i.e. voi, vei, and $\overline{v}_{i}$) for Parts 1 and 2 then can be estimated, where $\overline{n}_{i}$ in Eq. (6) is the arithmetic mean between density ni and ni+1, and can be obtained from the corresponding frequency. The quantities $\overline{n}_{1}$, $\Delta t_{1}$ and $\overline{n}_{2}$, $\Delta t_{2}$, for Parts 1 and 2 deduced from the observations are $1.042\times10^{11}/\rm cm^{3}$, 0.136 s and $1.169\times10^{11}/\rm cm^{3}$, 0.165 s, respectively. The local scale length $\lambda_{i}$ can be determined by expression (5)

\begin{displaymath}\lambda _{i}=h_{i}/ \ln \frac{n_{i+1}}{n_{i}},
\end{displaymath} (8)

where $h_{i}=\overline{v}_{i}\Delta t_{i}\cos\theta$. Based on the Eq. (1) the drift rate of Parts 1 and 2 can be determined by means of $\overline{v}_{i}$ and $\lambda_{i}$. The parameters and drift rates calculated are given in Table 1. It is evident that the drift rates calculated agree well with the observations.


  \begin{figure}
\par\includegraphics[width=8cm,clip]{3194.f4.eps}\end{figure} Figure 4: An exaggerating schematic of the 2nd V-shaped structure. The relevant parameters such as frequency f as well as the different part with differing drift rate are denoted by their subscriptions 1, 2, 3 and 4.
Open with DEXTER

3 Discussion

The frequency drift rate of type III or type III-like structure in the same drift direction looks the same in the dynamic spectrogram because of its small time scale and narrow bandwidth. Therefore, the beam velocity v and scale length $\lambda $ in Eq. (1) are usually considered as constant. However, the electron energy loss and hence the corresponding velocity variation, especially in the region with high electron density, must be considered for the estimation of electron lifetime corresponding to the deflection time $t_{\rm D}$ (usually replaced by duration $t_{\rm d}$). For a given E0 and n, the time interval $\Delta t$ (i.e. the life-time) can be determined by Eq. (4) i.e. $\Delta t=(E_{0}^{3/2}-E^{3/2})/(\frac{3}{2}K\sqrt{\frac{2}{m_{\rm e}}}n)$, where E corresponds to v, and $v_{0}\geq 0.1$ c has to be satisfied. To estimate the initial velocity v0 or initial energy E0 of an electron beam as well as its energy loss due to electron-electron collisions, the following approaches are adopted.

Setting $t_{\rm d}=t_{\rm D}$ in Eq. (7) then the initial velocity v0or energy E0 will be underestimated because $t_{\rm d}$ is shortened considerably by strong f-f absorption (Benz et al. 1992) and hence $t_{\rm d}<t_{\rm D}$. Therefore, a mean duration of 0.65 s averaged over the duration of typical type III bursts given by Isliker & Benz (1994) to determine the initial velocity v0or initial energy E0 is much more reasonable than the measurable duration. On the other hand, however, the initial electron velocity v0 in expression (1) will lead to an over-estimation of the scale length $\lambda $ if the drift rate is measurable. Therefore, in the course of derivation of the drift rates of Parts 1 and 2, an average velocity vi (i=1,2) instead of initial velocity v0 in Eq. (1) is suitable.

Considering an inhomogeneous region with two differing scale lengths $\lambda _{1}$ and $\lambda _{2}$, as well as different mean velocities $\overline{v}_{1}$ and $\overline{v}_{2}$, the estimated results of drift rates agree well with the observations. According to Eq. (8) $\lambda_{i}$ would be determined if $\theta $ is given. For example, $\theta =60^{0}$ and hence $\lambda _{1}$, $\lambda _{2}$ and h1h2 are about 5.5, $7.7\times10^{9}$ cm and 8.3, $6.3\times10^{8}$ cm, respectively. $\lambda
_{1}<\lambda _{2}$ represents the density gradient in Part 1 which is steeper than that in Part 2. $\lambda
_{1}<\lambda _{2}$ and $\overline{v}_{2}<\overline{v}_{1}$ will lead to a different drift rate of $({\rm d}f/{\rm d}t)_{1}>({\rm d}f/{\rm d}t)_{2}$. It is most likely that the jump in the drift rate is caused by the variations of both scale length and velocity, especially by velocity. The ratio of $\overline{v}_{2}$ to $\overline{v}_{1}$ has more effect on the jump of drift rate compared to the ratio of $\lambda _{1}$to $\lambda _{2}$ (see Table 1).

The V-shaped structure is formed by an electron beam moving downward along magnetic field lines to the mirror and then upward. Consequently, the normal drift rate $\left\vert ({\rm d}f/{\rm d}t)_{3}\right\vert $should be slower than the reverse drift rates $({\rm d}f/{\rm d}t)_{1}$ and $({\rm d}f/{\rm d}t)_{2}$ according to the electron velocity decrease with time due to energy loss. The observations, however, show that the normal drift rate $\left\vert ({\rm d}f/{\rm d}t)_{3}\right\vert $ is faster than the reverse drift rate $({\rm d}f/{\rm d}t)_{2}$. Either a small scale length or an increasing energy may cause this. According to the drift rate $({\rm d}f/{\rm d}t)_{3}=-3.53$ GHz/s and $v_{e2}\approx 0.1$ c given in Sect. 2.1 and Table 1, the scale length  $\lambda_{3}$ is about 109 cm which is much smaller than $\lambda _{1}$ and  $\lambda _{2}$ and is difficult to understand. Therefore, an increasing energy of the electron beam at the mirror should be taken into account. Considering an electron beam distribution with form of ${\rm d}F=f(E){\rm d}E$, then the mean energy of the electron beam may be expressed by

\begin{displaymath}E^{\prime }=\frac{\int Ef(E){\rm d}E}{\int f(E){\rm d}E},
\end{displaymath} (9)

where $E^{\prime
}$ is averaged over all electrons with energy E. For a power law distribution, for example, ${\rm d}F(E)=kE^{-\delta
}{\rm d}E$, the mean energy is

\begin{displaymath}E^{\prime }=\frac{\delta -1}{\delta -2}E_{\rm c},~~~~~~~2<\delta \leq
7
\end{displaymath} (10)

where $E_{\rm c}$ is a lower energy cut-off and $E^{\prime
}$ depends on $E_{\rm c}$. For a given $E^{\prime
}$, therefore, the lower energy cut-off can be determined by Eq. (10). For the sake of generality, setting $\delta=4.5$ and $E^{\prime }=E_{0}=50$ keV, then $E_{\rm c}\approx35$ keV. As seen from Table 1, the remaining energy Ee2 is only about 5 keV. Therefore, electrons with energy lower than 40 keV seem to be deflected, which is supported not only by Eq. (4) but also by observations. Based on Eq. (4), for a given density n the lift-time $\Delta t$ depends only on the initial energy E0, i.e. a lower energy E0 corresponds to a shorter lift-time $\Delta t$. Also as seen from Fig. 2 there is a strong diffusion within the range of turn-over frequency to $\sim $5.93 GHz, which is probably caused by the deflected electrons of lower energy. Therefore, the mean energy $E^{\prime
}$ of the electron beam at mirror would be 55 keV if $E_{\rm c}=40$ keV. It follows that the increasing energy is due to fewer electrons of lower energy. As a rough estimation, the normal drift rate $({\rm d}f/{\rm d}t)_{3}=-3.96$ GHz/s would be obtained by Eq. (1), where $\overline{v}_{3}$ is $\sim $0.43 c and $\lambda_{3}$ is determined by Eq. (8). It is larger than observations, -3.53 GHz/s, but such a trend of increasing drift rate in Part 3 is acceptable.

4 Conclusions

Fine structures in microwave emission are still poorly known. A spectrograph with high temporal and spectral resolution is needed to study the fine structures in the microwave band. The fine structures presented here consist of a descending branch followed by an ascending one to form a V-shaped structure, which is located between 5.57 and 6.26 GHz and lastsabout 0.4 s. The significant features of the V-shaped structure are: a distinct jump of drift rate from 3.22 to 1.53 GHz/s in the reverse drift component; a fastest drift rate of -3.53 GHz/s in the normal drift component.

Type III or type III-like events refer to the plasma-emitting process of moving electron beams. The most significant feature of type III bursts is the frequency drift. The drift rate can be determined by both excited velocity covering the range of $\beta
=0.1{-}0.6$ and scale length as a measure of the density gradient. As a bi-directional frequency drift event, the V-shaped structure cannot be explained by the common beam model due to different drift rates in the same drift component. The energy loss of the electron beam and the three-component solar atmospheric model with different scale lengths are proposed to explain the different drift rates. The variations of both local scale length $\lambda
_{2}>\lambda _{1}$ and the mean velocity $\overline{v}_{1}>\overline{v}_{2}$, especially $\overline{v}_{1}>\overline{v}_{2}$, are responsible for the jump of drift rate in the reverse drift component. While the fastest drift rate in the normal drift component is interpreted by means of the increase of the average energy of the electron beam due to the decrease of electrons lower than $\sim $40 keV caused by energy loss at the mirror. The local scale length $\lambda
_{1}<\lambda _{2}$ means an asymmetry of local inhomogeneity, which is probably more close to the real inhomegeneity.

Acknowledgements
The authors are grateful to the referee for helpful comments. This work was supported by the National Natural Science Foundation under grants Nos. 10333030 and 10273025 and by a key Project of the Chinese Academy of Science.

References

 

Copyright ESO 2004