A model of impulsive loop flares caused by anomalous heat conduction

Numerical simulation is made of the impulsive loop flare caused by transient heat conduction along the loop with an applied axial electric current.It is assumed that a segment near the top of the coronal loop is heated to above 10⁷ K by a heat input that is small compared with the total flare energy, which is given by the magnetic energy of the initial current. Due to the heat conduction, a hump appears in the velocity distribution of electrons, which may excite electron plasma waves with a sufficiently high intensity to cause an anomalous resistivity, as shown theoretically in a previous paper. In that paper, an effect of the plasma waves on the dynamics of electrons was taken into account consistently, but an anomalous heating due to an ohmic dissipation of the initial current under the anomalous resistivity was not taken into account.The aim of the present study is to study the subsequent dynamics of the heated gas caused by the anomalous heating, but in order to avoid an unpractically long computation time, the energy density of the plasma waves is estimated by the energy density of electrons in the velocity hump, without taking into account the effect of the plasma waves consistently in the dynamics of the electrons.The initial current starts to decay gradually by an ohmic dissipation under the anomalous resistivity occurring near the top of the loop to heat this region more. The enhanced heat conduction causes the velocity humps in a wider location. Consequently, the anomalous heating continues and spreads in a self-generating way even after the end of the initial minor heating. Thus the temperature near the loop top becomes above 10⁸ K and the high-temperature region spreads in both directions along the loop with such a high speed as (2–3) × 10⁴ km s^–1, which is nearly equal to the speed of flux-limited heat conduction. On the other hand, induced electric field estimated from the anomalous resistivity is 3.3 × 10⁷ V at the termination of the present simulation, under the modest initial current of 1.5 A m^–2.X-ray emissions expected from the present model loop, show three sources, two footpoints with unequal brightness and a coronal source expanding along the loop in both directions.     

Damping of slow magnetoacoustic waves in an inhomogeneous coronal plasma

We study the propagation and dissipation of slow magnetoacoustic waves in an inhomogeneous viscous coronal loop plasma permeated by uniform magnetic field. Only viscosity and thermal conductivity are taken into account as dissipative processes in the coronal loop. The damping length of slow-mode waves exhibit varying behaviour depending upon the physical parameters of the loop in an active region AR8270 observed by TRACE. The wave energy flux associated with slow magnetoacoustic waves turns out to be of the order of 10⁶ erg cm^?2 s^?1 which is high enough to replace the energy lost through optically thin coronal emission and the thermal conduction below to the transition region. It is also found that only those slow-mode waves which have periods more than 240s provide the required heating rate to balance the energy losses in the solar corona. Our calculated wave periods for slow-mode waves nearly match with the oscillation periods of loop observed by TRACE.     

Coronal loop heating by discrete Alfvén waves

We have modeled the solar coronal active loop heating by discrete Alfvén waves. Discrete Alfvén waves (DAW) are a new class of Alfvén waves which can be described by the two-fluid model with finite ion-cyclotron frequency, or the MHD model with plasma current along the magnetic field line as shown by Appert, Vaclavik, and Villar (1984). We have modeled the coronal loop as a semi-toroidal plasma with the major toroidal radius much larger than the plasma radius. We have shown that the absorption of discrete Alfvén waves by the plasma through viscosity can account for at least 30% of the coronal heating rate density of 10^–4 J m^–3 s^–1.     

A model of hot loops associated with solar flares

A dynamical model is proposed for the formation of soft X-ray emitting hot loops in solar flares. It is examined by numerical simulations how a solar model atmosphere in a magnetic loop changes its state and forms a hot loop when the flare energy is released in the form of heat liberation either at the top part or around the transition region in the loop.When the heat liberation takes place at the top part of the loop which arches in the corona, the plasma temperature around the loop apex rises rapidly and, as the result, the downward thermal conductive flux is increased along the magnetic tube of force. Soon after the thermal conduction front rushes into the upper chromosphere, a local peak of pressure is produced near the conduction front and the chromospheric material begins to expand into the corona to form a high-temperature (10⁷ K-3 × 10⁷ K at the loop apex) and high-density (10¹⁰ cm^–3-10¹¹ cm^–3 at the loop apex) loop. The velocity of the expanding material can reach a few hundred kilometres per second in the coronal part. The thermal conduction front also plays a role of piston pushing the chromospheric material downward and gives birth to a shock wave which propagates through the minimum temperature region into the photosphere. If, on the other hand, the heat source is placed around the transition region in the loop, the expansion of the material into the corona occurs from the beginning of the flare and the formation process of the hot loop differs somewhat from the case with the heat source at the top part of the loop.Thermal components of radiations emitted from flare regions, ranging from soft X-rays to radio wavelengths, are interpreted in a unified way by using physical quantities obtained as functions of time and position in our flare loop model as will be discussed in detail in a following paper.     

Steady heat conduction in coronal loop unstable against plasma instability

The Fokker-Planck equation is numerically solved to study the electron velocity distribution under steady heat conduction with an applied axial electric current in a model coronal loop.If the loop temperature is so high that the electron mean-free path is longer than the local temperature scale height along the loop, a velocity hump appears at about the local thermal electron velocity. The hump is attributed to cooler electrons moving up the temperature gradient to compensate for the runaway electrons moving down the gradient. If the ratio between the mean free path and temperature scale height is greater than about 2, negative absorption for the plasma waves can appear (waves grow). This effect is enhanced by the presence of axial electric current in the half of the coronal loop in which the electrons carrying the current are drifting up the temperature gradient. Thus, the plasma instability may occur in the coronal elementary magnetic flux tubes. Although the present paper is limited to show the critical condition and linear growth rate of the instability, the following scenarios may be inferred.If the flux tubes change from marginally stable to unstable against the plasma instability, due to an increase in the loop temperature, anomalous resistivity may suddenly appear because of the growth of plasma waves. Then a high axial electric field is induced that may accelerate particles. This could be the onset of impulsive loop flares.For a low electric current, if the loop temperature is sufficiently high to give the negative absorption for the plasma waves in a large part of the coronal loop, steady plasma turbulence may originate. This could be a source for the type I radio noise storm.     

Heat flux saturation in hydrodynamic soft X-ray solar flare plasmas

The role of heat flux limitation in soft X-ray emitting solar flare plasmas is considered. Simple analytic arguments suggest that flux limitation is likely to be important during the explosive heating phase, even for relatively modest coronal energy fluxes (say 10⁹ erg cm^-2 s^-1). This conclusion is reinforced by a detailed flare loop simulation of the heating phase. Since flux saturation effectively bottles up the coronal heat flux, mass motions now assume a dominant role in transferring energy from the coronal flare source to the lower transition region. The mass-energy exchange between the corona and chromosphere produces dramatic changes in the thermal structure of the plasma which are reflected in the differential emission measure profile of the flaring loop.     

Acceleration of the solar wind by whistler waves from coronal holes

We investigate the possibility of an additional acceleration of the high speed solar wind by whistler waves propagating outward from a coronal hole. We consider a stationary, spherically symmetric model and assume a radial wind flow as well as a radial magnetic field. The energy equation consists of (a) energy transfer of the electron beam which excites the whistler waves, and (b) energy transfer of the whistler waves described by conservation of wave action density. The momentum conservation equation includes the momentum transfer of two gases (a thermal gas and an electron beam). The variation of the temperature is described by a polytropic law. The variation of solar wind velocity with the radial distance is calculated for different values of energy density of the whistler waves. It is shown that the acceleration of high speed solar wind in the coronal hole due to the whistler waves is very important. We have calculated that the solar wind velocity at the earth's orbit is about equal to 670 km/sec (for wave energy density about 10^?4 erg cm^?3 at 1.1R⊙). It is in approximate agreement with the observed values.     

Steady models for the hard X-ray loops in the solar corona

In some gradual hard X-ray bursts with high intensity, hard X-ray source (15–40 keV) is steadily located in the corona along with softer X-ray source (5–10 keV).Two stationary models, high density and high temperature models, are proposed to solve the difficult problem of confinement of hot (or nonthermal) plasma in the direction of the magnetic field along the loops in the corona. In both models, an essential point is that the effective X-ray source is composed of fine dense filamentary loops imbeded in a larger rarefied coronal loop, and the electron number density in the filaments is so high as 10¹¹–10¹² cm^-3. If the density is so high heat conduction can be as reasonably small as of the order of 10²⁷ erg s ^-1 for the given emission measures of observed X-rays, since the required cross-sectional area is small and also classical conduction is valid. Collisional confinement of thermal tail, and nonthermal electrons if any, up to 50–60 keV in the filaments is also possible, so that the hard X-ray images can be loop like structure instead of double source (foot points).High density model is applicable to the coronal filamentary loops with temperature T _m < 5 × 10⁷ K at the loop summit. The heat flow from the summit downwards is lost almost completely by the radiation from the loop during the conduction to the foot points. A continuous energy release is assumed near the summit to maintain the stationary temperature T _m, and pressure balance is maintained along the loop. In this model, the number density at the summit is given by n _m - 10⁶ T _m ² /s_m, where s _m is the length of the loop from the summit to the foot point, and the distribution of temperature and density along the loop are given by T = T _m(s/s_m)^1/3 and n = n _m(s/s_m)^-1/3, respectively.High temperature model is applicable to the filamentary loops with higher temperature up to about 10^8.5 K and comparatively lower number density as 10¹¹ cm^-3 for the requirement of magnetic confinement of the hot plasma in radial direction. The radiation from the loop is negligibly small in this model so that the heat flux is nearly conserved down to the foot points. In this case, temperature gradient is smaller than that of the high density model, depending on the tapering of the magnetic bottle.In both models, the differential emission measure is maximum at the highest temperature T _m and the brightness distribution along the loop shows a maximum around the summit of the loop if some magnetic tapering is taken into account.     

Dynamics of electron beams in a coronal loop and the hard X-ray burst

A numerical simulation has been made for the dynamics of non-thermal electrons (> 10keV) injected with spatial, temporal and velocity distributions into a model coronal loop. The time variations of the spatial intensity distribution and the spectrum for the expected hard X-rays are computed for many models in order to find the important physical parameters for those characteristics.The most important one is the column density of plasma, CD, along the loop. If CD is smaller than 10²⁰ cm^–2, the expected X-rays behave like the solar impulsive hard X-ray bursts, that is the spatial maximum of X-rays shifts to the top of the loop in the later phase of the burst accompanying a spectral softening. On the other hand, if CD is greater than this value, quasi-steady decay appears in the later phase. In this case the intensity distribution of X-rays above about 20 keV along the loop shows a broad maximum away from the loop top giving an extended spatial distribution of hard X-rays, and spectral hardness is kept constant. These characteristics are similar to the solar gradual hard X-ray bursts (the so-called extended burst which is not a hot thermal gradual burst).     

Soft X-ray oscillations from AT Mic: Flare plasma diagnostics

We consider the flare oscillations from the active red dwarf AT Mic detected with the XMM-Newton space observatory in the soft X-ray energy range (0.2–12 keV). Following Mitra-Kraev et al. (2005a), we associate the observed oscillations with a period of ≈750 s with the excitation of a standing slow magnetoacoustic (SMA) wave in a coronal loop. The damping of flare loop SMA oscillations is shown to be governed by electron thermal conduction. We have estimated the plasma density (≈3 × 10¹⁰ cm^?3) and the minimum magnetic field strength (≈100 G) in the region of flare energy release. The adopted model is consistent with the results of a spectral analysis of the soft X-ray emission. The piston mechanism is assumed to be responsible for the excitation of loop SMA oscillations.     

Solar Wind and Chromospheric Network

A physical model of the transition region, including upflow of the plasma in magnetic field funnels that are open to the overlying corona, is presented. A numerical study of the effects of Alfvén waves on the heating and acceleration of the nascent solar wind originating in the chromospheric network is carried out within the framework of a two-fluid model for the plasma. It is shown that waves with reasonable amplitudes can, through their pressure gradient together with the thermal pressure gradient, cause a substantial initial acceleration of the wind (on scales of a few Mm) to locally supersonic flows in the rapidly expanding magnetic field trunks of the transition region network. The concurrent proton heating is due to the energy supplied by cyclotron damping of the high-frequency Alfvén waves, which are assumed to be created through small-scale magnetic activity. The wave energy flux of the model is given as a condition at the upper chromosphere boundary, located above the thin layer where the first ionization of hydrogen takes place.Among the new numerical results are the following: Alfvén waves with an assumed f ^-1 power spectrum in the frequency range from 1 to ⁴ Hz, and with an integrated mean amplitude ranging between 25 and 75 km s⁴, can produce very fast acceleration and also heating through wave dissipation. This can heat the lower corona to a temperature of 5× 10⁵ K at a height of h=12,000 km, starting from 5× 10⁴ K at h=3000 km. The resulting thermal and wave pressure gradients can accelerate the wind to speeds of up to 150 km s^-1 at h=12,000 km, starting from 20 km s^-1 at h=3000 km in a rapidly diverging flux tube. Thus the nascent solar wind becomes supersonic at heights well below the classical Parker-Type sonic point. This is a consequence of the fact that any large wave-energy flux, if it is to be conducted through the expanding funnel to the corona, implies the building-up of an associated wave-pressure gradient. Because of the diverging field geometry, this might lead to a strong initial acceleration of the flow. There is a multiplicity of solutions, depending mainly on the coronal pressure. Here we discuss two new (as compared with a static transition region model) possibilities, namely that either the flow remains supersonic or slows down abruptly by shock formation, which then yields substantial coronal heating up to the canonical 10⁶ K for the proton temperature.     

EUV Analysis of a Quasi-static Coronal Loop Structure

Decaying active region 10942 is investigated from 4:00?–?16:00 UT on 24 February 2007 using a suite of EUV observing instruments. Results from Hinode/EIS, STEREO and TRACE show that, although the active region has decayed and no sunspot is present, the physical mechanisms that produce distinguishable loop structures, spectral line broadening, and plasma flows still occur. A coronal loop that appears as a blue-shifted structure in Doppler maps is apparent in intensity images of log(T)=6.0?–?6.3 ions. The loop structure is found to be anti-correlated with spectral line broadening generally attributed to non-thermal velocities. This coronal loop structure is investigated physically (temperature, density, geometry) and temporally. Light curves created from imaging instruments show brightening and dimming of the loop structure on two different time scales; short pulses of 10?–?20?min and long duration dimming of two?–?four hours until its disappearance. The coronal loop structure, formed from relatively blue-shifted material that is anti-correlated with spectral line broadening, shows a density of 10¹⁰ to 10^9.3?cm^?3 and is visible for longer than characteristic cooling times. The maximum non-thermal spectral line broadenings are found to be adjacent to the footpoint of the coronal loop structure.     

Heating of coronal loops by fast mode MHD waves

A possible mechanism for the formation and heating of coronal loops through the propagation and damping of fast mode waves is proposed and studied in detail. Loop-like field structures are represented by a dipole field with the point dipole at a given distance below the solar surface. The density of the medium is determined by hydrostatic equilibrium along the field lines in an isothermal atmosphere. The fast mode waves propagating outward from the coronal base are refracted into regions with a low Alfvén speed and suffer collisionless damping when the gas pressure becomes comparable to the magnetic pressure. The propagation and damping of these waves are studied for three different cases: a uniform density at the coronal base, a density depletion within a given flux tube, and a density enhancement within a given flux tube. The fast mode waves are found to be important in the formation and heating of the loops if the wave energy flux density is of the order 10⁵ ergs cm^-2 s^-1 at the coronal base.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

Effect of Variable Background on an Oscillating Hot Coronal Loop

We investigate the effect of a variable, i.e. time-dependent, background on the standing acoustic (i.e. longitudinal) modes generated in a hot coronal loop. A theoretical model of 1D geometry describing the coronal loop is applied. The background temperature is allowed to change as a function of time and undergoes an exponential decay with characteristic cooling times typical for coronal loops. The magnetic field is assumed to be uniform. Thermal conduction is assumed to be the dominant mechanism for damping hot coronal oscillations in the presence of a physically unspecified thermodynamic source that maintains the initial equilibrium. The influence of the rapidly cooling background plasma on the behaviour of standing acoustic (longitudinal) waves is investigated analytically. The temporally evolving dispersion relation and wave amplitude are derived by using the Wenzel–Kramers–Brillouin theory. An analytic solution for the time-dependent amplitude that describes the influence of thermal conduction on the standing longitudinal (acoustic) wave is obtained by exploiting the properties of Sturm–Liouville problems. Next, numerical evaluations further illustrate the behaviour of the standing acoustic waves in a system with a variable, time-dependent background. The results are applied to a number of detected loop oscillations. We find a remarkable agreement between the theoretical predictions and the observations. Despite the emergence of the cooling background plasma in the medium, thermal conduction is found to cause a strong damping for the slow standing magneto–acoustic waves in hot coronal loops in general. In addition to this, the increase in the value of thermal conductivity leads to a strong decay in the amplitude of the longitudinal standing slow MHD waves.     

The temperature-density structure of coronal loops in hydrostatic equilibrium

The temperature and density structure are computed for a comprehensive set of coronal loops that are in hydrostatic and thermal equilibrium. The effect of gravity is to produce significant deviations from the usual uniform-pressure scaling law (T(pL) ^1/3) when the loops are taller than a scale height. For thermally isolated loops it lowers the pressure throughout the loop, which in turn lowers the density significantly and also the temperature slightly; this modifies the above scaling law considerably. For more general loops, where the base conductive flux does not vanish, gravity lowers the summit pressure and so makes the radiation decrease by more than the heating. This in turn raises the temperature above its uniform pressure value for loops of moderate length but lowers it for longer loops. A divergence in loop cross-section increases the summit temperature by typically a factor of 2, and decreases the density, while an increase in loop height (for constant loop length) changes the temperature very little but can halve the density.One feature of the results is a lack of equilibrium when the loop pressure becomes too large. This may explain the presence of cool cores in loops which originally had temperatures below 2 × 10⁶ K. Loops hotter than 2 × 10⁶ K are not expected to develop cool cores because the pressure necessary to produce non-equilibrium is larger than observed.     

Evolution of a flaring loop after injection of energetic electrons

For the November 5, 1980 flare it is investigated how the plasma in a large flaring loop responds to the injection of energetic electrons. Observations are compared with the results of a one-dimensional numerical simulation. For the simulation it is assumed that at the time the injection is started, the plasma is in an equilibrium state with a constant pressure along the loop and conductive heating compensated by radiative losses. Especially important for the evolution of the impulsively heated plasma is the penetration depth of the fast electrons compared to the depth of the transition layer. Both parameters are known from the observations. The injected energy is 2.6 × 10¹¹ ergs cm ^?2 in 30 s (as derived from the hard X-ray observations) and computations show that the high temperature plasma of the loop responds to it with upward motions of about 50 km s^?1, i.e. with velocities much smaller than the ion sound speed (≈ 500km s^?1). The heating of the plasma due to the absorption of beam energy can be understood using a constant density approximation. After the heating phase the plasma returns in about 5 min to its initial state by conductive cooling. The downward conducted energy is radiated away in the transition zone. The numerical simulation shows that impulsive heating by non-thermal electrons only does not explain the observed large increase in the density of the loop during the flare. It is therefore required that continuous energy and/or mass input occur after the impulsive phase.     

Initial phase of chromospheric evaporation in a solar flare

In this paper we discuss the initial phase of chromospheric evaporation during a solar flare observed with instruments on the Solar Maximum Mission on May 21, 1980 at 20:53 UT. Images of the flaring region taken with the Hard X-Ray Imaging Spectrometer in the energy bands from 3.5 to 8 keV and from 16 to 30 keV show that early in the event both the soft and hard X-ray emissions are localized near the footpoints, while they are weaker from the rest of the flaring loop system. This implies that there is no evidence for heating taking place at the top of the loops, but energy is deposited mainly at their base. The spectral analysis of the soft X-ray emission detected with the Bent Crystal Spectrometer evidences an initial phase of the flare, before the impulsive increase in hard X-ray emission, during which most of the thermal plasma at 10⁷ K was moving toward the observer with a mean velocity of about 80 km s^-1. At this time the plasma was highly turbulent. In a second phase, in coincidence with the impulsive rise in hard X-ray emission during the major burst, high-velocity (370 km s^-1) upward motions were observed. At this time, soft X-rays were still predominantly emitted near the loop footpoints. The energy deposition in the chromosphere by electrons accelerated in the flare region to energies above 25 keV, at the onset of the high-velocity upflows, was of the order of 4 × 10¹⁰ erg s^-1 cm^-2. These observations provide further support for interpreting the plasma upflows as the mechanism responsible for the formation of the soft X-ray flare, identified with chromospheric evaporation. Early in the flare soft X-rays are mainly from evaporating material close to the footpoints, while the magnetically confined coronal region is at lower density. The site where upflows originate is identified with the base of the loop system. Moreover, we can conclude that evaporation occurred in two regimes: an initial slow evaporation, observed as a motion of most of the thermal plasma, followed by a high-speed evaporation lasting as long as the soft X-ray emission of the flare was increasing, that is as long as plasma accumulation was observed in corona.     

Coronal loop heating by Alfvén waves

The excitation and dissipation of global and surface Alfvén waves and their conversion into kinetic Alfvén waves have been analyzed for solar coronal loops using a cylindrical model of a magnetized plasma. Also the optimal conditions for coronal loop heating regimes with density of dissipated power 10³ erg cm^–3 s^–1 by the new scheme named combined Alfvén wave resonance are found. Combined Alfvén wave heating regime appears when the global Alfvén wave is immersed into the Alfvén continuum with the condition of not-so-sharp distribution of axial current.Instituto de Matemática, Universidade Federal Fluminense, Niterói, RJ, Brazil