A Nonlinear Force-Free Magnetic Field Approximation Suitable for Fast Forward-Fitting to Coronal Loops. I. Theory

We derive an analytical approximation of nonlinear force-free magnetic field solutions (NLFFF) that can efficiently be used for fast forward-fitting to solar magnetic data, constrained either by observed line-of-sight magnetograms and stereoscopically triangulated coronal loops, or by 3D vector-magnetograph data. The derived NLFFF solutions provide the magnetic field components B _x (x), B _y (x), B _z (x), the force-free parameter α(x), the electric current density j(x), and are accurate to second-order (of the nonlinear force-free α-parameter). The explicit expressions of a force-free field can easily be applied to modeling or forward-fitting of many coronal phenomena.     

Alfvén waves in the solar atmosphere

We examine the propagation of Alfvén waves in the solar atmosphere. The principal theoretical virtues of this work are: (i) The full wave equation is solved without recourse to the small-wavelength eikonal approximation (ii) The background solar atmosphere is realistic, consisting of an HSRA/VAL representation of the photosphere and chromosphere, a 200 km thick transition region, a model for the upper transition region below a coronal hole (provided by R. Munro), and the Munro-Jackson model of a polar coronal hole. The principal results are:

1. If the wave source is taken to be near the top of the convection zone, where n _H = 5.2 × 10¹⁶ cm^?3, and if B _⊙ = 10.5 G, then the wave Poynting flux exhibits a series of strong resonant peaks at periods downwards from 1.6 hr. The resonant frequencies are in the ratios of the zeroes of J ₀, but depend on B _⊙, and on the density and scale height at the wave source. The longest period peaks may be the most important, because they are nearest to the supergranular periods and to the observed periods near 1 AU, and because they are the broadest in frequency.
2. The Poynting flux in the resonant peaks can be large enough, i.e. P _⊙ ≈ 10⁴–10⁵ erg cm^?2s^?1, to strongly affect the solar wind.
3. ¦δv¦ and ¦δB¦ also display resonant peaks.
4. In the chromosphere and low corona, ¦δv ≈ 7–25 kms^?1 and ¦δB¦ ≈0.3–1.0 G if P _⊙≈10⁴-10⁵ erg cm^?2s^?1.
5. The dependences of ¦δv¦ and ¦δB¦ on height are reduced by finite wavelength effects, except near the wave source where they are enhanced.
6. Near the base, ¦δB¦ ≈ 350–1200 G if P _⊙ ～- 10⁴–10⁵. This means that nonlinear effects may be important, and that some density and vertical velocity fluctuations may be associated with the Alfvén waves.
7. Below the low corona most wave energy is kinetic, except near the base where it becomes mostly magnetic at the resonances.
8. ?₀ < δv ² > v _A or < δB ² > v _A/4π are not good estimators of the energy flux.
9. The Alfvén wave pressure tensor will be important in the transition region only if the magnetic field diverges rapidly. But the Alfvén wave pressure can be important in the coronal hole.

     

Energy changes of cosmic rays in the interplanetary region

A clarification and discussion of the energy changes experienced by cosmic rays in the interplanetary region is presented. It is shown that the mean time rate of change of momentum of cosmic rays reckoned for a fixed volume in a reference frame fixed in the solar system is 〈p〉 =p V·G/3 (p=momentum,V is the solar wind velocity andG=cosmic-ray density gradient). This result is obtained in three ways:

1. by a rearrangement and reinterpretation of the cosmic-ray continuity equation;
2. by using a scattering analysis based on that of Gleeson and Axford (1967);
3. by using a special scattering model in which cosmic-rays are trapped in ‘magnetic boxes’ moving with the solar wind.

The third method also gives the rate of change of momentum of particles within a moving ‘magnetic box’ as 〈p〉_ad = ?p ?·V/3, which is the adiabatic deceleration rate of Parker (1965). We conclude that ‘turnaround’ energy change effects previously considered separately are already included in the equation of transport for cosmic rays.     

On the Origin of Pulsations of Sub-THz Emission from Solar Flares

We propose a model to explain fast pulsations in sub-THz emission from solar flares. The model is based on the approach of a flaring loop as an equivalent electric circuit and explains the pulse-repetition rate, the high-quality factor, Q≥10³, low modulation depth, pulse synchronism at different frequencies, and the dependence of the pulse-repetition rate on the emission flux, observed by Kaufmann et al. (Astrophys. J. 697, 420, 2009). We solved the nonlinear equation for electric current oscillations using a Van der Pol method and found the steady-state value for the amplitude of the current oscillations. Using the pulse rate variation during the flare on 4 November 2003, we found a decrease of the electric current from 1.7×10¹² A in the flare maximum to 4×10¹⁰ A just after the burst. Our model is consistent with the plasma mechanism of sub-THz emission suggested recently by Zaitsev, Stepanov, and Melnikov (Astron. Lett. 39, 650, 2013).     

The influence of fluctuation fields on the force-free field

In Paper I (Hu, 1982), we discussed the the influence of fluctuation fields on the force-free field for the case of conventional turbulence and demonstrated the general relationships. In the present paper, by using the approach of local expansion, the equation of average force-free field is obtained as (1+b)?×B ₀=(α#x002B;a)B ₀#x002B;a ⁽¹⁾×B ₀#x002B;K. The average coefficientsa,a ⁽¹⁾,b, andK show the influence of the fluctuation fields in small scale on the configurations of magnetic field in large scale. As the average magnetic field is no longer parallel to the average electric current, the average configurations of force-free fields are more general and complex than the usual ones. From the view point of physics, the energy and momentum of the turbulent structures should have influence on the equilibrium of the average fields. Several examples are discussed, and they show the basic features of the fluctuation fields and the influence of fluctuation fields on the average configurations of magnetic fields. The astrophysical environments are often in the turbulent state, the results of the present paper may be applied to the turbulent plasma where the magnetic field is strong.     

Drift theory of charged particles in electric and magnetic fields

In this paper we review the drift theory of charged particles in electric and magnetic fields. No new physical interpretations are added to this classical topic, but through an alternative, simplified derivation of the guiding centre velocity, several complexities are eliminated and possible misconceptions of the theory are clarified. It is shown that:

1. The curvature/gradient drift velocity in the magnetic field, averaged over a particle distribution function is to lowest order in the direction of?×B/B ², while the average particle velocity is in the direction ofB×? P withP the scalar particle pressure.
2. These drift directions are correct for first-order expansions of the particle distribution function, and only second-order or higher expansions change these directions.
3. The?×B/B ² drift, which is the standard gradient plus curvature drift, and which is usually considered as a ‘single particle’ drift, need not be ‘reconciled’ with theB×? P, or ‘macroscopic, collective’ drift, as is often asserted in the literature. They are in fact related per definition and we show how.
4. When viewed in fixed momentum intervals (p,p+dp), the so-called Compton-Getting factor enters into the electric field (E×B)/B ² drift term.
5. The results are independent of the scale length of variation ofE andB, in contrast to existing drift theory. We discuss the implications of this result for three important cases.

     

Dissipation of magnetic flux and magnetic energy in partially ionized gases

Macroscopic equations of motion are used to derive several forms of the generalized Ohm's law for partially ionized ternary gases in magnetic fields, and a conductivity σ is defined that is independent of the magnetic field. A flux theorem is derived using a velocityu _H that can be defined to be the velocity of magnetic field lines;u _H is only slightly different from the velocity of the electron component of the gas. It is shown that σ is the conductivity relevant to the decay of magnetic flux through any surface moving everywhere with velocityu _H . The rate of increase of the thermal energy density of the gas arising through collisions between particles of different species can be resolved into Joule heating at the ratej ²/σ, wherej is the current density, and heating associated with ambipolar drift. The latter, contrary to what has been claimed by some authors, is not necessarily fully compensated by a decrease in the energy of the electromagnetic field. In many applications such compensation does occur, but it may not in interstellar clouds where large amounts of gravitational energy can be made available by collapse, and then both heating and an increase in electromagnetic field energy may occur.     

A relation between Giorgilli-Galgani's and Deprit's algorithms

If \(T = \sum\nolimits_{i = 1}^\infty {\varepsilon ^i } T_i\) and \(W = \sum\nolimits_{n = 1}^\infty {n\varepsilon ^{n - 1} } W^{\left( n \right)}\) are respectively the generators of Giorgilli-Galgani's and Deprit's transformations, we show that the change of variables generated byT is the inverse of the one generated byW, ifT _i =W ⁽ⁱ⁾ for anyi. The method used is to show that the recurrence which defines the first algorithm can also be obtained with the second one.     

Some properties of finite energy constant-α force-free magnetic fields in a half-space

Some useful properties of a finite energy, constant-α, force-free magnetic field B _α occupying a half-space D are presented. In particular:

1. Fourier and Green representations of B _α are obtained and used to derive conditions for the existence and uniqueness of a B _α having a given normal component B _z on the boundary ?D.
2. The asymptotic behaviour of B _α at infinity as well as stability results against changes in the boundary condition on ?D and in the value of α are established.
3. The energy of B _α is shown to be smaller than the energy of the open field having the same B _z on ?D, thus confirming an earlier conjecture (Aly, 1984).
4. B _α is proved to not be a Taylor-Heyvaerts-Priest state, in spite of the fact that its relative helicity H is finite and that it is the only solution of the Lagrange-Euler equation associated with the problem of minimizing the energy among all the fields having the same value of H and the same B _z on ?D.

     

The two-dimensional structure of thin accretion disks

The two-dimensional structure of a thin accretion disk in the vicinity of a Schwarzschild black hole after passing a marginally stable orbit (r< 3r _g is discussed in terms of the Grad-Shafranov hydrodynamic equation. The accretion disk is shown to be sharply compressed as the sonic surface is approached, so the mass flow here is no longer radial. As a result, the dynamic forces ρ[(v ?)v] _θ , which are equal in magnitude to the pressure gradient ? _θ P on the sonic surface, become significant in vertical balance. Therefore, the disk thickness in the supersonic region (and, in particular, near the black-hole horizon) may be assumed to be determined not by the pressure gradient but by the shape of ballistic trajectories.     

A regularization of the three-body problem

Letr ₁,r ₂,r ₃ be arbitrary coordinates of the non-zero interacting mass-pointsm ₁,m ₂,m ₃ and define the distancesR ₁=|r ₁?r ₃|,R ₂=|r ₂?r ₃|,R=|r ₁?r ₂|. An eight-dimensional regularization of the general three-body problem is given which is based on Kustaanheimo-Stiefel regularization of a single binary and possesses the properties:

1. The equations of motion are regular for the two-body collisionsR ₁→0 orR ₂→0.
2. Provided thatR?R ₁ orR?R ₂, the equations of motion are numerically well behaved for close triple encounters.

Although the requirementR? min (R ₁,R ₂) may involve occasional transformations to physical variables in order to re-label the particles, all integrations are performed in regularized variables. Numerical comparisons with the standard Kustaanheimo-Stiefel regularization show that the new method gives improved accuracy per integration step at no extra computing time for a variety of examples. In addition, time reversal tests indicate that critical triple encounters may now be studied with confidence. The Hamiltonian formulation has been generalized to include the case of perturbed three-body motions and it is anticipated that this procedure will lead to further improvements ofN-body calculations.     

Electric fields in coronal magnetic loops

We analyze spectra taken with the 40 cm coronograph at Sacramento Peak Observatory, for evidence of Stark effect on Balmer lines formed in coronal magnetic structures. Several spectra taken near the apex of a bright post-flare loop prominence show significant broadening from H₁₀ to the limit of Balmer line visibility in these spectra, at about H₂₀ The most likely interpretation of the increasing width is Stark broadening, although unresolved blends of Balmer emissions with metallic lines could also contribute to the trend. Less significant broadening is seen in 3 other post-flare loops, and the data from 5 other active coronal condensations observed in this study show no broadening tendency at all, over this range of Balmer number. The trend clearly observed in one post-flare loop requires an ion density of n _i ? 2 × 10¹² cm^?3, if it is to be explained entirely as Stark effect caused by pressure broadening. But mean electron densities measured directly from the Thomson scattering at λ3875 in the same SPO spectra, yield n _e ? 3?7 × 10¹⁰ cm^?3 for the same condensations observed within that loop. Comparison of this evidence from electron scattering, with densities derived from emission measures and line-intensity ratios, argues against a volume filling factor small enough to reconcile the values of n _i and n _e derived in this study. This discrepancy leads us to suggest that the Stark effect observed in these loops, and possibly also in flares, could be caused by macroscopic electric fields, rather than by pressure broadening. The electric field required to explain the Stark broadening in the brightest post-flare loop observed here is approximately 170 V cm^?1. We suggest an origin for such an electric field and discuss its implications for coronal plasma dynamics.     

Multi-Wavelength Eclipse Observations of a Quiescent Prominence

We construct the maps of temperatures, geometrical thicknesses, electron densities and gas pressures in a quiescent prominence. For this we use the RGB signal of the prominence visible-light emission detected during the total solar eclipse of 1 August 2008 in Mongolia and quasi-simultaneous Hα spectra taken at Ond?ejov Observatory. The method of disentangling the electron density and geometrical (effective) thickness was described by Jej?i? and Heinzel (Solar Phys. 254, 89?–?100, 2009) and is used here for the first time to analyse the spatial variations of prominence parameters. For the studied prominence we obtained the following range of parameters: temperature 6000?–?15?000 K, effective thickness 200?–?15000 km, electron density 5×10⁹?–?10¹¹ cm^?3 and gas pressure 0.02?–?0.2 dyn?cm^?2 (assuming a fixed ionisation degree n _p/n _H=0.5). The electron density increases towards the bottom of the prominence, which we explain by an enhanced photoionisation due to the incident solar radiation. To confirm this, we construct a two-dimensional radiative-transfer model with realistic prominence illumination.     

On the velocity of jets powering hot spots similar to those in Cygnus A

Hot spots similar to those in the radio galaxy Cygnus A can be explained by the strong shock produced by a supersonic but classical jet \(\left( {u_{jet}< c/\sqrt 3 } \right)\) . The high integrated radio luminosity (L?2×10⁴⁴ erg s^?1) and the strength of mean magnetic field (B?2×10^?4 G) suggest the hot spots are the downstream flow of a very strong shock which generates the ultrarelativistic electrons of energy ?≥20 MeV. The fully-developed subsonic turbulence amplifies the magnetic field of the jet up to 1.6×10^?4 G by the dynamo effect. If we assume that the post-shock pressure is dominated by relativistic particles, the ratio between the magnetic energy density to the energy density in relativistic particles is found to be ?2×10^?2, showing that the generally accepted hypothesis of equipartition is not valid for hot spots. The current analysis allows the determination of physical parameters inside hot spots. It is found that:

1. The velocity of the upstream flow in the frame of reference of the shock isu ₁?0.2c. Radio observations indicate that the velocity of separation of hot spots isu _sep?0.05c, so that the velocity of the jet isu _jet=u ₁+u _sep?0.25c.
2. The density of the thermal electrons inside the hot spot isn ₂?5×10^?3 e ^? cm^?3 and the mass ejected per year to power the hot spot is ?4M ₀yr^?1.
3. The relativistic electron density is less than 20% of the thermal electron density inside the hot spot and the spectrum is a power law which continues to energies as low as 30 MeV.
4. The energy density of relativistic protons is lower than the energy density of relativistic electrons unlike the situation for cosmic rays in the Galaxy.

     

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.     

A Nonlinear Force-Free Magnetic Field Approximation Suitable for Fast Forward-Fitting to Coronal Loops. II. Numeric Code and Tests

Based on a second-order approximation of nonlinear force-free magnetic field solutions in terms of uniformly twisted field lines derived in Paper I, we develop here a numeric code that is capable to forward-fit such analytical solutions to arbitrary magnetogram (or vector magnetograph) data combined with (stereoscopically triangulated) coronal loop 3D coordinates. We test the code here by forward-fitting to six potential field and six nonpotential field cases simulated with our analytical model, as well as by forward-fitting to an exactly force-free solution of the Low and Lou (Astrophys. J. 352, 343, 1990) model. The forward-fitting tests demonstrate: i) a satisfactory convergence behavior (with typical misalignment angles of μ≈1^°?–?10^°), ii) relatively fast computation times (from seconds to a few minutes), and iii) the high fidelity of retrieved force-free α-parameters (α _fit/α _model≈0.9?–?1.0 for simulations and α _fit/α _model≈0.7±0.3 for the Low and Lou model). The salient feature of this numeric code is the relatively fast computation of a quasi-force-free magnetic field, which closely matches the geometry of coronal loops in active regions, and complements the existing nonlinear force-free field (NLFFF) codes based on photospheric magnetograms without coronal constraints.     

Magnetic field versus gas density,in different physical conditions

A search is made here for possible variations in the behavior of magnetic field valuesB at various gas density valuesn, when comparing low density gas versus high density gas, and when comparing compressed gas versus quiescent gas.

1. For thequiescent microturbulent interstellar gas (e.g., clouds, interclumps — see TableI), the statistical relationB ~n ^k yieldsk = 0.46 ± 0.07 forhigh gas densityn > 100 cm^?3, andk = 0.17 ± 0.03 forlow gas densityn < 100 cm^?3 (see Figure 1).
2. For thecompressed macroturbulent interstellar gas (e.g., masers, expanding shells — see Table II), the statistical relationB ~n ^K yieldsK = 0.61 ± 0.09 forhigh gas densityn > 100 cm^?3, andK = 0.37 ± 0.2 forlow gas densityn < 100 cm^?3 (see Figure 2).
3. The separation betweenlow density gas andhigh density gas is statistically significant. The 2 different physical behaviors (below and above the break at 100 cm^?3) are confirmed statistically (about 2 to 4 sigma away for the quiescent gas alone; about 3 to 6 sigma away for the combined quiescent plus compressed gas).
4. The separation betweencompressed gas andquiescent gas is not statistically significant now (see Figure 3). Atn > 100 cm^?3, a comparison of quiescent gas versus compressed gas shows no statistically significant differences in behavior (they are only about 1 sigma away). Atn < 100 cm^?3, a comparison of quiescent versus compressed gas also shows no statistically significant differences in behavior (less than 1 sigma away).
5. A relation between the densityn and the galactic-wide Star Formation Rate (SFR) can be made for galactic magnetic fields, i.e.: (SFR) ~n ⁿ . For galactic-wide parameters using quiescent, low densityn < 100 cm^?3, and the known relationshipsB ~n ^k/j withk = 0.17,B ~ (SFR) ^j withj = 0.13, then one gets here a lawSFR ~n ^k/j with an exponentk/j = 1.3. This is in rough accord with known data for the Milky Way and for NGC6946.

     

Frequency Fine Structures of Type III Bursts Due to Localized Medium-Scale Density Structures Along Paths of Type III Beams

Predictions from large-scale kinetic simulations are presented for the effects on coronal type III bursts of localized, medium-scale, enhanced density structures superposed on the coronal background along the paths of type III beams. The simulations show that these density structures can produce pronounced frequency fine structures in type III spectra. Flux intensifications and reductions of f _p and 2f _p emission relative to those for the unperturbed background corona occur at frequencies corresponding to the density structures, where f _p is the local electron plasma frequency. Frequency fine structures that are intense, slowly drifting, and narrowband, and thus resemble the characteristics of stria bursts, are predicted for the 2f _p emission. The 2f _p results are consistent with the qualitative proposal of Takakura and Yousef (Solar Phys. 40, 421, 1975) for the interpretation of stria/type IIIb bursts. However, the predicted f _p emission is much weaker than the 2f _p emission and generally below observable levels, and the predicted frequency fine structures do not always show stria characteristics. The predictions are thus inconsistent with the qualitative suggestion of Takakura and Yousef and the interpretations of many observers that stria bursts occur more often in f _p than in 2f _p emission. The significant discrepancies for f _p emission between our numerical calculations and the qualitative proposition of Takakura and Yousef (1975) are mainly caused by: i) differences in the detailed emission processes, ii) neglect of scattering of f _p emission off small-scale density fluctuations by Takakura and Yousef (1975), and iii) other simplifications made in both works. Possible improvements to the simulations are discussed, including improvements to the emission processes and the coronal and beam conditions (e.g., beam speed), in order to produce realistic stria/type IIIb bursts in f _p emission.