The rise of twisted magnetic flux tubes

Sunspots are caused by the eruption of magnetic flux tubes through the solar photosphere: current theories of the internal magnetic field of the Sun suggest that such tubes must rise relatively unscathed from the base of the convection zone. In order to understand how the structure of the magnetic field within a buoyant flux tube affects its stability as it rises, we have considered the quasi-two-dimensional rise of isolated magnetic flux tubes through an adiabatically stratified atmosphere. The magnetic field is initially helical; we have investigated a range of initial field configurations, varying the distribution and strength of the twist of the field.     

Radiative heating and the buoyant rise of magnetic flux tubes in the solar interior

We study the effect of radiative heating on the evolution of thin magnetic flux tubes in the solar interior and on the eruption of magnetic flux loops to the surface. Magnetic flux tubes experience radiative heating because (1) the mean temperature gradient in the lower convection zone and the overshoot region deviates substantially from that of radiative equilibrium, and hence there is a non-zero divergence of radiative heat flux; and (2) the magnetic pressure of the flux tube causes a small change of the thermodynamic properties within the tube relative to the surrounding field-free fluid, resulting in an additional divergence of radiative heat flux. Our calculations show that the former constitutes the dominant source of radiative heating experienced by the flux tube.In the overshoot region, the radiative heating is found to cause a quasi-static rising of the toroidal flux tubes with an upward drift velocity 10^-3|| cm s^-1, where _e – _ad < 0 describes the subadiabaticity in the overshoot layer. The upward drift velocity does not depend sensitively on the field strength of the flux tubes. Thus in order to store toroidal flux tubes in the overshoot region for a period comparable to the length of the solar cycle, the magnitude of the subadiabaticity (< 0) in the overshoot region must be as large as 3 × 10^–4. We discuss the possibilities for increasing the magnitude of and for reducing the rate of radiative heating of the flux tubes in the overshoot region.Using numerical simulations we study the formation of -shaped emerging loops from toroidal flux tubes in the overshoot region as a result of radiative heating. The initial toroidal tube is assumed to be non-uniform in its thermodynamic properties along the tube and lies at varying depths beneath the base of the convection zone. The tube is initially in a state of neutral buoyancy with the internal density of the tube plasma equal to the local external density. We find from our numerical simulations that such a toroidal tube rises quasi-statically due to radiative heating. The top portion of the nonuniform tube first enters the convection zone and may be brought to an unstable configuration which eventually leads to the eruption of an anchored flux loop to the surface. Assuming reasonable initial parameters, our numerical calculations yield fairly short rise times (2–4 months) for the development of the emerging flux loops. This suggests that radiative heating is an effective way of causing the eruption of magnetic flux loops, leading to the formation of active regions at the surface.The National Solar Observatory is one of the National Optical Astronomy Observatories by the Association of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation.     

Topological stability of finite-length magnetic flux tubes

It has been suggested that the activity of cosmical magnetic fields is a consequence of a general topological nonequilibrium in the neighbourhood of magnetostatic equilibria. Evidence for this suggestion can be obtained from the Kolmogorov-Arnold-Moser theorem of classical mechanics, applied to the magnetic field line flow as a Hamiltonian system. A finite-length magnetic flux tube, however, always possesses two independent sets of flux surfaces - or, equivalently, the corresponding Hamiltonian system two independent first integrals - and is topologically stable if in the volume occupied by the tube there are no singular (null) points of the magnetic field and the normal field component does not change its sign on the end faces of the tube. Therefore, the concept of nonequilibrium due to flux surface destruction is not applicable to solar atmospheric loops with each end situated in the interior of one polarity of the photospheric normal field component. Further, it seems unlikely that the tearing-mode mechanism can play a role in such loops.     

Thermodynamical properties of unresolved magnetic flux tubes

We propose a diagnostic method, based on the observation of circular polarization signals in line pairs, to derive the thermodynamical properties of unresolved magnetic elements in the solar atmosphere. The concept of response function for the ratio of circular polarization signals in two lines is introduced and its main properties are analyzed. Some detailed calculations for suitably selected line pairs are presented.     

Emergence of sheared magnetic flux tubes in an active region observed with the SVST and TRACE

The active region NOAA AR 8331 was a target of an international ground-based observational campaign in the Canaries and coordinated with space instruments (TRACE and Yohkoh). We focus our study on observations obtained with the SVST at LaPalma, and with TRACE. On 10 September 1998, arch-filament systems were observed with high spatial and temporal resolution, from the lower to the upper atmosphere of the Sun, during five hours. Flux tubes emerged with increasing shear, which apparently led to energy release and heating in the overlying corona. A model for filament formation by the emergence of U-shaped loops from the subphotosphere, as proposed by Rust and Kumar (1994), is supported by the present observations. The coronal response to these events is visualized by rising, medium-scale loop brightenings. The low-lying X-ray loops show short-lived, bright knots which are thought to result from interaction between different loop systems.     

Time-dependent interaction of granules with magnetic flux tubes

The time-dependent interaction of the granulation velocity field with a magnetic flux tube is investigated here. It is seen that when a magnetic field line is displaced normal to itself so as to simulate thebuffeting action of granules, a flow of gas is initiated along the field. By choosing a lateral velocity field which is consistent with observations of granules, it is found that the resulting gas motion is a downward flow with a velocity compatible with the observed downflow in isolated photospheric flux tubes. It is therefore proposed that the observed photospheric downflow is a manifestation of the interaction of granules with flux tubes.     

Active region magnetic fields

The tilt angles of sunspot groups are defined, using the Mount Wilson data set. It is shown that groups with tilt angles greater than or less than the average value (&amp;#x2248; 5 deg) show different latitude dependences. This effect is also seen in synoptic magnetic field data defining plages. The fraction of the total sunspot group area that is found in the leading spots is discussed as a parameter that can be useful in studying the dynamics of sunspot groups. This parameter is larger for low tilt angles, and small for extreme tilt angles in either direction. The daily variations of sunspot group tilt angles are discussed. The result that sunspot tilt angles tend to rotate toward the average value is reviewed. It is suggested that at some depth, perhaps 50 Mm, there is a flow relative to the surface that results from a rotation rate faster than the surface rate by about 60 m/sec and a meridional drift that is slower than the surface rate by about 5 m/sec. This results in a slanted relative flow at that depth that is in the direction of the average tilt angle and may be responsible for the tendency for sunspot groups (and plages) to rotate their magnetic axes in the direction of the average tilt angle.     

Twisting together magnetic flux tubes under helical symmetry

This paper investigates nonlinear interaction ofmagnetic flux tubes, as a possible cause of current-sheet formation. We focus attention on Gold–Hoyle tubes because of their simple analytic form, using a frictional magnetic relaxation code to find eventual equilibria. We assume that all fields possess helical symmetry so the problem becomes essentially two-dimensional and high resolution can be more easily achieved. When the tubes are not twisted together current sheets form only if the plasma pressure is zero. Twisting the tubes when the plasma pressure is small but finite results in curved current sheets and line currents. Current singularities are identified by performing calculations at increasing grid resolutions, and observing a regular increase in the maximum current.     

Leaky and non-leaky oscillations in magnetic flux tubes

An extensive analysis, both analytic and numerical, of waves in flux tubes imbedded in (possibly) magnetic surroundings is given. It is shown that any wave confined to the tube and its neighbourhood can be put into one of seven categories. Simple criteria for deciding the existence of each type in any particular case are derived. Many other (leaky) modes are found which excite waves in the external medium and thereby lose energy to the surroundings. A number of asymptotic analyses allow much information to be gained about these without the need for numerical solution of the complicated equations involved. Three particular cases, pertaining to photospheric flux tubes, H fibrils, and coronal loops, are considered in detail.     

The size dependence of contrasts and numbers of small magnetic flux tubes in an active region

Intensity contrasts and number densities of bright points, knots and pores ranging in size between 0.15 and 4 are studied using high resolution pictures in Mg b₁ of a young active region. On the average, the contrast in the wing of the line increases very strongly with decreasing size, while the continuum contrast increases more slowly. The ratio of contrast in the line to contrast in the continuum increases rapidly with decreasing size. The possibility is explored of using this contrast ratio as an indicator of size. The distribution of the contrast ratio in a part of the active region is used in this way to derive a size distribution of facular points. The resulting distribution has a limited accuracy, but is free from systematic distortion due to selection effects. Validity checks on the method are presented. We measure the size distribution of the pores in the same area, and combine the result with that for the facular points. The combined distribution shows that the surface area covered by magnetic elements with diameter has a maximum near = 0.8. It increases roughly proportional to for < 0.3 and falls off as for > 1.5. It is inferred that elements with 0.5 < < 1.6, which show no conspicuous contrast in the line wing or in the continuum, occupy as much area as the pores, and twice as much as the facular points. We suggest that the changing appearance of a facular area with increasing height of formation reflects at least as much the increasing weight of the small elements in the contrast as a real change in intrinsic properties (such as the diameter) of individual elements. A spatial resolution better than 0.1 may be needed to resolve the individual elements in plages and the chromospheric network.The observed variation of continuum contrast of facular points with size agrees with predictions based on magnetostatic flux tube models if a field strength of about 2000 G is assumed.     

Exact static equilibrium of vertically oriented magnetic flux tubes

A method is prescribed for generating exact solutions of magnetostatic equilibrium describing a cylindrically symmetric magnetic flux tube oriented vertically in a stratified medium. Given the geometric shape of the field lines, compact formulae are presented for the direct calculation of all the possible distributions of pressure, density, temperature and magnetic field strength compatible with these field lines under the condition of static equlibrium. The plasma satisfies the ideal gas law and gravity is uniform in space. A particular solution is obtained by this method for a medium sized sunspot whose magnetic field obeys the similarity law of Schlüter and Temesváry (1958). With this solution, it is possible for the first time to illustrate explicitly the confinement of the magnetic field of the cool sunspot by the hotter external plasma in an exact relationship involving both magnetic pressure and field tension as well as the support of the weight of the plasma by pressure gradients. It is found that the cool region of the sunspot is not likely to extend much more than a few density scale heights below the photosphere. The sunspot field approaches being potential in the neighbourhood of the photosphere so that the Lorentz force exerting on the photosphere is less than what the magnetic pressure would suggest. This accounts for how the sunspot field can be confined in the photosphere where its magnetic pressure is often observed to even exceed the normal photospheric pressure. The energy mechanism operating in the sunspot and the question of mechanical stability are not treated in this paper.Normally at Lau Kuei Huat (Singapore) Private Limited, 55 Shipyard Road, Singapore 22, Singapore.     

Twisted magnetic flux tubes: Effect of small twist

The shape of a magnetic flux tube is investigated when photospheric motion causes small twist at the magnetic footpoints. Using a Fourier-Bessel series expansion, the previous results of Zweibel and Boozer (1985) and Steinolfson and Tajima (1987), when the twist is small, are substantiated. A twisting motion that is restricted to a finite region is investigated. Inside the twisted region, the tube contracts, but in the outer region the field remains straight, except for a slight expansion at the outside of the loop near the footpoints. The amount of twist depends on the radial position and can in fact be larger in the contracted region with the twist decreasing as the tube expands. An axial boundary-layer region is present, as noted by the above authors, through which the field adjusts to the line-tied magnetic footpoint positions. An analysis of the boundary layer shows that the thickness remains constant as the loop-length is increased with the result that the main part of the loop has constant cross-sectional area and can be described by cylindrically-symmetric fields. This new 1-D model predicts the main behaviour of the loop without the need to solve the more complicated 2-D problem directly. It is speculated that the boundary layers will remain even when the twist becomes large and a simple example is presented. A detailed parametric study of different twist profiles shows how the central part of the loop responds.Using the result that the majority of the loop can be described by a constant cross-sectional area, a model for a toroidal loop is presented that models coronal loops in a more realistic manner. The main result from this section is that the coronal loops can only remain in equilibrium if they are confined by an external magnetic field (possibly potential in nature) and not by a gas pressure unless additional physical effects are included.     

Resonant behaviour of MHD waves on magnetic flux tubes

The resonances that appear in the linear compressible MHD formulation of waves are studied for equilibrium states with flow. The conservation laws and the jump conditions across the resonance point are determined for 1D cylindrical plasmas. For equilibrium states with straight magnetic field lines and flow along the field lines the conserved quantity is the Eulerian perturbation of total pressure. Curvature of the magnetic field lines and/or velocity field lines leads to more complicated conservation laws. Rewritten in terms of the displacement components in the magnetic surfaces parallel and perpendicular to the magnetic field lines, the conservation laws simply state that the waves are dominated by the parallel motions for the modified slow resonance and by the perpendicular motions for the modified Alfvén resonance.The conservation laws and the jump conditions are then used for studying surface waves in cylindrical plasmas. These waves are characterized by resonances and have complex eigenfrequencies when the classic true discontinuity is replaced by a nonuniform layer. A thin non-uniform layer is considered here in an attempt to obtain analytical results. An important result related to earlier work by Hollweg et al. (1990) for incompressible planar plasmas is found for equilibrium states with straight magnetic field lines and straight velocity field lines. For these equilibrium states the incompressible and compressible surface waves have the same frequencies at least in the long wavelength limit and there is an exact correspondence with the planar case. As a consequence, the conclusions formulated by Hollweg et al. still hold for the straight cylindrical case. The effects of curvature are subsequently considered.     

Resonant behaviour of MHD waves on magnetic flux tubes

Resonant absorption of MHD waves on a nonuniform flux tube is investigated as a driven problem for a 1D cylindrical equilibrium. The variation of the fractional absorption is studied as a function of the frequency and its relation to the eigenvalue problem of the MHD radiating eigenmodes of the nonuniform flux tube is established. The optimal frequencies producing maximal fractional absorption are determined and the condition for total absorption is obtained. This condition defines an impedance matching and is fulfilled for an equilibrium that is fine tuned with respect to the incoming wave. The variation of the spatial wave solutions with respect to the frequency is explained as due to the variation of the real and imaginary parts of the dispersion relation of the MHD radiating eigenmodes with respect to the real driving frequency.     

Resonant interactions of modes in coronal magnetic flux tubes

Nonlinear resonant interactions of different kinds of fast magnetosonic (FMS) waves trapped in the inhomogeneity of a low- plasma density, stretched along a magnetic field (as, for example, in coronal loops) are investigated. A set of equations describing the amplitudes of interactive modes is derived for an arbitrary density profile. The quantitative characteristics of such interactions are found. The decay instability of the wave with highest frequency is possible in the system. If amplitudes of interactive modes have close values, the long-period temporal and spatial oscillations are in the system.For a quantitative illustration, the parabolic approximation of the transverse density profile has been chosen. Dispersion relations of FMS waves trapped in a low- plasma slab with a parabolic transverse density profile are found. The transverse structure of the waves in this case can be expressed through Hermitian polynomials. The interaction of kink and sausage waves is investigated. The sausage wave, with a sufficiently large amplitude, may be unstable with respect to the decay into two kink waves, in particular. The spatial scale of a standing wave structure and the time spectrum of radiation are formed due to the nonlinear interactions of loop modes which contain information about the parameters of the plasma slab.     

Stability of magnetic flux tubes in close binary stars

Photometric and Doppler imaging observations of active binaries indicate the existence of starspots at preferred longitudes (position angles with respect to the companion star). We investigate the stability of magnetic flux tubes in the convection zone of close, fast‐rotating binary stars and explore whether the observed preferred longitudes could be caused by tidal forces and the deformation of the active star. We assume a synchronized binary system with spin axes perpendicular to the orbital plane and a rotation period of a few days. The tidal force and the deviation from spherical structure are considered in lowest‐order perturbation theory. The magnetic field is in the form of toroidal magnetic flux rings, which are stored in mechanical equilibrium within the stably stratified overshoot region beneath the convection zone until the field has grown sufficiently strong for the undulatory instability to initiate the formation of rising loops. Frequencies and geometry of stable as well as growth rates of unstable eigenmodes are determined by linear stability analysis. Particular consideration is given to the question whether the effects of tidal forces and perturbations of the stellar structure can force a rising flux loop to enter the convection zone at specific longitudes.     

Resonant behaviour of MHD waves on magnetic flux tubes

The absorption of solar five-min oscillations by sunspots is interpreted as the resonant absorption of sound waves by a magnetic cylinder. The absorption coefficient is calculated both analytically under certain simplifying assumptions, and numerically under more general conditions. The observed magnitude of the absorption coefficient, which is up to 0.5 or even more, can be explained for suitable ranges of parameters. Limitations in the present model are also discussed.