A&A 406, 345-356 (2003)
Department of Physics, University of Wales, Aberystwyth, SY23 3BZ, UK
Received 31 October 2002 / Accepted 28 April 2003
It is assumed that magnetic flux tubes are strongly concentrated at the boundaries of supergranule convection cells. A power law spectrum of high frequency Alfvén waves with a spectral index -1 originating from the sun is assumed to supply all the energy needed to energize the plasma flowing in such magnetic flux tubes. At the high frequency end, the waves are eroded by ions due to ion cyclotron resonance. The magnetic flux concentration is essential since it allows a sufficiently strong energy flux to be carried by high frequency ion cyclotron waves and these waves can be readily released at the coronal base by cyclotron resonance. The main results are: 1. The waves are capable of creating a steep transition region, a hot corona and a fast solar wind if both the wave frequency is high enough and the magnetic flux concentration is sufficiently strong in the boundaries of the supergranule convection zone. 2. By primarily heating alpha particles only, it is possible to produce a steep transition region, a hot corona and a fast solar wind. Coulomb coupling plays a key role in transferring the thermal energy of alpha particles to protons and electrons at the corona base. The electron thermal conduction then does the remaining job to create a sharp transition region. 3. Plasma species (even ions) may already partially lose thermal equilibrium in the transition region, and minor ions may already be faster than protons at the very base of the corona. 4. The model predicts high temperature alpha particles ( K) and low proton temperatures ( K) between 2 and 4 solar radii, suggesting that hydrogen Lyman lines observed by UVCS above coronal holes may be primarily broadened by Alfvén waves in this range.
Key words: Sun: corona - Sun: solar wind - Sun: transition region
It has been well known for several decades that the solar wind is the expansion of the hot corona along open magnetic field flux tubes. However, the exact mechanisms responsible for the creation of a hot corona and the acceleration of the fast solar wind continue to be an unresolved issue (Axford et al. 1999; Marsch 1999). Since the fast solar wind streams and the corona are connected, there has been a growing awareness that the two problems should be treated in a self-consistent way. It is also natural to seek a single theory to account for the coronal heating and the acceleration of the fast solar wind along open magnetic field fluxes. However, the Sun and the inner corona are far away, the environment near the Sun is harsh for the survival of space probes, so in situ measurements of the plasma parameters in the near Sun region are still not available. Without in situ measurements, it is difficult to tell what is responsible for the coronal heating and the acceleration of the fast solar wind. However, as the fast solar wind cruises into interplanetary space, manmade satellites have been able to detect the properties of the fast (and slow) solar wind streams over the past few decades. These direct measurements have shown that Alfvén waves are ubiquitous in fast solar wind streams (e.g., Belcher & Davis 1971; Smith et al. 1995). They provide an attractive mechanism for the solar wind heating and acceleration if their high frequency counterparts, whose energy can be directly absorbed by the solar wind plasmas, exist in the inner corona.
However, power spectra of magnetic field fluctuations in the solar wind, most likely of Alfvénic nature, are typically of the form , where k is the wavenumber and the spectral index. Since observations yield (Bavassano et al. 1982), most of the power in the Alfvénic fluctuations is in the frequency range below the ion cyclotron frequency. Hence no substantial influence of these waves on ions can be expected (see Isenberg & Hollweg 1983; Marsch et al. 1982). This difficulty can be overcome if the Sun launches waves in the ion cyclotron resonance range leading to the heating and acceleration of the solar wind (e.g., Axford & McKenzie 1995; McKenzie et al. 1995, 1997; Tu & Marsch 1997), or if there exists some mechanism which replenishes the resonant wave power as soon as it is absorbed by the particles.
To tap into the energy available in Alfvénic fluctuations, Hollweg (1986) assumed that a nonlinear cascade process, namely, the Kolmogorov dissipation of Alfvénic turbulence, transports energy from the low-frequency waves to the ion cyclotron resonant range, where it is picked up by the plasma through the cyclotron resonance. Hollweg & Johnson (1988) extended this idea to a two-fluid solar wind model. Their results yielded the right qualitative features: a steep temperature rise in the transition region to over 106 K in the corona and acceleration to supersonic speeds. These studies showed that the Kolmogorov dissipation rate is not capable of creating a transition region, a hot corona and a fast solar wind simultaneously, even though it seems that the mechanism is able to produce a fast solar wind (Isenberg 1990; Li et al. 1999).
Recent remote sensing observations from the Ultraviolet Coronagraph Spectrometer (UVCS) and Solar Ultraviolet Measurements of Emitted Radiation (SUMER) on board Solar and Heliospheric Observatory (SOHO) have enhanced our interests in the ion cyclotron resonance for coronal heating and solar wind acceleration. SUMER spectral line observations at the base of the corona in a southern coronal hole indicate that the temperature of minor ions decreases with increasing mass per charge (Tu et al. 1998). The temperature difference of various ions is more pronounced in the inner corona. UVCS observations reveal that ions not only have greater than mass-proportional temperatures in the inner corona, they are much hotter than electrons and are highly anisotropic as well (Kohl et al. 1998; Li et al. 1999; Cranmer et al. 1999). By recognizing the importance of the pumping by the second C II line in the interpretation of the UVCS observations, Li et al. (1998) concluded that the outflow speed of O+5 is about 400 at 3 RS, already faster than protons, and these oxygen ions also have a large temperature anisotropy at the same distance. These recent advances strongly indicate that ion cyclotron resonance may indeed play an important role in the coronal heating (at least along open field lines), thus enhancing the interest in ion cyclotron resonance. Various models have been published to explain observations and make predictions using ion cyclotron resonance (see a thorough review by Hollweg & Isenberg 2002 and the references).
One similar approach to tackle the lack of wave power at high frequencies in the solar wind is to include the non-linear interaction between outward propagating and inward propagating Alfvén waves. The non-linear interaction can result in a strong cascade. The idea was introduced by Tu et al. (1984). The cascade is able to explain the magnetic field spectrum evolution and to account for the proton heating beyond 10 solar radii (Tu 1987, 1988). In a series of papers, Hu and his co-authors applied this idea to the inner corona, and found that the cascade indeed is able to account for the heating and acceleration of fast solar wind (including minor ions) originating in coronal holes (Hu et al. 1999, 2000; Hu & Habbal 1999). In these studies, a hot corona was treated as a boundary condition, the creation of a steep transition region and a hot corona were not considered. Although these turbulence cascade models are promising, it is still not clear what is the source of the necessary inward propagating waves and whether these inward propagating waves exist at all in the inner corona. It is also not clear if the non-linearity of waves in the corona is strong enough to produce cascades since wave amplitudes are small.
High frequency Alfvén waves have also been suggested (Axford & McKenzie 1995), and been adopted to produce a steep transition region and a hot corona (Marsh & Tu 1997a,b; Hackenberg et al. 2000), or to accelerate the fast solar wind (Tu & Marsch 1997) by adopting a two fluid approach. Kinetic models have also been developed by assuming a relatively hot corona boundary (Vocks & Marsch 2001, 2002). It has been found that plasma species may already lose thermal equilibrium at the very base of the corona or the top of transition region. The two fluid models developed by Hackenberg et al. (2000) can only reach 2 solar radii. The kinetic models developed by Vocks & Marsch (2002) become numerically unstable at a few solar radii before reaching cruising speed and the preferential acceleration of minor ions is not accounted for. More recently, the above-mentioned two fluid transition region models were extended to a three fluid approach: the most abundant minor ions alpha particles were included in model calculations (Li 2002). It was shown that by heating alpha particles alone, a transition region, and a hot corona can be formed due to Coulomb coupling between alpha particles and protons, Coulomb coupling between protons and electrons, and the electron heat flux. Considering the strong evidence of ion cyclotron heating in the extended corona above polar coronal holes as mentioned above, one would naturally ask: can these high frequency Alfvén waves produce a self-consistent solution of a transition region, a hot corona and a fast wind if alpha particles are considered?
The goal of this paper is to develop a single coherent solution of a steep transition region, a hot corona and a fast solar wind stream using one single mechanism and including alpha particles. In this paper, it is assumed that an outward propagating power law spectrum of high frequency Alfvén waves (reaching ion cyclotron frequency) already exists below the transition region (Axford & McKenzie 1995; Tu & Marsch 1997). Physically, there are no "boundaries" between the corona and the solar wind. This one dimensional fluid model includes electrons, protons and alpha particles. A power law spectrum of Alfvén waves is assumed to supply all the energy needed to heat the corona and accelerate the solar wind. It is found that such a mechanism is capable of creating a steep transition region, a hot corona and a fast solar wind. Readers may also refer to the work by Hansteen & Leer (1995), Hansteen et al. (1997) and Lie-Svendsen et al. (2002) who also considered chromosphere and ionization in their models. In those models, different heating mechanisms were used. In particular, Hansteen et al. (1997) have also included alpha particles in their models.
Standard three-fluid plasma transport equations for electrons, protons and alpha particles are considered, including heat conduction and radiation losses. We consider an isotropic, charge-neutral plasma. The one-dimensional model studied here is the same as those in Li et al. (1997) and Hu & Habbal (1999), and the governing equations will not be repeated.
It is assumed that the magnetic flux tubes are strongly
concentrated at the boundaries of the supergranule convection
cells (Gabriel 1976; Marsch & Tu 1997; Hackenberg et al. 2000; Li
2002). The plasma moves along the open magnetic field lines,
and eventually builds up the solar wind. The flow tube cross section area is
a(r)=r2 f(r). We
assume that the magnetic flux tubes in this paper have a two-step
In this study the Coulomb logarithm
is taken to be 21,
it is assumed that plasmas in the transition region are collision dominated and
classical thermal conduction parallel to the magnetic field can be
used for the electrons, protons and alpha particles. However, the
heat flux of protons rapidly becomes negligible with increasing
heliocentric distances reflecting the collision-less nature of the
proton gas (Li 1999; Olsen & Leer 1999). To simulate this effect,
The radiative energy loss in the electron energy equation is assumed to have the form parameterized by Rosner et al. (1978) for an optically thin medium. It strongly depends on the density of electrons.
When high frequency Alfvén waves propagate along open coronal
magnetic field lines, ion cyclotron resonance may produce an upper
on the frequency of Alfvén waves by a sweeping
mechanism (Marsch & Tu 1997). Following Tu (1987) and Hu & Habbal
is taken as
The wave power spectrum density, P(f,r), is related to the
magnetic field variance
Alfvén wave pressure by
The treatment of the wave spectral evolution closely follows Hu & Habbal (1999), however non-linear cascade is ignored in this study. The plasma in a magnetic flux tube is only heated by a sweeping mechanism (Marsch & Tu 1997a). The total dissipation rate from this mechanism is divided between the acceleration and heating for different ion species determined by the microphysics of the resonant interaction. For a Maxwellian distribution of a thermally isotropic plasma interacting with a spectrum of parallel left-hand-polarized ion cyclotron waves, the resonant acceleration and heating rate of species i are given by Dusenberg & Hollweg (1981), Marsch et al. (1982) and also Hu & Habbal (1999). In the evaluation of the integrals involved in the quasi-linear terms, the wave power spectrum and the dispersion relation in the dissipation range must be known. Due to the strong magnetic field in the transition region, the plasma beta value is small. Hence cold plasma approximation for the dispersion of these waves can be adopted. Since ion cyclotron resonance is very efficient in dissipating wave power in small scales, is assumed in this study. It will be shown, in the transition region and near Sun region, that the outflow velocity difference between alpha particles and protons is small compared with Alfvén speed, as cold plasma dispersion effectively prohibits cyclotron resonance between waves and protons.
The time-dependent equations of the three-fluid plasma system are numerically integrated with variable step sizes in space and time using a fully implicit scheme described by Hu et al. (1997) and Li et al. (1997). Model calculations start from a point in the lower transition region where the temperature is K and extend to 50 . The radial step dr gradually increases from 35 m at the lower boundary, to about at the top of the computational domain, and 1100 grid points are used in calculations. The maximum error in any of the conserved quantities is less than 1% once a steady state is reached.
In this section, model calculations are presented. Calculations
starting temperature is chosen so that both hydrogen and helium
almost fully ionized, and complications involved with ionization
can be avoided (Hansteen et al. 1997). The corresponding electron
density at the boundary is
adopted electron thermal pressure is
This pressure is typical in
coronal holes according to (Warren & Hassler 1998). Warren &
Hassler (1998) found that thermal pressure is roughly a constant
in the transition region when the
electron temperature varies from
to 106 K. Our
computations end at 52.5 solar radii.
At 52.5 solar radii, the
fast solar wind should be almost fully developed, and beyond 10 solar radii, the cascading effect may not be ignored (Tu 1987).
Three models will be presented in this study, some common
parameters of these models at the base of the open magnetic flux
tube and the flux tube geometry are listed in Table 1.
We first present model calculations with strong magnetic flux concentration in the supergranule convection zone boundaries. Parameters at the lower computation boundary used in this study are listed in Table 1. The flux tube geometry parameters are as follows: is the maximum expansion factor of the cross-section of a flux tube in the transition region, is the heliocentric distance at which the major expansion takes place. The parameter represents the length over which the expansion occurs. This first super-radial expansion is complete by about 15 000 km. The parameters for the second super-radial expansion are: , and . This expansion is complete at about 3 , above 3 , the solar wind is almost radial.
The magnetic field at the bottom of the funnel is 123.5 G. This
means that we assume the coronal funnels considered here may
originate from deep in the photosphere (Dowdy et al. 1986). The
low-frequency end of the wave spectrum is fixed at
Hz. The wavelength of the waves below this frequency
is very large compared to the scale height in the transition
region and corona. Hence those waves are ignored. The
high-frequency end of the wave spectrum
treated as a free parameter. This parameter determines where ion
cyclotron resonance occurs in a flux tube, as can be seen from
|Figure 1: Gyro-frequency of alpha particles in a flux tube, the curve for (dashed line) is taken from model 1 in Fig. 2.|
Figure 1 shows the gyro-frequency of alpha particles and as a function of height h, which is defined by h=r-r0, and where r0 is the heliospheric distance at the lower boundary (assumed to be ). Below 105 km, flow velocities are much smaller than the Alfvén speed and is nearly identical to . As a result, , is roughly the frequency at which cyclotron resonance transfers the wave energy to alpha particles. However, can be Doppler shifted and become larger than at larger distances. If , no wave heating will occur but the waves can still transfer momentum to the plasma (via wave pressure gradient). The abundance of alpha particles is chosen at 6%. This abundance is slightly higher than its in situ measured value of about 5% in fast solar wind streams, but lower than the helium abundance at the photosphere. Previous solar wind models showed that this abundance at the coronal base can often lead to observed helium abundance in interplanetary space (Li et al. 1997; Hu & Habbal 1999).
|Figure 2: A coronal heating and fast solar wind model in which the plasma is heated by ion cyclotron waves, Model 1. a) Speed of protons (thin solid line) and alpha particles (thick solid line). b) Temperature of electrons (dotted line), protons (thin solid line) and alpha particles (thick solid line). The dashed-dotted line is the proton effective temperature calculated from model 1. The narrow and wide vertical bars represent the observed proton effective temperatures from UVCS observations reported by Esser et al. (1999) and Kohl et al. (1998), respectively. c) Density of protons (thin solid line) and alpha particles (thick solid line). d) Heating rates of protons (thin solid line) and alpha particles (thick solid line). Here Hz.|
|Figure 3: An enlarged version of Fig. 2 in the region just above the lower boundary.|
The model calculation predicts a rather low proton temperature
between 2 and 4 solar radii. The local minimum
K, only slightly higher
than the electron temperature. However, the Alfvén wave
amplitude is strong (as will be shown in Fig. 6) in the region.
Hence, the effective proton temperature
Note the wave amplitude at the base of the flux tube is 21 km s-1, much smaller than the value used in Li (2002). The wave amplitude is at the lower end observed by SUMER (Tu et al. 1998). However SUMER observations were done above the solar limb. The coronal funnels considered in this paper are deeply rooted in the lower transition region, a smaller wave amplitude can be justified due to the high density. By combining the solar wind acceleration and the transition region in a full treatment, the heat deposited below the subsonic point can increase mass fluxes substantially. To produce a particle flux constrained by in situ measurements, a smaller wave amplitude is necessary. This demonstrates that the coronal heating along open field lines and the solar wind acceleration must be treated as a single problem.
To show the transition region and the initial heating process more clearly, Fig. 3 shows a small region above the lower boundary of the solution in Fig. 2. When temperatures are still low, alpha particles and protons have virtually the same speed (h<30 km). At km, alpha particles are not only faster, but much hotter than protons as well. The proton velocity at the lower boundary is km s-1. The alpha particles reach a local maximum velocity of 28 km s-1 at a height of 600 km. This effect is totally due to the introduction of ion cyclotron resonance. Below this point, the wave frequency at the high end of the spectrum is still smaller than the frequency where alpha particles can resonate with the waves. Hence no wave dissipation occurs. The most interesting feature of Fig. 3 is that at the top of the transition region or at the bottom of the corona, alpha particles and protons (and electrons) can be significantly away from thermal equilibrium. In the region km, alpha particles are slower than protons. This is because the higher temperature gradient of alpha particles (compared with protons) actually produces a negative pressure gradient force for alpha particles (see Fig. 5), and alpha particles are slower than protons. Since is negligible at h<105 km, protons and electrons are heated by Coulomb coupling with protons and alpha particles, and by the electron heat flux. The evolution of speed and temperatures in Fig. 3 is similar to the results in Hackenberg et al. (2000). The loss of thermal equilibrium between protons and electrons at the top of transition region or the bottom of the corona has also been found by Marsch & Tu (1997a) and Vocks & Marsch (2002). Here it is shown that protons and alpha particles can also lose thermal equilibrium in the same region. The transition region of model 1 is very narrow: the jump of the electron temperature from to K takes place in 390 km.
A few comments can be made about the velocity profile of alpha particles in Fig. 3a. Transition region spectral line observations often find various velocities of ions. These velocities can vary from a few to about twenty kilometers per second. For instance, Hassler et al. (1999) found that in network boundaries an average outflow velocity is several kilometers per second when K, while Wilhelm et al. (2000) estimated that the outflow velocity of Ne7+ is about 14 km s-1 along open magnetic field lines from their SUMER data. The velocity of alpha particles in Fig. 3a is in the range 3 km s km s-1 at different heights. Even though no other minor ions have been included in this investigation, however, if those ions behave similar to the alpha particles, their velocity is expected to vary in a similar manner. Hence, if spectral lines are observed at different heights, observed outflow velocities may vary accordingly.
The fact that minor ions can flow much faster than the bulk plasma flow in
the upper transition region may have implications on the analysis of
transition region spectral line observations.
Traditional atomic physics calculations on which the diagnostics
of both temperatures and densities depend
often assume equal ion and electron temperatures. From Doppler
shift measurements of transition region spectral lines, Dere et al. (1989) found that the mass flux in the region is more than
enough to account for in situ measured solar wind mass flux. If
minor ions are indeed much faster than protons and electrons in
the transition region, one cannot safely determine the bulk
outflow velocity of plasmas in the region by using minor ion
|Figure 4: Ratio between alpha particle heating term due to waves and collisional heating terms in alpha particle energy energy equation (Hu & Habbal 1999).|
Needless to say, in the region where no direct external heating is
added, the strong heating of the three species is entirely due to
the heat fluxes of the species and Coulomb coupling between them.
The Coulomb collisional coupling between ions is much stronger
than the coupling between ions and electrons. Hence protons and
alpha particles have identical temperatures. When the electron
temperature gradient becomes very large, electrons become hotter
than protons and alpha particles (h<10 km, Fig. 3b). This can be
easily understood since the classical electron thermal
conductivity is more than one magnitude larger than the classical
ion thermal conductivities. Note that the strong Coulomb coupling
makes the temperature difference less than 20% (at h<10 km).
This seems to make our assumption that species have the same
K questionable. However,
the electron heat conductivity is proportional to
the heat flux is relatively small at low temperatures. And
most importantly, Coulomb collisions become extremely strong at
low temperatures. This small discrepancy does not influence the
model solution in a significant way. To avoid this discrepancy,
future models need to include the chromosphere and the detailed
ionization processes of hydrogen and helium, which is beyond the
scope of this paper.
|Figure 5: Model 1: thermal pressure of electrons (dotted line), protons (thin solid line), and alpha particles (dashed line). The thick solid line represents the total pressure .|
To show how the loss of thermal equilibrium between protons and
alpha particles can occur around h=1000 km in the model, R, the
and proton/alpha Coulomb coupling terms
in the energy equation of the alpha particles, is plotted in Fig. 4. Clearly only in the region where |R| is significantly larger
than 1 do alpha particles lose thermal equilibrium with protons.
Note the rate R is negative, suggesting that alpha particles
lose thermal energy to protons due to collisions. At h<105 km,
roughly speaking, only alpha particles are directly heated by
waves. Without any other external heat sources, protons are only
heated via Coulomb coupling with alpha particles. Above
h=300 000 km (or about 1.4
), this ratio is much larger
than 1, suggesting that Coulomb coupling is not important.
|Figure 6: Model 1: a) the Alfvén speed (solid line) and the wave amplitude (dashed line); b) ratio ; c) Proton beta value.|
Figure 5 shows the species thermal pressures from the model in Fig. 2. The total thermal pressure is nearly a constant between K and K. However, the individual species pressures are not as constant. Ion cyclotron resonance heating yields a large alpha thermal pressure bump at km.
The Alfvén speed
and Alfvén wave amplitude
are shown in Fig. 6a. Here
is the total plasma mass
The Alfvén wave amplitude reaches 35 km s-1 at
in Fig. 6b. To compare Fig. 6b with Fig. 2d,
the proton heating rate
becomes larger than at about 2.7
where the ratio
Note, if warm plasma dispersion relation is used, protons may get
more heating faster (Li & Habbal 1999). Hollweg & Isenberg
(2002) also discussed this issue. The "switch" from primarily
heating alpha particle to primarily heating protons depends on the
abundance of alpha particles, the plasma beta value and the
detailed dissipation range spectrum.
The proton beta value, the ratio between the proton thermal pressure and
|Figure 7: Contributions to the energy flux density, multiplied by for Model 1: Alfvén wave energy flux (dashed); total enthalpy (dotted); total kinetic energy flux (thin solid); total heat flux (dot-dashed); gravity (dash-triple-dotted); radiation (thick dashed). The thick solid line is the total energy flux density.|
Figure 7 shows the various contributions to the energy flux, scaled
to the inner boundary. The total energy flux is a sum of the
Alfvén wave energy flux; kinetic energy flux of protons and
alpha particles; gravitational energy flux; total enthalpy flux;
total heat flux of electrons, protons and alpha particles; and the
radiation loss (an integration from 1
to a given
distance). It is obvious that the radiation loss is a small
quantity in the energy budget. The heat flux is only important in
the transition region. The total energy flux is a straight
horizontal line, indicating an excellent energy conservation. In
fact, the maximum relative error in the total energy flux is
0.16%, which is too small to be seen in Fig. 7.
From Figs. 3 and 4, it is seen that alpha particles lose thermal
equilibrium with protons around h=1000 km when the heat gain that
the alpha particles obtained from waves is much more than the heat
loss through the Coulomb coupling with protons. In this region,
decreases most rapidly with height due to the fast
flux tube expansion. If the wave heating begins at an even higher
location in the coronal funnel, for instance, the wave heating
starts after the first flux tube expansion almost ends,
the heating rate
will be smaller since
now is smaller and there is less wave energy at the high
frequency end. One of course expects that the loss of an
equilibrium between alpha particles and protons is going to be
weaker. This is exactly the case when the wave frequency at the
high end of the spectrum
is 17 000 Hz. The overall
feature of this case is similar to those in Fig. 2 and will not be
shown here. The region just above the lower boundary is shown in
Fig. 8. The ion cyclotron heating is not introduced until at a
height 2070 km above the lower boundary (Fig. 8d). Even though the
wave heating still leads to different alpha particle and proton
km (Fig. 8d), the difference is far
less significant than in Fig. 3b. There is virtually no velocity
difference in this region. Plasma species are roughly in thermal
equilibrium at h<105 km (Figs. 8a and 8b). The temperatures of
Model 2 in the transition region below 5000 km are much lower than
in model 1 (see Figs. 3b and 7b). However, Model 2 has a slightly
lower maximum electron temperature (
K in model 1). Although electrons are not
directly heated, these maximum
observations in coronal holes (Wilhelm et al. 1998). Overall, the
difference between Model 1 and Model 2 is marginal except in the
|Figure 9: Proton particle flux density (solid line) at the lower computation boundary, proton speed (dashed line) at , and the transition region height (see text), as a function of .|
Figure 9 shows the proton particle flux density at 1
proton velocity at 52.5
and the height of the
transition region (defined as the distance over which the electron
temperature increases from
K), as a
function of ,
other parameters are the same as used in
Model 1. Once the wave amplitude at the lower boundary and the
magnetic field geometry are chosen, the total energy flux entering
the computational domain is fixed (the only energy losses are heat
flux and radiation loss). However, if
is smaller than
10 000 Hz, a steady state solution can not be obtained. The reason
is that as
is small, the electron temperature gradient
becomes smaller. The transition region becomes too wide so the
electron heat flux is unable to balance the local radiation
loss. Hence a steady state hot corona and wind cannot be found.
(Note, however, if the chromosphere is also included in the model
calculation, the position of transition region may change as the
position of the energy deposit changes. But this is beyond the
scope of the present study.) Of course, if the wave power at high
frequency is enhanced, as suggested by Tu & Marsch (1997), the
result may be different. Roughly speaking, the proton mass flux is
almost constant (the change is less than 10%), although the
plasma parameters at the coronal base can be very different as
shown in Figs. 3 and 8. The small variation reflects the different
heat flux losses to the chromosphere: about 8.6% of the total
Alfvén wave energy flux (deposited between 1 and 52.5 )
flows back to the chromosphere, while only less than 2.8% of
the Alfvén wave energy flux flows back to the chromosphere in
Model 2. Hence a smaller
means more energy flux will
be available for the solar wind acceleration, and a smaller
particle flux and a faster wind may result (Leer & Holzer 1980).
Note, the height of the transition region changes more
Finally we briefly present Model 3 with a much weaker magnetic flux concentration in the supergranule convection zone boundaries. In Model 3, and the magnetic field at the lower boundary is B0=36.1 G, so the magnetic field in interplanetary space is unchanged. Again the overall large scale velocity and temperature profiles (Model 3) are very close to Models 1 and 2 except a small region above the lower boundary, so only the region above the lower boundary is shown in Fig. 10. Because at the base of the corona, is much smaller than in model 1, and there is less wave energy available for alpha particles below 1.1 solar radii (h<105 km). This leads to a much colder corona below 1.1 solar radii and a less steep transition region (The maximum electron temperature is K.). Once again, species are in thermal equilibrium in the transition region (h<105 km). One noticeable feature of this model is that the flow speed at K level is quite small (less than 5 km s-1). The transition region is wide: the electron temperature increases from K to K in a distance of 3300 km.
The model also produces reasonable proton particle flux: . This model actually produces a faster wind than in Fig. 2. This can be understood since less wave energy is released just above the transition region, and the electron heat flux down to the region K is also smaller. Hence more wave energy flux becomes available for the solar wind acceleration. There is less energy at high frequencies suggesting less heat will be deposited in the subsonic region, leading to a faster but diluter wind (Leer & Holzer 1980). It is found that if is further decreased, with the same Alfvén wave amplitude, a steady state solution of a corona and wind can not be obtained.
Plasma heating in coronal funnels along open magnetic field lines by parallel propagating ion cyclotron waves is investigated in a three-fluid approach. It is shown that the Coulomb coupling between alpha particles and protons/electrons is strong enough to produce a sharp transition region and a hot corona by heating alpha particles alone. As shown by previous studies (Marsch & Tu 1997a,b; Tu & Marsch 1997; Hackeberg et al. 2000; Li 2002), a spectrum of ion cyclotron waves of wave spectrum index -1 is able to produce a transition region, a hot corona and to drive a fast solar wind if the wave spectrum extends to the ion cyclotron frequency range and the magnetic flux concentration in the super-granule convection zone boundaries is sufficiently strong. It is shown for the first time that by heating alpha particles alone, a sharp transition region, a hot corona and basically a fast solar wind can be produced. Hence, heavy ions may play a crucial role in the corona heating and the solar wind acceleration (Li et al. 1997).
It is found that the highest frequency of the waves, or the location of the heating determines whether a thermal equilibrium can be established between minor ions and protons. As the maximum frequency is increased, alpha particles and protons end up farther away from thermal equilibrium. As the maximum frequency becomes smaller, alpha particles and protons are more likely to be in thermal equilibrium. However, if the maximum frequency is below a critical value, ion cyclotron resonance occurs far away from the transition region and a hot corona cannot be produced. The fact that heating alpha particles alone can produce a hot corona and wind has interesting implications on the theory of coronal heating and expansion. This study suggests that a corona can be created by solely heating minor ions. Of course, the coronal heating physics may be far more complicated than this study shows. Since there are many ions heavier than the alpha particles in the transition region and their gyro-frequencies are smaller than the gyro-frequency of alpha particles, those heavy ions will be heated first if parallel propagating ion cyclotron waves indeed are responsible for the coronal heating. Since the abundances of those minor ions are so small in the transition region and corona, it remains an open question whether the Coulomb coupling is able to efficiently transfer the energy from these minor ions to the major species (protons and electrons) (see Cranmer 2000; Hollweg & Isenberg 2002).
For simplicity, fixed velocity distributions and spectra in the dissipation range are assumed. Tu & Marsch (2001) have pointed out that due to cyclotron resonance, it is impossible for the spectrum to maintain a rigid shape. Ion cyclotron resonance may make the wave spectrum to have a cutoff at higher wave-numbers (Li et al. 2001). Physically, a sufficiently high spectral index shall mimic this effect. Since we have already taken a quite large spectral index -5 in the dissipation range, a larger spectral index is not going to change model results presented in this study. The energy distribution between protons and alpha particles is mainly determined by the dispersion relation. The model results are not sensitive to the dissipation range spectral index in this study. However the assumption of a fixed spectrum and velocity distribution may give unrealistic energy distributions between species. To tackle the wave heating problem in the corona, self-consistent spectral evolution including wave dispersion and kinetic treatment of particles are needed. The issue is currently under study by many authors (Tam & Chang 1999; Galinsky & Shevchenko 2000; Isenberg et al. 2000, 2001; Vocks & Marsch 2001, 2002).
This study may have implications on the solar wind charge state investigations. By assuming plasma species in thermal equilibrium, coronal electron temperature can be calculated using charge fraction information from in situ charge state measurements. However, those calculations found that the coronal electron temperature above coronal holes is much higher than the electron temperature observed by SUMER (see, Ko et al. 1997; Esser et al. 1998). Since this investigation suggests that ions in the transition region may be already out of thermal equilibrium and have different velocities, the contradiction between SUMER and charge state calculations may be due to the loss of equilibrium in the transition region and corona. Recently, Esser & Edgar (2001) found that the velocity difference can significantly reduce the discrepancy between observed in situ ion fractions and low coronal temperatures.
We have neglected the temperature anisotropy of species in this study and have adopted a collision dominated standard transport equations and modified classical thermal conduction for the three-fluid modelling. One has to recognize that it is still an open question to what degree the treatment is valid. Certainly Coulomb collisions are playing a major role in the transition region. As the plasma density deceases rapidly with an increasing height, collisions become less important in the inner corona. During the transition from collision dominated regime to collisionless regime, a refined description of species heat flux is needed. The heat flux of ions may also play a role in shaping the temperature anisotropy of the ions (Olsen et al. 1998; Olsen & Leer 1999; Li 1999; Lie-Svendsen et al. 2001). A complete treatment to include species temperature anisotropies using high moment transport equations is under way.
This work is supported by PPARC rolling grant to University of Wales at Aberystwyth. Part of the work was supported by the grant NAG5-10873 to Smithsonian Astrophysical Observatory (SAO).